[問題] 用R跑segmented regression
想要請問一下 在R 統計軟體裡面的segmented函數裡
裡面要給定psi的值 可是我一直想不通要怎麼給定 好像找不到相關資料
我有210筆資料 這些資料的Time都是一樣的 可是有些資料跑得出來 有些不行
那些不行的試了好幾個範圍
結果都出現下列錯誤(附上參考程式碼)
Time<-seq(0,48,by=3)
d<-data.frame(Time,geneK[,1],low.coef.geneK[,1])
out.lm<-lm(low.coef.geneK[,1]~Time,data=d)
segmented(out.lm,seg.Z=~Time,psi=list(Time=c(3,17),
control=seg.control(display=FALSE)))
錯誤在segmented.lm(out.lm, seg.Z = ~Time, psi = list(Time =c(3,17)):
Length or names of Z and psi do not match
想要請問各位 是哪裡出了問題
謝謝~
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 118.171.161.120