## ---- include = FALSE---------------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

## ----eval=FALSE---------------------------------------------------------------
#  input <- read.table("methyrate", header = T)
#  e_value <- varevalue.single_general(methyrate=input, group1_name='g1', group2_name='g2', chr='chr21', start=9439679, end=9439679)
#  head(e_value)

## ----eval=FALSE---------------------------------------------------------------
#  input <- read.table("desq_out", header = T)
#  data_e <- metevalue.RNA_general(input, group1_name='treated', group2_name='untreated')
#  head(data_e)