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?

sas-innovate-2024.png

Available on demand!

Missed SAS Innovate Las Vegas? Watch all the action for free! View the keynotes, general sessions and 22 breakouts on demand.

 

Register now!

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
  • 829 views
  • 0 likes
  • 2 in conversation