CNVR output file contents

157 views
Skip to first unread message

Melody Rynerson Palmer

unread,
Feb 25, 2016, 12:45:35 PM2/25/16
to ParseCNV
Hello,
Is there a simple reason why there would be no entries in the {prefix}_PlinkP_CNVR_ALL_ReviewdCNVRs.txt file, but a few hundred in the {prefix}_PlinkP_CNVR_ALL_ReviewdCNVRs_verbose.txt file? There are also entries in the _ControlEnriched.txt files, just not the main ones. 

Also I'm having trouble importing the verbose file into R or Excel accurately. Many of the columns are empty, so it's not assigning them to the correct header regardless of whether I tell it it's tab or space delimited, and sometimes it's hard to figure out which values match which column names. Any tips for that?

My final problem is that the UCSC custom track files are getting rejected with the following error: Error Unrecognized binary data format in file {prefix}_PlinkP_CNV_UCSC_bedfile_ForCustomTrackReview.txt
I'm still looking into this one but if anyone has suggestions or success I'd love to see an example.

Thank you!
Melody Palmer

Joe Glessner

unread,
Feb 25, 2016, 1:51:31 PM2/25/16
to Melody Rynerson Palmer, ParseCNV
Hello,

It sounds like a process control problem with the code. The code is exiting without the file being written fully. The verbose.txt file is the main one while the others are generated as the last step of InsertPlinkPvalues.pl. To the code attached, I added a close statement for each of the filehandles which should prevent the code from exiting without completing the file writing.

I added a looser error handling for unmet dependencies (.stats,ParseCNV.log,Plink.log) for InsertPlinkPvalues.pl which may have been a mistake because it results in blank fields which can be confusing. Check the .log file for “WARNING:” lines which will tell you what unmet dependencies exist and try copying those files into the working directory. The reports should be tab delimited, not space delimited.

Could you run these couple of tests on your UCSC track to help narrow down the problem?
awk '{print NF}' {prefix}_PlinkP_CNV_UCSC_bedfile_ForCustomTrackReview.txt | sort | uniq -c
awk '{print $1}' {prefix}_PlinkP_CNV_UCSC_bedfile_ForCustomTrackReview.txt | sort | uniq -c

Joe Glessner

--
You received this message because you are subscribed to the Google Groups "ParseCNV" group.
To unsubscribe from this group and stop receiving emails from it, send an email to parsecnv+u...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

InsertPlinkPvalues.pl.txt
Reply all
Reply to author
Forward
0 new messages