t-test: p=2.79957e-11
Neutral set average density: 0.0132338
Different set average density: 0.00936181
motif_present gene_set absent present different 323 671 neutral 831 1122
Fisher's exact test: p=1.08694e-07
Binomial GLM: p=1.35291e-07
Call: glm(formula = motif_present ~ different, family = binomial, data = daf) Deviance Residuals: Min 1Q Median 3Q Max -1.4994 -1.3073 0.8865 1.0529 1.0529 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) 0.30024 0.04577 6.560 5.38e-11 *** differentTRUE 0.43088 0.08174 5.272 1.35e-07 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for binomial family taken to be 1) Null deviance: 3945.7 on 2946 degrees of freedom Residual deviance: 3917.4 on 2945 degrees of freedom AIC: 3921.4 Number of Fisher Scoring iterations: 4
Binomial GLM adjusted for length: p=1.46905e-05
Call: glm(formula = motif_present ~ length + different, family = binomial, data = daf) Deviance Residuals: Min 1Q Median 3Q Max -3.9338 -0.9894 0.3287 0.9809 1.7426 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) -1.1352012 0.0795665 -14.267 < 2e-16 *** length 0.0091278 0.0004602 19.833 < 2e-16 *** differentTRUE -0.4371164 0.1008743 -4.333 1.47e-05 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for binomial family taken to be 1) Null deviance: 3945.7 on 2946 degrees of freedom Residual deviance: 3210.1 on 2944 degrees of freedom AIC: 3216.1 Number of Fisher Scoring iterations: 5