similar to: Anova over a list of 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 works (although with warning): do.call(anova, models) # # after labeling the models, do.call gives an

I think that this is because `[<-` dispatches on S4 methods only if the first argument is S4. ?"[<-" says: "These operators are also implicit S4 generics, but as primitives, S4 methods will be dispatched only on S4 objects ?x?." Georgi Boshnakov -----Original Message----- Message: 19 Date: Tue, 8 Sep 2020 22:04:44 -0400 From: Ben Bolker <bbolker at

I am using numerical optimization to fit a 1 parameter model, in which the input parameter is bounded. I am currently using optimize(), however, the problem turns out to have local minima, and optimize does not always seem to find the global minimum. I could to write a wrapping function that tries multiple intervals or starting values, but I would prefer a package that has built-in methods to make

I am currently using the package 'adapt' for multivariate integration. However this package seems to be removed from CRAN (It is still referred to in the help file for integrate(stats) though). I assume it has been deprecated for a reason? Is there an alternative for multivariate numerical integration? ----- Jeroen Ooms * Dept. of Methodology and Statistics * Utrecht University Visit

I would like to print some tables and figures to a PDF device on a CentOS 5 vps. However, I cannot seem to get the latex function from Hmisc working. I followed the example, and got an error: sh: xdvi: command not found. I tried installing the 'tetex-xdvi' linux package, and now it returns: Error: Can't open display. I guess the reason for this is that the machine is a VPS terminal, so

I just grabbed the new EL5 binary's for my x64 VPS from CRAN. However, where updates usually go smoothly, I now get these errors: -bash-3.2# rpm -i R-2.9.1-1.el5.x86_64.rpm error: Failed dependencies: R-devel = 2.9.1-1.el5 is needed by R-2.9.1-1.el5.x86_64 -bash-3.2# rpm -i R-core-2.9.1-1.el5.x86_64.rpm error: Failed dependencies: perl(File::Copy::Recursive) is needed by

I would like to convert a dataframe to an array of lists, one for every record. A natural choide is apply as.list to the rows. However, as it seems, as.list() automatically converts all list elements to the same datatype. Eg: myData <- data.frame(a="foo",b=as.logical(rbinom(10,1,.5))); apply(myData,1,as.list); In this output, all boolean values have been converted to character

Hi all, I am having issues comparing models with lmer. As an example, when I run the code below the model summaries (AIC, BIC, loglik) differ between the summary() and anova() commands. Can anyone clear up what's wrong? Thank you! Darren Ward library(lme4) data(sleepstudy) fm1<-lmer(Reaction ~ Days + (1|Subject), sleepstudy) summary(fm1) fm2<-lmer(Reaction ~ Days +

We recently started noticing freezes that appear only on 32bit systems (both linux and windows) with a relatively recent versions of R 3.2.0, including the RC. It looks like the problem can be traced back to the use of R_ObjectTables (see R_ext/Callbacks.h) The problem is a bit difficult to reproduce because it does not appear on x64 and because the official R interface to this functionality, the

Full_Name: Dale Barr Version: 2.5.1 (patched) OS: Ubuntu linux x86_64 Submission from: (NULL) (138.23.70.108) When RMySQL is loaded in before lme4, the summary() function for lmer objects in the lme4 packages produces the following error: Error in printMer(object) : no slot of name "status" for this object of class "table" When RMySQL is loaded AFTER lme4, however, no such

On 9/9/20 8:48 AM, Hugh Parsonage wrote: > I am unable to set break or use gdb with any success when I use that version. > > On linux I would do R -d gdb but this gives "unknown option '-d' " > while gdb R.exe (in the same directory as the debug version) gives the > same output as before. > > I'm happy to help but I appreciate this list might not be the

On Tue, Sep 8, 2020 at 11:44 PM Jeroen Ooms <jeroenooms at gmail.com> wrote: > > On Tue, Sep 8, 2020 at 5:20 PM Tomas Kalibera <tomas.kalibera at gmail.com> wrote: > > > > On 9/8/20 4:48 PM, Hugh Parsonage wrote: > > > Unfortunately I only get > > > > > > [Thread 21752.0x4aa8 exited with code 3221225477] > > > [Thread 21752.0x4514

I am making a little GUI for lme4, and I was wondering if there is a function that automatically detects on which level every variable exists. Furtheremore I got kind of confused about what a random effects model actually calculates. I have some experience with commercial software packages for multilevel analysis, like HLM6, and I was surprised that lme4 does not require the user to specify the

I have setup a Solaris server to test packages before submitting to CRAN, in order to catch problems that might not reveal themselves on Fedora, Debian, OSX or Windows. The machine runs a Solaris 11.2 vm with Solaris Studio 12.3. I was able to compile current r-devel using the suggested environment variables from "R Installation and Administration" and: ./configure

For future reference: https://sourceforge.net/p/mingw-w64/mailman/message/35747206/ On Wed, Mar 22, 2017 at 2:12 PM, Jeroen Ooms <jeroenooms at gmail.com> wrote: > This looks like a bug in mingw-w64 CRT. The problem can be produced > with C++ without R: > > #include <iostream> > #include <cmath> > #include <complex> > > int main(){ >

On 9/9/20 9:30 AM, Hugh Parsonage wrote: > Thank you! > > I get > > Starting program: C:\R\R-devel-20200909\bin\x64\Rgui.exe > [New Thread 19940.0x638c] > [New Thread 19940.0x102c] > [New Thread 19940.0x329c] > [New Thread 19940.0x37dc] > warning: Invalid parameter passed to C runtime function. > > Program received signal SIGSEGV, Segmentation fault. >

On Tue, Feb 13, 2018 at 2:20 AM, Jeroen Ooms <jeroenooms at gmail.com> wrote: > Thanks for your question. Your logs show that 'gcc' is not found. As > of R 3.3 you need to set the path to the compiler using the BINPREF > variable. This is because we ship two separate versions of gcc, one > targeting win32 and one targeting win64. I am not sure what your > rtools

I'm trying to understand how to plot individual growth curve trajectories, with the overall mean trajectory superimposed (preferably in a slightly thicker line, maybe in black) over the individual trajectories. Using the sleepstudy data in lme4, here is the code I have so far: library(lme4) library(lattice) xyplot(Reaction ~ Days, data = sleepstudy, group = Subject, type = 'l')

I would like to get some idea of which R-packages are popular, and what R is used for in general. Are there any statistics available on which R packages are downloaded often, or is there something like a package-survey? Something similar to http://popcon.debian.org/ maybe? Any tips are welcome! ----- Jeroen Ooms * Dept. of Methodology and Statistics * Utrecht University Visit

When building R-devel with gcc 5.2.0 (mingw-w64 v4) on Windows, make check fails reg-tests-1b.R at the following check: x <- c(1:2, NA) sx <- sd(x) !is.nan(sx) Here 'sx' should be 'NA' but it is 'NaN'. It turns out this problem only appears when the function is byte compiled with optimization level 3: mysd <- function (x, na.rm = FALSE) sqrt(var(if