The scores matrix contains the gradient observations and are used to calculate the robust variance-covariance matrix. The function checks whether the scores matrix is present in the model object. The function is primarily used to control the flow in other functions.