% usefile C:\monographregression\computercode\prg\bayesx\forest.txt

% Create a dataset object and read the data
dataset d
d.infile using C:\monographregression\computercode\data\bayesx\beach_regmodels.raw

% Create a map object and read the data
map m
m.infile using C:\monographregression\computercode\data\bayesx\beach.bnd

% Create a regression object
remlreg r

% Set a new delimiter
delimiter=;

% Define the path where to save the results
r.outfile = C:\monographregression\computercode\results\bayesx\forest_additive;

% Estimation of the model without spatial effect
r.regress defol3 =  age(psplinerw2) + year(psplinerw2) + gradient(psplinerw2) + canopyd(psplinerw2) +  
alt(psplinerw2) + depth + ph + type  + fert + humus0 + humus2 + humus3 + humus4 + watermoisture1 + 
watermoisture3 + alkali1 + alkali3 + alkali4, family=cumlogit using d;

% Estimation of the model with spatial effect
r.outfile = C:\monographregression\computercode\results\bayesx\forest_spatial;
r.regress defol3 = id(geospline, map=m,gridsizex=20,gridsizey=20) + age(psplinerw2) + year(psplinerw2) + 
gradient(psplinerw2) + canopyd(psplinerw2) +  alt(psplinerw2) + depth + ph + type  + fert + humus0 + 
humus2 + humus3 + humus4 + watermoisture1 + watermoisture3 + alkali1 + alkali3 + alkali4, 
family=cumlogit using d;