Calculate degrees of freedom in r