Hi All,
Still hoping for some help to resolve a bootstrap error message. Here is the code I am trying to run with the error message:
Without the bootstrapping the competing risk regression runs fine. d2 is the failure variable with d2=1 as the main failure event and d2=2 being the competing risk event, with the following distribution:
. tab d2, miss
d2 | Freq. Percent Cum.
------------+-----------------------------------
0 | 47,948 94.98 94.98
1 | 158 0.31 95.30
2 | 2,374 4.70 100.00
------------+-----------------------------------
Total | 50,480 100.00
Could the error message be because the failure event is rare (i.e., only 0.31%)? I have also tried more reps (up to a million) and get the same error message.
I would be grateful for any help you may be able to offer.
Many thanks in advance,
Sumedha
P.S. Here is a dataex in case helpful.
Still hoping for some help to resolve a bootstrap error message. Here is the code I am trying to run with the error message:
Code:
. capture program drop myboot . program define myboot, rclass 1. . * first stage . reghdfe COVID Female age75plus Asian Black Hispanic lnPopDensity tavg lnprcp CumMonthsSAH lnnursing_visits deaths_rate cum_ > num_vacpct, absorb(grpatid monthlydate state, save) cluster(grpatid) residuals(resid) 2. predict double COVID_fe, r 3. . * second stage . eststo: stcrreg Female age75plus Asian Black Hispanic COVID COVID_fe, compete(d2=2) 4. estimates store twoSLS1 5. . end . xtset, clear . bootstrap _b, reps(1000) seed(123): myboot (running myboot on estimation sample) Bootstrap replications (1,000): xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx > xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx > xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx > xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx > xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx > xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx > xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx > xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx > xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx done x: Error occurred when bootstrap executed myboot. insufficient observations to compute bootstrap standard errors no results will be saved r(2000); end of do-file r(2000);
. tab d2, miss
d2 | Freq. Percent Cum.
------------+-----------------------------------
0 | 47,948 94.98 94.98
1 | 158 0.31 95.30
2 | 2,374 4.70 100.00
------------+-----------------------------------
Total | 50,480 100.00
Could the error message be because the failure event is rare (i.e., only 0.31%)? I have also tried more reps (up to a million) and get the same error message.
I would be grateful for any help you may be able to offer.
Many thanks in advance,
Sumedha
P.S. Here is a dataex in case helpful.
Code:
* Example generated by -dataex-. For more info, type help dataex clear input float(COVID Female age75plus) byte(Asian Black Hispanic) float(lnPopDensity tavg lnprcp CumMonthsSAH lnnursing_visits deaths_rate cum_num_vacpct grpatid monthlydate) str24 state float d2 byte(_st _d _t _t0) 0 0 0 0 0 0 4.597138 26.89387 4.625463 0 11.36908 0 0 7389284 714 "AL" 0 1 0 1 0 0 1 1 0 0 0 4.597138 26.89387 4.625463 0 11.36908 0 0 5707464 714 "AL" 0 1 0 1 0 0 1 1 1 0 0 4.597138 26.89387 4.625463 0 11.36908 0 0 7312555 714 "AL" 2 1 0 1 0 0 1 1 0 0 0 4.597138 26.89387 4.625463 0 11.36908 0 0 6168903 714 "AL" 2 1 0 1 0 0 1 0 0 1 0 4.597138 26.89387 4.625463 0 11.36908 0 0 7009023 714 "AL" 0 1 0 1 0 0 1 1 0 1 0 4.597138 26.89387 4.625463 0 11.36908 0 0 6340785 714 "AL" 0 1 0 1 0 0 0 0 0 1 0 4.597138 26.89387 4.625463 0 11.36908 0 0 861670 714 "AL" 2 1 0 1 0 0 1 0 0 1 0 4.597138 26.89387 4.625463 0 11.36908 0 0 1170181 714 "AL" 0 1 0 1 0 0 0 1 0 0 0 4.597138 26.89387 4.625463 0 11.36908 0 0 1921660 714 "AL" 0 1 0 1 0 0 1 0 0 0 0 4.597138 26.89387 4.625463 0 11.36908 0 0 2971145 714 "AL" 0 1 0 1 0 0 1 1 0 1 0 4.597138 26.89387 4.625463 0 11.36908 0 0 706786 714 "AL" 0 1 0 1 0 0 1 0 0 1 0 4.597138 26.89387 4.625463 0 11.36908 0 0 4296355 714 "AL" 0 1 0 1 0 0 0 0 0 1 0 4.597138 26.89387 4.625463 0 11.36908 0 0 7212486 714 "AL" 0 1 0 1 0 0 1 0 0 1 0 4.597138 26.89387 4.625463 0 11.36908 0 0 7042124 714 "AL" 2 1 0 1 0 0 0 0 0 0 0 4.597138 26.89387 4.625463 0 11.36908 0 0 6317893 714 "AL" 2 1 0 1 0 0 1 1 0 0 0 4.597138 26.89387 4.625463 0 11.36908 0 0 6502335 714 "AL" 0 1 0 1 0 0 1 0 0 1 0 4.597138 26.89387 4.625463 0 11.36908 0 0 4968653 714 "AL" 0 1 0 1 0 0 1 1 0 0 0 4.597138 26.89387 4.625463 0 11.36908 0 0 1135025 714 "AL" 0 1 0 1 0 0 1 0 0 0 0 4.597138 26.89387 4.625463 0 11.36908 0 0 7346079 714 "AL" 2 1 0 1 0 0 0 1 0 1 0 4.597138 26.89387 4.625463 0 11.36908 0 0 1927784 714 "AL" 0 1 0 1 0 0 0 0 0 1 0 4.597138 26.89387 4.625463 0 11.36908 0 0 5707029 714 "AL" 0 1 0 1 0 0 1 0 0 1 0 4.597138 26.89387 4.625463 0 11.36908 0 0 7474901 714 "AL" 0 1 0 1 0 0 0 0 0 0 0 4.597138 26.89387 4.625463 0 11.36908 0 0 7564318 714 "AL" 0 1 0 1 0 0 0 1 0 0 0 4.597138 26.844194 4.842059 0 11.405563 0 0 1921660 715 "AL" 0 1 0 2 1 0 1 1 0 0 0 4.597138 26.844194 4.842059 0 11.405563 0 0 1135025 715 "AL" 0 1 0 2 1 0 1 0 0 1 0 4.597138 26.844194 4.842059 0 11.405563 0 0 1170181 715 "AL" 0 1 0 2 1 0 1 1 0 0 0 4.597138 26.844194 4.842059 0 11.405563 0 0 5707464 715 "AL" 0 1 0 2 1 0 1 1 1 0 0 4.597138 26.844194 4.842059 0 11.405563 0 0 7312555 715 "AL" 2 1 0 2 1 0 1 0 0 1 0 4.597138 26.844194 4.842059 0 11.405563 0 0 4968653 715 "AL" 0 1 0 2 1 0 1 0 0 1 0 4.597138 26.844194 4.842059 0 11.405563 0 0 7009023 715 "AL" 0 1 0 2 1 0 1 1 0 1 0 4.597138 26.844194 4.842059 0 11.405563 0 0 706786 715 "AL" 0 1 0 2 1 0 1 0 0 1 0 4.597138 26.844194 4.842059 0 11.405563 0 0 4296355 715 "AL" 0 1 0 2 1 0 1 1 0 0 0 4.597138 26.844194 4.842059 0 11.405563 0 0 6502335 715 "AL" 0 1 0 2 1 0 1 1 0 1 0 4.597138 26.844194 4.842059 0 11.405563 0 0 6340785 715 "AL" 0 1 0 2 1 0 0 0 0 0 0 4.597138 26.844194 4.842059 0 11.405563 0 0 6317893 715 "AL" 2 1 0 2 1 0 0 0 0 1 0 4.597138 26.844194 4.842059 0 11.405563 0 0 861670 715 "AL" 2 1 0 2 1 0 1 0 0 0 0 4.597138 26.844194 4.842059 0 11.405563 0 0 2971145 715 "AL" 0 1 0 2 1 0 0 0 0 1 0 4.597138 26.844194 4.842059 0 11.405563 0 0 5707029 715 "AL" 0 1 0 2 1 0 0 0 0 0 0 4.597138 26.844194 4.842059 0 11.405563 0 0 7564318 715 "AL" 0 1 0 2 1 0 1 0 0 1 0 4.597138 26.844194 4.842059 0 11.405563 0 0 7042124 715 "AL" 2 1 0 2 1 0 1 0 0 0 0 4.597138 26.844194 4.842059 0 11.405563 0 0 7346079 715 "AL" 2 1 0 2 1 0 0 0 0 1 0 4.597138 26.844194 4.842059 0 11.405563 0 0 7212486 715 "AL" 0 1 0 2 1 0 1 1 0 0 0 4.597138 26.844194 4.842059 0 11.405563 0 0 6168903 715 "AL" 2 1 0 2 1 0 0 0 0 0 0 4.597138 26.844194 4.842059 0 11.405563 0 0 7389284 715 "AL" 0 1 0 2 1 0 0 1 0 1 0 4.597138 26.844194 4.842059 0 11.405563 0 0 1927784 715 "AL" 0 1 0 2 1 0 1 0 0 1 0 4.597138 26.844194 4.842059 0 11.405563 0 0 7474901 715 "AL" 0 1 0 2 1 0 1 1 1 0 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 7312555 716 "AL" 2 1 0 3 2 0 0 0 0 0 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 7564318 716 "AL" 0 1 0 3 2 0 0 1 0 1 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 1927784 716 "AL" 0 1 0 3 2 0 1 1 0 0 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 6502335 716 "AL" 0 1 0 3 2 0 0 0 0 1 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 861670 716 "AL" 2 1 0 3 2 0 0 0 0 0 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 7389284 716 "AL" 0 1 0 3 2 0 0 0 0 1 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 7212486 716 "AL" 0 1 0 3 2 0 1 0 0 1 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 7009023 716 "AL" 0 1 0 3 2 0 1 0 0 1 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 7042124 716 "AL" 2 1 0 3 2 0 1 1 0 1 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 6340785 716 "AL" 0 1 0 3 2 0 0 0 0 1 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 5707029 716 "AL" 0 1 0 3 2 0 1 1 0 0 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 5707464 716 "AL" 0 1 0 3 2 0 1 0 0 1 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 7474901 716 "AL" 0 1 0 3 2 0 1 0 0 1 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 1170181 716 "AL" 0 1 0 3 2 0 0 1 0 0 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 1921660 716 "AL" 0 1 0 3 2 0 1 0 0 1 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 4296355 716 "AL" 0 1 0 3 2 0 0 0 0 0 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 6317893 716 "AL" 2 1 0 3 2 0 1 1 0 0 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 6168903 716 "AL" 2 1 0 3 2 0 1 0 0 0 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 2971145 716 "AL" 0 1 0 3 2 0 1 0 0 0 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 7346079 716 "AL" 2 1 0 3 2 0 1 1 0 0 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 1135025 716 "AL" 0 1 0 3 2 0 1 0 0 1 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 4968653 716 "AL" 0 1 0 3 2 0 1 1 0 1 0 4.597138 26.732334 1.9600948 0 11.406064 0 0 706786 716 "AL" 0 1 0 3 2 0 1 0 0 1 0 4.597138 19.743225 5.197447 0 11.375776 0 0 1170181 717 "AL" 0 1 0 4 3 0 1 1 0 0 0 4.597138 19.743225 5.197447 0 11.375776 0 0 6502335 717 "AL" 0 1 0 4 3 0 1 0 0 1 0 4.597138 19.743225 5.197447 0 11.375776 0 0 7042124 717 "AL" 2 1 0 4 3 0 1 0 0 1 0 4.597138 19.743225 5.197447 0 11.375776 0 0 7474901 717 "AL" 0 1 0 4 3 0 0 1 0 0 0 4.597138 19.743225 5.197447 0 11.375776 0 0 1921660 717 "AL" 0 1 0 4 3 0 1 0 0 1 0 4.597138 19.743225 5.197447 0 11.375776 0 0 4296355 717 "AL" 0 1 0 4 3 0 1 1 0 0 0 4.597138 19.743225 5.197447 0 11.375776 0 0 1135025 717 "AL" 0 1 0 4 3 0 1 1 0 0 0 4.597138 19.743225 5.197447 0 11.375776 0 0 6168903 717 "AL" 2 1 0 4 3 0 1 0 0 1 0 4.597138 19.743225 5.197447 0 11.375776 0 0 4968653 717 "AL" 0 1 0 4 3 0 0 0 0 1 0 4.597138 19.743225 5.197447 0 11.375776 0 0 5707029 717 "AL" 0 1 0 4 3 0 0 0 0 0 0 4.597138 19.743225 5.197447 0 11.375776 0 0 7564318 717 "AL" 0 1 0 4 3 0 0 0 0 0 0 4.597138 19.743225 5.197447 0 11.375776 0 0 6317893 717 "AL" 2 1 0 4 3 0 1 0 0 0 0 4.597138 19.743225 5.197447 0 11.375776 0 0 2971145 717 "AL" 0 1 0 4 3 0 1 1 0 1 0 4.597138 19.743225 5.197447 0 11.375776 0 0 706786 717 "AL" 0 1 0 4 3 0 1 0 0 1 0 4.597138 19.743225 5.197447 0 11.375776 0 0 7009023 717 "AL" 0 1 0 4 3 0 1 1 0 1 0 4.597138 19.743225 5.197447 0 11.375776 0 0 6340785 717 "AL" 0 1 0 4 3 0 0 0 0 1 0 4.597138 19.743225 5.197447 0 11.375776 0 0 861670 717 "AL" 2 1 0 4 3 0 1 1 1 0 0 4.597138 19.743225 5.197447 0 11.375776 0 0 7312555 717 "AL" 2 1 0 4 3 0 1 0 0 0 0 4.597138 19.743225 5.197447 0 11.375776 0 0 7346079 717 "AL" 2 1 0 4 3 0 0 1 0 1 0 4.597138 19.743225 5.197447 0 11.375776 0 0 1927784 717 "AL" 0 1 0 4 3 0 1 1 0 0 0 4.597138 19.743225 5.197447 0 11.375776 0 0 5707464 717 "AL" 0 1 0 4 3 0 0 0 0 0 0 4.597138 19.743225 5.197447 0 11.375776 0 0 7389284 717 "AL" 0 1 0 4 3 0 0 0 0 1 0 4.597138 19.743225 5.197447 0 11.375776 0 0 7212486 717 "AL" 0 1 0 4 3 0 0 0 0 0 0 4.597138 10.332 4.4127984 0 11.393387 0 0 7389284 718 "AL" 0 1 0 5 4 0 1 0 0 1 0 4.597138 10.332 4.4127984 0 11.393387 0 0 4968653 718 "AL" 0 1 0 5 4 0 1 0 0 0 0 4.597138 10.332 4.4127984 0 11.393387 0 0 2971145 718 "AL" 0 1 0 5 4 0 1 0 0 1 0 4.597138 10.332 4.4127984 0 11.393387 0 0 7474901 718 "AL" 0 1 0 5 4 0 0 0 0 1 0 4.597138 10.332 4.4127984 0 11.393387 0 0 7212486 718 "AL" 0 1 0 5 4 0 1 0 0 1 0 4.597138 10.332 4.4127984 0 11.393387 0 0 4296355 718 "AL" 0 1 0 5 4 0 0 1 0 1 0 4.597138 10.332 4.4127984 0 11.393387 0 0 1927784 718 "AL" 0 1 0 5 4 0 1 1 0 1 0 4.597138 10.332 4.4127984 0 11.393387 0 0 706786 718 "AL" 0 1 0 5 4 end format %tm monthlydate label values COVID COVID label def COVID 0 "No confirmed COVID-19", modify
Comment