প্রেডিকশন ইন্টারভাল গণনা করা হচ্ছে

Question

আমার এখানে নিম্নলিখিত ডেটা রয়েছে । হাইড্রোকার্বন শতাংশের হার 1.0 হয় যখন আমি গড় বিশুদ্ধতার উপর 95% আত্মবিশ্বাসের ব্যবধান গণনা করার চেষ্টা করছি। আর এ, আমি নিম্নলিখিত লিখুন।

> predict(purity.lm, newdata=list(hydro=1.0), interval="confidence", level=.95)
   fit      lwr      upr
1 89.66431 87.51017 91.81845

যাইহোক, আমি কীভাবে এই ফলাফলটি নিজেই অর্জন করতে পারি? আমি নীচের সমীকরণটি ব্যবহার করার চেষ্টা করেছি।

$$s_{new}=\sqrt{s^2\left(1+\frac{1}{N}+\frac{(x_{new}-\bar x)^2}{\sum(x_i-\bar x)^2}\right)}$$

এবং আমি আর এ নিম্নলিখিত লিখুন।

> SSE_line = sum((purity - (77.863 + 11.801*hydro))^2)
> MSE = SSE_line/18
> t.quantiles <- qt(c(.025, .975), 18)
> prediction = B0 + B1*1
> SE_predict = sqrt(MSE)*sqrt(1+1/20+(mean(hydro)-1)^2/sum((hydro - mean(hydro))^2))
> prediction + SE_predict*t.quantiles
[1] 81.80716 97.52146

আমার ফলাফলগুলি আর এর পূর্বাভাস ফাংশন থেকে পৃথক। পূর্বাভাস অন্তরগুলি সম্পর্কে আমি কী ভুল বুঝাবুঝি করছি?

Comments

আপনি কীভাবে আপনার কোডটিতে এমএসই গণনা করছেন?

আমি পোস্টে গণনা যোগ করেছি।

— আদর্শস্টিক্জ

এমএমজে হিসাবে পরামর্শ দেওয়া হয়েছে আপনার চেষ্টা করা উচিতpredict(purity.lm, newdata=list(hydro=1.0), interval="prediction", level=.95)

— ভিনাক্স

Answer 1

আপনার predict.lmকোড লাগানো মানগুলির জন্য আত্মবিশ্বাসের ব্যবধানগুলি গণনা করছে। আপনার হাতের গণনা নতুন ডেটার জন্য পূর্বাভাস অন্তর গণনা করছে। আপনি যদি predict.lmহাতে গণনা থেকে পেয়েছেন একই ফলাফল পেতে চান তবে পরিবর্তন interval="confidence"করুন interval="prediction"

Answer 2

Dpel থেকে ভাল উত্তর। আমি যুক্ত করব যে আত্মবিশ্বাসের ব্যবধান এবং পূর্বাভাস ব্যবধানের মধ্যে পার্থক্যটি নীচের মত বর্ণনা করা যেতে পারে:

আস্থা ব্যবধান

$$s_{new}=\sqrt{s^2\left(\frac{1}{N}+\frac{(x_{new}-\bar x)^2}{\sum(x_i-\bar x)^2}\right)}$$

ভবিষ্যদ্বাণী ব্যবধান

$$s_{new}=\sqrt{s^2\left(1+\frac{1}{N}+\frac{(x_{new}-\bar x)^2}{\sum(x_i-\bar x)^2}\right)}$$

উত্স স্লাইড পৃষ্ঠা 5/17 এবং 11/17 দেখুন