Displaying 5 results from an estimated 5 matches for "crwn".
Did you mean:
cran
2011 May 13
1
graphs of gamma, normal fit to a histogram are about half as large as they should be
Hello,
I'm trying to compare the fit of two distributions, normal and gamma, to a
histogram of my response variable.
rate<-mean(na.omit(rwb$post.f.crwn.length))/var(na.omit(rwb$post.f.crwn.length))
shape<-rate*mean(na.omit(rwb$post.f.crwn.length))
hist((rwb$post.f.crwn.length), main="rwb$post.f.crwn.length")
lines(seq(0.01,70,0.01),length(rwb$post.f.crwn.length)*dgamma(seq(0.01,70,0.01),shape,rate))
lines(seq(0,70,0.1),length(na.omit(...
2011 Apr 28
1
using lme4 with three nested random effects
Hi all,
I'm trying to fit models for data with three levels of nested random
effects: site/transect/plot. For example,
modelincrBS<-glmer(l.ru.ba.incr~shigo.av+pre.f.crwn.length+bark.thick.bh+Date+slope.pos.num+dens.T+dbh+leaf.area+can.pos.num+(1|site/transect/plot),
data=rws30.UL, family=gaussian, na.action=na.omit)
but I get the following error:
Error: length(f1) == length(f2) is not TRUE
In addition: Warning messages:
1: In plot:(transect:site) :
numerical ex...
2011 May 17
0
hierarchical gamma model in lme4
...code I'm running---
>
> rws30.BL$site <- factor(rws30.BL$site)
> rws30.BL$transect <- interaction(rws30.BL$site, rws30.BL$transect, drop =
> TRUE)
> rws30.BL$plot <- interaction(rws30.BL$site, rws30.BL$transect,
> rws30.BL$plot, drop = TRUE)
> hist(rws30.BL$post.f.crwn.length)
> rws30.BL$gpost.f.crwn.length
>
> library("nlme")
> burnedmodel1.3<-lme(post.f.crwn.length~lg.shigo.av+dbh+leaf.area+
> bark.thick.bh+ht.any+ht.alive,
> random=(~1|site/transect/plot),na.action=na.omit, data=rws30.BL)
> Error: no valid set of coefficients...
2011 May 18
0
using hglm to fit a gamma GLMM with nested random effects?
...ample of one of those in the
package documentation. If it can, can anyone tell me what these errors are
trying to tell me?
If no, I promise I'll let this rest and just take B. Boker's advice to go
with a nice safe modified log transform of the data.
Best
test.gamma<-hglm(fixed=post.f.crwn.length~lg.shigo.av+dbh+leaf.area+
bark.thick.bh+ht.any, random=~1|site/transect/plot, family=Gamma(link=log),
data=rws30.BL)
Error in `contrasts<-`(`*tmp*`, value = "contr.treatment") :
contrasts can be applied only to factors with 2 or more levels
In addition: Warning messages:
1: I...
2011 May 13
1
using glmer to fit a mixed-effects model with gamma-distributed response variable
Sub: using glmer to fit a mixed-effects model with gamma-distributed
response variable
Hello,
I'm currently trying to fit a mixed effects model , i.e.:
> burnedmodel1.2<-glmer(gpost.f.crwn.length~lg.shigo.av+dbh+leaf.area+
bark.thick.bh+ht.any+ht.alive+(1|site/transect/plot), family=gaussian,
na.action=na.omit, data=rws30.BL)
If I run this code, I get the error below:
Error: length(f1) == length(f2) is not TRUE
In addition: Warning messages:
1: In plot:(transect:site) :
numerical...