Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • angel-sakura
    Member
    • Nov 2014
    • 23

    #1

    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~
  • SylvainL
    Senior Member
    • Feb 2012
    • 180

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

    Comment

    • angel-sakura
      Member
      • Nov 2014
      • 23

      #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

      • SylvainL
        Senior Member
        • Feb 2012
        • 180

        #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

        • blancha
          Senior Member
          • May 2013
          • 367

          #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

          • angel-sakura
            Member
            • Nov 2014
            • 23

            #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

            • angel-sakura
              Member
              • Nov 2014
              • 23

              #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

              • SylvainL
                Senior Member
                • Feb 2012
                • 180

                #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

                • angel-sakura
                  Member
                  • Nov 2014
                  • 23

                  #9
                  Yes! I understand it now!! lol

                  Thank you so much!

                  Comment

                  Latest Articles

                  Collapse

                  • SEQadmin2
                    Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
                    by SEQadmin2



                    Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

                    Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
                    ...
                    Today, 09:55 AM
                  • SEQadmin2
                    New Genomics Technologies Take Aim at Long-Standing Limits
                    by SEQadmin2


                    Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

                    We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
                    ...
                    09-28-2026, 10:25 AM

                  ad_right_rmr

                  Collapse

                  News

                  Collapse

                  Topics Statistics Last Post
                  Started by SEQadmin2, Today, 05:39 AM
                  0 responses
                  10 views
                  0 reactions
                  Last Post SEQadmin2  
                  Started by SEQadmin2, Yesterday, 09:56 AM
                  0 responses
                  13 views
                  0 reactions
                  Last Post SEQadmin2  
                  Started by SEQadmin2, 10-09-2026, 11:14 AM
                  0 responses
                  26 views
                  0 reactions
                  Last Post SEQadmin2  
                  Started by SEQadmin2, 09-29-2026, 09:51 AM
                  0 responses
                  77 views
                  0 reactions
                  Last Post SEQadmin2  
                  Working...