Displaying 4 results from an estimated 4 matches for "markov1".
Did you mean:
markov
2009 Feb 24
2
lmer, estimation of p-values and mcmcsamp
...+ (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[-1] > sqdist[1]/nrow(samp) }
markov1 <- mcmcsamp(lnmass, 10000)
HPDinterval(markov1)
mcmcpvalue(as.matrix(markov1[,1]))
mcmcpvalue(as.matrix(markov1[,2]))
mcmcpvalue(as.matrix(markov1[,3]))
mcmcpvalue(as.matrix(markov1[,4]))
mcmcpvalue(as.matrix(markov1[,5]))
mcmcpvalue(as.matrix(markov1[,6]))
The first time I tried it HPDinterva...
2007 Feb 13
1
lme4/lmer: P-Values from mcmc samples or chi2-tests?
...1 | site)
Df AIC BIC logLik Chisq Chi Df Pr(>Chisq)
m2c 4 89.036 96.173 -40.518
m2 5 84.834 93.755 -37.417 6.2023 1 0.01276 *
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
The p-values from mcmc are:
##
markov1=mcmcsamp(m2,5000)
HPDinterval(markov1)
lower upper
(Intercept) -1.394287660 0.6023229
logpatch 0.031154910 0.1906861
loghab 0.002961281 0.2165285
landscape_diversity 0.245623183 1.6442544
log(site.(In)) -41.156007604...
2007 Feb 12
1
lmer and estimation of p-values: error with mcmcpvalue()
...57 -0.5877
logpatch 0.1091 0.0491 2.2228
loghab 0.0875 0.0732 1.1954
landscape_diversity 1.0234 0.4850 2.1099
Correlation of Fixed Effects:
(Intr) lgptch loghab
logpatch 0.091
loghab -0.637 -0.121
lndscp_dvrs -0.483 -0.098 -0.348
markov1=mcmcsamp(m1,5000)
HPDinterval(markov1)
mcmcpvalue(as.matrix(markov1)[,1])
Error in colMeans(samp) : 'x' must be an array of at least two dimensions
2007 Mar 12
0
Pvalues and lme
...in the estimations of
> the p-values?
>
> Below is the corresponding R code with some output so that you can see:
>
> ##
> fit<-lmer(End~Treatment+offset(log(Area)+(1|Site/Treatment), family=poisson)
>
Summary
Intercept
> The p-values from mcmc are:
>
> ##
> markov1=mcmcsamp(m2,5000)
>
> HPDinterval(markov1)
> lower upper
> (Intercept) -1.394287660 0.6023229
> logpatch 0.031154910 0.1906861
> loghab 0.002961281 0.2165285
> landscape_diversity 0.245623183 1.6442544
> log(site.(In)) -41.156007604 -1.6993996
> attr(,"Probability&qu...