Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Help with RAPIDR package

    Hi all! I need your help with an R package named RAPIDR. I'm new to use this package, but I really need it, but I don't know how I can perform my analysis without errors. Now I explain my steps.

    1) I create a binned file with the following code:
    Code:
    > makeBinnedCountsFile(bam.file.list=c(bam1.bam, bam2.bam,bam3.bam,bam4.bam), sampleIDs=c("001", "003", "005", "013"), binned.counts.fname="output.csv", k=20000)
    2) I create a reference with the following code:
    Code:
    > rapidr.dir<-system.file(package="RAPIDR")
    > data(outcomes)
    > data(gcContent)
    > T21.pos<-which(outcomes$Dx=="T21")
    > chr.lens<-sapply(gcContent, length)
    > chr.names<-names(chr.lens)
    > header<-c("SampleID")
    > for(i in 1:length(chr.lens)){
    + header<-c(header,rep(chr.names[i], chr.lens[i]))}
    > nbins<-sum(chr.lens)
    > ncols<-nbins+1
    > binned.counts<-matrix(nrow=nrow(outcomes), ncol=ncols)
    > for(i in 1:nrow(binned.counts)){
    + binned.counts[i,]<-rpois(ncols, lambda=100)
    + if(i%in%T21.pos){
    + binned.counts[i,139087:141493]<-rpois(chr.lens[21], lambda=115)
    + }}
    > binned.counts[,1]<-outcomes$SampleID
    > colnames(binned.counts)<-header
    > t<-tempfile()
    > write.table(binned.counts, file=t, col.names=TRUE, row.names=FALSE, quote=FALSE, sep=",")
    > "output.csv"<-t
    > message(t)
    /tmp/RtmpwrkRt5/file33292c2c486
    > gcContent.fname<-paste(rapidr.dir, "/media/martina/VERBATIM HD/Altamedica/RAPIDR_0.1.1/RAPIDR/data/gcContent.rda", sep="")
    > head(outcomes)
      SampleID     Dx Gender
    1     1000 Normal Female
    2     1001 Normal Female
    3     1002    T21 Female
    4     1003 Normal   Male
    5     1004 Normal   Male
    6     1005 Normal   Male
    > ref.set<-createReferenceSetFromCounts("/media/martina/VERBATIM HD/risultati/output.csv", outcomes, gcCorrect=FALSE, PCA=FALSE, filterBin=FALSE, gcContentFile=gcContent.fname)
    But I receive this error:
    Code:
    Loading binned counts file
    Checking every sampleID has an outcome
    No outcomes for Sample 001
    No outcomes for Sample 003
    No outcomes for Sample 005
    No outcomes for Sample 013
    Error in `[.data.frame`(sampleIDs.with.outcomes, , "Gender") : 
      undefined columns selected

    The BAM files is from 2 male and 2 female, without aneuplodies and they come from ION Torrent, but I used the fastq file and re-alignd with bowtie2, and coverted into sorted-bam file with samtools.


    Why I receive this error? How can I perform the analysis? Is there anyone who have some experience with this package and can help me?

    Thank you
    Last edited by PandoraMid; 11-21-2017, 10:37 AM. Reason: I forgot some essential information

  • #2
    Sample Id in outcomes data frame should correspond to the sample file name . The reason for the error may be the change in sample Id .check whether both are having same name. If not, modifythe sample id or rename the sample file and do the process again.

    Comment


    • #3
      Hi! Thanks for your reply! I have a doubt: the only ID for sample that I create is in the creation of binned file. Where I can do this check? Is the second part of my code correct? Can I see all 3 aneuploidies with this code?

      I'm sorry for my question, but I'm new with this kind of analysis!

      Comment


      • #4
        Hi,
        No issue. Me too struggled a lot to understand the workflow.

        create a new text file with three columns like that in the outcomes data frame
        Use sampleid, Gender and Dx.Use your sample files id, gender and Dx details and save it as "outcome.txt".Use the Female/Male term as it is case sensitive.
        convert that into a data frame using
        outcome<-read.table("outcome.txt",header=T)

        head(outcome)

        should be similar to the outcomes dataframe.

        use this newly created dataframe in your script instead of outcomes.

        Comment


        • #5
          So, if I understood I have to create manually the outcomes, and then use this file with the function
          Code:
          createReferenceSetFromCounts
          with: my binned file, my outcome file (made by output.txt and your instuction), gcCorrect and the others things?

          Second question: is it ok if I create my outcomes file with csv extension?

          Third question: do you use RAPIDR-Plus too?

          Comment


          • #6
            Hi! I tryed to use your tips but I can't resolve the situation. This is my outcome.txt file:

            sampleID Gender Dx
            001 Female Normal
            003 Female Normal
            005 Male Normal
            013 Male Normal

            Here there is the code that I use

            Code:
            > makeBinnedCountsFile(bam.file.list=c(bam1, bam2, bam3, bam4), sampleIDs=c("1", "3", "5", "13"), binned.counts.fname="/media/martina/VERBATIM HD/Altamedica/risultati/output_uffa.csv", k=20000)
            Binning counts in bam files
            doing the binning
            Binning done in 42.556
            doing the binning
            Binning done in 31.375
            doing the binning
            Binning done in 49.443
            doing the binning
            Binning done in 35.978
            > outcome<-read.table("/media/martina/VERBATIM HD/Altamedica/risultati/outcomes.txt", header=T)
            > head(outcome)
              sampleID Gender     Dx
            1        1 Female Normal
            2        3 Female Normal
            3        5   Male Normal
            4       13   Male Normal
            > ref.set<-createReferenceSetFromCounts(output.binned, outcome)
            Loading binned counts file
            Checking every sampleID has an outcome
            No outcomes for Sample 1
            No outcomes for Sample 3
            No outcomes for Sample 5
            No outcomes for Sample 13
            Error in `[.data.frame`(sampleIDs.with.outcomes, , "Gender") : 
              undefined columns selected
            I tryed to use "sampleid" "SampleID", "sampleid" but with the same results. Make I some errors?

            Comment


            • #7
              I am getting the same error as above. Please anyone can help me regarding this issue

              Error in `[.data.frame`(sampleIDs.with.outcomes, , "Gender") :
              undefined columns selected
              Last edited by pooja.solanki2018; 05-23-2018, 09:17 PM.

              Comment

              Latest Articles

              Collapse

              • seqadmin
                Recent Advances in Sequencing Analysis Tools
                by seqadmin


                The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
                05-06-2024, 07:48 AM
              • 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

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, 05-14-2024, 07:03 AM
              0 responses
              19 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 05-10-2024, 06:35 AM
              0 responses
              44 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 05-09-2024, 02:46 PM
              0 responses
              54 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 05-07-2024, 06:57 AM
              0 responses
              42 views
              0 likes
              Last Post seqadmin  
              Working...
              X