Displaying 1 result from an estimated 1 matches for "lnmass".
2009 Feb 24
2
lmer, estimation of p-values and mcmcsamp
...been trying to follow the online guidance for estimating
p-values for parameter estimates on these and other help threads:
https://stat.ethz.ch/pipermail/r-help/2006-May/094765.html
http://pidgin.ucsd.edu/pipermail/r-lang/2007-August/000057.html
but have not yet been successful. This is my code:
lnmass <- lmer(log.mass ~ treatment + (1|block), data=exp1)
summary(lnmass)
samp <- rnorm(n=10000)
mcmcpvalue <- function(samp)
{std <- backsolve(chol(var(samp)),
cbind(0,t(samp)) - colMeans(samp),
transpose = TRUE)
sqdist <- colSums(std*std)
sum(sqdist[-...