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
- /
- odds ratio estimates in R/N logistic models

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
- Highlight
- Email to a Friend
- Report Inappropriate Content

09-07-2014 08:16 AM

I have a data set comparing animals that are treated with x versus y. The animals are in cages of varying number per cage (4-9 ). I count the number of sick animals at various time points in the cages.

My data lines include:

treatment day sick_animals animals_per_cage

I denote those R, and the rabbits per cage are N.

When I run the model

proc logistic data = ....;

class treatment (ref=control) day (ref=0):

sick_animals/animals_per_cage= treatment day/cl;

oddsratio treatment/diff=ref;

oddsratio day/diff=ref;

run;

I get different estimates for the parameters estimates (when I exponentiate these) in the 'Analysis of maximum likelihood estimates' compared to what the oddsratios that the model generates.

What is the true odds ratios, and where do I find level of significance (p-value) of the odds ratios. I can see the confidence intervals, but I need the p-value also.

Kind regards,

Anna Catharina

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

09-07-2014 01:59 PM

Your code looks like it may have been truncated as there's no model statement, and it looks like you have a colon instead of semi colon in the class statement.

try adding param=ref to the class statement. This tells SAS to use Referential coding rather than the default of Effect coding for the categorical variables.

SAS FAQ: In PROC LOGISTIC why aren't the coefficients consistent with the odds ratios?

proc logistic data = ....;

class treatment (ref=control) day (ref=0)/**Param=REF**;

sick_animals/animals_per_cage= treatment day/cl;

oddsratio treatment/diff=ref;

oddsratio day/diff=ref;

run;