/******************* plant97a.sps **************************/
set width=80  length=none.
title 'Logistic Regression on Died'

data list file = 'plantcat.dat'  free
   / obs plant mcg died whendied health7 girdling whengird progr .

variable labels
     mcg    'Mycelial Compatibility Group'
     plant  'Cultivar' .
value labels
     plant 1 'GP159'  2  'Hanna'   3  'Westar'/
     died  0 'No'     1  'Yes'.

crosstabs tables =  mcg plant by died
                   / mcg by died by plant
                   / cells=count row
                   / statistics=chisq .

logistic regression died with plant mcg plant by mcg
    / categorical plant mcg
    / enter plant
    / enter mcg
    / enter plant by mcg .