Displaying 4 results from an estimated 4 matches for "rws30".
2011 May 17
0
hierarchical gamma model in lme4
...oks like I could
> remove the values at the end of the distribution and get a decent fit. I'd
> still like to try a gamma model though, if that's possible. Is it possible
> in lme4 or another package I don't know about?
>
> ---This is the 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....
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 expression has 92 elements: only the first used
2: In plot:(transect:site) :
numerical expression has 92 ele...
2011 May 13
1
using glmer to fit a mixed-effects model with gamma-distributed response variable
...d-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 expression has 175 elements: only the first used
2: In plot:(transect:site) :
numerical expression has 175 elements: only the first...
2011 May 18
0
using hglm to fit a gamma GLMM with nested random effects?
...e?
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: In Ops.factor(site, transect) : / not meaningful for factors
2: In Ops.factor(site/transect, plot) : / not meaningful for...