cov1 <- c("NOx","ethnic","BMI")
> msm1<-msm(state ~ years, subject = f.eid, data = data8,qmatrix= matrixQ, death=4,
+ obstype= 1,gen.inits = TRUE,method ="BFGS",control=list(fnscale=500,maxit=40000),
+ covariates = as.formula(paste("~",paste0(cov1,collapse ="+"),sep="")))
警告信息:
In msm(state ~ years, subject = f.eid, data = data8, qmatrix = matrixQ, :
Optimisation has probably not converged to the maximum likelihood - Hessian is not positive definite.相似问题