similar to: error in lmLists in lme4 package (bug?)

Hi all, I'm trying to figure out the example about mixed models in the Pinheiro and Bates book (Mixed-Effects Models in S and S-Plus, 2000, pp. 135-137). One thing I don't understand is: When I run the command fm1Orth.lm <- lm( distance ~ age, Orthodont ) followed by fm2Orth.lm <- update( fm1Orth.lm, formula = distance ~ Sex*age ) and then do summary(fm2Orth.lm)

Dear R-devel members, I think I've found a minor bug in plot.intervals.lmList. ( The guide https://www.r-project.org/bugs.html suggests to report it here, as I do not have a bugzilla account.) Here is a minimal reproducible example to demonstrate the problem: fm1 <- lmList(distance ~ age | Subject, Orthodont) plot(intervals(fm1),ylab="a") This results in: "Error in

Dear List, I am plotting lmList objects using plot(intervals()) in nlme package. I want to make changes to the y-axis labels. When I try to change cex of y-axis labels using the following: fm1 <- lmList(distance ~ age | Subject, Orthodont) plot(intervals(fm1), scales=list(y=list(cex = .7))) I receive: Error in bwplot(formula = group ~ intervals | what, data = structure(list( : formal

Hello R-Helpers, I was trying to use the lmList function to get the lmList graphic similar to Pinheiro and Bates (pg 33). I did not have a problem creating the graphic when I used the Orthodont data frame or 2 other data sets when there are no missing values. My data has missing values. Do I need to remove the missing values before the lmList function will work? for a small example: > a

I'd like to fit a GLM to each of a number of subsets of some data. The `family' argument to `lmList' (in lme4) has given me cause for optimism, but so far I've only been able to achieve linear model fits. For example > df <- data.frame(gp = gp.temp <- factor(rep(1:3, each = 100)), x = x.temp <- rnorm(300), y = rpois(300, exp((-1:1)[gp.temp] + x.temp))) Then a call to

I am using rw0641. In my continuing quest to understand repeated measures analysis, I again return to lme. I exported the Potthoff and Roy data Orthodont.dat from S-PLUS 4.5 to avoid capture errors and ran the examples in the R help. I imported the data.frame with data <- read.table("Orthodont.dat",header=T) attach(data) and created the objects Orthodont.fit1 <-

Really old subject?, so, all my apologizes for digging up but, since I also ran into this? maybe this hack can be useful to someone I propose monkey patching here: library(ff) my.as.ffdf.data.frame <- function (x, vmode = NULL, col_args = list(), ...) { rnam <- attr(x, "row.names") if (is.integer(rnam)) { if (all(rnam == seq_along(rnam))) rnam <- NULL else

Dea-R community. I'd like to draw your attention to an issue I have recently encountered while doing my current data analysis. I've got an unexpected (to me) result from the command: > augPred(lmList(my.object)), 'my.object' being a grouped data frame of class: > class(my.object) [1] "nfnGroupedData" "nfGroupedData" "groupedData"

On Fri, 30 Oct 2020, Pages, Herve wrote: > > On 10/29/20 23:08, Pages, Herve wrote: > ... >> >> I can think of 2 ways to move forward: >> >> 1. Keep I()'s current implementation but suppress the warning. We'll >> make the necessary adjustments to DataFrame() to repair columns supplied >> as I(<S4>) objects. Note that we would still be in

Hi Martin, On 10/26/20 04:52, Martin Maechler wrote: >> >> Hi there, >> Is that change in R-devel intentional? >> >> library(Matrix) >> m <- as(matrix(c(0, 1)), "sparseMatrix") >> >> isS4(m) >> # [1] TRUE >> >> x <- I(m) >> # Warning message: >> # In `class<-`(x,

Dear Martin, Have you try to create a new variable for age squared, say agesq? If you fit the model using this new variable you should get the coefficients. So your new model is something like height~age+agesq I hope this helps, Saghir > -----Original Message----- > From: MARTIN Ludovic [SMTP:martinl@mathinfo.ens.univ-reims.fr] > Sent: Wednesday, 07 May, 2003 2:02 PM > To:

On 10/29/20 23:08, Pages, Herve wrote: ... > > I can think of 2 ways to move forward: > > 1. Keep I()'s current implementation but suppress the warning. We'll > make the necessary adjustments to DataFrame() to repair columns supplied > as I(<S4>) objects. Note that we would still be in the situation where > I(<S4>) objects break validObject() but

> > Hi there, > Is that change in R-devel intentional? > > library(Matrix) > m <- as(matrix(c(0, 1)), "sparseMatrix") > > isS4(m) > # [1] TRUE > > x <- I(m) > # Warning message: > # In `class<-`(x, unique.default(c("AsIs", oldClass(x)))) : > # Setting class(x) to multiple strings ("AsIs",

Dear r-helpers, After successfully running require(nlme) vfr.lmL <- lmList( estimate ~ (slant + respType + visField + hand)^2 | subject, vfr ) pairs(vfr.lmL, id = 0.01, adj = -0.5) # Pinheiro & Bates (p. 141) produces the following error: Error in sprintf(gettext(fmt, domain = domain), ...) : object "form" not found Any guesses as to what I may have done wrong?

Hi there, Is that change in R-devel intentional? library(Matrix) m <- as(matrix(c(0, 1)), "sparseMatrix") isS4(m) # [1] TRUE x <- I(m) # Warning message: # In `class<-`(x, unique.default(c("AsIs", oldClass(x)))) : # Setting class(x) to multiple strings ("AsIs", "dgCMatrix", ...); result will no longer be an S4 object

Dear R developers, It would be great if you could implement the two minor code changes suggested below, which would help processing large objects in R. Jens Oehlschl?gel # Wish no. 1: let [.AsIs return the class AFTER subsetting, not the class of the original object # Wish no. 2: adjust write.csv and write.csv2 for multiple calls in chunked writing # Rationale no. 1: a couple of packages

Dear colleges, I would like to use the lmList function. I have installed the lme4 library but when I try to use the lmList function I get this error message. Error: could not find function "lmList" Here you can see what kind of messages I am getting when installing lme4 package. Does anyone know how to solve this problem? Thanks Rosario > install.packages() --- Please select a

Hello, in R-2.1.0 I'm trying to prodice trellis plots from an lmList object as described in the help for plot.lmList. I can generate the plots from the help, but on my own data plotting fails with an error message that I cannot interpret (please see below). Any hints are greatly appreciapted. kind regards, Arne > dim(d) [1] 575 4 > d[1:3,] Level_of_Expression SSPos1 SSPos19

Hello, I installed the newest version of R (3.0.2) as well as the newest version of nlme (3.1-113) in order to use summary.lmList and other nlme functions. Once loading the new library, lmList and summary.lm can be found, but a number of additional functions cannot be found via command. Versions installed and loaded are correct. Any suggestions on how to allow for these functions of the new

Dear all, I've got a problem in including square variables in lme and lmlist functions. I've tried to work on Oxboys data of Pinheiro and Bates'book, which consist of the heights of 26 boys, each mesured on nine different occasions : > Oxboys Grouped Data: height ~ age | Subject Subject age height Occasion 1 1 -1.0000 140.50 1 2 1 -0.7479 143.40