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

% Open a log file
logopen using C:\monographregression\computercode\results\bayesx\orangejuice_additive.log

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

% Create a regression object
bayesreg m

% Set a new delimiter
delimiter=;

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

% Estimation of the model with nonlinear covariate effects and random effects
m.hregress y = outlet(random) +
logp(psplinerw2,monotone=decreasing) + 
logpc1(psplinerw2,monotone=increasing) +
logpc2(psplinerw2,monotone=increasing) +
logpc3(psplinerw2,monotone=increasing) +
logp*outlet(random,nofixed) + logpc1*outlet(random,nofixed) + 
logpc2*outlet(random,nofixed) + logpc3*outlet(random,nofixed), 
setseed=12345 family=gaussian iterations=102000 burnin=2000 step=100 predict using d;

% Close the log file
logclose