Use R to do the mean/variance fitting, use better fitting method