Re: [R] stepAIC problem

2013-08-22 Thread Prof Brian Ripley
On 22/08/2013 09:24, Sachinthaka Abeywardana wrote: Hi all, I get the expected behaviour of getting a useful model if I do the following This is an 'expectation problem', not a 'stepAIC problem'. And stepAIC is not part of R and you are not crediting the tools you use. fit<-lm( expressions

[R] stepAIC problem

2013-08-22 Thread Sachinthaka Abeywardana
Hi all, I get the expected behaviour of getting a useful model if I do the following fit<-lm( expressions[,i]~expressions[,pa_all[1]]+expressions[,pa_all[2]]+expressions[,pa_all[3]]+expressions[,pa_all[4]]+expressions[,pa_all[5]]) step<-stepAIC(fit, direction="both") Output: Step: AIC=-78.75 ex