t-test: p=0.650512
Neutral set average density: 0.000166193
Different set average density: 0.000187267
motif_present gene_set absent present different 953 41 neutral 1906 47
Fisher's exact test: p=0.0115953
Binomial GLM: p=0.010445
Call: glm(formula = motif_present ~ different, family = binomial, data = daf) Deviance Residuals: Min 1Q Median 3Q Max -0.2903 -0.2903 -0.2207 -0.2207 2.7302 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) -3.7026 0.1477 -25.077 <2e-16 *** differentTRUE 0.5566 0.2173 2.561 0.0104 * --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for binomial family taken to be 1) Null deviance: 791.32 on 2946 degrees of freedom Residual deviance: 784.91 on 2945 degrees of freedom AIC: 788.91 Number of Fisher Scoring iterations: 6
Binomial GLM adjusted for length: p=0.0597693
Call: glm(formula = motif_present ~ length + different, family = binomial, data = daf) Deviance Residuals: Min 1Q Median 3Q Max -0.4441 -0.2695 -0.2224 -0.2086 2.7937 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) -3.9152150 0.1877585 -20.852 <2e-16 *** length 0.0010417 0.0005248 1.985 0.0471 * differentTRUE 0.4283930 0.2275670 1.882 0.0598 . --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for binomial family taken to be 1) Null deviance: 791.32 on 2946 degrees of freedom Residual deviance: 781.28 on 2944 degrees of freedom AIC: 787.28 Number of Fisher Scoring iterations: 6