どのようにテストするかどうかを


9

平均して、私は、3つの独立したグループがあるとμ1, μ2, μ3それぞれ。

どのようにしてテストすることができるかどうかをμ1<μ2<μ3、または使用していないn1, n2, n3、各グループからのサンプル?

詳細な計算ではなく、一般的な方法論を知りたい。仮説H0H1設定方法がわかりませんでした。


1
これらは、順序が制限された統計的推論の場合です。あるトピックに関する本は
kjetil b halvorsen

1
注文制限(1973)の下での Barlow、Bartholemew、Bremner、およびBrunk Statistical推論による古い本もあります(それ以降、いくつかの進展がありました)。ノンパラメトリック検定に関する限り、Jonckheere-Terpstra検定(例:Conoverを参照)と1つの一致検定(NeaveとWorthingtonによる本を試してください)があります。通常、等価nullと順序付けられた代替を記述します。
Glen_b


ここで1は1つが持っている、ではないと言うべきでのグループからのサンプルを、私は 1の大きさのサンプルを持っていることが、N Iグループから私はnii,nii.
Michael Hardy

回答:


8

統計では、「Xが真かどうか」をテストすることはできません。帰無仮説が偽であるという証拠を見つけることのみを試みることができます。

さんがあなたの帰無仮説があるとしましょう

H01:μ1<μ2<μ3.
レッツ・また、あなたがベクトル推定する方法があることを前提とμ=(μ1,μ2,μ3)。物事を保つためには、単にあなたが推定持っていることを前提とし
xN(μ,Σ),
Σある3×3共変量行列を。帰無仮説を次のように書き直すことができます。
Aμ<0,
A=[110011].
このショーは、あなたの帰無仮説は、ベクター上の不等式制約として表現することができるというAμ。天然の推定Aμによって与えられ、
AxN(Aμ,AΣA).
これで、フレームワークを使用して、以下に示す法線ベクトルの不等式制約をテストできます

工藤昭夫(1963)「片側検定の多変量類似体」。で:Biometrika 50.3 / 4、ページ403–418。

n1,n2,n3Σ

(σ12/n1,σ22/n2,σ32/n3),
σk2k=1,2,3

H12:μ1<μ2<μ3
H02:NOT H1.
H1:Aμ<0
H02:there exists a k=1,2 such that (Aμ)k0.
H0,12:(Aμ)10.
H0,22:(Aμ)20.
H0H0H0Σ

H02

H02:maxk=1,2(Aμ)k0.
maxAx


結局のところ、私は自分の答えを編集しました。
Andreas Dzemski、

xμ^μ

1

@ andreas-dzemskiによって提供される答えは、データが通常に配布されていることがわかっている場合にのみ正しいです。

分布がわからない場合は、ノンパラメトリック検定を実行することをお勧めします。この場合、最も単純なものは順列テストを実行するようです。 これはトピックについての本であり、これは素晴らしいオンライン説明です。以下に、このテストを計算するためのRコードを含めます。

# some test data
D <- data.frame(group1=c(3,6,2,2,3,9,3,4,2,5), group2=c(5,3,10,1,10,2,4,4,2,2), group3=c(8,0,1,5,10,7,3,4,8,1))

# sample with replacement
resample <- function(X) sample(X, replace=TRUE)

# return true if mu1 < mu2 < mu3
test     <- function(mu1, mu2, mu3) (mu1 < mu2) & (mu2 < mu3)

# resampling test that returns the probability of observing the relationship
mean(replicate(1000, test(mean(resample(D$group1)), mean(resample(D$group2)), mean(resample(D$group3)))))
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.