[R] stepAIC problem
A.J. Rossini
rossini at blindglobe.net
Sun Oct 5 14:55:50 CEST 2003
"Hiroto Miyoshi" <h_m_ at po.harenet.ne.jp> writes:
> Could anyone help me please?
>
> --small example
>> library(MASS)
>> x1<-runif(100)
>> x2<-runif(100)
>> x3<-runif(100)
>> x4<-runif(100)
>> x5<-runif(100)
>> y<-x1+x2+x3+runif(100)
>> t<-data.frame(y=y,x1=x1,x2=x2,x3=x3,x4=x4,x5=x5)
It is going to be very unhappy if it tries to call the t() transpose
function. Please rename this variable and try again.
>> x<-lm(y~x1+x2+x3+x4+x5,data=t)
>> stepAIC(x)
> Start: AIC= -247.61
> y ~ x1 + x2 + x3 + x4 + x5
>
> Df Sum of Sq RSS AIC
> - x5 1 3.747e-06 7.456 -249.608
> - x4 1 0.026 7.483 -249.254
> <none> 7.456 -247.609
> - x1 1 4.866 12.322 -199.375
> - x2 1 8.182 15.639 -175.543
> - x3 1 8.597 16.054 -172.922
> Error in as.data.frame.default(data) : can't coerce function into a
> data.frame
best,
-tony
--
rossini at u.washington.edu http://www.analytics.washington.edu/
Biomedical and Health Informatics University of Washington
Biostatistics, SCHARP/HVTN Fred Hutchinson Cancer Research Center
UW (Tu/Th/F): 206-616-7630 FAX=206-543-3461 | Voicemail is unreliable
FHCRC (M/W): 206-667-7025 FAX=206-667-4812 | use Email
CONFIDENTIALITY NOTICE: This e-mail message and any attachme...{{dropped}}
More information about the R-help
mailing list