Added an expectation-propagataion fitting option and explained the difference in graf.Rd