《应用抽样技术》（第二版）附录数据、程序和PPTchapter5 系统抽样 程序部分电子版.docx

三、利用软件R实现上述过程

使用变量N、n和FWY代表总体单元数、样本量和总体数据，变量k表示抽样间隔。文件“chap5p4.csv”的第1行是每列数据的标题行，第1列“BH”给出了抽取的样本单元编号，第2列“CJDY”是村级单元的代码，第3列“XQBH”是村级单元包含的调查小区编号，第4列“YYSR”是每个调查小区的营业收入。

N-12779;

n-1575;

k-floor(N/n);

FWY-read.csv(chap5p4.csv,header=T);

attach(FWY);

采用简单随机抽样，随机抽取1个小于抽样间隔k的随机起点，存入变量r中，语句为

r-sample(c(1:k),1);

从随机起点开始，每隔抽样间隔k个编号抽取1个编号，相应的调查小区进入样本，保存在变量YB中，相应的语句为

YB-YYSR[(c(1:N)%%k)==(r%%k)][1:n];

对于每个调查小区规模以下服务业企业平均营业收入额的估计，相应的语句为

JZGJ-round(mean(YB),2);

JZFC-round(var(YB)*(1-n/N)/n,4);

JZBZC-round(sqrt(JZFC),2);

变量JZGJ为每个调查小区规模以下服务业企业平均营业收入额的估计值，变量JZFC为简单估计量的方差估计，变量JZBZC为简单估计量的标准差估计。对于某地区12779个调查小区规模以下服务业企业营业收入总额的估计，R语句为

ZZGJ-JZGJ*N;

ZZFC-JZFC*N^2;

ZZBZC-round(sqrt(ZZFC),2);

变量ZZGJ为某地区12779个调查小区规模以下服务业企业营业收入总额的估计值，变量ZZFC为估计量的方差估计，变量ZZBZC为估计量的标准差估计。