Lngmyers == Lngmyers [EMAIL PROTECTED]
on Sat, 04 Mar 2006 10:45:55 +0800 writes:
Lngmyers Thanks to Martin Maechler for the helpful
Lngmyers information below.
Lngmyers I guess I didn't make myself clear enough. I had
Lngmyers already discovered all them slots in lmer
I would like to write a function that runs GLMM using lmer
on a user-input model containing interactions, but if the
model doesn't produce significant results for the interaction,
a reduced model will be run without the interaction.
Doing this seems to require getting the p-values out of an
lmer
lngmyers == lngmyers [EMAIL PROTECTED]
on Fri, 3 Mar 2006 21:54:52 +0800 (CST) writes:
lngmyers I would like to write a function that runs GLMM using lmer
lngmyers on a user-input model containing interactions, but if the
lngmyers model doesn't produce significant results for
Thanks to Martin Maechler for the helpful information below.
I guess I didn't make myself clear enough. I had already discovered all
them slots in lmer objects, but these are apparently mostly input things
or things created in the course of the analysis, not the output itself
-- that is, it
Renaud == Renaud Lancelot [EMAIL PROTECTED]
on Tue, 6 Dec 2005 08:09:35 +0100 writes:
Renaud For example:
vc - vcov(m1, useScale = FALSE)
b - fixef(m1)
se - sqrt(diag(vc))
z - b / sqrt(diag(vc))
P - 2 * (1 - pnorm(abs(z)))
cbind(b, se, z, P)
Dear R users,
I've been struggling with the following problem: I want to extract the Wald
p-value
from an lmer() fit, i.e., consider
library(lme4)
n - 120
x1 - runif(n, -4, 4)
x2 - sample(0:1, n, TRUE)
z - rnorm(n)
id - 1:n
N - sample(20:200, n, TRUE)
y - rbinom(n, N, plogis(0.1 + 0.2 * x1 -
For example:
m1
Generalized linear mixed model fit using AGQ
Formula: cbind(y, N - y) ~ x1 + x2 + (1 | id)
Family: binomial(logit link)
AIC BIClogLik deviance
1137.308 1151.246 -563.6541 1127.308
Random effects:
GroupsNameVarianceStd.Dev.
id