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)