84699_zh-TW

liliconcom

是否可以計算矩陣中的分佈和置信區間

所有，
我期待計算一些數據的置信區間。我知道可以使用以下方法實現此目的： 

 

// Gets me what I want 
dt = Open( "$SAMPLE_DATA/Big Class.jmp" );
Obj = dt << Distribution( Column( :Age, :Weight ) );
Obj << Confidence Interval(0.95);

// How to achieve using matrix
Mat = dt << Get As Matrix; 
Obj = Distribution(Mat[0,1]) << Confidence Interval(0.95);

     如果我想通過小組進一步做到這一點，上面變得特別困難，讓我們假設我是否必須按年齡分組？任何指導表示讚賞 

Re：是否可以計算矩陣中的分佈和置信區間

我不會做矩陣操作，我會做這樣的事情，

 

dt = Open( "$SAMPLE_DATA/Big Class.jmp" );
Obj = dt << Distribution( Column( :height, :weight ), by(:sex), "Invisible" );
Obj << Confidence Interval(0.95);
Report(Obj[1])["Confidence Intervals"][TableBox(1)] << make combined data table;
Obj << close window;

您可以執行矩陣函數，但需要按組（年齡或性別）查找所有唯一，然后索引矩陣。  存在可能感興趣的矩陣CDF功能。  下面的腳本是一個FYI，不推薦，但這裡是公式。生成的矩陣沒有by組或哪個變量。我推薦上面的方法。

 

Mat = dt << Get As Matrix; 
//Obj = Distribution(Mat[0,1]) << Confidence Interval(0.95);



/* CDF function */
/* only one column at a time get weight*/
{Quant, CumProb} = CDF( Mat[0,3] ); 

mxb = V Mean(Mat);
msd = V Std(Mat);
mdf = V Sum( !IsMissing(Mat) ) -1;

xb_lci = mxb - msd:/sqrt(mdf+1) :* tquantile(0.975,mdf);
xb_uci = mxb + msd:/sqrt(mdf+1) :* tquantile(0.975,mdf);

sd_lci =  msd :* sqrt( mdf :/ ChiSquareQuantile(0.975,mdf)); 
sd_uci =  msd :* sqrt( mdf :/ ChiSquareQuantile(0.025, mdf));

show(Quant, CumProb, mxb[3], msd[3], xb_lci[3], xb_uci[3], sd_lci[3], sd_uci[3]);