polyA usage changes between groups?

WangJL2021 / scPolyA-seq

https://pubmed.ncbi.nlm.nih.gov/36454750/

2 stars 1 forks source link

polyA usage changes between groups? #2

Open DawnEve opened 1 year ago

DawnEve commented 1 year ago

set.seed(100)
getN=function(n=1000){
  #t1=round(rnorm(n, 100, 100))
  t1=round(runif(n, 0, 1000))
  abs(t1)
}

a1=getN()
a2=getN()
#
b1=getN()
b2=getN()

hist(a1/a2)
hist(b1/b2)

hist( log2((a1/a2) /(b1/b2)) )
plot( log2(a1/a2), log2(b1/b2), pch=19, col="red", cex=0.5)

APA_changes_ratio_group

DawnEve commented 1 year ago

https://pubmed.ncbi.nlm.nih.gov/30143597/ A compendium of conserved cleavage and polyadenylation events in mammalian genes Fig6

RED = Δlog2(distal PAS/proximal PAS), testis or brain versus heart