I have a marklist named "CACG.results" with 8 models that incorporate 2 covariates treated as factors and 1 continuous covariate. I want to produce 4 estimates of survival for the 4 possible combinations of the 2 binary covariates using the mean value for the continuous covariate based on the model average of the 8 models in the marklist.
However, when I run this script:
allModels85 <- covariate.predictions(CACG.results, data=data.frame(Relmed=mRelmed), indices=c(1,49,97,145), alpha = .075)
I get the error:
Error in model$results$beta.vcv[used.beta, used.beta] : incorrect number of dimensions
Interestingly, I do not get this error when estimating survival for only 1 model at a time as in the case below
MostComplex85 <- covariate.predictions(CACG.results$S.LayCamInit, data=data.frame(Relmed=mRelmed), indices=c(1,49,97,145), alpha = .075)
MostComplex85$estimates