inst/script/use_esophageal.R

In HiLDA: Conducting statistical inference on comparing the mutational exposures of mutational signatures by using hierarchical latent Dirichlet allocation

## load the library
library(HiLDA)

## read in the data
inputFile <- system.file("extdata/esophageal.mp.txt.gz", package="HiLDA")

## perform Shiraishi et al's method, pmsignature, to extract signatures
G <- hildaReadMPFile(inputFile, numBases=5, trDir=TRUE)