超酷的ggstatsplot包

整理好的数据格式

ggbetweenstats(

  data  = all.tpm.clean,

  x    = E.type_0h,

  centrality.type = 'nonparametric',

  type = 'nonparametric',

  y    = log2_0h,

  point.args = list(position = ggplot2::position_jitterdodge(dodge.width = 0.2), alpha

                    = 2, size = 1, stroke =0),

  title = "0h enhancer associated gene",

  xlab = 'Enhacner type',

  ylab = 'Log2(average(TPM) + 1)'

)

## 4h

ggbetweenstats(

  data  = all.tpm.clean,

  x    = E.type_4h,

  centrality.type = 'nonparametric',

  type = 'nonparametric',

  y    = log2_4h,

  point.args = list(position = ggplot2::position_jitterdodge(dodge.width = 0.2), alpha

                    = 2, size = 1, stroke =0),

  title = "4h enhancer associated gene",

  xlab = 'Enhacner type',

  ylab = 'Log2(average(TPM) + 1)'

)