Displaying 4 results from an estimated 4 matches for "mer_mcmcsamp".
2009 Feb 11
2
generalized mixed model + mcmcsamp
...using lmer
(library lme4), and the family = quasibinomial. I have tried to obtain a
MCMC sample, but on calling mcmcsamp(model1, 1000) I get the following
error which I don't understand at all:
Error in .local(object, n, verbose, ...) : Update not yet written
traceback() delivers:
4: .Call(mer_MCMCsamp, ans, object)
3: .local(object, n, verbose, ...)
2: mcmcsamp(model1, n = 1000, verbose = FALSE)
1: mcmcsamp(model1, n = 1000, verbose = FALSE)
which again doesn't particularly help me.
R is 2.8.1 under Windows, lme4 clean installed just today.
Before the model is fitted I just read in data,...
2006 Jan 28
1
yet another lmer question
...not for M2:
> mcmcsamp(M1)
(Intercept) log(sigma^2) log(grop.(In))
[1,] 9.099073 0.5711817 3.246981
attr(,"mcpar")
[1] 1 1 1
attr(,"class")
[1] "mcmc"
> mcmcsamp(M2)
Error: inconsistent degrees of freedom and dimension
Error in t(.Call("mer_MCMCsamp", object, saveb, n, trans, PACKAGE =
"Matrix")) :
unable to find the argument 'x' in selecting a method for
function 't'
--
Andrew Gelman
Professor, Department of Statistics
Professor, Department of Political Science
gelman at stat.columbia.edu
www.stat.c...
2007 Apr 27
1
Example of mcmcsamp() failing with lmer() output
...dapting to lmer:
<<2>>=
require(lme4)
gl.lmer <- lmer(rounds ~ method * physique + (1 | physique/team),
data = gunload)
mcmcsamp(gl.lmer)
@
lmer runs as expected, but mcmcsamp() gives:
Error: Leading minor of order 1 in downdated X'X is not positive
definite
Error in t(.Call(mer_MCMCsamp, object, saveb, n, trans, verbose,
deviance)) :
error in evaluating the argument 'x' in selecting a method for
function 't'
I tried:
<<3>>=
require(nlme)
gl.lme <- lme(rounds ~ method * physique, random = ~1 | physique/
team, data = gunload)
@
which gave all t...
2006 Dec 11
2
How to write a two-way interaction as a random effect in a lmer model?
Dear All,
I am working with linear mixed-effects models using the lme4 package in
R. I created a model with the lmer function including some main effects,
a two-way interaction and a random effect. Now I am searching how I
could incorporate an interaction between the random effect and one of
the fixed effects.
I tried to express the interaction in: