Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Problem with the tol CoNiFer

    Hi,

    Does anyone here has ever worked with conifer.
    I have this strange kind of error.
    I am running 240 samples but it shows some problems with probes on cgr24.
    I used grep -v '24' from the probes.txt file as it says to do so but it was of no use.
    Before this, I used conifer on ~500 samples and it's working perfectly fine.
    Can anyone give me some idea how to fix this issue.
    Thanks in advance
    Vishal


    Error while analysing:
    [RUNNING: chr23] Masking 92 probes with median RPKM < 1.000000
    [RUNNING: chr23] Calculating ZRPKM scores...
    [RUNNING: chr23] SVD decomposition...
    [RUNNING: chr23] Saving SVD-ZRPKM values
    [RUNNING: chr24] Now on: chrY
    [RUNNING: chr24] Found 465 probes; probeID range is [193615-194080]
    [RUNNING: chr24] Calculating median RPKM
    [RUNNING: chr24] Masking 432 probes with median RPKM < 1.000000
    [ERROR] This chromosome has fewer informative probes than there are samples in the analysis! There are probably no mappings on this chromosome. Please remove these probes from the probes.txt file
    Closing remaining open files: /vishal/conifer/ige/output/analysis.hd5... done

    Error while calling :
    [INIT] Initializing caller at threshold = 1.500000
    Traceback (most recent call last):
    File "conifer.py", line 682, in <module>
    args.func(args)
    File "conifer.py", line 338, in CF_call
    r = cf.rpkm_reader(h5file_in_fn)
    File "/sw/conifer_v0.2.2/conifer_functions.py", line 256, in __init__
    self.sample_table = self.h5file.root.samples.samples
    File "/opt/rrzk/software/python/python-2.7.5/lib/python2.7/site-packages/tables/group.py", line 813, in __getattr__
    return self._f_get_child(name)
    File "/opt/rrzk/software/python/python-2.7.5/lib/python2.7/site-packages/tables/group.py", line 683, in _f_get_child
    self._g_check_has_child(childname)
    File "/opt/rrzk/software/python/python-2.7.5/lib/python2.7/site-packages/tables/group.py", line 407, in _g_check_has_child
    % (self._v_pathname, name))
    tables.exceptions.NoSuchNodeError: group ``/`` does not have a child named ``samples``
    Closing remaining open files: /scratch/ccg-ngs/tmp/vishal/conifer/ige/output/analysis.hd5... done

  • #2
    no one who has worked with conifer?

    Comment


    • #3
      Hi Vishal,

      I am facing the similar problem.
      Were you able to fix this issue?

      Thanks, Braj

      Comment


      • #4
        Hi Braj,

        I think that there is a bug in the tool, I wrote to the author.

        Go to conifer folder and edit conifer.py
        Go to line 153, 154 and 155 and comment them (put # in front of them) and then save it and run conifer.

        #if num_chr_probes <= len(samples):
        #print "[ERROR] This chromosome has fewer informative probes than there are samples in the analysis! There are probably no mappings on this chromosome. Please remove these probes from the probes.txt file"
        # sys.exit(0)


        Write to me again if it doesn't work.
        Regards
        Vishal
        Last edited by vishal.rossi; 02-04-2014, 01:57 AM.

        Comment


        • #5
          Hi Vishal,

          Thanks for your reply.

          I tried the way you suggested, but this time I am getting some new error at analysis step. The error which I am getting this time is :

          .
          .
          .
          [RUNNING: chr22] Now on: chr22
          [RUNNING: chr22] Found 3961 probes; probeID range is [114481-118442]
          [RUNNING: chr22] Calculating median RPKM
          [RUNNING: chr22] Masking 49 probes with median RPKM < 1.000000
          [RUNNING: chr22] Calculating ZRPKM scores...
          [RUNNING: chr22] SVD decomposition...
          [RUNNING: chr22] Saving SVD-ZRPKM values
          [RUNNING: chr23] Now on: chrX
          [RUNNING: chr23] Found 6157 probes; probeID range is [179253-185410]
          [RUNNING: chr23] Calculating median RPKM
          [RUNNING: chr23] Masking 87 probes with median RPKM < 1.000000
          [RUNNING: chr23] Calculating ZRPKM scores...
          [RUNNING: chr23] SVD decomposition...
          [RUNNING: chr23] Saving SVD-ZRPKM values
          [RUNNING: chr24] Now on: chrY
          [RUNNING: chr24] Found 226 probes; probeID range is [185410-185636]
          [RUNNING: chr24] Calculating median RPKM
          [RUNNING: chr24] Masking 221 probes with median RPKM < 1.000000
          [RUNNING: chr24] Calculating ZRPKM scores...
          [RUNNING: chr24] SVD decomposition...
          Traceback (most recent call last):
          File "../conifer_v0.2.2/conifer.py", line 682, in <module>
          args.func(args)
          File "../conifer_v0.2.2/conifer.py", line 189, in CF_analyze
          rpkm = np.dot(U, np.dot(new_S, Vt))
          ValueError: objects are not aligned
          Closing remaining open files: analysis.hdf5... done

          -----------------------------
          I want to give a try to troubleshoot the issue...Let me look at the script, if I could understand.

          Please do share in case you got it fixed, it will be a great help.

          Thanks & Regards,

          Braj

          Comment


          • #6
            Hi Vishal,

            I tried to understand this issue what I understand is :

            Since you have total 240 samples.

            Now if you see the error which you get

            .
            .
            .
            [RUNNING: chr24] Found 465 probes; probeID range is [193615-194080]
            [RUNNING: chr24] Calculating median RPKM
            [RUNNING: chr24] Masking 432 probes with median RPKM < 1.000000
            [ERROR] This chromosome has fewer informative probes than there are samples in the analysis! There are probably no mappings on this chromosome. Please remove these probes from the probes.txt file
            Closing remaining open files: /vishal/conifer/ige/output/analysis.hd5... done
            -----------------------------
            In your "Y' chromosome out of 465 probes 432 are masked and remaining is only 33 probes. Since this probe count i.e 33 is lower than the number of total samples i.e 240 hence statistical calculation can't be performed uniformly for this chromosome "Y'.

            Therefore no need to comment the lines in conifer.py, it is correct only according to the statistical requirement.

            In my probes.txt file , I removed the targets related to chr 'Y' and repeating my analysis. Hope this will be the correct approach according to the error obtained.

            Let me know if you don't agree with this explanation.

            Thanks & Regards

            Braj

            Comment


            • #7
              Hi,

              I removed the probes and tried to run conifer but it didn't go through. Maybe, split your bam file by chromosomes, remove the Y chromosome alignment and merge the files again and try it.
              I didn't do this step, as I run some simulations and found that I am getting the same results even after making comments in the code.


              [And did you comment out the 155th line.... # sys.exit(0)]

              Regards,
              Vishal
              Last edited by vishal.rossi; 02-04-2014, 05:23 AM.

              Comment


              • #8
                Hi Vishal,

                Once you remove your target coordinates related to 'Y' chromosome from "probes.txt", then rerun the first step i.e RPKM calculation, and follow the next steps as instructed in the tutorial. You don't need to split your bam file, as tool will predict RPKM only for targets in the probes.txt.

                No, I didn't commented any line in the script "conifer.py'

                I got my result by doing so.


                Regards, Braj
                Last edited by Brajbio; 02-05-2014, 07:14 AM.

                Comment


                • #9
                  Hey Vishal & Braj,

                  Looks like you guys are experts in using conifer by now.

                  I am newbie to bioinfo as well to using NGS & conifer.

                  Some how, i am learning & working on the stuff.

                  I was able to succesfully run the commands to the sample data of conifer.
                  But, when I give commands to my data, errors come.

                  What is the first step in using conifer analysis.
                  Hope it's creating RPKM files.
                  I am not able to that properly.

                  Here is my commands & error.

                  python conifer.py rpkm --probes sampledata/probes.txt --input AD1cord-rmdup-realigned.bam --output RPKM/AD1.rpkm.txt

                  [INIT] Successfully read in 194080 probes from sampledata/probes.txt
                  [ERROR] Cannot open rpkm file for writing: ['../../odity/IndelTargetCreator/RealignedbamFiles/RPKM/AD1.rpkm.txt']

                  Help me with proper commands & solving the issue.

                  Comment


                  • #10
                    Hi I also am running into problems getting Conifer up and running. I have successfully run through the quick start so I believe all is correctly installed. I also generated 10 rpkm files from bam files utilizing the conifer.py rpkm command. When I run the confier.py analyze command I get the following error:
                    [RUNNING: chr1] Now on: chr1
                    [RUNNING: chr1] Found 45973 probes; probeID range is [0-224692]
                    [RUNNING: chr1] Calculating median RPKM
                    [RUNNING: chr1] Masking 3125 probes with median RPKM < 1.000000
                    Traceback (most recent call last):
                    File "/usr/local/conifer_v0.2.2/conifer.py", line 682, in <module>
                    args.func(args)
                    File "/usr/local/conifer_v0.2.2/conifer.py", line 157, in CF_analyze
                    probeIDs = np.array(map(operator.itemgetter("probeID"),chr_probes))[probe_mask]
                    ValueError: too many boolean indices
                    Closing remaining open files:analysis.hdf5...done
                    My call is: python conifer.py analyze --probes AnnotatedSeqCap_EZ_Exome_v2nochr.bed --rpkm_dir RPKM/ --output analysis.hdf5 --svd 2 --write_svals singular_values.txt

                    Any suggestions how I can get this off the ground? Thanks

                    Comment


                    • #11
                      Hi,

                      This error might be due to old numpy version. Can you check which numpy version is installed?
                      Make sure it is numpy1.8 or above.

                      Vishal

                      Comment


                      • #12
                        Thank you Vishal. I am running numpy 1.8. I updated to 1.9 and this still resulted in the same error.
                        :
                        [RUNNING: chr1] Now on: chr1
                        [RUNNING: chr1] Found 45973 probes; probeID range is [0-224692]
                        [RUNNING: chr1] Calculating median RPKM
                        [RUNNING: chr1] Masking 3125 probes with median RPKM < 1.000000
                        Traceback (most recent call last):
                        File "/usr/local/conifer_v0.2.2/conifer.py", line 682, in <module>
                        args.func(args)
                        File "/usr/local/conifer_v0.2.2/conifer.py", line 157, in CF_analyze
                        probeIDs = np.array(map(operator.itemgetter("probeID"),chr_probes))[probe_mask]
                        IndexError: index 45973 is out of bounds for axis 1 with size 45973
                        Closing remaining open files:analysis.hdf5...done
                        Last edited by wdemos; 10-31-2014, 11:06 AM. Reason: additional info

                        Comment


                        • #13
                          Originally posted by wdemos View Post
                          Thank you Vishal. I am running numpy 1.8. I updated to 1.9 and this still resulted in the same error.
                          :
                          [RUNNING: chr1] Now on: chr1
                          [RUNNING: chr1] Found 45973 probes; probeID range is [0-224692]
                          [RUNNING: chr1] Calculating median RPKM
                          [RUNNING: chr1] Masking 3125 probes with median RPKM < 1.000000
                          Traceback (most recent call last):
                          File "/usr/local/conifer_v0.2.2/conifer.py", line 682, in <module>
                          args.func(args)
                          File "/usr/local/conifer_v0.2.2/conifer.py", line 157, in CF_analyze
                          probeIDs = np.array(map(operator.itemgetter("probeID"),chr_probes))[probe_mask]
                          IndexError: index 45973 is out of bounds for axis 1 with size 45973
                          Closing remaining open files:analysis.hdf5...done


                          Was this issue ever solved? I've gotten this error too.

                          Comment


                          • #14
                            No, I have not resolved this issue

                            Comment


                            • #15
                              Hello,

                              I am trying to run Conifer on my exome seq data and I am not able to get the step that calculates RPKM values to work. An rpkm.txt file gets created but it is empty. Does anyone have an idea about what is going wrong ?

                              Comment

                              Latest Articles

                              Collapse

                              • 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 on Modified Bases...
                                Yesterday, 07:01 AM
                              • seqadmin
                                Current Approaches to Protein Sequencing
                                by seqadmin


                                Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
                                04-04-2024, 04:25 PM

                              ad_right_rmr

                              Collapse

                              News

                              Collapse

                              Topics Statistics Last Post
                              Started by seqadmin, 04-11-2024, 12:08 PM
                              0 responses
                              55 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 04-10-2024, 10:19 PM
                              0 responses
                              52 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 04-10-2024, 09:21 AM
                              0 responses
                              45 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 04-04-2024, 09:00 AM
                              0 responses
                              55 views
                              0 likes
                              Last Post seqadmin  
                              Working...
                              X