# NOT RUN {
#File path for output files
wd <- tempdir()
#wd <-getwd()
file <- file.path(wd,"test_seq2")
#Importing MSA file
msf <- system.file("msa/toy_align.msf", package = "Bios2cor")
align <- import.msf(msf)
#Creating OMES correlation object
omes <- omes(align, gap_ratio = 0.2)
#Creating entropy object
#entropy <- entropy(align)
#Selecting correlation matrix
omes <-omes$Zscore
#Creating PCA object for selected correlation matrix and saving eigen values
pca <- centered_pca(omes, filepathroot= file, pc= NULL, dec_val= 5, filter= NULL)
#Creating 2D plot of elements projected on selected [abs,ord] plane
pca_2d(pca, abs = 1, ord = 3, file)
#Creating 2D plot of elements projected on default [1,2] plane
#pca_2d(pca)
# }
Run the code above in your browser using DataLab