Displaying 2 results from an estimated 2 matches for "fitsi".
Did you mean:
fits
2007 Mar 12
2
Lmer Mcmc Summary and p values
...ects
Intra T2 T3
T2 -0.989
T3 -0.745 0.737
T4 -0.577 0.570 0.430
> The p-values from mcmc are:
>
mcmcpvalue<-function(samp)
{
std<-backsolve(chol(var(samp)),
cbind(0,t(samp))-colMeans(samp),
transpose=TRUE)
sqdist<-colSums(std*std)
sum(sqdist[-1]>sqdist[1]/nrow(samp)
}
fitSI<-mcmcsamp(fit,50000)
library(coda)
HPDinterval(fitSI)
lower upper
Intercept -4.0778905 -3.1366836
Treatment2 3.4455972 4.3196598
Treatment 3 0.399302 1.287747
Treatment 4 -1.7898933 -0.2980325
log(Treatment*Site.(in)) -22.2198233 -19.7342530
log(Site.(In)) -28.7857...
2007 Mar 13
1
lme4 and mcmcamp
...effects
Intra T2 T3
T2 -0.989
T3 -0.745 0.737
T4 -0.577 0.570 0.430
> The p-values from mcmc are:
>
mcmcpvalue<-function(samp)
{
std<-backsolve(chol(var(samp)),
cbind(0,t(samp))-colMeans(samp),
transpose=TRUE)
sqdist<-colSums(std*std)
sum(sqdist[-1]>sqdist[1]/nrow(samp)
}
fitSI<-mcmcsamp(fit,50000)
library(coda)
HPDinterval(fitSI)
lower upper
Intercept -4.0778905 -3.1366836
Treatment2 3.4455972 4.3196598
Treatment 3 0.399302 1.287747
Treatment 4 -1.7898933 -0.2980325
log(Treatment*Site.(in)) -22.2198233 -19.7342530
log(Site.(In)) -28.78576...