Andrea.Capodicasa@albacom.it
2003-Aug-27 10:39 UTC
[Rd] Problem in step() and stepAIC() when a name of a regressors has b (PR#3991)
Hi all, I've experienced this problem using step() and stepAIC() when a name of a regressors has blanks in between (R:R1.7.0, os: w2ksp4). Please look at the following code: "x" <- c(14.122739306734, 14.4831100207131, 14.5556459667089, 14.5777151911177, 14.5285815352327, 14.0217803203846, 14.0732571632964, 14.7801310180502, 14.7839362960477, 14.7862217992577) "reg" <- structure(list(feriale = c(0, 1, 1, 1, 1, 0, 0, 1, 1, 1), pasquetta = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0), giorno = structure(c(2, 3, 4, 5, 6, 7, 1, 2, 3, 4), class = "factor", .Label = c("1-DOM", "2-LUN", "3-MAR", "4-MER", "5-GIO", "6-VEN", "7-SAB"))), .Names c("feriale", "pasquetta", "giorno"), row.names = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10"), class = "data.frame") #this works lm(x~., data=reg) #this works step(lm(x~., data=reg)) #this works library(MASS) stepAIC(lm(x~., data=reg)) #when you change a name of a regressors inserting blanks step and stepAIC won't work names(reg)[1]="with blank" names(reg) #this still works lm(x~., data=reg) #this doesn't work step(lm(x~., data=reg)) #this is the error #"Error in parse(file, n, text, prompt) : parse error" #this doesn't work stepAIC(lm(x~., data=reg)) # the error is the same as step #"Error in parse(file, n, text, prompt) : parse error" These are my R.version() and Sys.info() relevant infos.> R.Version()$platform [1] "i386-pc-mingw32" $arch [1] "i386" $os [1] "mingw32" $system [1] "i386, mingw32" $status [1] "" $major [1] "1" $minor [1] "7.0" $year [1] "2003" $month [1] "04" $day [1] "16" $language [1] "R"> Sys.info()sysname release version "Windows" "NT 5.0" "(build 2195) Service Pack 4" machine "x86"
Prof Brian D Ripley
2003-Aug-27 10:50 UTC
(PR#3991) Re: [Rd] Problem in step() and stepAIC() when a name of a regressors has b
Terms in formulae which are not syntactic names are not supported in R at present, and nothing that I know of says they are. We are working towards that as a goal, but this is not a bug. On Wed, 27 Aug 2003 Andrea.Capodicasa@albacom.it wrote:> Hi all, > I've experienced this problem using step() and stepAIC() when a name of a > regressors has blanks in between (R:R1.7.0, os: w2ksp4). > > Please look at the following code: > > "x" <- > c(14.122739306734, 14.4831100207131, 14.5556459667089, > 14.5777151911177, > 14.5285815352327, 14.0217803203846, 14.0732571632964, > 14.7801310180502, > 14.7839362960477, 14.7862217992577) > > "reg" <- > structure(list(feriale = c(0, 1, 1, 1, 1, 0, 0, 1, 1, 1), > pasquetta = c(0, > 0, 0, 0, 0, 0, 0, 0, 0, 0), giorno = structure(c(2, 3, 4, 5, > > 6, 7, 1, 2, 3, 4), class = "factor", .Label = c("1-DOM", > "2-LUN", > "3-MAR", "4-MER", "5-GIO", "6-VEN", "7-SAB"))), .Names > c("feriale", > "pasquetta", "giorno"), row.names = c("1", "2", "3", "4", > "5", > "6", "7", "8", "9", "10"), class = "data.frame") > > #this works > lm(x~., data=reg) > > #this works > step(lm(x~., data=reg)) > > #this works > library(MASS) > stepAIC(lm(x~., data=reg)) > > #when you change a name of a regressors inserting blanks > step and stepAIC won't work > names(reg)[1]="with blank" > names(reg) > > #this still works > lm(x~., data=reg) > > #this doesn't work > step(lm(x~., data=reg)) > > #this is the error > #"Error in parse(file, n, text, prompt) : parse error" > > > #this doesn't work > stepAIC(lm(x~., data=reg)) > > # the error is the same as step > #"Error in parse(file, n, text, prompt) : parse error" > > > These are my R.version() and Sys.info() relevant infos. > > R.Version() > $platform > [1] "i386-pc-mingw32" > > $arch > [1] "i386" > > $os > [1] "mingw32" > > $system > [1] "i386, mingw32" > > $status > [1] "" > > $major > [1] "1" > > $minor > [1] "7.0" > > $year > [1] "2003" > > $month > [1] "04" > > $day > [1] "16" > > $language > [1] "R" > > > Sys.info() > sysname release > version > "Windows" "NT 5.0" "(build 2195) > Service Pack 4" > machine > "x86" > > ______________________________________________ > R-devel@stat.math.ethz.ch mailing list > https://www.stat.math.ethz.ch/mailman/listinfo/r-devel > >-- Brian D. Ripley, ripley@stats.ox.ac.uk Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UK Fax: +44 1865 272595