[R-sig-eco] pollination experiment with missing value
Roth Romain
romain.roth at cemagref.fr
Tue Dec 8 08:16:37 CET 2009
Dear Humberto,
Have you tried to replace the missing values?
Some process does exist like multiple imputation or NIPALS algorythm.
You can see in Legendre and Legendre, Numerical Ecology p47 for further
informations.
I don't know if those solutions could solve your problem, as they are
usually used in multivariate analysis...
Hope this helps,
Romain
Humberto Dutra a écrit :
> Dear subscribers
> I performed a quite simple pollination experiment:
> 12 plants
> 12 pollination treatments per plant,
> plant 1 had treatments a, b, c,
> d... etc plant 2 had treatments a, b, c, d...
> The idea was to run a simple anova using plants as blocks to see the effects of treatments on fruit production, but I lost some a few treatments in some plants, thus I have to deal with missing values
> First I tried the obvious, arc transformed fruits
>
>> fruset.aov<-aov(arcfruit~treatment+Error(Plant),
>>
> na.action=na.omit, fruset)
>
> but I get the wrong DF because of the missing values, and I cannot perform the multiple comparisons test. Is there a better way to deal with this unbalanced design?
>
> I did a little bit more research and I decided to run a GLM for mixed effects using a binomial distribution
>
>> y<-cbind(fruits,nofruits)
>> model<-glmer(y ~ treatment+ (1|Plant), binomial,data=fruset)
>>
>
> but them I cannot get the anova table
>
>> anova(model,test="F")
>>
> Error in anova(model, test = "F") :
> single argument anova for GLMMs not yet implemented
>
> Some people suggested using Anova function on the package car, but I don't see how can I get it to work with a mixed effects model like this.
>
> Any suggestions are appreciated. Are there other straight forward ways to analyze such data, given the missing values, and multiple comparisons follow up?
> Thank you
>
> Humberto
>
>
>
>
>
> ==========================================================
> 'Discipline - Success doesn't just happen. You have to be intentional about it, and that takes discipline.' - John Maxwell
>
>
>
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-ecology mailing list
> R-sig-ecology at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-ecology
>
More information about the R-sig-ecology
mailing list