আরে আনোভা () ফাংশন ব্যবহার করে দুটি মডেলের তুলনা করা


32

এর জন্য ডকুমেন্টেশন থেকে anova():

যখন অবজেক্টের ক্রম দেওয়া হয়, তখন 'আনোভা' নির্দিষ্ট ক্রমে মডেলগুলির একে অপরের বিরুদ্ধে পরীক্ষা করে ...

একে অপরের বিরুদ্ধে মডেলগুলি পরীক্ষা করার অর্থ কী? এবং আদেশ কেন ব্যাপার?

জেনাবেল টিউটোরিয়াল থেকে এখানে একটি উদাহরণ দেওয়া হয়েছে :

    >  modelAdd = lm(qt~as.numeric(snp1))
    >  modelDom = lm(qt~I(as.numeric(snp1)>=2))
    >  modelRec = lm(qt~I(as.numeric(snp1)>=3))
     anova(modelAdd, modelGen, test="Chisq")
    Analysis of Variance Table

    Model 1: qt ~ as.numeric(snp1)
    Model 2: qt ~ snp1
      Res.Df  RSS Df Sum of Sq Pr(>Chi)
    1   2372 2320                      
    2   2371 2320  1    0.0489     0.82
     anova(modelDom, modelGen, test="Chisq")
    Analysis of Variance Table

    Model 1: qt ~ I(as.numeric(snp1) >= 2)
    Model 2: qt ~ snp1
      Res.Df  RSS Df Sum of Sq Pr(>Chi)
    1   2372 2322                      
    2   2371 2320  1      1.77     0.18
     anova(modelRec, modelGen, test="Chisq")
    Analysis of Variance Table

    Model 1: qt ~ I(as.numeric(snp1) >= 3)
    Model 2: qt ~ snp1
      Res.Df  RSS Df Sum of Sq Pr(>Chi)  
    1   2372 2324                        
    2   2371 2320  1      3.53    0.057 .
    ---
    Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1   1 

আমি এই আউটপুটটি কীভাবে ব্যাখ্যা করব?

উত্তর:


31

আপনি যখন ব্যবহার করবেন anova(lm.1,lm.2,test="Chisq")এটি তুলনা করার জন্য চি-বর্গ পরীক্ষা করে lm.1এবং lm.2(যেমন এটি স্কোয়ারের অবশিষ্টাংশের যোগফল হ্রাসের পরিসংখ্যানগতভাবে গুরুত্বপূর্ণ কিনা তা পরীক্ষা করে)। মনে রাখবেন যে এটি কেবলমাত্র lm.1এবং lm.2নেস্টেড মডেলগুলির দ্বারা উপলব্ধি করে।

উদাহরণস্বরূপ, আপনি যে 1 তম আনোভা ব্যবহার করেছেন তাতে পরীক্ষার পি-মান 0.82 0. এর অর্থ the এর স্তরে মডেলজেন থেকে লাগানো মডেল "মডেলএড" উল্লেখযোগ্যভাবে আলাদা নয় । তবে, তৃতীয় আনোভাতে পি-মানটি ব্যবহার করে, "মডেলআরসি" মডেলটি significantly এ উল্লেখযোগ্যভাবে আলাদা ফর্ম মডেল "মডেলজেন" ।α=0.05α=0.1

লিনিয়ার মডেল ফিটগুলির জন্যও এনওভা পরীক্ষা করে দেখুন ।


1
কিন্তু এটি কি দু'জনের মধ্যে একটির চেয়ে ভাল কিনা সে সম্পর্কে কিছু বোঝায়? ধন্যবাদ!
কিড

5
এটি আপনি "আরও ভাল" শব্দটি কীভাবে সংজ্ঞায়িত করেন তার উপর নির্ভর করে। আপনি যদি এটিকে এমন মডেল হিসাবে সংজ্ঞায়িত করেন যা স্কোয়ারের কম অবশিষ্টাংশ যোগ করে দেয় তবে উত্তরটি হ্যাঁ। এটি কারণ, এই পরীক্ষাটি স্কোয়ারের অবশিষ্টাংশের যোগফলকে হ্রাসের সাথে তুলনা করে।
স্টেট

3
অন্যদিকে, দুটি মডেল যদি উল্লেখযোগ্যভাবে আলাদা না হয়, তবে কেউ কি যুক্তি দিতে পারে যে সহজ মডেলটি "আরও ভাল"? আমি এখানে পার্সিমনি নিয়ে ভাবছি।
সিনিনহো

1
আনোভা (মোড 1, মোড 2, পরীক্ষা = "এলআরটি") এর ফলে কী পার্থক্য আসে?
এলেরিএল
আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.