% den Datensatz einlesen dataset d d.infile using c:\monographieregression\homepage\download\buche.raw % trivariate Response generieren d.generate buche3 = 0 d.replace buche3 = 1 if buche >= 12.5 d.replace buche3 = 2 if buche >= 50 % Dummy-Variablen fuer Humus d.generate humus0 = 0 d.replace humus0 = 1 if humus=0 d.replace humus0 = -1 if humus=1 d.generate humus2 = 0 d.replace humus2 = 1 if humus=2 d.replace humus2 = -1 if humus=1 d.generate humus3 = 0 d.replace humus3 = 1 if humus=3 d.replace humus3 = -1 if humus=1 d.generate humus4 = 0 d.replace humus4 = 1 if humus>3 d.replace humus4 = -1 if humus=1 % Dummy-Variablen fuer Frischegrad d.generate frische1 = 0 d.replace frische1 = 1 if frische=1 d.replace frische1 = -1 if frische=2 d.generate frische3 = 0 d.replace frische3 = 1 if frische=3 d.replace frische3 = -1 if frische=2 % Dummy-Variablen fuer Alkali d.generate alkali1 = 0 d.replace alkali1 = 1 if alkali=1 d.replace alkali1 = -1 if alkali=2 d.generate alkali3 = 0 d.replace alkali3 = 1 if alkali=3 d.replace alkali3 = -1 if alkali=2 d.generate alkali4 = 0 d.replace alkali4 = 1 if alkali=4 d.replace alkali4 = -1 if alkali=2 % Kubische Polynome fuer Alter und Jahr d.generate alter2=alter*alter d.generate alter3=alter2*alter d.replace jahr=jahr-1982 d.generate jahr2=jahr*jahr d.generate jahr3=jahr2*jahr % Anlegen eines Regressionsobjekts remlreg r % Aufruf der Schaetzung % die Option family=cumlogit bewirkt die Schaetzung eines kumulativen Logit-Modells r.regress buche3 = hang + schirm + alter + alter2 + alter3 + jahr + jahr2 + jahr3 + hoehe + grund + ph + artkat + dueng + humus0 + humus2 + humus3 + humus4 + frische1 + frische3 + alkali1 + alkali3 + alkali4, family=cumlogit using d