search for: lmlist

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 `glm...

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 version of nlme to be found and u...

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

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

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

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 they only w...

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

...---Original Message----- > From: MARTIN Ludovic [SMTP:martinl@mathinfo.ens.univ-reims.fr] > Sent: Wednesday, 07 May, 2003 2:02 PM > To: r-help@stat.math.ethz.ch > Subject: [R] [R ] Query : problems with the arithmetic operator "^" > with function "lme" and "lmList" > > 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 : > &...

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 had a look around the web and in R-devel, but not been able to find any references to these issues. -----------------------------------------------------------------------...

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),...

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

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

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 7.106531e-01 0.4...

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=mydat...

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

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

...;-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 #coefficients and store it (lmList is from nlme package) b1 <- lmList(Y~X1+X2| year, na.action=na.omit, data, subset=industry==1) b2<-summary(b1) b3<-as.data.frame(b2$coefficients) b4<-round(b3,3) b4 #Instead of changing the industry in the subset I want to do it at once, so I?ve thinking in doing a loop. Before I...

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: &qu...

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(for...

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 originalmodel <- lm(log(v.dep) ~ log(v2) +log(v4) + v3 + v5+ v6, data=data) (Th...