求问,如果想在gWQS中的logistic回归模型中矫正协变量应该怎么写代码,我按照正常回归模型中矫正协变量的方法写:OW~wqs+age+sex,但是总是报错说:
Can't subset columns that dosen't exist.但是明明变量age,sex已经包含在data中了。
代码如下:
datawqs<-read_xlsx("D:/SAS/Rwqs.xlsx")
a = names(datawqs)[8:14]
results = gwqs(OW ~ wqs, mix_name = a, data = datawqs, q = 10, validation = 0.6,
b = 1000, b1_pos = FALSE, b1_constr = TRUE, family = "binomial", seed = 2020)
summary(results)