R ဖြင့် nested anova ကို မည်သို့လုပ်ဆောင်ရမည်နည်း။

nested ANOVA သည် ANOVA (“ ကွဲလွဲမှုကို ခွဲခြမ်းစိတ်ဖြာခြင်း” ) အမျိုးအစားတစ်ခုဖြစ်ပြီး အနည်းဆုံးအချက်တစ်ခုသည် အခြားအချက်တစ်ခုအတွင်းတွင် အသိုက်အမြုံ ရှိသည်။

ဥပမာအားဖြင့်၊ သုတေသီတစ်ဦးသည် မတူညီသောမြေသြဇာသုံးမျိုးမှ အပင်ကြီးထွားမှုအဆင့် မတူညီသည်ကို သိရှိလိုသည်ဆိုပါစို့။

၎င်းကို စမ်းသပ်ရန်အတွက် မတူညီသော ပညာရှင် သုံးဦးစီက အပင်လေးမျိုးတွင် မြေသြဇာ A ဖြန်းကြပြီး အခြားနည်းပညာရှင် သုံးဦးက အပင်လေးမျိုးတွင် ဓာတ် B ဖြန်းကြပြီး အခြားပညာရှင် သုံးဦးက အပင်လေးမျိုးတွင် ဓာတ် C ဖြန်းကြသည်။

ဤအခြေအနေတွင်၊ တုံ့ပြန်မှုပြောင်းလဲမှု မှာ အပင်ကြီးထွားမှုဖြစ်ပြီး အချက်နှစ်ချက်မှာ နည်းပညာရှင်နှင့် ဓာတ်မြေသြဇာဖြစ်သည်။ နည်းပညာရှင်သည် ဓာတ်မြေသြဇာတွင် ချုပ်နှောင်ထား ကြောင်း တွေ့ရှိရပါသည်။

အောက်ဖော်ပြပါ အဆင့်ဆင့် ဥပမာသည် R တွင် ဤ nested ANOVA ကို မည်သို့လုပ်ဆောင်ရမည်ကို ပြသထားသည်။

အဆင့် 1: ဒေတာကိုဖန်တီးပါ။

ဦးစွာ၊ ကျွန်ုပ်တို့၏ဒေတာကို R တွင်ထိန်းထားရန် ဒေတာဘောင်တစ်ခု ဖန်တီးကြပါစို့။

 #create data
df <- data. frame (growth=c(13, 16, 16, 12, 15, 16, 19, 16, 15, 15, 12, 15,
                          19, 19, 20, 22, 23, 18, 16, 18, 19, 20, 21, 21,
                          21, 23, 24, 22, 25, 20, 20, 22, 24, 22, 25, 26),
                 fertilizer=c(rep(c(' A ', ' B ', ' C '), each= 12 )),
                 tech=c(rep(1:9, each= 4 )))

#view first six rows of data
head(df)

  growth fertilizer tech
1 13 A 1
2 16 A 1
3 16 A 1
4 12 A 1
5 15 A 2
6 16 A 2

အဆင့် 2- Nested ANOVA ကို ချိန်ညှိပါ။

R တွင် nested ANOVA နှင့် အံဝင်ခွင်ကျဖြစ်စေရန် အောက်ပါ syntax ကို သုံးနိုင်သည်။

aov (အဖြေ ~ အချက် A / အချက် B)

ရွှေ-

• တုံ့ပြန်မှု- တုံ့ပြန်မှု ပြောင်းလဲနိုင်သော
• factor A: ပထမအချက်
• factor B: ပထမအချက်တွင် ဒုတိယအချက်

အောက်ဖော်ပြပါကုဒ်သည် ကျွန်ုပ်တို့၏ဒေတာအတွဲအတွက် nested ANOVA နှင့် မည်သို့ကိုက်ညီရမည်ကို ပြသသည်-

 #fit nested ANOVA
nest <- aov(df$growth ~ df$fertilizer / factor(df$tech))

#view summary of nested ANOVA
summary(nest)

                              Df Sum Sq Mean Sq F value Pr(>F)    
df$fertilizer 2 372.7 186.33 53.238 4.27e-10 ***
df$fertilizer:factor(df$tech) 6 31.8 5.31 1.516 0.211    
Residuals 27 94.5 3.50                     
---
Significant. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

အဆင့် 3- ရလဒ်ကို ဘာသာပြန်ပါ။

အချက်တစ်ခုစီသည် အပင်ကြီးထွားမှုအပေါ် ကိန်းဂဏန်းအရ သိသာထင်ရှားသော သက်ရောက်မှုရှိမရှိ ဆုံးဖြတ်ရန် p-value ကော်လံကို ကြည့်ရှုနိုင်ပါသည်။

အထက်ဖော်ပြပါဇယားမှ ဓာတ်မြေသြဇာသည် အပင်ကြီးထွားမှုအပေါ် ကိန်းဂဏန်းအရ သိသာထင်ရှားသော အကျိုးသက်ရောက်မှုရှိသည်ကို တွေ့နိုင်သည် (p-value < 0.05) ဖြစ်သော်လည်း နည်းပညာရှင် (p-value = 0.211) မရှိပါ။

ကျွန်ုပ်တို့သည် အပင်ကြီးထွားမှု တိုးလိုပါက ဓာတ်မြေသြဇာကို အသုံးပြုသည့် နည်းပညာရှင် တစ်ဦးချင်းထက် အသုံးပြုသည့် မြေသြဇာကို အာရုံစိုက်ရန် လိုအပ်ကြောင်း ၎င်းကဆိုသည်။

အဆင့် 4- ရလဒ်များကို မြင်ယောင်ကြည့်ပါ။

နောက်ဆုံးတွင်၊ ဓာတ်မြေဩဇာနှင့် နည်းပညာရှင်အားဖြင့် အပင်ကြီးထွားမှုဖြန့်ဝေမှုကို မြင်သာစေရန် ဘောက်ကွက်များကို အသုံးပြုနိုင်သည်။

 #load ggplot2 data visualization package
library (ggplot2)

#create boxplots to visualize plant growth
ggplot(df, aes (x=factor(tech), y=growth, fill=fertilizer)) +
  geom_boxplot()

ကွဲပြားသောမြေသြဇာသုံးမျိုးကြားတွင် ကြီးထွားမှုမှာ သိသာထင်ရှားစွာ ကွဲပြားမှုရှိသည်ကို ဂရပ်ဖ်တွင်ပြသထားသော်လည်း ဓာတ်မြေသြဇာအုပ်စုတစ်ခုစီရှိ နည်းပညာရှင်များကြားတွင် များစွာကွဲပြားမှုမရှိပါ။

၎င်းသည် nested ANOVA ၏ရလဒ်များနှင့် ကိုက်ညီပုံပေါ်ပြီး ဓာတ်မြေသြဇာများသည် အပင်ကြီးထွားမှုအပေါ် သိသာထင်ရှားသောအကျိုးသက်ရောက်မှုရှိသည်ကို အတည်ပြုပေးသော်လည်း နည်းပညာရှင်တစ်ဦးချင်းမှ မလုပ်ပါ။

ထပ်လောင်းအရင်းအမြစ်များ

R ဖြင့် တစ်လမ်းမောင်း ANOVA ကို မည်သို့လုပ်ဆောင်ရမည်နည်း
R ဖြင့် နှစ်လမ်းသွား ANOVA ကို မည်သို့လုပ်ဆောင်ရမည်နည်း။
R ဖြင့် ANOVA ထပ်ခါတလဲလဲတိုင်းတာနည်း