BookmarkSubscribeRSS Feed
ANKH1
Pyrite | Level 9

Hello, 

We are running the following code:

 

proc surveyreg data=sample1 varmethod=brr;
CLASS ANIMALCAT;
MODEL boxcox_FSDDFI1=ANIMALCAT SMK_202 /ANOVA;
LSMEANS ANIMALCAT/DIFF=CONTROL('50-74.99') MEANS CL ADJUST=DUNNETT PDIFF;
WEIGHT WTS_P;
REPWEIGHTS BSW1-BSW500;
RUN;

 

And we get the following output:

 

Capture.PNG

 

But we don't get an adjusted p-value for each comparison with control. What are we missing?

 

Thank you in advance!

 

4 REPLIES 4
Reeza
Super User
What does the log show? Notes and Warnings.
ANKH1
Pyrite | Level 9

It says the following:

WARNING: The ADJUST=DUNNETT p-values and confidence limits cannot be computed because the
Dunnett-Hsu approximation did not converge. Try ADJUST=SIMULATE.
NOTE: PROCEDURE SURVEYREG used (Total process time):
real time 27.27 seconds
cpu time 23.83 seconds

 

 

I used ADJUST= SIMULATE and got this

 

Capture.PNG

 

Is SIMULATE the same as DUNNETT? Can I report that p-value and say post-hoc DUNNETT was used? 

Reeza
Super User
I don't think they're the same. Check the docs for the differences and for why it may not converge.
ANKH1
Pyrite | Level 9

I checked support SAS website and found this: 

 

Problem Note 36073: ADJUST=DUNNETT DIFF=CONTROL multiple comparisons may not converge in PROC LIFETEST

 

"The Dunnett-Hsu approximation depends on an iterative algorithm to approximate the covariance matrix. In the k-sample tests scenario, the covariance matrix is singular, which causes problems with the iterative algorithm.

To circumvent the problem, use ADJUST=SIMULATE."

 

Should it be ok then to use SIMULATE? Do you recommend any documents that I can look at?

hackathon24-white-horiz.png

The 2025 SAS Hackathon Kicks Off on June 11!

Watch the live Hackathon Kickoff to get all the essential information about the SAS Hackathon—including how to join, how to participate, and expert tips for success.

YouTube LinkedIn

What is ANOVA?

ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.

Find more tutorials on the SAS Users YouTube channel.

Discussion stats
  • 4 replies
  • 1646 views
  • 0 likes
  • 2 in conversation