Search Archives

Dear All, I just tried to fit a mixed-effects logistic regression model with version 0.99999911-5 of lme4 (installed from github). The model includes a random effect for clusters and a random group/treatment effect. I received the...

Thanks to all who responded so quickly! Yes, I totally overlooked par(mpg). Exactly what I was looking for. -- Wolfgang Viechtbauer ?Department of Methodology and Statistics ?University of Maastricht, The Netherlands ?http://www.wvbauer.com/ > -----Original Message----- > From: John Kane...

Indeed, the original post leaves some room for interpretation. In any case, I hope the OP has enough information now to figure out what approach is best for his data. Best, Wolfgang > -----Original Message----- > From: Bert Gunter [mailto:gunter.berton...

Yes, making blup() and ranef() work with 'rma.mv' objects is on my to-do list. Best, Wolfgang > -----Original Message----- > From: Gustaf Granath [mailto:gustaf.granath at gmail.com] > Sent: Thursday, December 15, 2016 10:03 > To: r-sig-mixed...

Trying to make the fastest referal to FAQ 7.31 is becoming quite a sport around here. David beat us all to it on this round. Chances are there will be many more rounds to come. Best, Wolfgang ----Original Message...

Hi all, I'm trying to fit a marginal (longitudinal) model with an exponential serial correlation function to the Orange tree data set. However, R crashes frequently when using the gnls() function. With the following simple example, I was able...

Dear All, The metafor package currently fails CRAN checks on Sparc Solaris: http://cran.r-project.org/web/checks/check_results_metafor.html The problem is probably due to an unintended (= stupid) use of identical() in a couple tests. I...

This is not an odds ratio. You simply calculated 27/131 = .206, which is the odds of having the disease given <40%. You can calculate a confidence interval for this with: exp(log(27/131) - 1.96 * sqrt(1/27...

Hi Naike, As far as I am concerned, stuff like this is more than appropriate for this mailing list. Best, Wolfgang -----Original Message----- From: R-sig-meta-analysis [mailto:r-sig-meta-analysis-bounces at r-project.org] On Behalf...

Dear All, I have a question about formulas and model.matrix(). If one specifies a model via a formula, the corresponding design matrix can be obtained with the model.matrix() function. For example: x1 <- c(1,4,2,3,5...

Dear All, I just noticed this paper, which may be of interest to some of the readers of this list: Li, B., Lingsma, H. F., Steyerberg, E. W., & Lesaffre, E. (2011). Logistic random effects regression models: A comparison of statistical...

Dear All, I am working on a function that has an optional data argument, just like lm(). If the user sets the data argument equal to some dataframe, then the function should look inside the dataframe for the variables given...

Hello All, I just stumbled across this article and I figured this may be of interest to some of you: Austin, P. C. (2010). Estimating Multilevel Logistic Regression Models When the Number of Clusters is Low: A Comparison of Different...

Dear Roberto, It depends on what parts of the output you want to put in a table. If you just want the "Model Results" part, then this will work: xtable(coef(summary(res))) where 'res' is an object returned by...

Hello All, This is actually more of a question directly for the lme4 maintainers, but I think others may be interested as well, so I am sending this to the mailing list. I noticed that there are already versions 1...

Dear All, An article was just published in PeerJ that may be of interest to some readers of this mailing list: Harrison, X. A. (2015). A comparison of observation-level random effect and beta-binomial models for modelling overdispersion in...

Dear All, Suppose I want to write a method for the generic function confint(): > args(confint) function (object, parm, level = 0.95, ...) So, it looks like the second and third argument have been predefined in the generic function. Suppose one...

As far as I can tell, you can't. But what is that 0.07 SD supposed to be a standard deviation of anyway? And n=3 for each group or was there only a single group that was measured...

Hello, In help(clusterApply) of the 'parallel' package, it says: "clusterApply calls fun on the first node with arguments seq[[1]] and ..., on the second node with seq[[2]] and ..., and so on, recycling nodes as needed." But the argument...

Dear Yalem, Naike Wang posted a video tutorial on this: https://stat.ethz.ch/pipermail/r-sig-meta-analysis/2017-September/000181.html The link to the video is: https://youtu.be/2wbXTFvaRnM There are also some notes on the...