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
                    Strategies for Sequencing Challenging Samples
                    by seqadmin


                    Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                    03-22-2024, 06:39 AM
                  • seqadmin
                    Techniques and Challenges in Conservation Genomics
                    by seqadmin



                    The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                    Avian Conservation
                    Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                    03-08-2024, 10:41 AM

                  ad_right_rmr

                  Collapse

                  News

                  Collapse

                  Topics Statistics Last Post
                  Started by seqadmin, Yesterday, 06:37 PM
                  0 responses
                  11 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, Yesterday, 06:07 PM
                  0 responses
                  10 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 03-22-2024, 10:03 AM
                  0 responses
                  51 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 03-21-2024, 07:32 AM
                  0 responses
                  68 views
                  0 likes
                  Last Post seqadmin  
                  Working...
                  X