Я пытаюсь воспроизвести график пример кривых отношения рисков, в которых использовались данные NHANES и из:
«Низкое воздействие свинца и смертность среди взрослых в США: население- на основе когортного исследования "Брюс П. Ланфер, Стивен Раух, Пегги Ауингер, Райан В. Аллен, Ричард В. Хорнунг
Но каким-то образом использование pspline() дает мне следующее сообщение об ошибке:
Error in solve.default(g$var, coef(g)) : system is computationally singular: reciprocal condition number = 1.23319e-20
Когда Я использую HE_pb (уровни опережения) вместо pspline(HE_pb,df=10), это работает и дает мне приличный график. Из других сообщений в большинстве ответов говорится: «Это означает, что ваша матрица дизайна необратима и, следовательно, не может использоваться для разработки модели». Но почему эта проблема возникает, когда я использую одну независимую переменную pspline ()?
library('survival')
library(survey)
library('pspline')
bpdsn<-svydesign(id=~psu, strata=~kstrata, weights=~wt_hm, data=knh_final, nest=T)
svycox.HE_pb.sp <- svycoxph(Surv(person_y,total)~pspline(HE_pb,df=10)+
age+factor(sex)+factor(sm_status)+sm_pack_yr+
factor(ho_incm)+HE_BMI, design=bpdsn)