Dear Stata Forum,
I have data with a small number of clusters but with a large number of individuals within a class.
I have four household surveys over 4 countries and 26 regions (Pooled-Crosssection).

To overcome the overstated standard errors I´d like to bootstrap my results.
However, a regional dummy variable is causing me trouble.
I already know that sometimes a region is included in a bootstrap sample and sometimes left out.
That is alright since I am not interested in the coefficients of the dummy variable.

I am using Stata 16.

Below you can see my data:
Code:
 Stata 16
input double rem_edu float(handling_edu1 ln_exp_total ln_remittance) double(rem_freq oecd femmigrant childprim childsec childtert elderlyageabove62 womenabove15) float(hh_edu head_gender) byte(hagriland urban) float(incomepoverty school) byte region
                  0 .03949346  6.050091  4.544163  5 0 1 0 0 0 1 1 2 0 1 1 .5724485 .9207876 11
                  0 .03949346   6.18593  5.642776  3 1 1 0 0 0 1 1 1 0 1 1 .5724485 .9207876 11
   .800000011920929 .03949346  6.777756  6.510276 10 0 0 0 0 0 0 4 3 1 1 1 .5724485 .9207876 11
                  1 .03949346  6.255384  6.510276  5 0 1 2 0 0 0 1 1 1 1 1 .5724485 .9207876 11
                  0 .03949346  5.062107  5.370842  5 1 1 3 2 0 0 2 1 1 1 1 .5724485 .9207876 11
                  0 .03949346  5.837912  7.818609 16 1 1 2 0 0 0 3 3 1 0 1 .5724485 .9207876 11
                .75 .03949346  7.687905  8.685028  3 1 0 2 1 0 0 3 3 0 0 1 .5724485 .9207876 11
                  0 .03949346  6.927889  5.163203  4 1 1 0 0 0 1 3 2 0 1 1 .5724485 .9207876 11
                  0         0  4.999174  5.999451  2 1 1 1 0 0 0 1 2 0 1 0 .6428978 .7905115 12
                  0         0  4.515176         .  3 1 1 0 0 0 0 1 2 0 1 1 .6428978 .7905115 12
  .4000000059604645         0  8.181123  7.203424  2 1 1 0 0 0 0 3 2 0 1 0 .6428978 .7905115 12
  .5882353186607361         0  6.092485  7.916373 12 0 0 3 0 0 0 3 1 1 1 0 .6767302        1 12
 .44117647409439087         0  4.240481  6.181772 12 0 0 5 2 0 0 1 1 1 1 0 .6767302        1 12
and this is the regression command I use:
[CODE][xtset country
bootstrap, nodrop cluster(country) idcluster(countryb): xtreg rem_edu handling_edu1 ln_exp_total ln_remittance rem_freq oecd femmigrant childprim childsec childtert elderlyageabove62 womenabove15 i.hh_edu head_gender i.hagriland i.urban incomepoverty school i.region, re vce(cluster country)/CODE]

Is there any solution to the dummy problem?

Kind regards
Alexander