- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
I'm trying to get SAS to spit out t-values for my contrasts in LSMEANS. Is there any way to achieve this? It just reports p-values for each comparison.
My code is
LSMEANS species*wax / pdiff stderr;
Accepted Solutions
- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
I think you should switch to PROC MIXED (or GLIMMIX). Much better functionality, and you get the t values and df (using many different possible df methods).
- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
In general, if p is a probability for a t-distributed statistic, then you can recover the t value by using the QUANTILE function. You can use the OUT= option on the LSMEANS statement to produce an output data set, or you can use ODS OUTPUT.
However, you have to know the degrees of freedom for the pairwise comparison and the significance level of the test. A discussion with formulas is in the doc: SAS/STAT(R) 13.1 User's Guide
I am not an expert in multiple comparisons, but I know that they can be tricky, which might be why the procedure does not output the degrees of freedom.
- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
I think you should switch to PROC MIXED (or GLIMMIX). Much better functionality, and you get the t values and df (using many different possible df methods).
- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
Thanks all.