Re: [R] data frame question

Hi Andras, here is an other solution which also works if b contains missing values: a <-seq(0,10,by=1) b <-c(NA, 11:20) f <-16 # a[which.max(b[b If it's not homework, then I'm happy to provide more help: > > > a <-seq(0,10,by=1) > b <-c(10:20) > d <-data.frame(a=a,b=b) > f <-16 > > subset(d, b <

Re: [R] Repeated measures ANOVA for unbalanc

Hi, The new lme4 package has a bootMer function, so it is fairly easy to compute bootstrap statistics. The following packages could be what you are looking for: 1) lmerTest (see ?lmerTest:::lmer) 2) car (see Anova) 3) afex (see mixed) 4) LMERConvenienceFunctions (see pamer.fnc) HTH, Denes >

Re: [R] Bug in predict.lm?

The same problem appears on a 64-bit linux: platform x86_64-pc-linux-gnu arch x86_64 os linux-gnu system x86_64, linux-gnu status major 3 minor 0.2 year 2013 month 09 day25 svn rev63987 language R ve

Re: [R] lapply?

Hi, the output of lapply() is a list; see ?lapply and ?sapply. # if you know the length of your list in advance, # this definition is better: uu <- vector("list", 2) # list elements uu[[1]] <- c(1,2,3) uu[[2]] <- c(3,4,5) # some options to achieve what you want: matrix(unlist(uu), 2, 3, T) do

Re: [R] named lmer.models in do.call(anova,models)

Sorry, I missed to attach the sessionInfo, here it goes: > sessionInfo() R version 3.0.1 (2013-05-16) Platform: x86_64-pc-linux-gnu (64-bit) locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C LC_TIME=en_US.UTF-8LC_COLLATE=en_US.UTF-8 [5] LC_MONETARY=en_US.UTF-8LC_MES

[R] named lmer.models in do.call(anova,models)

Hi, For some reason do.call on anova fails if the models are named lmer objects. Consider the following example: library(lme4) models <- list( lmer(Reaction ~ Days + (1| Subject), sleepstudy), lmer(Reaction ~ Days + (Days | Subject), sleepstudy)) # # models is an unnamed list, do.call wor