একটি "সাধারণ" পরিমাপ ত্রুটি মডেল ফিট করার জন্য পদ্ধতি


13

আমি এমন পদ্ধতিগুলির সন্ধান করছি যা "ওএলএস" পরিমাপ ত্রুটির মডেলটি অনুমান করতে ব্যবহার করা যেতে পারে।

x i = X i + e x , i Y i = α + β X i

yi=Yi+ey,i
xi=Xi+ex,i
Yi=α+βXi

কোথায় ত্রুটি স্বাধীন স্বাভাবিক অজানা সঙ্গে ভেরিয়ানস এবং σ 2 এক্স । "স্ট্যান্ডার্ড" ওএলএস এই ক্ষেত্রে কাজ করবে না।σy2σx2

উইকিপিডিয়ায় কিছু অপ্রীতিকর সমাধান রয়েছে - দুটি প্রদত্ত দুটি আপনাকে ধরে নিতে বাধ্য করে যে হয় "বৈকল্পিক অনুপাত" বা "নির্ভরযোগ্যতা অনুপাত"λ=σ 2 এক্সδ=σy2σx2 , পরিচিত হয় যেখানেσ 2 এক্স সত্য regressor ভ্যারিয়েন্স হয়এক্সআমি। আমি এতে সন্তুষ্ট নই, কারণ যে রূপটি জানেন না এমন ব্যক্তি কীভাবে তাদের অনুপাতটি জানতে পারে?λ=σX2σx2+σX2σX2Xi

যাইহোক, এই দুটি ছাড়াও কি অন্য কোনও সমাধান রয়েছে যা আমাকে প্যারামিটারগুলি সম্পর্কে কিছু "জানতে" প্রয়োজন হয় না?

কেবল ইন্টারসেপ্ট এবং opeালের জন্য সমাধানগুলি ভাল।


উইকিপিডিয়া নিবন্ধ নিজেই আপনাকে এই প্রশ্নের উত্তর সরবরাহ করে। যদি আপনি "সত্য" রেজিস্ট্রারের স্বাভাবিকতা অনুমান করেন, তবে ত্রুটিগুলি বিতরণ করার জন্য আপনার আরও শর্ত প্রয়োজন। সত্যিকারের রেজিস্ট্রার যদি গাউসিয়ান না হন তবে আপনার কিছুটা আশা আছে। দেখুন Reiersol (1950)
কার্ডিনাল

এছাড়াও, "কেবলমাত্র ইন্টারসেপ্ট এবং opeালু সমাধানের সমাধানগুলি" বলতে আপনার অর্থ কী? এগুলি আপনার মাত্র দুটি পরামিতি! অথবা আপনিও "সত্য" রেজিস্ট্রারকেও ব্যাকআপ করার চেষ্টা করবেন বলে আশা করছেন?
কার্ডিনাল

@ কার্ডিনাল - আমি বোঝাতে চেয়েছিলাম যে আমি বিশেষত দুটি স্কেল পরামিতি সম্পর্কে যত্ন নিই না, এবং আপনি যেমনটি বলেছিলেন, "সত্য" রেজিস্ট্রার Xi
সম্ভাব্যতা ব্লগ

আমি দেখি. এটা বোধগম্য.
কার্ডিনাল

উত্তর:


7

উভয় ভেরিয়েবলের ত্রুটিগুলি সহ লিনিয়ার রিগ্রেশন সম্পর্কিত একটি Overতিহাসিক ওভারভিউতে জে ডাব্লু গিলার্ড দ্বারা বর্ণিত বিভিন্ন সম্ভাবনার রয়েছে

(x¯,y¯)β^=sy/sxxyyα^=y¯β^x¯.

এই বিশেষ পদ্ধতির গুণাবলী হয়

  1. xyyx
  2. এটি স্কেল-ইনগ্রেন্টেন্ট তাই আপনাকে ইউনিটগুলির বিষয়ে চিন্তা করার দরকার নেই,
  3. এটি দুটি সাধারণ লিনিয়ার রিগ্রেশন লাইনের মধ্যে রয়েছে
  4. এটি পর্যবেক্ষণের সেন্ট্রয়েডে এবং একে অপরকে যেখান থেকে পেরিয়ে যায় সেগুলি এটি অতিক্রম করে
  5. এটি গণনা করা খুব সহজ।

xy

YXXY

X0 <- 1600:3600
Y0 <- 5*X0 + 700
X1 <- X0 + 400*rnorm(2001)
Y1 <- Y0 + 2000*rnorm(2001)
slopeOLSXY  <- lm(Y1 ~ X1)$coefficients[2]     #OLS slope of Y on X
slopeOLSYX  <- 1/lm(X1 ~ Y1)$coefficients[2]   #Inverse of OLS slope of X on Y
slopesimple <- sd(Y1)/sd(X1) *sign(cov(X1,Y1)) #Simple slope
c(slopeOLSXY, slopeOLSYX, slopesimple)         #Show the three slopes
plot(Y1~X1)
abline(mean(Y1) - slopeOLSXY  * mean(X1), slopeOLSXY,  col="red")
abline(mean(Y1) - slopeOLSYX  * mean(X1), slopeOLSYX,  col="blue")
abline(mean(Y1) - slopesimple * mean(X1), slopesimple, col="green")

β^

{yi}{xi}σs

β^xyyxβ^yxxyβ^xy=ρ^sy/sxβ^yx=ρ^sx/syρ^xyρ^

x=by+c1/by=x/bc/byxρ^sy/sxsy/ρ^sxsy/sxyx

YX
আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.