topic Re: Code for Logistic Regression in SAS Programming

Code for Logistic Regression

prjadhav00 — Fri, 18 Jan 2019 01:06:33 GMT

Hello everyone,

Can anyone suggest me a code for logistic regression

Question: Perform a logistic regression with covariates age/sex/race on the outcome of having at least 1 visit to a Primary Care Provider. (PCP = Primary Care Provider)

ID PCP SEX RACE AGE

1 8   1   1   71
2 8   1   1   65
3   8   2   1   81
4   8   1   1   75
5   8   2   1   83
6   8   2   1   79
7   8   1   1   66
8   8   2   1   71
9 8   1   1   86
10 8   2   1   70
11 8   2   1   76
12 8   1   1   75
13 8   2   1   66
14 8   2   1   91
15 8   2   1   67
    Any help will be appreciated.

Re: Code for Logistic Regression

PeterClemmensen — Fri, 18 Jan 2019 07:46:33 GMT

The first thing you need to do is to get your data in a form that suits the question. If you want to model the outcome of having at least 1 visit to a Primary Care Provider, then PCP should be a binary variable, created like below. The example code for the logistic regression can then easily be found in the Examples Section of the PROC LOGISTIC Documentation.

data have;
input ID PCP SEX RACE AGE;
datalines;
1 8 1 1 71
2 8 1 1 65
3 8 2 1 81
4 8 1 1 75
5 8 2 1 83
6 8 2 1 79
7 8 1 1 66
8 8 2 1 71
9 8 1 1 86
10 8 2 1 70
11 8 2 1 76
12 8 1 1 75
13 8 2 1 66
14 8 2 1 91
15 8 2 1 67
;

data RegData;
   set have;
   PCP=(PCP>0);
run;

Re: Code for Logistic Regression

PaigeMiller — Fri, 18 Jan 2019 14:22:13 GMT

If I am understanding your problem properly, this data doesn't seem to match the problem description, nor does the addition by @PeterClemmensen. You would be doing a logistic regression where the responses are 1 for every observation, which doesn't work in logistic regression. You would need some zeros in the Y variable to fit a logistic regression.

Re: Code for Logistic Regression

prjadhav00 — Fri, 18 Jan 2019 14:40:55 GMT

Thanks for the suggestion

Re: Code for Logistic Regression

prjadhav00 — Fri, 18 Jan 2019 15:37:42 GMT

Could you tell me what I am supposed to put in place of 'have'.

data RegData;
   set have..........NOT CLEAR;
   PCP=(PCP>0);
run;

Re: Code for Logistic Regression

Reeza — Fri, 18 Jan 2019 16:25:28 GMT

SET is telling SAS what your input data set is, so what's your input data set.