In this model, the between subject variability is characterized by the RESIDUAL variance, since you are doing an OLS solution (due to the method=type3 option). The model does not have a within-subject structure, so that is not presented.
So, when you divide the residual standard deviation by the treatment means, you get different CV's for each treatment. That is not surprising, and in fact, exactly what I would expect, as the model assumes a constant variance by group, and you have one value divided successively by 3 different values. Consequently, I think your macro is doing exactly what the data tells us. You will have to explain why you think something different should occur, and then we could perhaps come up with an alternative. And by the way, the variable se_exp is not the standard error of estim_exp. Consider this: estim = 2, se =0.2. Applying your formulas, I get estim_exp = 7.389 and se_exp=1.221. Note that se_exp is bounded below by 1, whereas the standard error of data with a mean of 2 is bounded below by 0, as n goes to infinity. Check some references for a proper conversion - the formula will be almost identical to the one you have for CV_inter_geometric (divided by sqrt(N))..
Finally, since this is in log space, you are implicitly making the assumption that the CV is identical for all treatment groups - it is a multiplicative model, and the standard deviation is in a constant ratio to the mean.
You may want to use PROC GENMOD on your untransformed data to get the information you want.
SteveDenham
... View more