我想运行调整后的细分/断点/变更点分析。感兴趣的主要预测因子是年龄。我对估算种族,性别和教育程度时的断点很感兴趣。下面是我的代码:
o <- lm(ab42~race+sex+educ, data=merged)
o.seg <- segmented(o, seg.Z = age+race+sex+educ, psi=list(age=NA,race=NA,sex=NA,educ=NA))
当我运行上面的代码时,它为psi选项中列出的所有变量提供了断点。我只想按年龄细分。
如何在分段函数中指定?
o <- lm(ab42~race+sex+educ, data=merged)
o.seg <- segmented(o, seg.Z = age+race+sex+educ, psi=list(age=NA,race=NA,sex=NA,educ=NA))