tests/testthat/setup.R

mFile <- system.file("extdata/pwm1", package = "seqLogo")
dfDNA <- read.table(mFile)

## DNA
mDNA <- as.matrix(dfDNA)
rownames(mDNA) <- c("A", "C", "G", "T")
colnames(mDNA) <- seq_len(ncol(mDNA))

## RNA
mRNA <- mDNA
rownames(mRNA) <- c("A", "C", "G", "U")

## Amino Acid
mAA <- rbind(mDNA, matrix(0, nrow=16, ncol=ncol(mDNA)))
rownames(mAA) <- c("A", "C", "D", "E", "F", "G", "H", "I", "K", "L", 
                   "M", "N", "P", "Q", "R", "S", "T", "V", "W", "Y")

# default fills
fills = c(A = "#61D04F", C = "#2297E6", G = "#F5C710", T = "#DF536B")

Try the seqLogo package in your browser

Any scripts or data that you put into this service are public.

seqLogo documentation built on Nov. 8, 2020, 8:20 p.m.