// import the data into a dataset called Data and generate a constant term
Data = webimport("www.bayeconsoft.com/datasets/CCES2017.csv");
Data.constant = 1;
// run the model
myMVProbit = mvprobit( {
investor ~ constant age educ male faminc,
savings ~ constant age educ male faminc,
homeowner ~ constant age educ male faminc
} );
// calculate marginal effects marginally with respect to the remaining
// response variables
mfx( "model"=myMVProbit );
// calculate marginal effects conditionally on the remaining response variables
// being equal to 1
mfx( "model"=myMVProbit, "type"=2, "opt" = [1,1,1] );