Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- >>TTBK2_Hsap ,(CK1/TTBK)
- MSGGGEQLDILSVGILVKERWKVLRKIGGGGFGEIYDALDMLTRENVALKVESAQQPKQVLKMEVAVLKKLQGKDHVCRFIGCGRNDRFNYVVMQLQGRNLADLRRSQSRGTFT
- >>TTBK2_Hsap ,(CK1/TTBK)
- MSGGGEQLDILSVGILVKERWKVLRKIGGGGFGEIYDALDMLTRENVALKVESAQQPKQVLKMEVAVLKKLQGKDHVCRFIGCGRNDRFNYVVMQLQGRNLADLRRSQSRGTFT
- extractAPAAC(x, props = c("Hydrophobicity", "Hydrophilicity"), lambda = 30,
- w = 0.05, customprops = NULL)
- read_and_extract <- function(fasta) {
- seq <- readFASTA(fasta)[[1]]
- return(extractAPAAC(seq, props = c("Hydrophobicity", "Hydrophilicity"), lambda = 30,
- w = 0.05, customprops = NULL))
- }
- setwd("H:\CC")
- fasta_files <- dir(pattern = "[.]fasta$")
- aa_comp <- vapply(fasta_files, read_and_extract, rep(pi, 80))
- write.csv(aa_comp, file = "C:\Users\PAAC.csv")
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement