topic Re: Print more decimals for p-value in proc freq in New SAS User

Print more decimals for p-value in proc freq

mlensing — Mon, 30 Nov 2020 21:01:40 GMT

Hi SAS communities,

I am looking for a solution to print more decimals for my p-values generated via a proc freq chisq procedure. I would ideally like to have 10-12 decimals printed as I am performing bivariate analyses and need to compare highly significant values (currently all my p-values are printing <0.0001). I'm definitely open to other ways to get a p-value other than a proc freq chisq procedure, I just need to be able to test the association between two variables.

Thank you!

Re: Print more decimals for p-value in proc freq

FreelanceReinh — Mon, 30 Nov 2020 21:51:26 GMT

Hi @mlensing,

I mostly use an ODS output dataset if I want to see more decimals of a statistic. In your example add

ods output chisq=csq;

before (or in) the PROC FREQ step -- csq is just an arbitrary dataset name. Then use, e.g., PROC PRINT to show the values with a suitable format:

proc print data=csq;
format prob 14.12;
run;

To change the display format directly in PROC FREQ output, you would need to modify the ODS template. I recently suggested this in another thread: https://communities.sas.com/t5/Statistical-Procedures/Proc-mixed-formats-of-the-output-dataset/m-p/699850/highlight/true#M33761

Re: Print more decimals for p-value in proc freq

ballardw — Mon, 30 Nov 2020 21:55:48 GMT

You will likely get some push back on "more significant" (I'm one) but just specify a larger value for the format such as

format variablename f32.30; to display up to 30 decimals. 32 is the maximum number of digits the format supports.

So send the output of the test to data set and use another procedure to print the result.

Note that there is likely to be some machine precision issues with anything showing 15 or more decimal places.

If you didn't pick your rejection regions/ significance level before the test the values are "not more significant".

Re: Print more decimals for p-value in proc freq

mkeintz — Tue, 01 Dec 2020 01:01:41 GMT

Are you aware of statistical analyses or reports that find justification in ranking/comparing numerous p-values, all less than <.0001? I'd be very interested in knowing the context.

Re: Print more decimals for p-value in proc freq

mlensing — Tue, 01 Dec 2020 02:25:50 GMT

@ballardw @mkeintz

Hi all, thank you for your clarifying comments. I am needing to "compare significance" among my bivariate interactions to determine the order in which terms should be introduced into my logistic regression model to compare model fit through -2LL tests. The instructors for the course I am in that suggest this step tend to be against automated variable selection tools and so this is part of our intentional/purposeful modeling.

Re: Print more decimals for p-value in proc freq

mlensing — Tue, 01 Dec 2020 04:54:07 GMT

Thank you so much!