Displaying 3 results from an estimated 3 matches for "resumry".
Did you mean:
resume
2005 May 09
1
bootstap and lme4
Hi,
I am trying to get bootstrap confidence intervals on variance
components and related statistics. To calculate the variance components
I use the package lme4.
> off.fun <- function(data, i){
d <- data[i,]
lme1<- lmer(y ~ trt + (trt-1|group), d)
VarCorr(lme1)@reSumry$group[2,1] #just as an example
}
> off.boot <- boot(data=data.sim, statistic=off.fun, R=100)
If I choose small values of R (<10) then I get very reasonable results
but for large R I get the following error massage:
*** malloc: vm_allocate(size=1069056) failed (error code=3)...
2005 Jun 24
1
lme4 extracting individual variance components
Hi,
For further calculations I need to extract indivdual Variances of
different random effects from a fitted model.
I found out how to extract the correlations
(VarCorr(m1)@reSumry$group1) but I was not able to find a way to
extract the other components individually.
To extract the Residuals I tried: (ranef(m1)@ stdErr) which
unfortunately did not work.
Thank you very much for your help!
Thomas
My model:
> m1<- lmer(y ~ trtt + (trtt-1|group3) + (trtt-1|group2) +...
2004 Jul 06
2
lme: extract variance estimate
For a Monte Carlo study I need to extract from an lme model
the estimated standard deviation of a random effect
and store it in a vector. If I do a print() or summary()
on the model, the number I need is displayed in the Console
[it's the 0.1590195 in the output below]
>print(fit)
>Linear mixed-effects model fit by maximum likelihood
> Data: datag2
> Log-likelihood: