Dear Mike,
When I call summary (random2a, R=5000), z-values and standard errors of some estimates are significantly changed. The estimates and goodness-of-fit indices remain the same.
I didn't try (R=10000), as "R=5000" already took 7 hours to finish;)
However, if the number of replications dramatically influences the z-values and std.errors, are there any criteria for choosing one result for reporting?
Thanks!
Ryan