Output from logidiab.sas

Source
4 Graphs

Diabetes data: logistic regression 1

                            The LOGISTIC Procedure
    Data Set: WORK.DIABETES
    Response Variable: GROUP     Clinical Group
    Response Levels: 3
    Number of Observations: 145
    Link Function: Logit

                               Response Profile

                    Ordered
                      Value  GROUP                   Count
                          1  (3) Overt Diabetic         33
                          2  (2) Chem. Diabetic         36
                          3  (1) Normal                 76
WARNING: 55 observation(s) were deleted due to missing values for the
         response or explanatory variables.

               Score Test for the Proportional Odds Assumption
                  Chi-Square = 16.0881 with 1 DF (p=0.0001)

     Model Fitting Information and Testing Global Null Hypothesis BETA=0

                              Intercept
                Intercept        and
  Criterion       Only       Covariates    Chi-Square for Covariates
  AIC             300.198        48.358         .
  SC              306.152        57.289         .
  -2 LOG L        296.198        42.358      253.840 with 1 DF (p=0.0001)
  Score              .             .          92.434 with 1 DF (p=0.0001)

                  Analysis of Maximum Likelihood Estimates

             Parameter Standard    Wald       Pr >    Standardized     Odds
 Variable DF  Estimate   Error  Chi-Square Chi-Square   Estimate      Ratio
 INTERCP1 1   -41.0995   8.3690    24.1172     0.0001            .     .
 INTERCP2 1   -29.1971   5.9300    24.2425     0.0001            .     .
 GLUTEST  1     0.0691   0.0142    23.5179     0.0001    12.070608    1.072

Diabetes data: logistic regression 2

                            The LOGISTIC Procedure

        Association of Predicted Probabilities and Observed Responses
                  Concordant = 72.6%          Somers' D = 0.723
                  Discordant =  0.2%          Gamma     = 0.993
                  Tied       = 27.2%          Tau-a     = 0.446
                  (6432 pairs)                c         = 0.862

Printout of first 10 observations of PROBS 3

  OBS    PATIENT      GROUP       GLUTEST         _LEVEL_            PROB
    1       1       (1) Normal      356      (3) Overt Diabetic    .0000001
    2       1       (1) Normal      356      (2) Chem. Diabetic    .0098922
    3       2       (1) Normal      289      (3) Overt Diabetic    .0000000
    4       2       (1) Normal      289      (2) Chem. Diabetic    .0000976
    5       3       (1) Normal      319      (3) Overt Diabetic    .0000000
    6       3       (1) Normal      319      (2) Chem. Diabetic    .0007750
    7       4       (1) Normal      356      (3) Overt Diabetic    .0000001
    8       4       (1) Normal      356      (2) Chem. Diabetic    .0098922
    9       5       (1) Normal      323      (3) Overt Diabetic    .0000000
   10       5       (1) Normal      323      (2) Chem. Diabetic    .0010214

4 Individual Density Plots

                            The LOGISTIC Procedure
    Data Set: WORK.DIABETES
    Response Variable: GROUP     Clinical Group
    Response Levels: 3
    Number of Observations: 145
    Link Function: Logit

                               Response Profile

                    Ordered
                      Value  GROUP                   Count
                          1  (3) Overt Diabetic         33
                          2  (2) Chem. Diabetic         36
                          3  (1) Normal                 76
WARNING: 55 observation(s) were deleted due to missing values for the
         response or explanatory variables.

               Score Test for the Proportional Odds Assumption
                  Chi-Square = 16.0881 with 1 DF (p=0.0001)

     Model Fitting Information and Testing Global Null Hypothesis BETA=0

                              Intercept
                Intercept        and
  Criterion       Only       Covariates    Chi-Square for Covariates
  AIC             300.198        48.358         .
  SC              306.152        57.289         .
  -2 LOG L        296.198        42.358      253.840 with 1 DF (p=0.0001)
  Score              .             .          92.434 with 1 DF (p=0.0001)

                  Analysis of Maximum Likelihood Estimates

             Parameter Standard    Wald       Pr >    Standardized     Odds
 Variable DF  Estimate   Error  Chi-Square Chi-Square   Estimate      Ratio
 INTERCP1 1   -41.0995   8.3690    24.1172     0.0001            .     .
 INTERCP2 1   -29.1971   5.9300    24.2425     0.0001            .     .
 GLUTEST  1     0.0691   0.0142    23.5179     0.0001    12.070608    1.072

5 Individual Density Plots

                            The LOGISTIC Procedure

        Association of Predicted Probabilities and Observed Responses
                  Concordant = 72.6%          Somers' D = 0.723
                  Discordant =  0.2%          Gamma     = 0.993
                  Tied       = 27.2%          Tau-a     = 0.446
                  (6432 pairs)                c         = 0.862

6 Individual Density Plots

                               CATMOD PROCEDURE
       Response: GROUP                       Response Levels (R)=     3
       Weight Variable: None                 Populations     (S)=   118
       Data Set: DIABETES                    Total Frequency (N)=   145
       Frequency Missing: 55                 Observations  (Obs)=   145



7 Individual Density Plots

                             POPULATION PROFILES
                                             Sample
                          Sample  GLUTEST     Size
                          -------------------------
                              1     269           1
                              2     277           1
                              3     289           1
                              4     290           1
                              5     291           2
                              6     296           2
                              7     301           1
                              8     304           1
                              9     306           2
                             10     312           1
                             11     313           1
                             12     318           2
                             13     319           2
                             14     321           1
                             15     323           3
                             16     325           1
                             17     327           1
                             18     328           1
                             19     332           1
                             20     333           1
                             21     334           1
                             22     335           1
                             23     336           1
                             24     345           1
                             25     347           1
                             26     349           1
                             27     350           1
                             28     351           1
                             29     352           2
                             30     353           2
                             31     356           4
                             32     359           1
                             33     360           3
                             34     364           2
                             35     365           2
                             36     367           1
                             37     371           1
                             38     373           1
                             39     376           2
                             40     378           2
                             41     379           1
                             42     381           1
                             43     385           1
                             44     386           1
                             45     390           1
                             46     391           2
                             47     393           3
                             48     396           1

8 Individual Density Plots

                             POPULATION PROFILES
                                             Sample
                          Sample  GLUTEST     Size
                          -------------------------
                             49     398           2
                             50     403           1
                             51     413           1
                             52     414           1
                             53     416           1
                             54     418           1
                             55     423           1
                             56     425           1
                             57     426           2
                             58     429           1
                             59     436           1
                             60     439           1
                             61     442           1
                             62     455           1
                             63     456           1
                             64     465           1
                             65     466           1
                             66     469           1
                             67     472           4
                             68     476           1
                             69     477           1
                             70     478           1
                             71     486           1
                             72     503           2
                             73     517           1
                             74     522           1
                             75     527           1
                             76     533           1
                             77     537           1
                             78     538           1
                             79     540           1
                             80     541           1
                             81     557           1
                             82     558           1
                             83     562           1
                             84     568           1
                             85     580           1
                             86     599           1
                             87     613           1
                             88     636           1
                             89     643           1
                             90     670           1
                             91     714           1
                             92     741           1
                             93     832           1
                             94     847           1
                             95     849           1
                             96     854           1

9 Individual Density Plots

                             POPULATION PROFILES
                                             Sample
                          Sample  GLUTEST     Size
                          -------------------------
                             97     857           1
                             98     920           1
                             99     923           1
                            100     958           1
                            101     967           1
                            102     972           1
                            103    1001           1
                            104    1025           1
                            105    1113           1
                            106    1133           1
                            107    1183           1
                            108    1246           1
                            109    1263           1
                            110    1354           1
                            111    1364           1
                            112    1373           1
                            113    1428           1
                            114    1468           1
                            115    1470           1
                            116    1487           1
                            117    1520           1
                            118    1568           1



                              RESPONSE PROFILES

                               Response  GROUP
                               ---------------
                                    1      1
                                    2      2
                                    3      3



10 Individual Density Plots

                          MAXIMUM-LIKELIHOOD ANALYSIS

                              Sub        -2 Log     Convergence
               Iteration   Iteration   Likelihood    Criterion
               ------------------------------------------------
                    0           0       318.59756       1.0000
                    1           0       183.83628       0.4230
                    2           0       132.96952       0.2767
                    3           0       85.735908       0.3552
                    4           0       57.297412       0.3317
                    5           0        42.24083       0.2628
                    6           0       34.720224       0.1780
                    7           0        31.32446       0.0978
                    8           0       29.850327       0.0471
                    9           0        29.33869       0.0171
                   10           0       29.248154     0.003086
                   11           0       29.244598     0.000122
                   12           0       29.244592    2.1775E-7
                   13           0       29.244592    7.316E-13


                                   Parameter Estimates
          Iteration        1           2           3           4
          ---------------------------------------------------------
               0              0           0           0           0
               1         4.4095      2.1243   -0.006475   -0.003794
               2         8.2124      4.0249     -0.0139   -0.005714
               3        15.1736      6.3116     -0.0288   -0.009138
               4        24.6085      9.3020     -0.0495     -0.0140
               5        36.3701     12.7753     -0.0757     -0.0201
               6        50.6706     16.5016     -0.1078     -0.0267
               7        66.3692     18.8878     -0.1440     -0.0309
               8        82.5671     19.4327     -0.1823     -0.0318
               9        97.7806     19.4527     -0.2185     -0.0319
              10       107.3411     19.4528     -0.2412     -0.0319
              11       109.7202     19.4528     -0.2469     -0.0319
              12       109.8264     19.4528     -0.2471     -0.0319
              13       109.8266     19.4528     -0.2471     -0.0319



                MAXIMUM-LIKELIHOOD ANALYSIS-OF-VARIANCE TABLE

              Source                   DF   Chi-Square      Prob
              --------------------------------------------------
              INTERCEPT                 2        13.57    0.0011
              GLUTEST                   2        12.70    0.0017

              LIKELIHOOD RATIO        232        26.47    1.0000

11 Individual Density Plots

                   ANALYSIS OF MAXIMUM-LIKELIHOOD ESTIMATES

                                              Standard    Chi-
       Effect            Parameter  Estimate    Error    Square   Prob
       ----------------------------------------------------------------
       INTERCEPT                 1     109.8   43.0212     6.52  0.0107
                                 2   19.4528    6.4615     9.06  0.0026
       GLUTEST                   3   -0.2471    0.1018     5.90  0.0152
                                 4   -0.0319    0.0111     8.18  0.0042



12

Graphic output from logidiab.sas

PLOT OF LOGIT * GLUTEST = _LEVEL_
[GPLOT 09/27/99 ]
PLOT OF PROB * GLUTEST = _LEVEL_
[GPLOT1 09/27/99 ]
PLOT OF PROB_3 * GLUTEST
[GPLOT2 09/27/99 ]
PLOT OF PROB_3 * GLUTEST
[GPLOT3 09/27/99 ]