[R-sig-ME] LME and TukeyHSD
Christoph Scherber
Christoph.Scherber at agr.uni-goettingen.de
Fri Nov 19 10:58:34 CET 2010
Dear Sibylle,
You could also try to remove the NA´s like this:
mydata<-mydata[is.na(mydata$sinmortality)==F,]
or:
mydata<-na.exclude(mydata) # this may be a bit too strong, especially if many other variables contain NAs
Best wishes
Christoph
Sibylle Stöckli wrote:
> Dear Christoph
>
> Thanks a lot! I still have some trouble with the is.na function
>
> #
> mydata<-subset(Kaltenborn,ADDspecies!=1)
> mydata$sinmortality<-asin(sqrt(Kaltenborn$PropMortality))
> mydata<-mydata[!is.na(mydata$sinmortality),]
> error: is.na() applied to non-(list or vector)of type "NULL".
>
> I tried to use as.list
> #
> mydata<-mydata[!is.na(as.list(mydata$sinmortality)),]
>
> but this function produced some error within the lme model (varialbe lenght of sinmortality).
>
> Do you have some idea, how to change the !is.na function for non-list data?
>
> Thanks Sibylle
>
>
>
> Dear Sibylle,
>
> The reason for the error message may be that you transform your response
> variable inside the model formula. This frequently causes errors when
> working with components of lme models. In addition, you should try "ML"
> instead of "REML" (when working with the fixed effects).
>
> You may want to try:
>
> ##
>
> mydata<-subset(Kaltenborn,ADDspecies!=1)
> mydata$sinmortality<-asin(sqrt(Kaltenborn$PropMortality))
> mydata<-mydata[!is.na(mydata$sinmortality),]
>
> modelF<-lme(sinmortality~Diversity+Management+Species
> +Height+Year+Diversity*Year, data=mydata, random=~1|Plot/SubPlot,
> weights=varPower(form=~Diversity), method="ML")
>
> summary(glht(modelF,linfct=mcp(Species="Tukey")))
>
> ##
>
> Be reminded that it may not be sensible to test for Species main effects
> when interactions are present; have you tested for interactions between
> Species and other terms in the model?
>
> Best wishes,
> Christoph
>
>
>
>
>
> Dear R users
>
> I used lme to fit a linear mixed model inlcuding weights=varPower()
> and subset. Additionally I wanted to use glht to calculate Tukey-
> Kramer multiple comparisions.
>
> I think that the error in glht() is probably caused by the subset and
> varPower function within the lme model? I would very much appreciate
> an input.
>
> Thanks
> Sibylle
>
> LME MODEL
> modelF<-lme(asin(sqrt(PropMortality))~Diversity+Management+Species
> +Height+Year+Diversity*Year, data=Kaltenborn, random=~1|Plot/SubPlot,
> na.action=na.omit, weights=varPower(form=~Diversity), subset=Kaltenborn
> $ADDspecies!=1, method="REML")
>
> anova(modelF, type="marginal")
> numDF denDF F-value p-value
> (Intercept) 1 162 7.12789 0.0084
> Diversity 1 14 12.89284 0.0030
> Management 2 30 5.52544 0.0091
> Species 3 162 41.81003 <.0001
> Height 1 162 2.59134 0.1094
> Year 1 162 7.07660 0.0086
> Diversity:Year 1 162 12.88095 0.0004
>
> TukeyHSD
> library(multcomp)
> Loading required package: mvtnorm
> Warning messages:
> 1: package 'multcomp' was built under R version 2.7.2
> 2: package 'mvtnorm' was built under R version 2.7.2
> summary(glht(modelF,linfct=mcp(Species="Tukey")))
>
> error
> Error in glht.matrix(model = list(modelStruct = list(reStruct =
> list(SubPlot = -0.305856275920955, :
> ‘ncol(linfct)’ is not equal to ‘length(coef(model))’
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>
>
>
>
More information about the R-sig-mixed-models
mailing list