I have repeated measures data that contains subjects' performance on a task over many trials. I also have a group variable as well (three groups, and this variable is categorical- there is no order to the groups).
In NLMIXED, I get the parameters for an exponential function that describes the learning curves. Specifically, I have a parameter that reflects the asymptote and another that reflects the learning rate (ie, a time constant). I also find that there is a significant effect of my group variable on the asymptote parameter, but not the learning rate. This all looks good and the fit is nice.
**How can I test for differences between pairs of groups?** Ie, the hypothesis that Group A does not differ from Group B, but does differ from Group C. Here are some avenues I've tried, with little luck:
1. Use the CONTRAST statement. This really seems like the what I'm supposed to do, but I just can't figure out the syntax, and how to code the data to use it! I read that it's *not* executed the same was as the contrast statement in PROC MIXED, but I can't seem to find the right way.
2. Get each individual subjects' parameter estimate, then determine (in a separate analysis) whether there is a signif effect of group on the estimate. The output gives me the mean and variance of the estimate, but can I get the estimate for each individual? Or, should I run NLMIXED on each subject (without the group variable, of course), and get the parameters from each output?
Would either of these work? Thanks for any insight!
First of all, I don't think you have coded your group effects appropriately. Since you have three unordered groups, you should have three models for the asymptote and also three models for the learning rate. The three asymptote models can all be written employing a single statement of the form:
Just be sure to substitute specific values for everywhere in the above equation. Similarly, you can construct contrasts between expected values for groups 1 and 3 and also between expected values for groups 2 and 3. I'll let you work out the code for these latter two contrasts. It should be readily apparent from everything that has been shown above.
Right, of course! I felt like my version was off but I didn't know what the solution was. This makes complete sense, thanks very much for the clear and thorough reply. And thanks also for pointing out the trial contrast.