% usefile C:\monographregression\computercode\prg\bayesx\rent_geoadditive.txt % Create a dataset object and read the data dataset d d.infile, maxobs=5000 using C:\monographregression\computercode\data\bayesx\rent99.raw % Create a map object and read the data map m m.infile using C:\monographregression\computercode\data\bayesx\munich_withoutparks.bnd % Set a new delimiter delimiter=; % Open a log file logopen , replace using C:\monographregression\computercode\results\bayesx\reml_geoadditive_kriging.log; % Create a regression object remlreg r; % Define the path where to save the results r.outfile = C:\monographregression\computercode\results\bayesx\reml_geoadditive_without; % Estimation of the geoadditive model without location-dummies r.regress rentsqm = area(psplinerw2) + yearc(psplinerw2) + district(geokriging, map=m), family=gaussian using d; % Generate dummy variables d.generate glocation=1*(location=2); d.generate tlocation=1*(location=3); % Estimation of the geoadditive model with location-dummies r.outfile = C:\monographregression\computercode\results\bayesx\reml_geoadditive_with; r.regress rentsqm = glocation+tlocation+area(psplinerw2) + yearc(psplinerw2) + district(geokriging, map=m), family=gaussian using d; % Close the log file logclose;