Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • How can I extract sequence name from fast file in R

    Here is my code:
    Code:
    setwd("/Users/xshang/Dropbox/project-ing/peptidome")
    library("Biostrings")
    human_protein <- readAAStringSet("HUMAN.fasta",use.names=TRUE)
    When I use
    Code:
    names(human_protein[[1]])
    , it returns "NULL". How can I extract sequence name from fast file in R?

    Thanks!

  • #2
    Try
    Code:
    names(human_protein)[1]
    The names() function is an accessor for the names in the AAStringSet object and returns a character vector. Alternatively, you can also use
    Code:
    names(human_protein[1])
    The [ function also returns an AAStringSet, rather than an AAString (returned by [[).

    Comment


    • #3
      Originally posted by dpryan View Post
      Try
      Code:
      names(human_protein)[1]
      The names() function is an accessor for the names in the AAStringSet object and returns a character vector. Alternatively, you can also use
      Code:
      names(human_protein[1])
      The [ function also returns an AAStringSet, rather than an AAString (returned by [[).
      It works! Many thanks!

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Exploring the Dynamics of the Tumor Microenvironment
        by seqadmin




        The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
        07-08-2024, 03:19 PM
      • seqadmin
        Exploring Human Diversity Through Large-Scale Omics
        by seqadmin


        In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
        06-25-2024, 06:43 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, 07-19-2024, 07:20 AM
      0 responses
      28 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 07-16-2024, 05:49 AM
      0 responses
      41 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 07-15-2024, 06:53 AM
      0 responses
      46 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 07-10-2024, 07:30 AM
      0 responses
      43 views
      0 likes
      Last Post seqadmin  
      Working...
      X