Try this or search help for probit model.

For the probit model, it’s

glm (formula , family=binomial(link=”probit”))

where ‘formula’ stands for the description of your model to befitted (e.g. for an binary outcome Y dependent on two variables X andZ, replace formula by Y ~ X + Z). ‘family” describes the errordistribution and the link function to be used in the model (the R helpprovides more detailed information).