Hi all, I analysed my data with lme and after that I spent a lot of time for mean separation of treatments (post hoc). But still I couldn’t make through it. This is my data set and R scripts I tried. replication fertilizer variety plot height 1 level1 var1 1504 52 1 level1 var3 1506 59 1 level1 var4 1509 54 1 level1 var2 1510 48 2 level1 var1 2604 47 2 level1 var4 2606 51 2 level1 var3 2607 55 2 level1 var2 2609 44 3 level1 var2 3401 46 3 level1 var3 3402 64 3 level1 var4 3403 64 3 level1 var1 3404 50 1 level2 var3 1601 59 1 level2 var1 1605 56 1 level2 var2 1610 53 1 level2 var4 1611 53 2 level2 var2 2403 56 2 level2 var1 2405 61 2 level2 var4 2407 69 2 level2 var3 2413 70 3 level2 var3 3508 67 3 level2 var4 3511 73 3 level2 var2 3512 62 3 level2 var1 3513 67 1 level3 var4 1406 77 1 level3 var3 1408 74 1 level3 var1 1409 71 1 level3 var2 1410 69 2 level3 var4 2501 62 2 level3 var3 2507 58 2 level3 var2 2508 56 2 level3 var1 2513 63 3 level3 var3 3601 73 3 level3 var2 3603 59 3 level3 var1 3609 56 3 level3 var4 3612 61 modela<-lme(height~variety*fertilizer, random=~1|replication) summary(modela) anova(modela) library(multcomp) hgt <- glht(modela,linfct=mcp(fertilizer="Tukey")) summary(hgt) Any body can help me to proceed tukey (or lsd) with lme that would be highly appreciated. Prabhath University of Saskatchewan [[alternative HTML version deleted]]
-- View this message in context: http://r.789695.n4.nabble.com/lme-post-hoc-tp3224436p3224652.html Sent from the R help mailing list archive at Nabble.com.
plsc wrote:> > > I analysed my data with lme and after that I spent a lot of time for > mean separation of treatments (post hoc). But still I couldn???t make > through it. This is my data set and R scripts I tried. > > .... > > 3 level3 var4 3612 61 > > modela<-lme(height~variety*fertilizer, random=~1|replication) > library(multcomp) > hgt <- glht(modela,linfct=mcp(fertilizer="Tukey")) > Any body can help me to proceed tukey (or lsd) with lme that would be > highly appreciated. > >Try a search on tukey and lme, and make sure that you have a current version of R installed. The combination was broken on some earlier versions. http://r-project.markmail.org/search/?q=tukey+lme Dieter -- View this message in context: http://r.789695.n4.nabble.com/lme-post-hoc-tp3224436p3224655.html Sent from the R help mailing list archive at Nabble.com.
plsc wrote:> > > I analysed my data with lme and after that I spent a lot of time for > mean separation of treatments (post hoc). But still I couldn???t make > through it. This is my data set and R scripts I tried. > > .... > > 3 level3 var4 3612 61 > > modela<-lme(height~variety*fertilizer, random=~1|replication) > library(multcomp) > hgt <- glht(modela,linfct=mcp(fertilizer="Tukey")) > Any body can help me to proceed tukey (or lsd) with lme that would be > highly appreciated. > >Try a search on tukey and lme, and make sure that you have a current version of R installed. The combination was broken on some earlier versions. http://r-project.markmail.org/search/?q=tukey+lme Dieter -- View this message in context: http://r.789695.n4.nabble.com/lme-post-hoc-tp3224436p3224656.html Sent from the R help mailing list archive at Nabble.com.
Seemingly Similar Threads
- [LLVMdev] Closures, newbie question
- Improvement: SiteMapper - working ideas as a possible RoR''s routing replacement
- multiline string continuation
- [LLVMdev] "Recursive compilation detected" and signals
- Samba 3.0.x access rights issue with secondary groups or Unix rights