Search Archives

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-sig-mixed-models/attachments/20100701/0443bad1/attachment.pl>

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-sig-mixed-models/attachments/20100701/15de19c9/attachment.pl>

On Thu, Jul 28, 2011 at 3:37 PM, Ben Bolker <bbolker at gmail.com> wrote: > I would try drop1() [works with version ...-40 of lme4, I believe] for > more convenient output. Just spotted this and can confirm it does...

Chuck Cleland wrote: > > Allison, P.D. (2005). Fixed Effects Regression Methods for Longitudinal > Data Using SAS. Cary, NC: SAS Institute. > > covers some of this territory. As someone else pointed out, a > limitation is that including "Subject" as a factor precludes...

Hi, Is there any way to get fitted values out of lmer so that any NAs in the predictors/outcome pop out in the correct position? From what I can tell, the "fitted" function doesn't have access to missingness...

Hi, Is the following function built in somewhere? concat = function(v) { res = "" for (i in 1:length(v)) res = paste(res,v[i],sep="") res } e.g. > concat(c("12","3","45")) [1] "12345" Cheers, Andy -- Andy Fugard, Postgraduate Research...

Douglas Bates wrote: > What I was referring to is the na.action attribute of the model frame Aha, I didn't wade that deeply. > > At this point one can get the fitted values with the NA's in the > correct...

Hello all, There's a new mailing list for researchers in psychology who are learning and using R. New users of R are especially welcome. To join, venture to http://www.jiscmail.ac.uk/archives/psych-r.html and click...

Hi there, Am I correct to believe that tetrachoric correlation is a special case of polychoric correlation when there are only two levels to the ordered factor? Thus it should be okay to use hetcor from the polycor package to...

On 25 Feb 2007, at 06:51, Anup Nandialath wrote: > > Is there an R function to generate random integers? Thanks in > advance. The package Random does what you want, but requires a net connection. http://cran.r-project.org/src...

Adam, See Baayen, R.H., Davidson, D.J. and Bates, D.M. (2008). Mixed-effects modeling with crossed random effects for subjects and items. Journal of Memory and Language 59, 390--412. beginning p. 396. Cheers, Andy On Mon, Jun...

Hello, I was wondering how to get residuals out of binomial lmers as the "residuals" function isn't implemented. I've noticed a few other people ask the question too but get no response. (Or at least I haven't...

Greetings, Assuming you're using lmer, here's an example which does what you need: > (fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy)) Linear mixed model fit by REML Formula: Reaction ~ Days + (Days | Subject) Data: sleepstudy AIC BIC logLik deviance REMLdev...

Dear John, Yes, that's great - thanks! Andy John Fox wrote: > Dear Andy, > > Yes, the tetrachoric correlation is a special case of the polychoric > correlation when both factors are dichotomous. > > The 95-percent confidence interval that you suggest might be...

Sorry, that works modulo variable names... Second try. med = median(sleepstudy$Reaction) sleepstudy$bin = (sleepstudy$Reaction > med) + 0 mod = lmer(bin ~ Days + (1|Subject) + (0+Days|Subject), data = sleepstudy, family = binomial) ilog = function(x) { 1/(1 + exp(-x)) } boxplot(ilog...

Ken Beath wrote: > On 15/01/2008, at 9:42 AM, Zhong, Xiao wrote: > >> Thanks, Ken. >> >> I tried your advice on my model: >> >> model2.growth.mcas5 <- lmer(response ~ bs(monthElapsed) + skills + >> (1|studentID), data= mcas5, family=binomial(link="logit"), control...

On 28 Jun 2008, at 16:20, Peng Jiang wrote: > Hi, Andy > I am a little confused, why don't you just use paste() directly? > > > paste("12","3","45",sep="") > > produce the same result with your concat. The problem is...

I made (more) sense of mixed effects models when I went back to standard regression and ANOVAs with the philosophy: Everything is a Comparison. So for instance noting different ways of getting the magical numbers that result from doing regression...

Dear all, Are there any good reasons not to graph data points from fixed effects estimates + observation-level residuals in models with random effects at levels above residuals? Some code below to demonstrate what I mean. I can see arguments...

Peter Dixon wrote: > On Sep 11, 2008, at 1:15 PM, Douglas Bates wrote: > >> I should definitely add a check on p to the validate method. (In some >> ways I'm surprised that it got as far as mer_finalize...