2016年8月9日火曜日

R:区分回帰

区分回帰分析のまとめ(できたところまで)


データは
data.frame( )
などでまとめてから使用した.

test.lm <- lm(y~x, data)
test.seg <- segmented(test.lm, ~x)

プロットするなら,
plot(data)
par(new=T)
plot(test.seg)

近似直線は
test.seg$psi
で出たけど一個分だけのような.

まだまだウヌヌな状態.