def test_rmcorr(self):
"""Test function rm_corr"""
df = read_dataset('rm_corr')
# Test again rmcorr R package.
stats = rm_corr(data=df, x='pH', y='PacO2', subject='Subject').round(3)
assert stats.at["rm_corr", "r"] == -0.507
assert stats.at["rm_corr", "dof"] == 38
assert np.allclose(np.round(stats.at["rm_corr", "CI95%"], 2),
[-0.71, -0.23])
assert stats.at["rm_corr", "pval"] == 0.001
# Test with less than 3 subjects (same behavior as R package)
with pytest.raises(ValueError):
rm_corr(data=df[df['Subject'].isin([1, 2])], x='pH', y='PacO2',
subject='Subject')
def test_rmcorr(self):
"""Test function rm_corr"""
df = read_dataset('rm_corr')
# Test again rmcorr R package.
r, p, dof = rm_corr(data=df, x='pH', y='PacO2', subject='Subject')
assert r == -0.507
assert dof == 38
assert np.round(p, 3) == 0.001
