# NOT RUN {
#File path for output files
wd <- tempdir()
#wd <-getwd()
file <- file.path(wd,"test_seq3")
#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)
#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)
#Plotting scree plot
pca_screeplot(pca, file)
#pca_screeplot(pca)
# }
Run the code above in your browser using DataLab