t-test: p=0.00768647
Neutral set average density: 0.00384848
Different set average density: 0.00438948
motif_present gene_set absent present different 360 634 neutral 1112 841
Fisher's exact test: p=1.26738e-26
Binomial GLM: p=6.26938e-26
Call:
glm(formula = motif_present ~ different, family = binomial, data = daf)
Deviance Residuals:
Min 1Q Median 3Q Max
-1.4252 -1.0613 0.9484 1.2981 1.2981
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -0.27932 0.04570 -6.112 9.82e-10 ***
differentTRUE 0.84527 0.08027 10.530 < 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: 4085.4 on 2946 degrees of freedom
Residual deviance: 3971.2 on 2945 degrees of freedom
AIC: 3975.2
Number of Fisher Scoring iterations: 4
Binomial GLM adjusted for length: p=0.0489928
Call:
glm(formula = motif_present ~ length + different, family = binomial,
data = daf)
Deviance Residuals:
Min 1Q Median 3Q Max
-3.2997 -0.8697 0.0660 0.9849 1.8160
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -1.6421776 0.0791511 -20.747 <2e-16 ***
length 0.0076538 0.0003745 20.438 <2e-16 ***
differentTRUE 0.1859130 0.0944366 1.969 0.049 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for binomial family taken to be 1)
Null deviance: 4085.4 on 2946 degrees of freedom
Residual deviance: 3287.6 on 2944 degrees of freedom
AIC: 3293.6
Number of Fisher Scoring iterations: 5