Hello R-help I am having trouble getting gls to find the R objects that comprise a linear model when the data=named.object option(option!) is not specified. In the gls() help it states data is "an optional data frame containing the variables named in model, correlation, weights, and subset. By default the variables are taken from the environment from which gls is called".
An example: > temp <- data.frame(x=1:10,y=11:20+rnorm(10)) > temp x y 1 1 11.52458 2 2 10.77643 3 3 12.56845 4 4 14.48822 5 5 13.58116 6 6 16.26223 7 7 17.89619 8 8 19.40359 9 9 18.56699 10 10 21.05374 > gls(temp$y~temp$x) Error in eval(expr, envir, enclos) : object "y" not found > gls(y~x,data=temp) Generalized least squares fit by REML Model: y ~ x Data: temp Log-restricted-likelihood: -14.00387 Coefficients: (Intercept) x 9.366256 1.135619 Degrees of freedom: 10 total; 8 residual Residual standard error: 0.9156084 I'm trying hard to avoid having to specify the data option if at all possible. Paul ------------------------------- R version 2.8.0 (2008-10-20) (yes I know its not 2.8.1 but nothing in the 2.8.1 news seemed to be relevant, I also tried on a different PC with R 2.7.0 but same problem) i386-pc-mingw32 locale: LC_COLLATE=English_Australia.1252;LC_CTYPE=English_Australia.1252;LC_MONETARY=English_Australia.1252;LC_NUMERIC=C;LC_TIME=English_Australia.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] nlme_3.1-89 loaded via a namespace (and not attached): [1] grid_2.8.0 lattice_0.17-17 ----------------------------------------------- Information on package 'nlme' Description: Package: nlme Version: 3.1-89 Date: 2008-06-07 Priority: recommended Title: Linear and Nonlinear Mixed Effects Models Author: Jose Pinheiro <jose.pinhe...@pharma.novartis.com>, Douglas Bates <ba...@stat.wisc.edu>, Saikat DebRoy <sai...@stat.wisc.edu>, Deepayan Sarkar <deepayan.sar...@r-project.org>, the R Core team. Maintainer: R-core <r-c...@r-project.org> Description: Fit and compare Gaussian linear and nonlinear mixed-effects models. Paul Rustomji Rivers and Estuaries CSIRO Land and Water GPO Box 1666 Canberra ACT 2601 ph +61 2 6246 5810 mobile 0406 375 739 ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.