95% 置信區間大家好,
我正在嘗試從數據表中獲取摘要,我想要變量 l10t 的平均值、標準誤差和 95% 置信區間。 這一切看起來相當簡單,但我在 95% 置信區間方面遇到了麻煩,嘗試了不同的方法但無濟於事,目前有
dt2<< Summary(
Group( :Sample ),
Mean( :dilution factor, :l10t ),
Freq( "None" ),
Weight( "None" ),
Output Table( "log10t" ),
Std Err( :l10t ),
Confidence Interval( :l10t, 0.95),
Link to original data table(0)
我顯然做錯了什麼,任何人都可以提供一些建議。
謝謝堆
回复:95% 置信區間嗨姆茲瓦爾德,
我已經使用了以下並且效果很好
New Column( "Lower 95% CI", Numeric, Continuous, Formula ( Round(:Name ("Mean(l10t)") - (:Name("Std Err(l10t)") * t Quantile (0.975,23)), 4)));如果我使用以下命令創建變量
dt2<< Summary(
Group( :Sample ),
Mean( :dilution factor, :l10t ),
Freq( "None" ),
Weight( "None" ),
Output Table( "log10t" ),
Std Err( :l10t ),
Link to original data table(0)
);
我想將變量 Mean(l10t) 四捨五入到小數點後四位,你能從上面的語法塊中做到這一點嗎???
謝謝
米奇男孩
回复:95% 置信區間嗯,看起來 Summary 目前不支持 CI,但您始終可以自己為它們制定一個公式。
95% 置信區間(假設正態分佈)的公式為:
Upper95CI = 平均值 + 1.96*stdev
Lower95CI = 平均值 - 1.96*stdev
您可以根據摘要平台返回的字段生成這兩列。 這些公式將使您非常接近您正在尋找的東西。
回复:95% 置信區間嗨姆茲瓦爾德,
非常感謝您的幫助和建議,我嘗試了這段代碼
platform = dt2 << Oneway(
Y( :l10t ),
X( :fold dilution ),
Means and Std Dev( 1 ),
Box Plots( 1 ),
Mean Diamonds( 1 ),
Mean Error Bars( 1 ),
Std Dev Lines( 1 ),
Connect Means( 1 ),
Points Jittered( 1 ),
Grand Mean( 0 )
);
Wait( 0 );
Report( platform )[Outline Box( "Oneway Analysis of height By sex" )][
Outline Box( "Means and Std Deviations" )][Table Box( 1 )] << Make Into Data Table;
Report( platform ) << Close Window;
並收到以下錯誤消息
Cannot subscript Display Box in access or evaluation of 'Report(platform)[Outline Box("Oneway Analysis of height By sex")]' , Report( platform )[/*###*/Outline Box( "Oneway Analysis of height By sex" )]然而,我幾乎已經知道我想用這個來為我的數據表實現什麼,
dt2<< Summary(
Group( :Sample ),
Mean( :dilution factor, :l10t ),
Freq( "None" ),
Weight( "None" ),
Output Table( "log10t" ),
Std Err( :l10t ),
Confidence Interval( :Name ("Mean(l10t)")),
Link to original data table(0)
);
Loq = Data Table( "log10t" );
Loq << New Column( "Back T'sform Mean", Numeric, Continuous, Formula( Round(10^:Name ("Mean(l10t)"),0)));
Loq << New Column( "Back T'sform Std error", Numeric, Continuous, Formula(10^:Name ("Std Err(l10t)"), Format( "fixed dec", 6, 0 ) ) );給我這個,
所以我就快到了,我希望置信區間是一行簡單的代碼。
謝謝
米奇男孩
回复:95% 置信區間您可以完成類似的事情:
Names Default To Here( 1 );
dt = Open( "$SAMPLE_DATA/Big Class.jmp" );
platform = dt << Oneway(
Y( :height ),
X( :sex ),
Means and Std Dev( 1 ),
Box Plots( 1 ),
Mean Diamonds( 1 ),
Mean Error Bars( 1 ),
Std Dev Lines( 1 ),
Connect Means( 1 ),
Points Jittered( 1 ),
Grand Mean( 0 )
);
Wait( 0 );
Report( platform )[Outline Box( "Oneway Analysis of height By sex" )][
Outline Box( "Means and Std Deviations" )][Table Box( 1 )] << Make Into Data Table;
Report( platform ) << Close Window;
將您的組變量放入 X() 角色中。