Блог

Не можете определить эквивалентную функцию для команды «lincom» в R? Следовательно, невозможно реально вычислить p-значения. Если вообще кто-нибудь может мне помочь, какой stat-тест использовать для вычисления p-значений?

 Variables -
disease = numeric =0,1
race= char= European, Indian, Asian
Age1=num= 0,1
sex =factor= M,F
 

Код Stata —

 glm disease i.raceage1 sex, link(identity) family(binomial) 
testparm X 
di "Omnibus test of age1race interaction p = " %8.6f r(p) 
lincom age1 di %6.4f r(p) " | | " 
lincom age1   _IraceXage1__i' 
di %6.4f r(p) " | " _continue 
lincom _IraceXage__i' 
di %6.4f r(p) " | "
 

Пока что в R дошли до этого, любая помощь будет оценена по достоинству —

R-

 # install.packages("remotes")
# remotes::install_github("bmckuw/UWbe536")
library(UWbe536)

model1 <- glm(formula = disease ~ race * age1   sex, family = binomial(link = "identity"), data = MORT)
summary(model1)
fit2 <- augment(model1, newdata = MORT, type.predict = "response")
fit2$.fitted <- fit2$.fitted*100
lincom(model1, lc = "age1 == 0")
lincom(model1, lc = "race age1 = 0")