Can the monte carlo utilities be used to estimate power for covariate effect in a DTSA model? I was trying to use the corresponding monte carlo code for example 6.20 as a guide to develop a model to estimate survival across 7 periods, but having some difficulty.
Specifically: 1) not clear exactly how to modify the population statement to set-up a given hazard model with low frequency event 2) the %sig for x on f (at least as I specified the model) jumps from 0.000 at one sample size (e.g., 225) to the next increment (e.g., 226).
Sorry, just realized that #2 (above) should have read: he %sig for x on f (at least as I specified the model) jumps from 0.000 at one sample size (e.g., 225) to 1.000 at the next increment (e.g., 226).
I discovered the error that was leading to my sample size finding (had another model statement later in the code that was incorrect). However, now have an issue where the covariate effect power jumps the same way (0.000 at .3 vs. 1.000 at .31) for a given n. Will send outputs/license as indicated.
1. By SS method, I assume you refer to the population covariance structure approach - that does not apply to categorical outcomes.
2. Chapter 12 doesn't cover MCMC (i.e. Bayesian analysis), but Monte Carlo simulations using various estimators (you can go to a UG example with growth modeling of categorical outcomes and then find the corresponding Monte Carlo setup on our website to see how that is done). If you want MCMC estimation you have to look for example with Estimator = Bayes.
Thank you for your quick reply and apologies for the confusion: I indeed meant Monte Carlo simulation, not Markov Chain Monte Carlo.
To clarify: I have fitted a rather basic LGC model with categorical outcomes and covariates using robust maximum likelihood, exported the estimates into EST1.DAT. and used these as population values in a Monte Carlo simulation to get estimates of the power of the regression coefficient for x10, following the broad logic of example 12.6.
I managed to make it successfully work with another LGC model, which had continuous outcomes, but I get these error messages when I declare the outcome variable to be categorical as you can see above.
Does that mean that such a 2 stage MC simulation study cannot run with categorical outcomes and covariates? If so, are there any alternatives for me to estimate the power of this regression coefficient?
It sounds like you want to do Step 2, not Step 1 of UG ex 12.6. In which case you should look at the input on page 430. Here you see Data Type=MonteCarlo, not MonteCarlo: as you have it. And Step 2 does not have the Model Population statements.
If this doesn't help, send to Support with license number.