Я пытался выполнить взвешивание обратной вероятности с помощью следующего сценария:
temp1<-ipwpoint(exposure=bmi_ly,
family="gaussian", numerator=~1,
denominator=~ age+sex+pack_years+edu_cat+drinkstatus,data=devdata)
where
bmi_ly=continuous (range 13-78),
age=continuous (range 18-79),
sex=dichotomous (1,2),
pack_years= categorical (1,2,3) and
drinkstatus=categorical (1,2,3).
but it returns with an error msg
Er
ror in dnorm(tempdat$exposure, predict(mod1), sd(mod1$residuals))
Нечисловой аргумент математической функции
Я пытался выяснить проблему, но не смог.Ребята, вы можете мне помочь.Спасибо.