|
setup diary | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Rの直線での回帰は、次のように簡単にできる。
x<-c(1,2,3,4,5) y<-c(1,4,8,11,18) plot(x,y) fit1<-lm(y~x) abline(fit1)
しかし、多項式fitは面倒である。以前は、次のようにしていた。
fit2<-lm(y~x+I(x^2))
当たり前なのだが、つぎのようにしてもできるらしい。
x2<-x^2 fit3<-lm(y~x+x2)
実際に使うときには、どちらが良いのだろう。メモリを食わないという意味では、前者の方が良いような気もするし、後者の方が計算コストが少ない気もする。