कैसे टेस्ट करना है या नहीं


9

मान लीजिए कि मेरे पास तीन स्वतंत्र समूह हैं, मतलब के साथ μ1, μ2, μ3 क्रमशः।

मैं कैसे परीक्षण कर सकता हूं μ1<μ2<μ3 या उपयोग नहीं कर रहा है n1, n2, n3 प्रत्येक समूह से नमूने?

मैं कुछ सामान्य कार्यप्रणाली जानना चाहता हूं, विस्तृत गणना नहीं। मैं समझ नहीं पा रहा था कि अपनी परिकल्पना को कैसे स्थापित किया जाएH0 तथा H1


1
यह आदेश प्रतिबंधित सांख्यिकीय निष्कर्ष का मामला है । कर रहे हैं विषय पर किताबें
kjetil b halvorsen

1
ऑर्डर प्रतिबंध (1973) के तहत बार्लो, बार्थोलेमेव, ब्रेमर और ब्रंक सांख्यिकीय निष्कर्ष द्वारा पुरानी किताब भी है (हालांकि तब से कुछ विकास हुए हैं); जहां तक ​​अप्रमाणिक परीक्षणों की बात है, वहां जोन्केरे-टेरपस्ट्रा टेस्ट (जैसे कॉनओवर देखें) और एक मैच टेस्ट (नीव और वर्थिंगटन द्वारा पुस्तक की कोशिश करें)। आप आमतौर पर एक समानता शून्य और एक आदेशित विकल्प लिखेंगे।
Glen_b -Reinstate मोनिका


यहां किसी को कहना चाहिए, ऐसा नहीं है ni समूह से नमूने i, लेकिन उस आकार का एक नमूना है ni समूह से i.
माइकल हार्डी

जवाबों:


8

आंकड़ों में आप यह जांच नहीं कर सकते कि "X सही है या नहीं"। आप केवल सबूत खोजने की कोशिश कर सकते हैं कि एक अशक्त परिकल्पना झूठी है।

मान लीजिए कि आपकी अशक्त परिकल्पना है

H01:μ1<μ2<μ3.
आइए यह भी मान लें कि आपके पास वेक्टर का आकलन करने का एक तरीका है μ=(μ1,μ2,μ3)। चीजों को बस रखने के लिए मान लें कि आपके पास एक अनुमानक है
xN(μ,Σ),
कहाँ पे Σ है 3×3कोवरिएट मैट्रिक्स। हम शून्य परिकल्पना को फिर से लिख सकते हैं
Aμ<0,
कहाँ पे
A=[110011].
यह दर्शाता है कि आपकी अशक्त परिकल्पना को वेक्टर पर असमानता प्रतिबंध के रूप में व्यक्त किया जा सकता है Aμ। का एक प्राकृतिक अनुमानकAμ द्वारा दिया गया है
AxN(Aμ,AΣA).
अब आप दिए गए सामान्य वैक्टर पर असमानता की कसौटी के लिए रूपरेखा का उपयोग कर सकते हैं:

कुडो, एकियो (1963)। "एक तरफा परीक्षण का एक बहुभिन्नरूपी एनालॉग"। इन: बायोमेट्रिक 50.3 / 4, पीपी 403–418।

यह परीक्षण भी काम करेगा यदि सामान्य धारणा केवल (लगभग "asymptotically") रखती है। उदाहरण के लिए, यह तब काम करेगा जब आप समूहों से नमूना साधन बना सकते हैं। यदि आप आकार के नमूने आकर्षित करते हैंn1,n2,n3 और यदि आप समूहों से स्वतंत्र रूप से आकर्षित कर सकते हैं Σ विकर्ण के साथ एक विकर्ण मैट्रिक्स है

(σ12/n1,σ22/n2,σ32/n3),
कहाँ पे σk2 समूह में विचरण है k=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.
यदि आप दोनों बार अस्वीकार करते हैं तो आपको ऐसे प्रमाण मिले हैं H0 गलत है और आप अस्वीकार करते हैं H0। यदि आप नहीं करते हैं, तो आप अस्वीकार नहीं करते हैंH0। चूंकि आप कई बार परीक्षण कर रहे हैं, इसलिए आपको उप-स्तर के नाममात्र स्तर को समायोजित करना होगा। आप एक बोन्फ्रॉनी सुधार का उपयोग कर सकते हैं या एक सटीक सुधार का पता लगा सकते हैं (जब से आप जानते हैंΣ)।

के लिए एक परीक्षण के निर्माण का दूसरा तरीका H02 उस पर ध्यान देना है

H02:maxk=1,2(Aμ)k0.
इसका तात्पर्य परीक्षण आँकड़ा के रूप में का उपयोग करना है । परीक्षण में अशक्त के तहत एक गैर-मानक वितरण होगा, लेकिन उचित महत्वपूर्ण मान अभी भी गणना करने के लिए काफी आसान होना चाहिए।maxAx


पर्याप्त रूप से, मैंने अपना उत्तर संपादित किया।
एंड्रियास डेजीमेस्की

अच्छा उत्तर (+1)। बस इसे थोड़ा और बेहतर बनाने के लिए, क्या मैं को साथ बदलने की सलाह दे सकता हूं, ताकि संकेतन इस आशय को दर्शाता है कि यह वस्तु लिए एक अनुमानक है । xμ^μ
बेन -

1

@ Andreas-dzemski द्वारा प्रदान किया गया उत्तर केवल सही है अगर हम जानते हैं कि डेटा सामान्य रूप से वितरित किया गया है।

अगर हम वितरण को नहीं जानते हैं, तो मेरा मानना ​​है कि एक गैर-परीक्षण पद्धति को चलाना बेहतर होगा। इस मामले में, सबसे सरल एक क्रमपरिवर्तन परीक्षण चलाने के लिए लगता है। यह विषय के बारे में एक पुस्तक है और यह एक अच्छा ऑनलाइन स्पष्टीकरण है। नीचे मैं इस परीक्षण की गणना करने के लिए आर कोड शामिल करता हूं।

# 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)))))
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.