t-test: p=4.9392e-14
Neutral set average density: 0.0252905
Different set average density: 0.0298549
motif_present gene_set absent present different 36 958 neutral 281 1672
Fisher's exact test: p=9.07633e-22
Binomial GLM: p=1.60517e-16
Call: glm(formula = motif_present ~ different, family = binomial, data = daf) Deviance Residuals: Min 1Q Median 3Q Max -2.5761 0.2716 0.5574 0.5574 0.5574 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) 1.78342 0.06447 27.661 <2e-16 *** differentTRUE 1.49791 0.18160 8.248 <2e-16 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for binomial family taken to be 1) Null deviance: 2012.2 on 2946 degrees of freedom Residual deviance: 1918.7 on 2945 degrees of freedom AIC: 1922.7 Number of Fisher Scoring iterations: 6
Binomial GLM adjusted for length: p=0.314652
Call: glm(formula = motif_present ~ length + different, family = binomial, data = daf) Deviance Residuals: Min 1Q Median 3Q Max -3.06538 0.00179 0.06804 0.39087 1.58894 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) -1.636469 0.180526 -9.065 <2e-16 *** length 0.035339 0.002263 15.616 <2e-16 *** differentTRUE 0.216339 0.215154 1.006 0.315 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for binomial family taken to be 1) Null deviance: 2012.2 on 2946 degrees of freedom Residual deviance: 1260.3 on 2944 degrees of freedom AIC: 1266.3 Number of Fisher Scoring iterations: 8