関数 Strata では、各層から描画されるユニットの数を宣言し、それを実行すると奇妙な結果が得られます。ユニットの数が私が申告したものと異なります。
次に例を示します。
library(sampling)
data(swissmunicipalities)
st=strata(swissmunicipalities,stratanames=c("REG"),size=c(30,20,45,15,20,11,44),method="srswor")
結果:
> table(st$REG)
1 2 3 4 5 6 7
20 15 45 30 20 11 44
次のようにする必要があります。
30,20,45,15,20,11,44
アンサーにとても感謝しています。