Иногда я использую $ в формуле линейных моделей, которая работает, но дает странные ответы.Рассмотрим следующее
library (ISLR)
agelims<-range(Wage$age)
age.grid<-seq(from=agelims[1],to=agelims[2])
fitX<-lm(wage~poly(age,4),data=Wage)
predsX<-predict(fitX,newdata=list(age=age.grid),se=TRUE)
fit<-lm(Wage$wage~poly(Wage$age,4),data=Wage)
preds<-predict(fit,newdata=list(age=age.grid),se=TRUE)
, поэтому в качестве подгонки я называю заработную плату Wage$wage, которая кажется правильной, но излишней.R, очевидно, не согласен, потому что если вы запустите preds$fit и predsX$fit, то один даст вам тысячи предсказаний, а другой даст вам только 63.
Какого черта здесь происходит?