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

- Home
- /
- Analytics
- /
- Stat Procs
- /
- PROC GENMOD statistics

Topic Options

- 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
- RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

02-22-2018 04:38 PM

Good afternoon all, I hope everyone is having a great day. I am writing because I am working with cross-sectional complex survey data and generating a log-binomial model (for prevalence ratios). The default statistical test used for each regression term is a Z statistic. Since this is a sample (without the assumption of population normality) I want to generate t statistics (or F statistics) for each term instead. How do I accomplish this? Below is an example of my code.

PROC GENMOD data=Cancer;

CLASS Cluster Stratum;

MODEL Kidney (event='1') = BMI Gluc

/ d=binomial link=log;

TITLE 'Model 1';

WEIGHT Weight;

REPEATED SUBJECT = Cluster (Stratum);

RUN;

I appreciate any help you can provide. Thanks, and have a great rest of the day.

Accepted Solutions

Solution

03-01-2018
12:25 AM

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

Posted in reply to RTurpin

02-28-2018 02:10 PM

The parameter tests produced for a GEE model in PROC GENMOD are large-sample Wald (Z) tests. If you particularly want to use t-tests instead, then you can simply use the PROBT function in the DATA step with the statistic value and the error degrees of freedom as arguments. For example, the following produces the t-test p-value for the Smoke predictor in the GEE example in the Getting Started section of the GENMOD documentation.

```
data _null_;
p = 1-probt(-0.2036,60);
put p=;
run;
```

The result is shown in the log.

```
30 data _null_;
31 p = probt(-0.2036,60);
32 put p=;
33 run;
p=0.4196775261
NOTE: DATA statement used (Total process time):
real time 0.01 seconds
cpu time 0.01 seconds
```

All Replies

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

Posted in reply to RTurpin

02-28-2018 10:09 AM

If you have survey data you should not use PROC GENMOD. Only the survey procedures have the variance estimators necessary for a proper analysis of survey data. Unfortunately, there is no survey procedure that will fit a log binomial model, but you could use PROC SURVEYLOGISTIC for fitting a logistic model to binary response data.

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

Posted in reply to StatDave_sas

02-28-2018 10:31 AM - edited 02-28-2018 10:32 AM

Thank you so much for the information. If I were using unweighted, unstratified, unclustered data (which has none of the complex survey data parameters) how would I go about reporting a t or F test for regression terms in a sample in PROC GENMOD (instead of the default Z test)?

Thanks,

Rod Turpin

Solution

03-01-2018
12:25 AM

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

Posted in reply to RTurpin

02-28-2018 02:10 PM

The parameter tests produced for a GEE model in PROC GENMOD are large-sample Wald (Z) tests. If you particularly want to use t-tests instead, then you can simply use the PROBT function in the DATA step with the statistic value and the error degrees of freedom as arguments. For example, the following produces the t-test p-value for the Smoke predictor in the GEE example in the Getting Started section of the GENMOD documentation.

```
data _null_;
p = 1-probt(-0.2036,60);
put p=;
run;
```

The result is shown in the log.

```
30 data _null_;
31 p = probt(-0.2036,60);
32 put p=;
33 run;
p=0.4196775261
NOTE: DATA statement used (Total process time):
real time 0.01 seconds
cpu time 0.01 seconds
```

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

Posted in reply to StatDave_sas

03-01-2018 12:28 AM

Thanks! This was very helpful.