monitoring of autocorrelated linear profiles using generalised estimation equation


I am Charles Okon a research student from Nigeria.

How can I compute the ARL comparisons under different shifts of intercept, slope and standard deviation of Hotelling T^2, MEWMA and MCUSUM at autocorrelated levels of 0.1, 0.5 and 0.9?

I will like to know the statistical package that can run this analysis and the corresponding code.

A sample of these computations is seen in the work of "A, Narvand and P, Soleimani" on "Phase 11 monitoring of autocorrelated linear profiles using linear mixed model" under the simulation studies section.

I will be eternally grateful if you can get me out of this quagmire; I have been stuck in it for too long and can't make progress without the knowledge.