从已知数据列中进行放回(或不放回)抽样,再统计每个样本的均值及分布,这个脚本如何编辑?例如,有2000个数据,每次随机抽取20个数据计算其均值,一共抽取100次,计算这100个样本的均值以及分布状况;
以下是一个用于抽取 100 个子集并计算其均值的脚本:
dt = Open( "$SAMPLE_DATA/Aircraft Incidents.jmp" );
:Latitude << Format("Best");
dt << Select Where( Is Missing( :Latitude ) ) << Delete Rows;<BR /><BR />runs = 100;<BR />//runs= 100000;
results = new table("results", add rows (runs), new column ("mean"));
rows = index (1, nrows(dt));
for each({i},index(1,runs),
myRows = Random Shuffle( rows)[1::20];
data = dt[myRows, "Latitude"];
results[i,"mean"]=Mean(data);
);
my Mean = Col Mean(:mean);
Show(myMean);
……以及均值的均值。
使用此脚本,您可以使用 Distribution 平台将 100 个平均值与“所有”可用分布进行拟合:
dist = results << Distribution(Continuous Distribution(Column( :mean )); );
dist << Fit All
