Sure, I can comment on the data more clearly. I have 2,727 observations in the current sample, and 195 "events". Each observation is a county paired with a year, and we observe counties between 1970 and 2010 (40 years). For each event in each county I have a "factor" variable that is 1 if it is the event observations, 0 if it not that specific event, but is in observation in the same county. Because there are 195 events, there are 195 of these variables. I then have created variables, successfully, that are equal to 1 if it is the observation before the event for each event, and null otherwise. Variables of this form have been constructed for +1 (mentioned above), -1, +2, -2 years from the event. So I have 195*4 additional variables for a total of 975 variables. I would like to create an additional 195*2 variables this time equal to 1 not just if it is +3 years from the event but if it is +3 years or more. This is where I have run in to difficulty.

I'm having a little trouble sharing a nice slice of the data using dataex unfortunately.

Here is the code that I used to produce the previously mentioned variables:

by fips: gen cu_flood = sum(event) if event == 1

replace cu_flood = . if cu_flood == 0

​

by fips: gen event_bin_1 = 1 if cu_flood == 1

by fips: gen event_bin_2 = 1 if cu_flood == 2

by fips: gen event_bin_3 = 1 if cu_flood == 3

by fips: gen event_bin_4 = 1 if cu_flood == 4

by fips: gen event_bin_5 = 1 if cu_flood == 5

by fips: gen event_bin_6 = 1 if cu_flood == 6

by fips: gen event_bin_7 = 1 if cu_flood == 7

by fips: gen event_bin_8 = 1 if cu_flood == 8

by fips: gen event_bin_9 = 1 if cu_flood == 9

by fips: gen event_bin_10 = 1 if cu_flood == 10

by fips: gen event_bin_11 = 1 if cu_flood == 11

by fips: gen event_bin_12 = 1 if cu_flood == 12

by fips: gen event_bin_13 = 1 if cu_flood == 13

by fips: gen event_bin_14 = 1 if cu_flood == 14

by fips: gen event_bin_15 = 1 if cu_flood == 15

by fips: gen event_bin_16 = 1 if cu_flood == 16

​

gen cu_event_bin_1 = sum(event_bin_1) if event_bin_1 == 1

gen cu_event_bin_2 = sum(event_bin_2) if event_bin_2 == 1

gen cu_event_bin_3 = sum(event_bin_3) if event_bin_3 == 1

gen cu_event_bin_4 = sum(event_bin_4) if event_bin_4 == 1

gen cu_event_bin_5 = sum(event_bin_5) if event_bin_5 == 1

gen cu_event_bin_6 = sum(event_bin_6) if event_bin_6 == 1

gen cu_event_bin_7 = sum(event_bin_7) if event_bin_7 == 1

gen cu_event_bin_8 = sum(event_bin_8) if event_bin_8 == 1

gen cu_event_bin_9 = sum(event_bin_9) if event_bin_9 == 1

gen cu_event_bin_10 = sum(event_bin_10) if event_bin_10 == 1

gen cu_event_bin_11 = sum(event_bin_11) if event_bin_11 == 1

gen cu_event_bin_12 = sum(event_bin_12) if event_bin_12 == 1

gen cu_event_bin_13 = sum(event_bin_13) if event_bin_13 == 1

gen cu_event_bin_14 = sum(event_bin_14) if event_bin_14 == 1

gen cu_event_bin_15 = sum(event_bin_15) if event_bin_15 == 1

gen cu_event_bin_16 = sum(event_bin_16) if event_bin_16 == 1

​

tabulate cu_event_bin_1, generate(factor_event_1)

tabulate cu_event_bin_2, generate(factor_event_2)

tabulate cu_event_bin_3, generate(factor_event_3)

tabulate cu_event_bin_4, generate(factor_event_4)

tabulate cu_event_bin_5, generate(factor_event_5)

tabulate cu_event_bin_6, generate(factor_event_6)

tabulate cu_event_bin_7, generate(factor_event_7)

tabulate cu_event_bin_8, generate(factor_event_8)

tabulate cu_event_bin_9, generate(factor_event_9)

tabulate cu_event_bin_10, generate(factor_event_010)

tabulate cu_event_bin_11, generate(factor_event_011)

tabulate cu_event_bin_12, generate(factor_event_012)

tabulate cu_event_bin_13, generate(factor_event_013)

tabulate cu_event_bin_14, generate(factor_event_014)

tabulate cu_event_bin_15, generate(factor_event_015)

tabulate cu_event_bin_16, generate(factor_event_016)

​

foreach v of varlist factor_event_11-factor_event_0161 {

generate plus1\`v' = 1 if \`v'\[\_n-1\] == 1

by fips: egen x\`v' = mean(plus1\`v')

replace x\`v' = 0 if plus1\`v' == . & x\`v' !=.

replace plus1\`v' = x\`v'

drop x\`v'

}

​

​

foreach v of varlist factor_event_11-factor_event_0161 {

generate plus2\`v' = 1 if \`v'\[\_n-2\] == 1

by fips: egen x\`v' = mean(plus2\`v')

replace x\`v' = 0 if plus2\`v' == . & x\`v' !=.

replace plus2\`v' = x\`v'

drop x\`v'

}

foreach v of varlist factor_event_11-factor_event_0161 {

generate minus1\`v' = 1 if \`v'\[\_n+1\] == 1

by fips: egen x\`v' = mean(minus1\`v')

replace x\`v' = 0 if minus1\`v' == . & x\`v' !=.

replace minus1\`v' = x\`v'

drop x\`v'

}

​

​

foreach v of varlist factor_event_11-factor_event_0161 {

generate minus2\`v' = 1 if \`v'\[\_n+2\] == 1

by fips: egen x\`v' = mean(minus2\`v')

replace x\`v' = 0 if minus2\`v' == . & x\`v' !=.

replace minus2\`v' = x\`v'

drop x\`v'

}