Displaying 8 results from an estimated 8 matches for "ergostool".
2015 Mar 02
5
Import data set from another package?
I've moved nlme from Depends to Imports in my coxme package. However, a few of the
examples for lmekin use one of the data sets from nlme. This is on purpose, to show how
the results are the same and how they differ.
If I use data(nlme::ergoStool) the data is not found, data(nlme:::ergoStool) does no
better.
If I add importFrom(nlme, "ergoStool") the error message is that ergoStool is not exported.
There likely is a simple way, but I currently don't see it.
Terry T.
2008 Jan 29
3
on trellis.par.set/get (reproducing figures from Pinheiro & Bates)
Dear R users,
I would like to exactly reproduce a figure like the 1.5 or 1.9 or 4.13
from the book
Mixed effects models in S and S-Plus.
Not for the sake of it, but because I have my own data I would like to
plot in that fashion
(no colors)
If I write
plot(ergoStool)
I can get a good informative plot with colors, but I would like to have
a B&W one instead.
I've played a little with trellis.par.set
managing to change some aspect of the appearance of the plot,
for example with
prove.theme<-
list(
dot.line = list(col = 1, lty = 3, lwd = 1)...
2006 Mar 23
2
lme plot
...ve a questions regarding mixed effects models:
I'm trying to plot e.g. fitted vs residuals for each level of the
random effects, and i'm getting the same error.
I guess this might be a problem of the graphic capabilities of R.
Is there any way to obtain those plots?
library(nlme)
attach(ergoStool)
names(ergoStool)
[1] "effort" "Type" "Subject"
m1<-lme(effort~Type,random=~1|Subject)
plot(m1,form=resid(.,type="p")~fitted(.)|Subject,abline=0)#resid and
fitted for each level of subject
Error in as.vector(x, "list") : cannot coerce to vec...
2008 Apr 16
2
Post hoc tests with lme
Using the "ergoStool" data cited in Mixed-Effects Models in S and
S-PLUS by Pinheiro and Bates as an example, we have
========
> library(nlme)
> fm <- lme(effort~Type-1, data=ergoStool, random=~1|Subject)
> summary(fm)
Linear mixed-effects model fit by REML
Data: ergoStool
AIC BIC log...
2001 Dec 09
1
plot.design()
Greetings-
I'm working through Pinheiro and Bates' _Mixed Effects Models in S and
S-Plus_ using R (1.3.1 for linux). On page 13 (okay, so I haven't got that
far :)) is:
plot.design( ergoStool)
which returns on my system:
> plot.design(ergoStool)
Error: couldn't find function "plot.design"
any ideas?
Thanks.
---------------------------------------------------------
Andrew J. Perrin - Assistant Professor of Sociology
University of North Carolina, Chapel Hi...
2008 Jan 10
1
general linear hypothesis glht() to work with lme()
Hi,
I am trying to test some contrasts, using glht() in
multcomp package on fixed effects in a linear mixed
model fitted with lme() in nlme package. The command I
used is:
## a simple randomized block design,
## type is fixed effect
## batch is random effect
## model with interaction
dat.lme<-lme(info.index~type, random=~1|batch/type,
data=dat)
glht(dat.lme, linfct = mcp(type
2004 Jan 21
0
intervals in lme() and ill-defined models
...the random effects term without appropriate replication in
the data, the model will appear to fit but the confidence intervals for the
variance components will be very large. They suggest using intervals() as a
check that the model is appropriately defined:
> test.1 <- lme(effort~Type, data=ergoStool, random=~1|Subject)
> test.2 <- lme(effort~Type, data=ergoStool, random=~1|Subject/Type)
> intervals(test.2)
Random Effects:
Within-group standard error:
lower est. upper
1.054760e-07 4.599834e-01 2.005999e+06
In fact, using anova() to compare these two models s...
2004 Dec 19
1
PBIB datataset
I'm looking at Pinheiro & Bates "Mixed-Effects Models in
S and S-PLUS" at the moment. Several datasets are used,
one of which is called "PBIB" (a partially balanced
incomplete block design).
All the other datasets can be found somewhere or other in R.
However, I cannot locate PBIB, and it does not seem to
be mentioned in the latest edition of the R Full Reference