Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • get average gene length file

    Hi everyone,

    I have a very quick question...probably super easy for experienced programmers!
    I need to generate an average gene length file from the .gtf file I downloaded from Ensembl. In detail, I will need a 2 column file with ENSG id on the first and gene length on the second column (I need for HTseq normalization...)

    Can anyone help me with scripting?

    Thanks!!!
    Manu

  • #2
    In R:
    Code:
    library(GenomicFeatures)
    txdb <- makeTranscriptDbFromGFF("foo.gtf", type="gtf")
    trans <- transcripts(txdb, columns=c("GENEID"))
    df <- data.frame(gene=trans$GENEID, len=width(trans))
    You probably need to specify something else for "GENEID". You can see available options with columns(txdb).

    Comment


    • #3
      Thank you dpryan!

      I get this error message
      txdb <- makeTranscriptDbFromGFF("foo.gtf", type="gtf")
      Error in makeTranscriptDbFromGFF("foo.gtf", type = "gtf") :
      unused argument (type = "gtf")

      what is it?

      Comment


      • #4
        It should have been format= rather than type=. Mea culpa.

        Comment


        • #5

          thank you!
          it's running now...i'll update the result soon!!!

          Comment


          • #6
            Hi dpryan!


            after 10 minutes of running this command:
            Originally posted by dpryan View Post
            In R:
            Code:
            df <- data.frame(gene=trans$GENEID, len=width(trans))
            R Studio simply crashes... any idea why???
            It takes forever to run txdb <- makeTranscriptDbFromGFF("genes.gtf", format="gtf"), but it ultimately does it in 15 min or so...

            Thanks!
            Manu

            Comment


            • #7
              Presumably you're running out of memory and that's causing the crash. R isn't always the best when it comes to memory. You might quickly do:
              Code:
              head(trans$GENEID)
              head(width(trans))
              just to ensure that there's nothing strange that happened while making the "trans" object. The odds of an error there are crazy low, but it couldn't hurt to double check.

              BTW, the txdb object can be saved to a file, should you ever need it again. Just use something like
              Code:
              saveDb(txdb, file="myOrganism.sqlite")
              and you can simply load it again later with
              Code:
              loadDb("myOrganism.sqlite")
              This ends up saving a lot of time if you work with the same organism often.

              Comment


              • #8
                Originally posted by dpryan View Post
                Presumably you're running out of memory and that's causing the crash. R isn't always the best when it comes to memory. You might quickly do:
                Code:
                head(trans$GENEID)
                head(width(trans))
                just to ensure that there's nothing strange that happened while making the "trans" object. The odds of an error there are crazy low, but it couldn't hurt to double check.

                BTW, the txdb object can be saved to a file, should you ever need it again. Just use something like
                Code:
                saveDb(txdb, file="myOrganism.sqlite")
                and you can simply load it again later with
                Code:
                loadDb("myOrganism.sqlite")
                This ends up saving a lot of time if you work with the same organism often.
                Thanks dpryan!
                as a matter of fact my computer was a bit cranky...it worked perfect after i rebooted it!!!
                Manu

                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, 03-27-2024, 06:37 PM
                0 responses
                12 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-27-2024, 06:07 PM
                0 responses
                11 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-22-2024, 10:03 AM
                0 responses
                53 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-21-2024, 07:32 AM
                0 responses
                69 views
                0 likes
                Last Post seqadmin  
                Working...
                X