topic Re: PROC MIXED and non normal data in Statistical Procedures

PROC MIXED and non normal data

AnalytX — Wed, 07 Nov 2012 11:25:58 GMT

Hello everybody,

I have non normal data to analyse using a repeated measures model. I schedule to use PROC MIXED on rank transformed data. Is it a good methodology?

Thanking you in advance for your answers.

Re: PROC MIXED and non normal data

SteveDenham — Wed, 07 Nov 2012 13:44:57 GMT

I have done this. Generate the ranks separately by time point. But, if you are analyzing the ranks using PROC MIXED, beware of unbalanced data. If you have 20 observations at the first time point, and only 10 at some later time point, there is no way the mean ranks at each time could be equal. Look up Friedman's Test on Google as an example of repeated measures on ranks.

Remember, the assumption in PROC MIXED about normality applies to the normality of the residuals, not of the data itself. For example, if you plotted the data and saw a bi-modal distribution, you would assume that the distribution is non-normal. Drill a little deeper, and it might be that you have an unaccounted for covariate (say gender) that leads to two peaks. If you can identify the process that generates the errors/residuals, you could specify the distribution and use PROC GLIMMIX. This is especially the case for errors from various exponential family distributions.

Steve Denham

Re: PROC MIXED and non normal data

lvm — Wed, 07 Nov 2012 14:00:13 GMT

There are many issues to consider when analyzing factorials (including repeated measures) using ranks. By definition, the ranks will have unequal variances, and an unstructured covariance matrix. Full details, including SAS code, can be found in the book:

Brunner, E., Domhof, S., Langer, F. 2002. Nonparametric Analysis of Longitudinal Data in Factorial Experiments. Wiley Publ.

Without taking some precautions in the analysis, you will get grossly inflated type I errors.

Re: PROC MIXED and non normal data

SteveDenham — Wed, 07 Nov 2012 14:12:51 GMT

I definitely agree with this, and it is easy to see if you transform all observations to ranks. However, with equal observations at each time point, and computing ranks separately at each time point, at least some of the concerns regarding unequal variances is addressed (i.e., for 20 observations per time point, the variance of the numbers 1 to 20 is asymptotically constant under a null of no group effect). I'm not so sure about the covariances, though, and this gives me room to think about some of the analyses we have been doing. Generally, autoregressive models have yielded lower information criterion values than the unstructured models, given the by-timepoint ranking.

Anyway, I am willing to wager that the OP's original question could be better answered by using GLIMMIX with a proper selection of error distribution. I have never been a big fan of nonparametrics, going back to a FORTRAN program I wrote as an undergraduate to do Spearman correlation.

Steve Denham

Re: PROC MIXED and non normal data

AnalytX — Tue, 13 Nov 2012 15:36:51 GMT

Thanks a lot for your answers which give me helpful inormation.

Regards,

Re: PROC MIXED and non normal data

plf515 — Fri, 16 Nov 2012 11:44:17 GMT

I agree with Steve's idea to use GLIMMIX or possibly NLMIXED, after checking out the residuals from the models generated by MIXED. GLIMMIX has a lot of distributions and with NLMIXED - well, there's almost nothing you can't do. There was an interesting paper at the most recent NESUG on fitting W shaped distributions.

Re: PROC MIXED and non normal data

Edgar_Brunner — Mon, 01 Jun 2015 12:15:19 GMT

Please find some answers to your question in the following papers:

Akritas, M. G., Arnold, S. F. and Brunner, E. (1997). Nonparametric hypotheses and rank statistics for unbalanced factorial designs. Journal of the American Statistical Association 92, 258-265.

Brunner, E. and Puri, M.L. (2001). Nonparametric Methods in Factorial Designs. Statistical Papers 42, 1-52.

Shah, D.A. and Madden, L.V. (2004). Nonparametric Analysis of Ordinal Data in Designed Factorial Experiments. Phytopathology 94, 33-43. Electronic Appendix: Instructions on the Use of Software and Applications (e-extra).

Re: PROC MIXED and non normal data

Kip1 — Tue, 24 Jan 2023 22:32:41 GMT

I want to resurrect this topic if you don't mind. According to SAS documentation Overview PROC MIXED

"The primary assumptions underlying the analyses performed by PROC MIXED are as follows:

The data are normally distributed (Gaussian).
The means (expected values) of the data are linear in terms of a certain set of parameters.
The variances and covariances of the data are in terms of a different set of parameters, and they exhibit a structure matching one of those available in PROC MIXED." Here it states explicitly "the data" not "the residuals". Anyone care to elaborate?