この例では、 segmented.lm (またはその他の関数) が明らかなブレークポイントを見つけられるように助けていただければ幸いです。
data = list(x=c(50,60,70,80,90) , y= c(703.786,705.857,708.153,711.056,709.257))
plot(data, type='b')
require(segmented)
model.lm = segmented(lm(y~x,data = data),seg.Z = ~x, psi = NA)
次のエラーで返されます。
solve.default(crossprod(x1), crossprod(x1, y1)) のエラー: システムは計算上特異です: 逆条件数 = 1.51417e-20
K を変更すると:
model.lm = segmented(lm(y~x,data = data),seg.Z = ~x, psi = NA, control = seg.control(K=1))
別のエラーが発生します：
segmented.lm(lm(y ~ x, data = data), seg.Z = ~x, psi = NA, control = seg.control(K = 1)) のエラー: 間隔内の 1 つのデータのみ: ブレークポイント(s ) 境界または近すぎる