custom Rsp_3_CAGGAC 3' UTR

[CSV file]

Density of motif

t-test: p=0.434136

Neutral set average density: 1.19197e-05

Different set average density: 2.64189e-05

Presence of motif

           motif_present
gene_set    absent present
  different    990       4
  neutral     1947       6

Fisher's exact test: p=0.741115

Binomial GLM: p=0.675304

Call:
glm(formula = motif_present ~ different, family = binomial, data = daf)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-0.0898  -0.0898  -0.0784  -0.0784   3.4016  

Coefficients:
              Estimate Std. Error z value Pr(>|z|)    
(Intercept)    -5.7823     0.4089 -14.142   <2e-16 ***
differentTRUE   0.2709     0.6467   0.419    0.675    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 133.69  on 2946  degrees of freedom
Residual deviance: 133.51  on 2945  degrees of freedom
AIC: 137.51

Number of Fisher Scoring iterations: 8

Binomial GLM adjusted for length: p=0.860495

Call:
glm(formula = motif_present ~ length + different, family = binomial, 
    data = daf)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-0.2828  -0.0837  -0.0688  -0.0623   3.5879  

Coefficients:
               Estimate Std. Error z value Pr(>|z|)    
(Intercept)   -6.494536   0.565005 -11.495   <2e-16 ***
length         0.002915   0.001227   2.375   0.0175 *  
differentTRUE -0.118038   0.671648  -0.176   0.8605    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 133.69  on 2946  degrees of freedom
Residual deviance: 128.91  on 2944  degrees of freedom
AIC: 134.91

Number of Fisher Scoring iterations: 8