[R] lmer, MCMCsamp and ranef samples?

2011-02-19 Thread Jukka Koskela
I really hope sombody could help me with the following, I'm having problems accessing the random effect samples following the example on MCMCsamp: (fm1 - lmer(Reaction ~ Days + (1|Subject) + (0+Days|Subject), sleepstudy)) set.seed(101); samp0 - mcmcsamp(fm1, n = 1000, saveb=TRUE)

[R] BMS package and leading minor problems

2011-02-13 Thread Jukka Koskela
Hello, Could somebody please tell me what the following is about? I try to run BMS function for quite a small dataset (400 rows and 50 columns) to get model suggestions, but I end up getting following error message: bms-bms(dat) Error in chol.default(XtX.start) : the leading minor of

Re: [R] Class coef.mer into a data.frame?

2011-01-02 Thread Jukka Koskela
Thank you! This was exactly what I was looking for. Jukka Lainaus David Winsemius dwinsem...@comcast.net: On Dec 31, 2010, at 2:49 AM, Jukka Koskela wrote: Hello, Could somebody please tell me what am I doing wrong in following? I try extract coefficients (using arm-package) from

[R] Class coef.mer into a data.frame?

2010-12-30 Thread Jukka Koskela
Hello, Could somebody please tell me what am I doing wrong in following? I try extract coefficients (using arm-package) from the lmer frunction, but I get the following warning: a-data.frame(coef(res)) Error in as.data.frame.default(x[[i]], optional = TRUE, stringsAsFactors =

[R] Plotting longitudinal data

2010-09-07 Thread Jukka Koskela
Hello, Hope that someone could help me plotting longitudinal data below: 7213330001 0.8300 13.05.091 1 330001 0.8700 09.02.05NULL 4797330001 0.7700 21.03.07NULL 2399330001 0.7800 12.04.06NULL 2400