% usefile c:\oupbuch\foresthealth\foresthealth_ch2.prg dataset d d.infile using c:\oupbuch\foresthealth\foresthealth.dat remlreg r % Trend models without any further covariates % global trend common to all categories r.outfile = c:\oupbuch\foresthealth\ch2\trendsglobal logopen using c:\oupbuch\foresthealth\ch2\trendsglobal.log r.regress bu3 = time(psplinerw2), family=cumlogit using d logclose % category-specific trend functions r.outfile = c:\oupbuch\foresthealth\ch2\trendscatspec logopen using c:\oupbuch\foresthealth\ch2\trendscatspec.log r.regress bu3 = time(psplinerw2, catspecific), family=cumlogit using d logclose % Additive models r.outfile = c:\oupbuch\foresthealth\ch2\additiveglobal logopen using c:\oupbuch\foresthealth\ch2\trendscatspec.log r.regress bu3 = time(psplinerw2) + age(psplinerw2) + elevation + inclination + soil + ph + canopy + standcat + fertilisation + humus0 + humus2 + humus3 + humus4 + moisture1 + moisture3 + alkali1 + alkali3 + alkali4, family=cumlogit using d logclose r.outfile = c:\oupbuch\foresthealth\ch2\additivecatspec logopen using c:\oupbuch\foresthealth\ch2\trendscatspec.log r.regress bu3 = time(psplinerw2, catspecific) + age(psplinerw2) + elevation + inclination + soil + ph + canopy + standcat + fertilisation + humus0 + humus2 + humus3 + humus4 + moisture1 + moisture3 + alkali1 + alkali3 + alkali4, family=cumlogit using d logclose drop d r