// import the data into a dataset and generate a constant term
Data = webimport("www.bayeconsoft.com/datasets/UnionMembership.csv");
Data.constant = 1;
// run the Probit and Logit models
Probit = probit( union ~ constant hours married black hisp health );
Logit = logit( union ~ constant hours married black hisp health );
// calculate marginal effects from the two models at the means of the data
mfx( "model"=Probit, "point"="mean" );
mfx( "model"=Logit, "point"="mean" );