मैं अपने मिश्रित मॉडल विश्लेषण को चलाने के लिए साहित्य में खुदाई करने का प्रयास करने के बाद अभिभूत महसूस करता हूं कि एआईसी का उपयोग करने के साथ सबसे अच्छा मॉडल या मॉडल का चयन करने के बाद। मुझे नहीं लगता कि मेरा डेटा इतना जटिल है, लेकिन मैं इस बात की पुष्टि कर रहा हूं कि मैंने जो किया है वह सही है, और फिर आगे बढ़ने की सलाह दें। मैं अनिश्चित हूं कि मुझे lme या lmer का उपयोग करना चाहिए और फिर उन दोनों के साथ, यदि मुझे REML या ML का उपयोग करना चाहिए।
मेरे पास चयन का एक मूल्य है और मैं जानना चाहता हूं कि कौन से कोवरिएट्स उस मूल्य को सबसे अच्छा प्रभावित करते हैं और भविष्यवाणियों के लिए अनुमति देते हैं। यहाँ मेरे परीक्षण के लिए कुछ उदाहरण दिए गए डेटा और मेरे कोड हैं जिनके साथ मैं काम कर रहा हूं:
ID=as.character(rep(1:5,3))
season=c("s","w","w","s","s","s","s","w","w","w","s","w","s","w","w")
time=c("n","d","d","n","d","d","n","n","n","n","n","n","d","d","d")
repro=as.character(rep(1:3,5))
risk=runif(15, min=0, max=1.1)
comp1=rnorm(15, mean = 0, sd = 1)
mydata=data.frame(ID, season, time, repro, risk, comp1)
c1.mod1<-lmer(comp1~1+(1|ID),REML=T,data=mydata)
c1.mod2<-lmer(comp1~risk+(1|ID),REML=T,data=mydata)
c1.mod3<-lmer(comp1~season+(1|ID),REML=T,data=mydata)
c1.mod4<-lmer(comp1~repro+(1|ID),REML=T,data=mydata)
c1.mod5<-lmer(comp1~time+(1|ID),REML=T,data=mydata)
c1.mod6<-lmer(comp1~season+repro+time+(1|ID),REML=T,data=mydata)
c1.mod7<-lmer(comp1~risk+season+season*time+(1|ID),REML=T,data=mydata)मेरे पास ~ 19 मॉडल हैं जो इस डेटा को विभिन्न संयोजनों के साथ और एक 2 तरह से बातचीत की शर्तों तक खोजते हैं, लेकिन हमेशा एक यादृच्छिक प्रभाव के रूप में आईडी के साथ और मेरे आश्रित चर के रूप में comp1।
- Q1। कौन सा उपयोग करें? लम या लमर? फर्क पड़ता है क्या?
इन दोनों में, मेरे पास ML या REML का उपयोग करने का विकल्प है - और मुझे बहुत अलग-अलग उत्तर मिलते हैं - AIC द्वारा अनुसरण किए गए ML का उपयोग करते हुए मैं सभी समान AIC मानों के साथ 6 मॉडल समाप्त करता हूं और मॉडल संयोजन का कोई मतलब नहीं है, जबकि REML 2 सबसे अधिक संभावना मॉडल में परिणाम सबसे अच्छा है। हालाँकि, जब REML चल रहा है मैं किसी भी लंबे समय तक aova का उपयोग नहीं कर सकता।
- Q2। एवीओएए के साथ एमएमएल का उपयोग करने का मुख्य कारण क्या है? यह मेरे लिए स्पष्ट नहीं है।
मैं अभी भी stepAIC को चलाने में सक्षम नहीं हूं या मुझे उन 19 मॉडलों को संकीर्ण करने का एक और तरीका नहीं पता है।
- Q3। क्या इस बिंदु पर stepAIC का उपयोग करने का कोई तरीका है?