কীভাবে পরীক্ষা করবেন


9

ধরুন আমার কাছে তিনটি স্বতন্ত্র গ্রুপ রয়েছে, যার সাথে আছে μ1, μ2, μ3 যথাক্রমে।

আমি কীভাবে পরীক্ষা করতে পারি μ1<μ2<μ3 বা ব্যবহার না n1, n2, n3 প্রতিটি গ্রুপ থেকে নমুনা?

আমি বিশদ গণনা নয়, কিছু সাধারণ পদ্ধতি জানতে চাই। কীভাবে আমার অনুমান সেট করা যায় তা আমি বুঝতে পারি নাH0 এবং H1


1
এটি অর্ডার সংক্রান্ত একটি পরিসংখ্যানগত অনুমিতিবিষয়ের উপর বই আছে ।
কেজেটিল বি হালওয়ারসেন

1
অর্ডার বিধিনিষেধের আওতায় (১৯ 197৩) বার্লো, বার্থলেমিউ, ব্রামনার এবং ব্রাঙ্ক স্ট্যাটিস্টিকাল ইনফারেন্সের পুরানো বইও রয়েছে (যদিও এরপরে কিছু উন্নয়ন হয়েছে); ননপ্যারমেট্রিক পরীক্ষাগুলি যতদূর আছে, জোনকিয়ার-টেরপস্ট্রার পরীক্ষা আছে (যেমন কনভারের দেখুন) এবং একটি ম্যাচ টেস্টের (নিভ এবং ওয়ারথিংটনের বইটি চেষ্টা করুন)। আপনি সাধারণত একটি সাম্য নাল এবং একটি বিকল্প বিকল্প লিখবেন।
গ্লেন_বি -রিনস্টেট মনিকা


এখানে একজনের উচিত বলা উচিত নয়, যার রয়েছে ni দল থেকে নমুনা i, তবে সেটির আকারের একটি নমুনা রয়েছে ni গ্রুপ থেকে i.
মাইকেল হার্ডি

উত্তর:


8

পরিসংখ্যানগুলিতে আপনি "এক্স সত্য কিনা" তা পরীক্ষা করতে পারবেন না। নাল অনুমানটি মিথ্যা যে আপনি কেবল তার প্রমাণ সন্ধানের চেষ্টা করতে পারেন।

আপনার নাল হাইপোথিসিসটি বলা যাক

H01:μ1<μ2<μ3.
আসুন ধরে নেওয়া যাক আপনার কাছে ভেক্টরটি অনুমান করার একটি উপায় রয়েছে μ=(μ1,μ2,μ3)। জিনিসগুলি সহজেই ধরে রাখুন যে আপনার অনুমানকারী রয়েছে
xN(μ,Σ),
কোথায় Σ হয় 3×3covariate ম্যাট্রিক্স। আমরা নাল অনুমান হিসাবে আবার লিখতে পারেন
Aμ<0,
কোথায়
A=[110011].
এটি দেখায় যে আপনার নাল হাইপোথিসিসটি ভেক্টরের উপর একটি অসমতা সীমাবদ্ধতা হিসাবে প্রকাশ করা যেতে পারে Aμ। একটি প্রাকৃতিক অনুমানকারীAμ দেওয়া হয়
AxN(Aμ,AΣA).
আপনি এখন দেওয়া সাধারণ ভেক্টরগুলিতে অসমতার সীমাবদ্ধতা পরীক্ষা করার জন্য ফ্রেমওয়ার্কটি ব্যবহার করতে পারেন:

কুডো, আকিও (1963)। "একতরফা পরীক্ষার একটি মাল্টিভিয়ারেট অ্যানালগ"। ইন: বায়োমেট্রিকা 50.3 / 4, পৃষ্ঠা 403–418।

যদি স্বাভাবিকতা অনুমানটি প্রায় আনুমানিক ("asympotically") ধারণ করে তবে এই পরীক্ষাটিও কাজ করবে। উদাহরণস্বরূপ, আপনি গ্রুপগুলি থেকে নমুনা মাধ্যম আঁকতে পারলে এটি কাজ করবে। আপনি যদি আকারের নমুনা আঁকেন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.
If you reject both times then you have found evidence that H0 is false and you reject H0. If you don't, then you don't reject H0. Since you are testing multiple times you have to adjust the nominal level of the subtest. You can use a Bonferroni correction or figure out an exact correction (since you know Σ).

Another way of constructing a test for H02 is to note that

H02:maxk=1,2(Aμ)k0.
This implies using maxAx as a test statistic. The test will have a non-standard distribution under the null, but the appropriate critical value should still be fairly easy to compute.


Fair enough, I edited my answer.
Andreas Dzemski

Good answer (+1). Just to improve it a bit more, may I recommend replacing x with μ^ so that the notation reflects the intention that this object is an estimator for μ.
Ben - Reinstate Monica

1

The answer provided by @andreas-dzemski is correct only if we know that the data is normally distributed.

If we do not know the distribution, I believe it would be better to run a nonparametric test. In this case, the simplest seems to run a permutation test. This is a book about the topic and this is a nice online explanation. Below I include R code to compute this test.

# 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.