I have data on fire incidences:
structure(list(Season = c("Winter", "Winter", "Winter", "Winter", "Winter", "Winter", "Winter", "Winter", "Winter", "Winter", "Winter", "Winter", "Winter", "Winter", "Winter", "Winter", "Winter", "Winter", "Winter", "Winter", "Winter", "Winter", "Winter", "Winter", "Summer", "Summer", "Summer", "Summer", "Summer", "Summer", "Summer", "Summer", "Summer", "Summer", "Summer", "Summer", "Summer", "Summer", "Summer", "Summer", "Summer", "Summer", "Summer", "Summer", "Summer", "Summer", "Summer", "Summer", "Summer"), Time = c("Early", "Early", "Early", "Early", "Early", "Early", "Early", "Early", "Early", "Early", "Early", "Early", "Late", "Late", "Late", "Late", "Late", "Late", "Late", "Late", "Late", "Late", "Late", "Late", "Early", "Early", "Early", "Early", "Early", "Early", "Early", "Early", "Early", "Early", "Early", "Early", "Early", "Early", "Early", "Early", "Late", "Late", "Late", "Late", "Late", "Late", "Late", "Late", "Late"), Area = c("A1", "A1", "A1", "A1", "A1", "A2", "A2", "A2", "A2", "A2", "A3", "A3", "A3", "A3", "A2", "A2", "A2", "A2", "A2", "A1", "A1", "A1", "A1", "A1", "A1", "A1", "A1", "A1", "A1", "A2", "A2", "A2", "A2", "A2", "A3", "A3", "A3", "A3", "A2", "A2", "A2", "A2", "A2", "A1", "A1", "A1", "A1", "A1", "A1"), Small_fire = c(0L, 0L, 1L, 0L, 1L, 0L, 0L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 0L, 1L, 0L, 1L, 1L, 1L ), Big_fire = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L)), class = "data.frame", row.names = c(NA, -49L))
Here, we have variables - seasons (summer & winter), time taken by fire truck to arrive (Early (<30 mins) and Late (>30 mins)), area (area code), small fire incident (1=Yes, 0=No), and big fire incident (1=Yes, 0=No).
I have a few questions in mind:
- Whether the frequency of fire incidence increase in the summer compared to winter? (I guess, I could use logistic regression for this).
- Whether the frequency of small fire converting to big fire incidence increase in summer while adjusting for the arrival time of arrival truck and area?
I would be grateful if someone could suggest a way to analyze this data in R or an R function that I can use.
Thanks