t-test: p=1.55453e-48
Neutral set average density: 0.0685891
Different set average density: 0.0838143
motif_present gene_set absent present different 4 990 neutral 60 1893
Fisher's exact test: p=2.7891e-07
Binomial GLM: p=6.96676e-05
Call: glm(formula = motif_present ~ different, family = binomial, data = daf) Deviance Residuals: Min 1Q Median 3Q Max -3.3213 0.0898 0.2498 0.2498 0.2498 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) 3.4516 0.1311 26.322 < 2e-16 *** differentTRUE 2.0598 0.5179 3.977 6.97e-05 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for binomial family taken to be 1) Null deviance: 616.80 on 2946 degrees of freedom Residual deviance: 588.18 on 2945 degrees of freedom AIC: 592.18 Number of Fisher Scoring iterations: 8
Binomial GLM adjusted for length: p=0.127486
Call: glm(formula = motif_present ~ length + different, family = binomial, data = daf) Deviance Residuals: Min 1Q Median 3Q Max -3.2781 0.0001 0.0131 0.1246 1.0557 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) -0.847385 0.380352 -2.228 0.0259 * length 0.057026 0.006894 8.272 <2e-16 *** differentTRUE 0.830053 0.544622 1.524 0.1275 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for binomial family taken to be 1) Null deviance: 616.8 on 2946 degrees of freedom Residual deviance: 395.3 on 2944 degrees of freedom AIC: 401.3 Number of Fisher Scoring iterations: 11