Displaying 2 results from an estimated 2 matches for "meanvar".
Did you mean:
meanval
2012 Feb 21
0
mvabund package: errors using manyglm() and meanvar.plot()
...rep, depth and
lit.
I then created a formula:
>
mvform<-mvformula(faunadat~harv*rep*lit+harv*rep*depth+harv*lit*depth+rep*lit*depth)
When I try to create a mean-variance plot, I get the following error. I
have checked and all the factors are the same length as the dependent
variables.
> meanvar.plot(mvform)
Error in default.meanvar.plot.mvformula(x = x, n.vars = n.vars, var.subset
= var.subset, :
the dimensions of the variables do not match
When I fit a multivariate linear model, there is no problem. All the
following statements give the desired result.
> lm<-manylm(mvform)
>...
2013 Jun 07
0
error running mvabund package
...mvabund(cplflorabund)
##Fitting predictive models
transect=as.factor(transect$transect)
flormat=as.matrix(florabund)
abund.nb <-manyglm(flormat~transect*vegtype, family="nagative.binomial")
predict(abund.nb, type= "response")
##Checking Model Assumptions
plot(abund.nb)
meanvar.plot(flormat~transect, col=as.numeric (vegtype)) #direct plot of
mean-variace relationship
anova(abund.nb, p.uni="adjusted")
It's after running the last code line [>anova(abund.nb, p.uni="adjusted")],
when I get thiis error:
Error in XvarIn[nterms - i, varseq > i...