这次用RStudio的BQt包跑数据的时候,突然提示有错误
qrjm1 <- qrjm.BQt(formFixed = TC ~ OBSTIME+SEX+age+TG+HDL+LDL+BMI+SBP+DBP,
formRandom = ~ OBSTIME,
formGroup = ~ ID,
formSurv = Surv(TIME,STATUS2) ~ SEX+age ,
survMod = "weibull",
param = "value",
timeVar= "OBSTIME",
data = data1,
tau = 0.05,
n.adapt=20000)
数据格式如下:
> str(data1)
'data.frame': 1581 obs. of 14 variables:
$ ID : Factor w/ 234 levels "1","2","3","4",..: 1 1 1 1 1 1 1 2 2 2 ...
$ OBSTIME: int 1 2 3 4 5 6 7 0 1 2 ...
$ TIME : num 7 7 7 7 7 7 7 8 8 8 ...
$ STATUS : int 0 0 0 0 0 0 0 0 0 0 ...
$ STATUS2: int 0 0 0 0 0 0 0 0 0 0 ...
$ SEX : Factor w/ 2 levels "1","2": 1 1 1 1 1 1 1 1 1 1 ...
$ age : int 51 51 51 51 51 51 51 62 62 62 ...
$ TC : num 6.12 5.09 5.91 5.84 5.27 ...
$ TG : num 4.18 3.55 3.62 4.88 3.71 ...
$ HDL : num 2.1 1.7 1.44 1.53 1.17 ...
$ LDL : num 5.04 4.24 4.65 3.77 3.52 ...
$ BMI : num 29.4 23.8 30 25.5 28.7 ...
$ SBP : num 120 110 136 135 124 ...
$ DBP : num 90 79.3 87 79 86
报错如下:
Error in [[<-.data.frame
(*tmp*
, timeVar, value = c(7, 8, 4, 8, 8, :
替换数据里有242行,但数据有234
尝试过增加空白行,但是会直接提示有缺失值