Я хочу выполнить скорректированный анализ по сегментам / точкам останова / изменения. Основным предиктором интереса является возраст. Я заинтересован в оценке контрольной точки по возрасту с учетом расы, пола и образования. Ниже мой код:
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))