在分段包中使用davies.test需要幫助

Question

對不起，提出這樣一個愚蠢的問題，但我真的很困惑！ 我有一個斷點的分段線性迴歸，現在我想在包segmented要測試與davies.test斜坡顯著差異：

現在哪一個

lin.mod <- lm(y~x) 
segmented.mod <- segmented(lin.mod, seg.Z=~x) 

是戴維斯測試是否正確？

davies.test(lin.mod, seg.Z=~x) 

或

davies.test(segmented.mod, seg.Z=~x) 

Answer 1

davies.test(lin.mod, seg.Z=~x)是正確的作爲davies.test（）只接受LM（）或GLM（）對象作爲合適的輸入。該功能在內部計算斷點。