Displaying 5 results from an estimated 5 matches for "shigo".
Did you mean:
saigo
2011 Apr 28
3
subset without removing NAs
Hi folks,
I'm dealing with a dataset with a lot of NAs, and want to use subset on the
data without removing the NAs from the the dataframe; e.g.
rws50 <- subset(rw.fire.RW,shigo.av<50)
This removes instances with NA for the value in addition to anything <50.
Any suggestions much appreciated.
*Ben Caldwell*
[[alternative HTML version deleted]]
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) :...
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 expression has...
2011 May 18
0
using hglm to fit a gamma GLMM with nested random effects?
...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: In Ops.factor(sit...
2011 May 17
0
hierarchical gamma model in lme4
...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 has been found: please supply starting
> values
> In addition: Warning message:
> In log(ifelse(y == 0, 1, y/mu)) : NaNs pro...