t-test: p=2.72468e-26
Neutral set average density: 0.0181745
Different set average density: 0.0230244
motif_present gene_set absent present different 45 949 neutral 340 1613
Fisher's exact test: p=7.1071e-26
Binomial GLM: p=8.43429e-20
Call:
glm(formula = motif_present ~ different, family = binomial, data = daf)
Deviance Residuals:
Min 1Q Median 3Q Max
-2.4880 0.3044 0.6185 0.6185 0.6185
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 1.55691 0.05968 26.090 <2e-16 ***
differentTRUE 1.49184 0.16380 9.107 <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: 2284.5 on 2946 degrees of freedom
Residual deviance: 2172.3 on 2945 degrees of freedom
AIC: 2176.3
Number of Fisher Scoring iterations: 5
Binomial GLM adjusted for length: p=0.19446
Call:
glm(formula = motif_present ~ length + different, family = binomial,
data = daf)
Deviance Residuals:
Min 1Q Median 3Q Max
-4.2838 0.0026 0.0909 0.4444 1.5297
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -1.642561 0.163094 -10.071 <2e-16 ***
length 0.031266 0.001894 16.505 <2e-16 ***
differentTRUE 0.252447 0.194564 1.297 0.194
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for binomial family taken to be 1)
Null deviance: 2284.5 on 2946 degrees of freedom
Residual deviance: 1450.9 on 2944 degrees of freedom
AIC: 1456.9
Number of Fisher Scoring iterations: 8