Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • problem with biomart---give 0 results

    Hi everyone:

    I'm trying to convert ensembl-peptide-id to GO-id with biomart. However, when I was using the following command to run it, I got 0 results.

    Here are commands I was using:
    ensembl = useMart("ensembl",dataset="hsapiens_gene_ensembl")
    id_list=c("ENSP00000298510","ENSP00000474109")
    trans_id=getBM(attributes=c("ensembl_peptide_id","entrezgene","go_id","name_1006","definition_1006","namespace_1003"), filter="ensembl_peptide_id", values="id_list", mart=ensembl)

    Could someone tell me what cause the problem please?

    Thanks a lot~

  • #2
    In my code, I use filters instead of filter, don't know which version you are using....

    Comment


    • #3
      Originally posted by SylvainL View Post
      In my code, I use filters instead of filter, don't know which version you are using....
      Thank you very much for your reply.
      I'm pretty sure it is not the problem of typo, as it was working with the example data with the same commands.

      Comment


      • #4
        Ah...

        did you update your biomaRt package? Then you have to check whether "ensembl_peptide_id" is still the proper filter name... I had this problem once with biomaRt, with "ensembl_gene_id"

        Comment


        • #5
          Remove the commas around id_list.

          Incorrect.
          trans_id=getBM(attributes=c("ensembl_peptide_id","entrezgene","go_id","name_1006","definition_1006","namespace_1003"), filter="ensembl_peptide_id", values="id_list", mart=ensembl)

          Correct.
          trans_id=getBM(attributes=c("ensembl_peptide_id","entrezgene","go_id","name_1006","definition_1006","namespace_1003"), filter="ensembl_peptide_id", values=id_list, mart=ensembl)

          Comment


          • #6
            Originally posted by blancha View Post
            Remove the commas around id_list.

            Incorrect.
            trans_id=getBM(attributes=c("ensembl_peptide_id","entrezgene","go_id","name_1006","definition_1006","namespace_1003"), filter="ensembl_peptide_id", values="id_list", mart=ensembl)

            Correct.
            trans_id=getBM(attributes=c("ensembl_peptide_id","entrezgene","go_id","name_1006","definition_1006","namespace_1003"), filter="ensembl_peptide_id", values=id_list, mart=ensembl)
            Thank you very much!!!
            It is working now!!

            Comment


            • #7
              Originally posted by SylvainL View Post
              Ah...

              did you update your biomaRt package? Then you have to check whether "ensembl_peptide_id" is still the proper filter name... I had this problem once with biomaRt, with "ensembl_gene_id"
              The problem is solved!!
              The "value" doesn't need the quotation mark!

              Comment


              • #8
                Can't believe I didn't see it...

                So for your future scripts:
                "object" means it's not a variable but a value
                object (without ") refers to an existing object.

                In your case, it was looking to "id_list" as an "ensembl_peptide_id", meaning if you replace "id_list" by "ENSP00000298510" in your first command, it will give you the results for this protein...

                Hope I'm clear enough

                Comment


                • #9
                  Yes! I understand it now!! lol

                  Thank you so much!

                  Comment

                  Latest Articles

                  Collapse

                  • seqadmin
                    Essential Discoveries and Tools in Epitranscriptomics
                    by seqadmin




                    The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
                    04-22-2024, 07:01 AM
                  • seqadmin
                    Current Approaches to Protein Sequencing
                    by seqadmin


                    Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
                    04-04-2024, 04:25 PM

                  ad_right_rmr

                  Collapse

                  News

                  Collapse

                  Topics Statistics Last Post
                  Started by seqadmin, Today, 08:47 AM
                  0 responses
                  12 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 04-11-2024, 12:08 PM
                  0 responses
                  60 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 04-10-2024, 10:19 PM
                  0 responses
                  59 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 04-10-2024, 09:21 AM
                  0 responses
                  54 views
                  0 likes
                  Last Post seqadmin  
                  Working...
                  X