topic Re: Chi square in Statistical Procedures

Chi square

L_L — Tue, 10 May 2011 11:56:48 GMT

Dear all,

I would like to add in a contingency table (created by proc tabulate) the corrispondence p-value: is it possible?
Anyway, how can I do it?

Thanks in advance
Regards

Re: Chi square

Paige — Tue, 10 May 2011 11:57:39 GMT

I don't think PROC TABULATE will do this.

I believe you have to create the table in PROC FREQ and use the CHISQ option.

Re: Chi square

L_L — Tue, 10 May 2011 12:42:31 GMT

Ok. The PROC FREQ output (with the chi square option) is
1) contingency table
2) Chi square test
I would like to have in the contingency table the p-value, for example:
Group1 Group2 Group3 *
Male xxx xxx xxx
xx% xx% xx%
Female xxx xxx xxx
xx% xx% xx%
*p=x.xxx

Any ideas?
thx

Re: Chi square

Paige — Tue, 10 May 2011 15:14:44 GMT

It shouldn't be hard.

proc freq data=whatever;
table gender*group/chisq;
run;

Re: Chi square

L_L — Tue, 10 May 2011 17:18:23 GMT

Thank you. Unfortunaley it's not exactly what I mean.
The aim is to improve the output and to have all the information in one table.
Maybe I have to use a PUT on a file and a PRINTTO...

Re: Chi square

SPR — Tue, 10 May 2011 18:19:17 GMT

Hello L&L,

Is this what you need?
[pre]
ODS output CHISQ=Chi (keep=Statistic Prob where=(Statistic="Chi-Square"));
proc freq data=SASHELP.CLASS;
table sex*age/out=freq chisq;
run;
proc transpose data=freq out=r(drop=_:) prefix=Age_;
var percent;
id age;
by sex;
run;
data r;
length Sex $10;
set r CHI(in=c rename=(Statistic=Sex Prob=Age_11));
if c then Sex="p-value";
run;
[/pre]
Sincerely,
SPR Message was edited by: SPR

Re: Chi square

L_L — Tue, 10 May 2011 20:23:32 GMT

Yes!! Thank you so much :-))