Advertisement
Guest User

Untitled

a guest
Apr 16th, 2014
56
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.75 KB | None | 0 0
  1. >>TTBK2_Hsap ,(CK1/TTBK)
  2. MSGGGEQLDILSVGILVKERWKVLRKIGGGGFGEIYDALDMLTRENVALKVESAQQPKQVLKMEVAVLKKLQGKDHVCRFIGCGRNDRFNYVVMQLQGRNLADLRRSQSRGTFT
  3.  
  4. >>TTBK2_Hsap ,(CK1/TTBK)
  5. MSGGGEQLDILSVGILVKERWKVLRKIGGGGFGEIYDALDMLTRENVALKVESAQQPKQVLKMEVAVLKKLQGKDHVCRFIGCGRNDRFNYVVMQLQGRNLADLRRSQSRGTFT
  6.  
  7. extractAPAAC(x, props = c("Hydrophobicity", "Hydrophilicity"), lambda = 30,
  8. w = 0.05, customprops = NULL)
  9.  
  10. read_and_extract <- function(fasta) {
  11. seq <- readFASTA(fasta)[[1]]
  12. return(extractAPAAC(seq, props = c("Hydrophobicity", "Hydrophilicity"), lambda = 30,
  13. w = 0.05, customprops = NULL))
  14.  
  15. }
  16.  
  17. setwd("H:\CC")
  18. fasta_files <- dir(pattern = "[.]fasta$")
  19. aa_comp <- vapply(fasta_files, read_and_extract, rep(pi, 80))
  20. write.csv(aa_comp, file = "C:\Users\PAAC.csv")
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement