I just found a most astounding blog via Gelman's blog:
http://emotion.inrialpes.fr/~dangauthier/blog/
This blog is a repository of cool things relating to statistical computing, simulation and stochastic modeling.
Search
Saturday, April 21, 2007
Tuesday, April 17, 2007
How to extract SEs from lmer fixed effects estimates
Extracting fixed effects coefficients from lmer is easy:
fixef(lmer.fit)
But extracting SEs of those coefficients is, well, trivial, but you have to know what to do. It's not obvious:
Vcov <- vcov(lmer.fit, useScale = FALSE)
se <- sqrt(diag(Vcov))
fixef(lmer.fit)
But extracting SEs of those coefficients is, well, trivial, but you have to know what to do. It's not obvious:
Vcov <- vcov(lmer.fit, useScale = FALSE)
se <- sqrt(diag(Vcov))
Saturday, February 17, 2007
Hmisc: how to increase magnification
One non-obvious thing (at least to me) about Hmisc's xYplot function is that to increase magnification or other parameters of a graph component, you have to do the following.
xlab=list("Condition",cex=2)
I.e., you have to make a list out of the parameter, and add whatever information you need. This works generally for any of the xYplot parameters.
xlab=list("Condition",cex=2)
I.e., you have to make a list out of the parameter, and add whatever information you need. This works generally for any of the xYplot parameters.
Subscribe to:
Posts (Atom)