search for: lmer4

I seem to come to this problem alot, and I can find my way out of it with a loop, but I wish, and wonder if there is a better way. Here's an example (lmer1-5 are a series of lmer objects): bs=data.frame(bic=BIC(lmer1,lmer2,lmer3,lmer4,lmer5)$BIC) rownames(bs)=c('lmer1','lmer2','lmer3','lmer4','lmer5') best=rownames(bs)[bs==min(bs)] > best [1] "lmer5" This tells me that lmer5 is the model with the lowest BIC. I want to start working with lmer5 as the best model, such as fi...

...ion when it comes to lmer. I've read about the mcmcsamp()-function, but I'm not entirely sure on how to use it or on how to intrepret the output. Any advice would be most appreciated. Kind regards, Andreas Nord -- View this message in context: http://www.nabble.com/lmer4-and-variable-selection-tp19146850p19146850.html Sent from the R help mailing list archive at Nabble.com.

...null] equation like this: L1 - Yij = b0 + e L2 - b0 = B00 + u0 The following command in R provides results that are identical to the HLM program. results <- lmer( Y ~ 1 |id , PanelData4) I can't seem to find any examples on-line nor in the help about how to write the lmer4 formula that contains two predictor variables at level 1 with fixed slopes. L1 - Yij = b0 + b1(x) + b2(z) + e L2 - b0 = B00 + u0 b1 = B10 b2 = B20 Can someone give me an example? Thanks Jeff

Dear r-helpers, I would very much appreciate help with the following problem: The following command (in a .Rnw file) latex(anova(e7.lmer3, e7.lmer4), file = 'e7lmer34.tex', rowname = c ('nonlinear', 'linear'), longtable = FALSE, dcolumn = T, booktabs = T, table.env = F) produces the following output after running Sweave: % latex.default(anova(e7.lmer1, e7.lmer2), file = "e7lmer12.tex", rowname = c(...

I'm trying to figure out how to carry out a Poisson regression fit to longitudinal data with a gamma distribution with unknown shape and scale parameters. I've tried the 'lmer4' package's glmer() function, which fits the Poisson regression using: library('lme4') fit5<- glmer(seizures ~ time + progabide + timeXprog + offset(lnPeriod) + (1|id), data=pdata, nAGQ=1, family=poisson) #note: can't use nAGQ>1, not yet implemented summary(fit5) He...

Dear R users I am trying to load the package lmer4 but it seems the Matrix library is failing to load and giving me the following error message : >library(Matrix) Error in importIntoEnv(impenv, impnames, ns, impvars) : object 'Logic' is not exported by 'namespace:methods' Error: package/namespace load failed for 'Ma...

I have recently started learning R. I converted a Stata datafile into an R image file, but I am unable to do anything with the data. For example, I am unable to calculate the mean of any variable or create a scatterplot of some of the data. I always get a message saying, "Error in plot(rgnpc, incmean) : object "rgnpc" not found" when I attempt to make a simple plot of the

...replacement has length zero In addition: Warning message: In matrix(rep(q, m), nrow = n, ncol = sum(m), byrow = TRUE) : data length exceeds size of matrix I get the same error message when I try with the full model or a null model! I don't understand what this means. The model works ok in lmer4 with a poisson distribution (althought it's very overdispersed), so it's not as though it wont work at all. Any help with this would be greatly appreciated cheers Rachel -- View this message in context: http://r.789695.n4.nabble.com/glmmADMB-tp4616701.html Sent from the R help mailing l...