def test_odof_with_without_media():
gdsc = ANOVA(ic50_test)
_res = gdsc.anova_one_drug_one_feature('Drug_1047_IC50', 'TP53_mut')
dd1 = gdsc._get_anova_summary(gdsc.data_lm, output='dict')
assert_list_almost_equal([dd1['feature'], dd1['msi'], dd1['tissue']],
[1.5750735472022118e-58, 0.025902887791637515,
5.541879283763767e-44])
gdsc = ANOVA(ic50_test, set_media_factor=True)
_res = gdsc.anova_one_drug_one_feature('Drug_1047_IC50', 'TP53_mut')
dd2 = gdsc._get_anova_summary(gdsc.data_lm, output='dict')
assert_list_almost_equal([dd2['feature'], dd2['media'], dd2['msi'],
dd2['tissue']], [ 2.9236500715529455e-58, 0.7762487502315283,
0.023777744527686766, 1.5729157319290974e-44])
def test_odof_with_without_media():
gdsc = ANOVA(ic50_test)
_res = gdsc.anova_one_drug_one_feature('Drug_1047_IC50', 'TP53_mut')
dd1 = gdsc._get_anova_summary(gdsc.data_lm, output='dict')
assert_list_almost_equal(
[dd1['feature'], dd1['msi'], dd1['tissue']],
[1.5750735472022118e-58, 0.025902887791637515, 5.541879283763767e-44])
gdsc = ANOVA(ic50_test, set_media_factor=True)
_res = gdsc.anova_one_drug_one_feature('Drug_1047_IC50', 'TP53_mut')
dd2 = gdsc._get_anova_summary(gdsc.data_lm, output='dict')
assert_list_almost_equal(
[dd2['feature'], dd2['media'], dd2['msi'], dd2['tissue']], [
2.9236500715529455e-58, 0.7762487502315283, 0.023777744527686766,
1.5729157319290974e-44
])
def test_odof_with_without_media():
gdsc = ANOVA(ic50_test)
_res = gdsc.anova_one_drug_one_feature(1047, 'TP53_mut')
odof = gdsc._get_one_drug_one_feature_data(1047, 'TP53_mut')
dd1 = gdsc._get_anova_summary(gdsc.data_lm, output='dict', odof=odof)
assert_almost_equal(dd1['feature'] - 1.5750735472022118e-58, delta=1e-64)
assert_almost_equal(dd1['msi'] - 0.025902887791637515, delta=1e-10)
# Change in Nov 2016
#assert_almost_equal(dd1['tissue'], 5.541879283763767e-44, delta=1e-50)
assert_almost_equal(dd1['tissue'] - 1.0258702741509e-44, delta=1e-50)
gdsc = ANOVA(ic50_test, set_media_factor=True)
_res = gdsc.anova_one_drug_one_feature(1047, 'TP53_mut')
odof = gdsc._get_one_drug_one_feature_data(1047, 'TP53_mut')
dd2 = gdsc._get_anova_summary(gdsc.data_lm, output='dict', odof=odof)
assert_almost_equal(dd2['feature'] - 2.9236500715529455e-58, delta=1e-64)
assert_almost_equal(dd2['media'] - 0.7762487502315283, delta=1e-10)
assert_almost_equal(dd2['msi'] - 0.023777744527686766, delta=1e-10)
assert_almost_equal(dd2['tissue'] - 1.5729157319290974e-44, delta=1e-50)
def test_odof_with_without_media():
gdsc = ANOVA(ic50_test)
_res = gdsc.anova_one_drug_one_feature(1047, "TP53_mut")
odof = gdsc._get_one_drug_one_feature_data(1047, "TP53_mut")
dd1 = gdsc._get_anova_summary(gdsc.data_lm, output="dict", odof=odof)
assert_almost_equal(dd1["feature"], 1.5750735472022118e-58, delta=1e-64)
assert_almost_equal(dd1["msi"], 0.025902887791637515, delta=1e-10)
# Change in Nov 2016
# assert_almost_equal(dd1['tissue'], 5.541879283763767e-44, delta=1e-50)
assert_almost_equal(dd1["tissue"], 1.0258702741509e-44, delta=1e-50)
gdsc = ANOVA(ic50_test, set_media_factor=True)
_res = gdsc.anova_one_drug_one_feature(1047, "TP53_mut")
odof = gdsc._get_one_drug_one_feature_data(1047, "TP53_mut")
dd2 = gdsc._get_anova_summary(gdsc.data_lm, output="dict", odof=odof)
assert_almost_equal(dd2["feature"], 2.9236500715529455e-58, delta=1e-64)
assert_almost_equal(dd2["media"], 0.7762487502315283, delta=1e-10)
assert_almost_equal(dd2["msi"], 0.023777744527686766, delta=1e-10)
assert_almost_equal(dd2["tissue"], 1.5729157319290974e-44, delta=1e-50)
