মিশ্র-প্রভাবগুলির মডেলটির সাথে ইন্টারঅ্যাকশন টার্মে পোস্ট-হক তুলনা কীভাবে করবেন?

Question

পলির মাইক্রোবায়াল ক্রিয়াকলাপগুলিতে শুকানোর প্রভাবটি মূল্যায়নের জন্য আমি একটি ডেটা সেটে কাজ করছি। শুকানোর প্রভাব পলির ধরণের এবং / অথবা পলির মধ্যে গভীরতার মধ্যে পরিবর্তিত হয় কিনা তা নির্ধারণ করা হয়।

পরীক্ষামূলক নকশাটি নিম্নরূপ:

• প্রথম ফ্যাক্টর সেডিমেন্টটি তিন ধরনের পলির সাথে মিলিত হয় (কোডড শেড 1, সেড 2, সেড 3)। প্রতিটি ধরণের সেডিমেন্টের জন্য, তিনটি সাইটে (সেড 1 এর জন্য 3 সাইট, শেড 2 এর জন্য 3 সাইট, শেড 3 এর জন্য 3 সাইট) নমুনা তৈরি করা হয়েছিল।
• সাইট কোড করা হয়েছে: সাইট 1 , সাইট 2, ..., সাইট 9।
• পরবর্তী বিষয় হাইড্রোলজি : প্রতিটি সাইটের মধ্যে, একটি শুকনো প্লটে এবং একটি ভিজা প্লটে (শুকনো / ভিজা কোডযুক্ত) নমুনা দেওয়া হয়।

পূর্ববর্তী প্রতিটি প্লটের মধ্যে, ত্রিভুজের মধ্যে দুটি গভীরতায় (ডি 1, ডি 2) নমুনা দেওয়া হয় ।

মোট এন = 108 নমুনা = 3 টি পলি * 3 টি সাইট * 2 হাইড্রোলজি * 2 গভীরতা * 3 প্রতিলিপি রয়েছে।

আমি নীচে lme()আর ( এনএলএম প্যাকেজ) এ ফাংশনটি ব্যবহার করি :

Sediment <- as.factor(rep(c("Sed1","Sed2","Sed3"),each=36))
Site <- as.factor(rep(c("Site1","Site2","Site3","Site4","Site5",            
                        "Site6","Site7","Site8","Site9"),each=12))
Hydrology <- as.factor(rep(rep(c("Dry","Wet"),each=6),9))
Depth <- as.factor(rep(rep(c("D1","D2"),each=3),18))
Variable <- rnorm(108)

mydata <- data.frame(Sediment,Site,Hydrology,Depth,Variable)

mod1 <- lme(Variable ~ Sediment*Hydrology*Depth, data=mydata, 
             random=~1|Site/Hydrology/Depth)
anova(mod1)

আমি পদটি উল্লেখযোগ্য কিনা তা পরীক্ষা করার জন্য পোস্ট-হক তুলনা চালাতে চাই ।

আমি এটি একটি সাধারণ মূল প্রভাবের জন্য করতে সক্ষম হয়েছি (উদাঃ, পলি ):

summary(glht(mod1,linfct=mcp(Sediment="Tukey")))

তবে glht()ফাংশন মিথস্ক্রিয়া শর্তাবলী জন্য কাজ করে না।

আমি দেখতে পেয়েছি যে নিম্নলিখিতগুলি 2-উপায় আনোভার জন্য কাজ করতে পারে:

mod1 <- lme(Variable~Sediment*Hydrology, data=mydata, 
            random=~1|Site/Hydrology)
mydata$SH <- interaction(mydata$Sediment, mydata$Hydrology)
mod2 <- lme(Variable ~ -1 + SH, data=mydata, random=~1|Site/Hydrology)
summary(glht(mod2, linfct=mcp(SH="Tukey")))

3-উপায় অ্যানোভা ক্ষেত্রে একই পন্থাটি ব্যবহার করা কি সম্ভব? এক্ষেত্রে ইন্টারঅ্যাকশন শর্তাদির সাথে পোস্ট-হক তুলনা করার পথে যে কোনও সহায়তা করা প্রশংসিত হবে।

Answer 1

আমি "lsmeans" প্যাকেজটি বেশ দরকারী বলে মনে করেছি বিশেষত যখন কুড়াল * z * ভি ইন্টারঅ্যাকশন থাকে। তবে, প্যাকেজটি কেবলমাত্র আর এর নতুন সংস্করণগুলির জন্য উপলব্ধ

http://cran.r-project.org/web/packages/lsmeans/vignettes/using-lsmeans.pdf

Answer 2

আপনি কি বলতে চান যে তিনটি কারণের জন্য আপনি সমস্ত জুটিযুক্ত তুলনা করতে চান?

mod1<-lme(Variable~Sediment*Hydrology*Depth, data=mydata, random=~1|Site/Hydrology/Depth)
mydata$SHD<-interaction(mydata$Sediment,mydata$Hydrology,mydata$Depth)
mod2<-lme(Variable~-1+SHD, data=mydata, random=~1|Site/Hydrology/Depth)
summary(glht(mod2,linfct=mcp(SHD="Tukey")))

Comments

আমি দ্বি-মুখী মিথস্ক্রিয়া শর্তগুলির জন্য যুগলভাবে তুলনা করতে চাই, উদাহরণস্বরূপ "পলি জলবিদ্যুৎ"। যদি আনোভা টেবিলটি ইঙ্গিত দেয় যে পুরো মডেলের জন্য সেডিমেন্ট হাইড্রোলজি ইন্টারেক্টিশনটি উল্লেখযোগ্য তবে mod1<-lme(Variable~Sediment*Hydrology*Depth, data=mydata, random=~1|Site/Hydrology/Depth)আমি মনে করি না যে নিম্নলিখিত পদ্ধতিতে (সিডিমেন্ট * হাইড্রোলজি মিথস্ক্রিয়াতে) নিম্নলিখিত পদ্ধতিটি করা সঠিক:mydata$SH<-interaction(mydata$Sediment,mydata$Hydrology) mod2<-lme(Variable~-1+SH, data=mydata, random=~1|Site/Hydrology/Depth) summary(glht(mod2,linfct=mcp(SH="Tukey")))

— জন স্মিথ