turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Find a Community

Topic Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

12-18-2015 04:05 PM

I am looking for SAS syntax or SAS macro with Code that will perform Goodness of Fits on discrete categorical data.

Specifically one that shows how to do the math described by N. Cressie, T.R.C. Read, in their paper entitled Multinomial goodness-of-fit tests, J. Roy. Statist. Soc. B 46 (1984) 440–464.

I do have code now for 4 of the 5 test statistics I need.

- proc freq gives the Pearson's Chi Square,
- Proc Univarite produces the Komogorov Smirnov test statistic,
- Proc Probit produces the Log Likelihood Ratio as does Proc GenMod,
- Proc Multtest calculates the Freeman-Tukey test statistic,

but I have not been able to find the Power divergence test show in the image attached here. I found some R code so maybe I can run it inside SAS enterprise guide? But I do NOT have a SAS/IML or SAS Enterprise Miner license so I can not run the R code inside a node that way. Please advise.

Message 1 of 2
(449 Views)

Accepted Solutions

Solution

12-28-2015
03:33 PM

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

Posted in reply to marygc

12-19-2015 06:35 AM

Looks like @PGStatsalready gave you an answer in the Statistics Comunity: https://communities.sas.com/t5/SAS-Statistical-Procedures/Power-Divergence-Test-Statistic/m-p/240111

These computations (all five of them!) can be done with SQL or IML.

```
proc freq data=sashelp.heart;
where weight_status is not missing;
table sex*weight_status / out=cellcounts outexpect;
run;
proc iml;
use cellCounts;
read all var {count expected};
close cellCounts;
lambda = 2/3;
powerDivergence = 2/(lambda*(1+lambda))*sum(COUNT#(COUNT/EXPECTED)##lambda-1);
print lambda powerDivergence;
```

All Replies

Solution

12-28-2015
03:33 PM

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

Posted in reply to marygc

12-19-2015 06:35 AM

Looks like @PGStatsalready gave you an answer in the Statistics Comunity: https://communities.sas.com/t5/SAS-Statistical-Procedures/Power-Divergence-Test-Statistic/m-p/240111

These computations (all five of them!) can be done with SQL or IML.

```
proc freq data=sashelp.heart;
where weight_status is not missing;
table sex*weight_status / out=cellcounts outexpect;
run;
proc iml;
use cellCounts;
read all var {count expected};
close cellCounts;
lambda = 2/3;
powerDivergence = 2/(lambda*(1+lambda))*sum(COUNT#(COUNT/EXPECTED)##lambda-1);
print lambda powerDivergence;
```