List the output of the probit model in an organized fashion.
Calculate the mean of the data. Sum all the data points and divide by the number of data points. Call this value “m.”
Subtract the mean from every individual data point, creating a new dataset. Rewrite this new dataset in an organized fashion.
Square the new data, creating again a new dataset.
Sum the data in the new dataset.
Calculate the variance of the new dataset. Divide the sum just calculated by the number of data points. This is the estimation for the variance for the original dataset.
Assign a variable to the output of the probit model. For example, in R, decide on a variable (e.g. x), and write “x <- probitmodel,” where “probitmodel” is the output of the probit model.
Call the standard deviation command for the created variable. In R, type “sd(x).” The result will be a single number, the standard deviation of the model’s output.
Square the standard deviation. The square of the standard deviation is the estimated variance for the probit model.