[R] Problem with pairs() in nlme
Michael Kubovy
kubovy at virginia.edu
Tue Nov 28 20:45:02 CET 2006
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?
vfr is not a grouped object because:
vfrg <- groupedData(estimate ~ slant | subject, data = vfr)
vfrg.lmL <- lmList(
estimate ~ (slant + respType + visField + hand)^2 | subject,
vfrg)
Error in function (classes, fdef, mtable) :
unable to find an inherited method for function "lmList", for
signature "formula", "nfnGroupedData"
The data are at
http://www.people.virginia.edu/~mk9y/LongFormat.csv
To get the working data.frame:
vf <- read.table("http://www.people.virginia.edu/~mk9y/
LongFormat.csv", header = TRUE, sep=',')
names(vf) <- c(
'slant', 'estimate', 'respType', 'visField', 'hand',
'order', 'sex', 'subject', 'handPref'
)
require(reshape)
vf <- sort_df(vf, vars = 'slant')
vf$subject <- factor(vf$subject)
summary(vf)
vfr <- subset(vf, handPref == 'R', select = -handPref)
vfr <- subset(vfr, slant %in% c(10, 20, 30))
_____________________________
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400 Charlottesville, VA 22904-4400
Parcels: Room 102 Gilmer Hall
McCormick Road Charlottesville, VA 22903
Office: B011 +1-434-982-4729
Lab: B019 +1-434-982-4751
Fax: +1-434-982-4766
WWW: http://www.people.virginia.edu/~mk9y/
More information about the R-help
mailing list