分布プラットフォームからの幾何平均の信頼区間こんにちは、
分散プラットフォームから幾何平均の信頼区間を取得するショートカットはありますか?
基本的なデータ分析とモデリングデータの探索と視覚化Re: 分布プラットフォームからの幾何平均の信頼区間こんにちは、イアンさん
それは最高です。 ご提案ありがとうございます。
Re: 分布プラットフォームからの幾何平均の信頼区間また、間隔だけを気にする場合は、この種のアプローチを使用して間隔を取得できます。ブートストラップの汎用性に必要なすべてのプラットフォーム起動のオーバーヘッドがないため、はるかに高速です。
NamesDefaultToHere(1);
// Get some data
dt = Open("$SAMPLE_DATA/Big Class.jmp");
n = NRow(dt);
dataVec = dt[0, "height"]; // Get the requisite data as a vector
nb = 100000; // Number of bootstrsp samples
gms = J(nb, 1, .); // Vector to hold the geometric mean of each bootstrap sample
for(s=1, s<=nb, s++, // Loop over bootstrap samples
pickMe = J(n, 1, RandomInteger(1, n)); // Indices of elements to compose this bootstrap sample
bootstrapSampleVec = dataVec[pickMe]; // Values in this bootstrap sample
gm = exp(Mean(ln(bootstrapSampleVec))); // Geometric mean of this bootstrap sample
gms[s] = gm;
);
gm95 = Quantile(0.95, gms);
gm05 = Quantile(0.05, gms);
Print("The Geometric Mean lies betweeen "||Char(Round(gm05, 3))||" and "||Char(Round(gm95, 3))||".");
おそらく、このコードを並列化してさらに高速化することもできますが、読みやすさは多少失われます。
Re: 分布プラットフォームからの幾何平均の信頼区間詳細をガイダンスと比較しませんでした。 @Mark_Bailey引用符で囲まれていますが、ブートストラップも (対話的にまたは JSL 経由で) 試すことができます。 例えば：
NamesDefaultToHere(1);
// Get some data
dt = Open("$SAMPLE_DATA/Big Class.jmp");
// Use Distribution
dist = dt << Distribution(
Continuous Distribution(
Column( :height ),
Quantiles( 0 ),
Vertical( 0 ),
Outlier Box Plot( 0 ),
Customize Summary Statistics(
Mean( 0 ),
Std Dev( 0 ),
Std Err Mean( 0 ),
Upper Mean Confidence Interval( 0 ),
Lower Mean Confidence Interval( 0 ),
N( 0 ),
N Missing( 0 ),
Geometric Mean( 1 )
)
)
);
// Bootstrp the statistic in question (note 'dt2' is a list)
dt2 = Report(dist)[TableBox(1)] << Bootstrap(1000);
// Inspect the bootstrp distribution
dist2 = dt2[2] << runScript("Distribution");
// Get the confidence interval for the statistic by inspection or further JSL . . .Re: 分布プラットフォームからの幾何平均の信頼区間別の方法では、JSL を使用せず、対数正規分布について同じ仮定を立てます。配信プラットフォームを利用してデータの自然対数を解析します。 取得した平均値と信頼限界は、Exp を介して逆変換できます。 結果として得られる区間は、GM の信頼区間です。
Re: 分布プラットフォームからの幾何平均の信頼区間ハイマーク、
素早いフィードバックをいただきまして誠にありがとうございます。
Re: 分布プラットフォームからの幾何平均の信頼区間その結果は、配布プラットフォームまたは JMP の他の場所からは入手できません。 計算は、このガイドのガイダンスを使用してスクリプト化できます。 DO ループのブログ投稿。 SAS/IML は JMP JSL に似ているため、変換はそれほど悪くありません。
