আসল পরীক্ষার পরিসংখ্যানগুলির জন্য আপনার সূত্রটি একটি বিশেষ ক্ষেত্রে যার সাধারণ ক্ষেত্রে ফলাফলটি দেখান। সাধারণভাবে, আমাদের যাচাই করতে হবে যে পরিসংখ্যানগুলি F বন্টনের বৈশিষ্ট্যএফ অনুসারে স্বাধীন χ 2 এর অনুপাত হিসাবে লেখা যেতে পারেχ2 , স্বাধীনতার ডিগ্রি দ্বারা বিভক্ত ।
যাক এইচ0: আর'β= আর সঙ্গে ও পরিচিত, nonrandom এবং পূর্ণ কলাম র্যাঙ্ক হয়েছে । এই প্রতিনিধিত্ব করে কুই (অসদৃশ Ops স্বরলিপি) জন্য রৈখিক বিধিনিষেধ ট ধ্রুবক মেয়াদ সহ regressors। সুতরাং, @ user1627466 এর উদাহরণে, P - 1 সাথে সঙ্গতিপূর্ণ কুই = ট - 1 শুন্যতে সব ঢাল কোফিসিয়েন্টস স্থাপনের সীমাবদ্ধতা।আরRআর : কে × কিউকুইকুইটপি - 1কুই= কে - 1
পরিপ্রেক্ষিতে ভীa r ( β^OLS) =σ2( এক্স'এক্স)- 1 , আমরা
R′(β^ols−β)∼N(0,σ2R′(X′X)−1R),
যাতে (সঙ্গে B−1/2={R′(X′X)−1R}−1/2 একটি "ম্যাট্রিক্স বর্গমূল" হচ্ছেB−1={R′(X′X)−1R}−1 মাধ্যমে, যেমন, একটি Cholesky পচানি)
n:=B−1/2σR′(β^ols−β)∼N(0,Iq),
যেমন
Var(n)==B−1/2σR′Var(β^ols)RB−1/2σB−1/2σσ2BB−1/2σ=I
যেখানে দ্বিতীয় লাইন ওলএসইয়ের বৈকল্পিকতা ব্যবহার করে।
এই যেমন দেখানো উত্তর আপনি এর প্রতি সংযোগ আছে (এছাড়াও দেখুন এখানে ), স্বাধীন d:=(n−k)σ^2σ2∼χ2n−k,
যেখানে σ 2=Y'এমএক্সY/(এন-ট)চলিত পক্ষপাতিত্বহীন ত্রুটি ভ্যারিয়েন্স অনুমান, সাথে আছেনএমএক্স=আমি-এক্স(এক্স'এক্স)-1এক্স'হয়এক্স এ পুনরায়চাপ দেওয়া থেকে "অবশিষ্টাংশ প্রস্তুতকারক ম্যাট্রিক্স"।σ^2=y′MXy/(n−k)MX=I−X(X′X)−1X′X
সুতরাং, যেমন ′ n′n স্বাভাবিক হিসাবে একটি চতুর্ভুজ রূপ,
n′n∼χ2q/qd/(n−k)=(β^ols−β)′R{R′(X′X)−1R}−1R′(β^ols−β)/qσ^2∼Fq,n−k.
বিশেষ করে, অধীনH0:R′β=r, এই পরিসংখ্যাত করার হ্রাস
এফ= ( আর'β^OLS- আর )'। আর'( এক্স'এক্স)- 1আর }- 1( আর'β^OLS- র ) / কিউσ^2। চকুই, এন - কে।
চিত্রণ জন্য, বিশেষ ক্ষেত্রে বিবেচনা আর'= আমি , r=0 , q=2 , σ 2 = 1 এবং এক্স ' এক্স = আমি । তারপর,
এফ = β ' OLS β OLS / 2 = β 2 OLS , 1 + + β 2 OLS , 2σ^2=1X′X=IF=β^′olsβ^ols/2=β^2ols,1+β^2ols,22,
OLS ঔজ্জ্বল্যের প্রেক্ষাপটে এর স্কোয়ারড ইউক্লিডিয় দূরত্ব উৎপত্তি উপাদানের সংখ্যা দ্বারা প্রমিত থেকে অনুমান - যে হাইলাইট, যেহেতুβ2OLS,2মান লম্ব ছক করা হয় এবং অত: পরχ21,এফবন্টন "গড়পড়তা একটি হিসেবে দেখা যেতে পারেχ2বন্টন।β^2ols,2χ21Fχ2
আপনি যদি একটি সামান্য সিমুলেশন পছন্দ করেন (যা অবশ্যই প্রমাণ হিসাবে নয়!), নাল পরীক্ষা করা হয় যে k রেজিস্ট্রারগুলির কোনওটাই বিবেচনা করে না - যা তারা প্রকৃতপক্ষে না, যাতে আমরা নাল বন্টন অনুকরণ করি।
আমরা তাত্ত্বিক ঘনত্ব এবং মন্টে কার্লো পরীক্ষার পরিসংখ্যানের হিস্টোগ্রামের মধ্যে খুব ভাল চুক্তি দেখতে পাই।
library(lmtest)
n <- 100
reps <- 20000
sloperegs <- 5 # number of slope regressors, q or k-1 (minus the constant) in the above notation
critical.value <- qf(p = .95, df1 = sloperegs, df2 = n-sloperegs-1)
# for the null that none of the slope regrssors matter
Fstat <- rep(NA,reps)
for (i in 1:reps){
y <- rnorm(n)
X <- matrix(rnorm(n*sloperegs), ncol=sloperegs)
reg <- lm(y~X)
Fstat[i] <- waldtest(reg, test="F")$F[2]
}
mean(Fstat>critical.value) # very close to 0.05
hist(Fstat, breaks = 60, col="lightblue", freq = F, xlim=c(0,4))
x <- seq(0,6,by=.1)
lines(x, df(x, df1 = sloperegs, df2 = n-sloperegs-1), lwd=2, col="purple")
প্রশ্ন এবং উত্তরের পরীক্ষার পরিসংখ্যানগুলির সংস্করণগুলি প্রকৃতপক্ষে সমতুল্য তা দেখতে নোটটি বিধিনিষেধের সাথে মিলে যায় তা লক্ষ করুন R′=[0I] এবংr=0 ।
যাক X=[X1X2] ভাগ করা হবে যার অনুসারে গুণফলগুলি শূন্যের নীচে শূন্যের মধ্যে সীমাবদ্ধ রয়েছে (আপনার ক্ষেত্রে, ধ্রুবক ব্যতীত সমস্ত কিছুই রয়েছে, তবে অনুসরণযোগ্য ব্যয়টি সাধারণ। এছাড়াও, দিন β OLS = ( β ' OLS , 1 , β ' OLS , 2 ) 'β^ols=(β^′ols,1,β^′ols,2)′ হতে উপযুক্ত বিভক্ত OLS ঔজ্জ্বল্যের প্রেক্ষাপটে অনুমান।
R′β^ols=β^ols,2
আর ' β OLS = β OLS , 2
এবং
R′(X′X)−1R≡D~,
নীচে ডান দিকে ব্লক
( এক্সটিএক্স)- 1= ( এক্স'1এক্স1এক্স'2এক্স1এক্স'1এক্স2এক্স'2এক্স2)- 1≡ ( এ।)~সি~বি~ডি~)
এখন, ˜ ডি =(এক্স ′ 2 এক্স2-এক্স ′ 2 এক্স1(এক্স ′ 1 এক্স1)-1এক্স ′ 1 এক্স2)প্রাপ্ত
বিভাজনযুক্ত বিপরীতগুলির ফলাফলগুলি ব্যবহার করুন-1=(এক্স ′ 2 এম এক্স 1 এক্স2)ডি~= ( এক্স'2এক্স2- এক্স'2এক্স1( এক্স'1এক্স1)- 1এক্স'1এক্স2)- 1= ( এক্স'2এমএক্স1এক্স2)- 1
যেখানেএমএক্স1=আমি- এক্স1( এক্স'1এক্স1)- 1এক্স'1 ।
এফকুইএফn আপনি মি= β^'ol , 2( এক্স'2এমএক্স1এক্স2) β^ol , 2
β^ol , 2= ( এক্স'2এমএক্স1এক্স2)- 1এক্স'2এমএক্স1Y
এফn আপনি মি= y'এমএক্স1এক্স2( এক্স'2এমএক্স1এক্স2)- 1( এক্স'2এমএক্স1এক্স2) ( এক্স'2এমএক্স1এক্স2)- 1এক্স'2এমএক্স1Y= y'এমএক্স1এক্স2( এক্স'2এমএক্স1এক্স2)- 1এক্স'2এমএক্স1Y
ইউএসএসআর - আরএসএসআর
আরএসএসআর = y'এমএক্স1Y
Yএক্স1এইচ0টিএসএস= ∑আমি( y)আমি- y¯)2
ইউএসএসআরএমএক্স1Yচালুএমএক্স1এক্স2
ইউএসএসআর====Y'এম'এক্স1এমএমএক্স1এক্স2এমএক্স1YY'এম'এক্স1( আমি- পিএমএক্স1এক্স2) এমএক্স1YY'এমএক্স1Y- y'এমএক্স1এমএক্স1এক্স2( ( এমএক্স1এক্স2)'এমএক্স1এক্স2)- 1( এমএক্স1এক্স2)'এমএক্স1YY'এমএক্স1Y- y'এমএক্স1এক্স2( এক্স'2এমএক্স1এক্স2)- 1এক্স'2এমএক্স1Y
সুতরাং,
আরএসএসআর - ইউএসএসআর==Y'এমএক্স1Y- ( y)'এমএক্স1Y- y'এমএক্স1এক্স2( এক্স'2এমএক্স1এক্স2)- 1এক্স'2এমএক্স1Y)Y'এমএক্স1এক্স2( এক্স'2এমএক্স1এক্স2)- 1এক্স'2এমএক্স1Y