Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • RAPIDR package for the NIPT

    Hello,

    I am a young student. I would like to know where we can integrate the RAPIDR package in a workflow for the NIPT.
    Does someone use this package ?

    This package uses the human genome reference hg19. Can we use RAPIDR with an other ?

    What can be a substitute to RAPIDR in the workflow for NIPT ?

    Thanking you in advance.

  • #3
    Rapidr

    Hi Laquais,
    Have you tried RAPIDR?
    I have some doubts regarding the the data 'outcomes'
    Thanks in advance.

    Comment


    • #4
      Hey Manonathan,

      Yes, I integrated RAPIDR in a pipeline and I compared its results with a homemade pipeline. The results were convincing.
      Today, the RAPIDR pipeline is used routinely to detect T21, T18, T13. Last I noticed, the results are reliable (I have now a new job but I continue to discuss with my old colleagues).

      What are your doubts ?

      Best regards.

      Comment


      • #5
        Hi Laquais,

        Thanks for your reply.
        I have finished makeBinnedCountsFile(c("file1.bam", "file2.bam"), c("sample1", "sample2"), 'output.fname',mask=NULL, k = 20000)
        I got an output file 'output.fname'.
        file1.bam and file2.bam are two bam files corresponding to aneuploidy samples.
        I don't understand 'outcomes' dataframe.
        I couldn't run createReferenceSetFromCounts('87.fname','a',method='zscore',gcCorrect = FALSE,gcContentFile = 'gcContent.fname')

        I am getting error like this
        "Loading binned counts file
        Checking every sampleID has an outcome
        Error in outcomes$SampleID : $ operator is invalid for atomic vectors"

        Thanks and Regards,
        Manonathan

        Comment


        • #6
          RAPIDR transforms bam files into binned counts files to compare binned counts files with a reference set.


          Reference set

          That's why you have to gather some diploidy samples (more you have better is) and to use makeBinnedCountsFile with all your samples of the reference set. You can also add some aneuploidy samples (optionnal) in your reference set to have positive samples in your reference set (that's why you have a column Dx in outcomes). makeBinnedCountsFile with all your samples of the reference set will create an referenceSet.fname that you will use with createReferenceSetFromCounts.
          You have to make an outcomes dataframe to define what is the sexe of the sample and if the sample is an aneuploidy sample (RAPIDR uses outcomes to calculate its zscores).
          I advise you to create a file outcomes.txt and convert it in a dataframe with a R script.


          Unknown samples

          After that, you use makeBinnedCountsFile with your unknown samples (you don't know if these samples are diploide or not) to create unknown.fname
          You use testUnknowns(ref.set, unknown.fname) to have the results.


          "Error in outcomes$SampleID : $ operator is invalid for atomic vectors"
          ==> I think your column sampleID of outcomes is wrong. The SampleID must match with the sampleId you wrote in makeBinnedCountsFile

          If you have :
          makeBinnedCountsFile(c("file1.bam", "file2.bam",...), c("sample1", "sample2",...), 'reference.fname')

          outcomes must be :
          SampleID Dx Gender
          1 sample1 Normal Female
          2 sample2 Normal Female
          3 sample3 T21 Female
          ......
          ......


          I hope I am understable... My english is not very good. ^^

          Comment


          • #7
            Rapidr

            Hi Laquais,

            Thanks for your reply.
            Now I got some idea about how to proceed with this package.
            Thanks a lot. I have been struggling since December'15 to get output.
            Thanks and Regards,
            Manonathan

            Comment


            • #8
              You're welcome,

              I hope you will manage to obtain your results

              Best regards.

              Comment


              • #9
                HI Manonathan,
                I am using the RAPIDR, could you send me your complete scripts to perform it ?, from bam files to get the output ..... I am a new guy about NIPT
                My email is : [email protected]
                Thank you very much!!

                Comment


                • #10
                  Hi, I know that this is an old post, but I really need your help! I'm trying to use RAPIDR, but I have some problems with the reference creation. I explained my problems here: http://seqanswers.com/forums/showthr...ghlight=rapidr, can you help me? It is very important, I don't found anything on the web!

                  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
                  12 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
                  52 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