search for: lmlists

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

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 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

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

Hi all, Forgive me if this is an obvious one.... I want to make a plot of confidence intervals from an lmList object with a collection of simple linear models (lm(y~x)) using: plot(intervals(mylmList)) and sort the plot by increasing mean values for the intercept. Is there an easy way of doing this? I've tried the "order()" and "sort.list()" functions, but I suspect

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

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:

Hello R-Help, ----------------------------------------------------------------------------------------------------------------------------------------- Issues (there are 2): 1) Possible bug when using lmtest::encomptest() with a linear model created using nlme::lmList() 2) Possible modification to lmtest::encomptest() to fix confusing fail when models provided are, in fact, nested. I have

I'm having trouble understanding how functions in the subset argument for lmList search for the objects they need. This trivial example produces "Error in rownames(fakedf) : Object "fakedf" not found": library(nlme) fitbyID <- function() { fakedf <- data.frame(ID = gl(5, 10, 50), A = sample(1:100, 50), B = rnorm(50))

Howdee, *** I know that the lmList() function exists, yet I don't want to use it. *** Would anyone be kind enough to tell how I can apply the function lm() to each level of a given factor so to obtain the intercept and slope for each factor level within a matrix? For instance, suppose a dataframe containing 3 variables: id, x and y. I want to compute the function lm() for each level

Hello, I am trying to use lmLists in the lme4 package and copying over very standard code from the nlme package given in 'Mixed-Effects Models in S and S-Plus'. It appears to not accept an 'I(age-11)' in the formula, though it will accept the formula with out the subtraction of 11 from age. This seems like it wo...

Hi r-helpers! I need to save the output of summary() function that I?ve runned like this: z<- lmList(y~x1+x2| x3, na.action=na.omit,data1,subset=year==1999) w<-summary(z) The output (w) is something like this: Call: Model: y ~ x1 + x2 | x3 Data: data1 Coefficients: (Intercept) Estimate Std. Error t value Pr(>|t|) 1 0.081110514 1.141352e-01

Dear all, I have a question on handling of missing values in lmList. My data set have continuous predictor and response, x and y, and a grouping variable group.id. All these variables have NAs and the data set also has several other variables that also contains NAs. To create the lmList-object seems to work fine: y.list <- lmList(y ~ x | group.id, data=mydata, na.action=na.omit) However,

In R, we currently have the functions anovalist.lm and anova.glmlist S / S-plus has anova.lmlist anova.glmlist On the other hand, the [n]lme package (library) of Doug Bates and Jose Pinheiro has an "lmList" class and an anova.lmList(.) method for that. We are considering to use anovalist.lm and anovalist.glm instead of the S/S-plus names mentioned above. These functions

in lmList ,lmList.formula and lmList.groupedData functions in help pages there is an argument called "groups" that doesn't exist in the code of the functions. <<insert bug report here>> --please do not edit the information below-- Version: platform = i686-unknown-linux arch = i686 os = linux system = i686, linux status = Patched (unreleased snapshot) status.rev

Hi R-helpers. #I start with the reproducible example: firm<-c(rep(1,10),rep(2,10),rep(3,10),rep(4,10),rep(5,10)) year<-c(rep(1998:2007,5)) industry<-c(rep(1,20),rep(5,10),rep(7,10),rep(9,10)) X1<-rnorm(50) X2<-rnorm(50,mean=0.5,sd=0.1) Y<-rnorm(50,mean=0,sd=0.5) data<-data.frame(firm, industry,year,X1,X2,Y) data #I need to calculate for all the industries the following

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

Dear community, I'm trying to examine my grouped data following page 6 http://cran.r-project.org/doc/contrib/Fox-Companion/appendix-mixed-models.pdf. I'm trying lmList this way: model.list <- lmList(log(v.dep) ~ log(v2) +log(v4) + v3 | v5, subset = v6=0, data=data) and obtain this error message: In Ops.factor(v4, v4) : | not meaningful for factors My original model is