Displaying 20 results from an estimated 8000 matches similar to: "summary(lmList) issues"
2003 May 07
1
[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 :
> Oxboys
Grouped Data: height ~ age | Subject
Subject age height Occasion
1 1 -1.0000 140.50 1
2 1 -0.7479 143.40
2003 May 07
1
[R ] Query : problems with the arithmetic operator "^" wi th function "lme" and "lmList"
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:
2006 Sep 14
1
plotting all subgroups with augPred
All,
I have a question RE plotting the prediction lines of a random effects
model via augPred. I'll illustrate via the Pixel dataset within the
nlme package:
library(nlme)
attach(Pixel)
fm1Pixel = lme(pixel ~ day + I(day^2), data = Pixel, random = list(Dog =
~ 1))
plot(augPred(fm1Pixel)) ### 10 fitted lines since there are 10 dogs
fm2Pixel = update(fm1Pixel, . ~ . + Side)
2004 Oct 11
0
lmList - strange problem
I am using lmList, which takes 'Data' partitioned according to the
levels of a grouping factor (gf) and individual 'lm' fits are obtained
for each 'data'partition, using a model defined as in "lm". So my call
to lmList looks something like:
> mg.lis <- lmList(strmbf ~ age1c + gender1 + sysbp.clinic +
+ diabp.clinic + ldl1 + mets.total + bmi1c |
2010 Feb 25
1
error in lmLists in lme4 package (bug?)
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 would be a bug, since this is standard formula
syntax, unless
2013 Nov 20
1
nlme function summary.lmList cannot be found with new versions
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
2008 Oct 08
1
Observed responses in 'augPred' data frame - Wrong order ?
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"
2009 Mar 17
1
coefficient graph
Dear R list members,
I'd like to make a graph of coefficients of the intercept, variable 1, and
variable 2 (and possibly the interaction between variable 1 and variable
2). When I use the lmList function as attached below, it shows a nice
coefficient graph.
> PACRP.lis <- lmList(PAffect ~ CRPC + CRPT + CINT | ID, redinteract)
> coef(PACRP.lis)
> PACRPlis.coef <-
2009 Jun 16
2
save the output of summary(lmList(x)) into a dataframe
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
2003 Sep 08
0
lmList with NAs
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
2010 Oct 25
1
Error: could not find function "lmList"
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
2006 Aug 18
0
lmList and missing values
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,
1998 Oct 21
0
anovalist.lm or anova.lmlist?
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
2005 May 13
1
error in plot.lmList
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
2012 Jan 18
0
examine grouped data lmList
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
2006 Feb 08
2
lme syntax for P&B examples
Hi helpeRs,
I've been working through some examples in Pinhiero & Bates( 2000)
trying to understand how to translate to the new Lme4 syntax but without
much luck.
Below is what I think I should do, but either the answers don't come out
the same or I get errors.
In the Oxide problems I'm particularly interested in obtaining the
levels coeficients but this options no longer seems
2006 Jan 17
2
multiple GLMs with lmList in lme4
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
2009 Jun 30
0
bug in nlme package function predict.lmList (PR#13788)
Full_Name: Vaidotas Zemlys
Version: 2.9.0
OS: Ubuntu 8.10
Submission from: (NULL) (213.197.173.50)
Steps to reproduce the bug:
library(nlme)
data(Oxboys)
qm=lmList(height~age|Subject,data=Oxboys)
grid=with(Oxboys,expand.grid(age=seq(min(age),max(age),length=50),Subject=levels(Subject)))
res <- predict(qm,grid,se=TRUE)
Erreur dans if (pool) { : l'argument est de longueur nulle
res <-
2009 Apr 16
0
incorrect handling of NAs by na.action with lmList (package nlme) (PR#13658)
Greetings,
I just found out a bug in the function lmList of the package nlme with R
2.8.1 running under windows XP 32-bits. I have a data table with various
columns corresponding to continuous variables as well as treatment
variables taken on several years and several sites. Here is an example :
Id year treatment A treatment B variable1
variable2 variable3
1
2009 Apr 14
0
Fitted values and MSE of individual fits in lmList
Dear useRs,
I am working on a series of field experiments (159 in total) carried out in different years in several locations. The cultivars in each experiment are not always the same, in fact they change over
time. I would like to get the fitted values and MSE of the individual fits from the following lmList object, so I can use them to fit a mixed model using the fitted values and weight each