LISTSERV at the University of Georgia
Menubar Imagemap
Home Browse Manage Request Manuals Register
Previous messageNext messagePrevious in topicNext in topicPrevious by same authorNext by same authorPrevious page (August 2008, week 2)Back to main SAS-L pageJoin or leave SAS-L (or change settings)ReplyPost a new messageSearchProportional fontNon-proportional font
Date:   Tue, 12 Aug 2008 12:35:55 -0400
Reply-To:   cooch17@NOSPAMVERIZON.NET
Sender:   "SAS(r) Discussion" <SAS-L@LISTSERV.UGA.EDU>
From:   cooch17@NOSPAMVERIZON.NET
Organization:   Cornell University
Subject:   interpreting ESTIMATE output GLM
Comments:   To: sas-l@uga.edu
Content-Type:   text/plain; charset=ISO-8859-1; format=flowed

Suppose I have a single-classification ANOVA, where the sample means show ordinal trend among successive levels. I construct a CONTRAST to test for trend, in the usual manner. Fine - but I'm puzzled how to interpret the parameter output from ESTIMATE (for the same contrast).

For example, consider the following code snippet:

data test;

do trt=1 to 5;

do i=1 to 5000;

x=normal(0)*1.25+trt;

output;

end;

end;

proc glm;

class trt;

model x=trt / ss3;

lsmeans trt;

contrast 'linear' trt 2 1 0 -1 -2;

estimate 'linear' trt 2 1 0 -1 -2;

estimate '1 v 2' trt 1 -1 0 0 0;

estimate '2 v 3' trt 0 1 -1 0 0;

estimate '3 v 4' trt 0 0 1 -1 0;

estimate '4 v 5' trt 0 0 0 1 -1;

run;

So, 5 treatment levels, linear increase of the treatment mean by 1 for each level of trt. So, in fact, the 'slope' of the increase in the treatment effect is 1 (confirmed by each of the paired estimate statement contrasting successive treatment means). However, the output from the linear ESTIMATE is (for a given simulated sample) something like

estimate: -9.9968 SE: 0.056

So, rounding off the estimate to -10 (which is what it should be...

(2)(1)+(1)(2)+(0)(3)+(-1)(4)+(-2)(5) = 10

how do I interpret it? It clearly isn't the slope of the change in treatment effect (even if I take it and divide it by 5 - the number of levels of the treatment - it still isn't). While I *know* this is an acceptable approach to testing significance of a trend, I also need to estimate the slope of the trend. So, how to go from what ESTIMATE gives me, to what I actually want/need.

Pointers to the obvious?

Thanks in advance.


Back to: Top of message | Previous page | Main SAS-L page