I have a dataset that consists of binary output variables (0 or 1) across a dependent variable range. I am trying to plot a probit regression model with that dataset using the stat_smooth function to plot a sigmoidal trendline. However, the stat_smooth function is not recognizing the arguments I am using to make the distribution binomial.