從已知資料列中進行放回(或不放回)抽樣,再統計每個樣本的平均數及分佈,這個腳本如何編輯?例如,有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
