Code Monkey home page Code Monkey logo

Comments (12)

leepc12 avatar leepc12 commented on July 21, 2024

from chipseq_pipeline.

kdemuren avatar kdemuren commented on July 21, 2024

I am also having this problem - the pipeline stops at the jsd step. I've done the above to no avail. Here's what the error looks like for me:

# SYS command. line 813

 if [[ -f $(which conda) && $(conda env list | grep aquas_chipseq | wc -l) != "0" ]]; then source activate aquas_chipseq; sleep 5; fi;  export PATH=/net/bmc-pub9/data/boyerlab/users/kdemuren/TF_chipseq_pipeline/.:/net/bmc-pub9/data/boyerlab/users/kdemuren/TF_chipseq_pipeline/modules:/net/bmc-pub9/data/boyerlab/users/kdemuren/TF_chipseq_pipeline/utils:${PATH}:/bin:/usr/bin:/usr/local/bin:${HOME}/.bds; set -o pipefail; STARTTIME=$(date +%s)

# SYS command. line 814

 plotFingerprint -b /home/kdemuren/ChIP_pipeline/CM_CTCF/align/rep1/D10_CTCF_rep4.nodup.bam /home/kdemuren/ChIP_pipeline/CM_CTCF/align/rep2/D10_CTCF_rep5.nodup.bam /home/kdemuren/ChIP_pipeline/CM_CTCF/align/ctl1/D10_TF_WCE.nodup.bam --JSDsample /home/kdemuren/ChIP_pipeline/CM_CTCF/align/ctl1/D10_TF_WCE.nodup.bam \
			--labels rep1 rep2 ctl1 \
			--outQualityMetrics /home/kdemuren/ChIP_pipeline/CM_CTCF/qc/CM_CTCF_jsd.qc \
			--minMappingQuality 30 \
			-T "Fingerprints of different samples" \
			--blackListFileName /net/bmc-pub9/data/boyerlab/users/kdemuren/genome_data/mm9/mm9-blacklist.bed.gz \
			--numberOfProcessors 4 \
			--plotFile /home/kdemuren/ChIP_pipeline/CM_CTCF/qc/CM_CTCF_jsd.png

# SYS command. line 823

 TASKTIME=$[$(date +%s)-${STARTTIME}]; if [ ${TASKTIME} -lt 60 ]; then echo "Waiting for $[60-${TASKTIME}] seconds."; sleep $[60-${TASKTIME}]; sleep 0; fi
 
     
--------------------Post mortem info--------------------
==============================================================
job_number:                 2951688
exec_file:                  job_scripts/2951688
submission_time:            Fri Feb 16 11:32:35 2018
owner:                      kdemuren
uid:                        103035
group:                      bmc-boyer-nfs
gid:                        102080
sge_o_home:                 /home/kdemuren
sge_o_log_name:             kdemuren
sge_o_path:                 /home/kdemuren/miniconda3/bin:/home/software/r/r-3.0.2/bin:/home/software/jre/jre-1.8.0-77/bin:/home/software/bowtie2/bowtie2-2.2.0/scripts:/home/software/bowtie2/bowtie2-2.2.0:/home/software/bedtools/bedtools-2.20.1/bin:/home/sge/sge6.2u5/bin/lx24-amd64:/usr/local/sbin:/usr/local/bin:/usr/sbin:/usr/bin:/sbin:/bin:/usr/games:/opt/454/bin:/home/kdemuren/weblogo/:/home/kdemuren/HOMER//bin/:/home/kdemuren/.bds:/home/kdemuren/bin
sge_o_shell:                /bin/bash
sge_o_workdir:              /net/bmc-pub9/data/boyerlab/users/kdemuren/TF_chipseq_pipeline
sge_o_host:                 rous
account:                    sge
stderr_path_list:           NONE:NONE:/net/bmc-pub9/data/boyerlab/users/kdemuren/TF_chipseq_pipeline/chipseq.bds.20180216_113223_255/task.postalign_bam.jsd.line_812.id_10.stderr.cluster
hard resource_list:         h_vmem=7516192768,h_rt=21120,s_rt=21120
mail_list:                  kdemuren@rous
notify:                     FALSE
job_name:                   STDIN
stdout_path_list:           NONE:NONE:/net/bmc-pub9/data/boyerlab/users/kdemuren/TF_chipseq_pipeline/chipseq.bds.20180216_113223_255/task.postalign_bam.jsd.line_812.id_10.stdout.cluster
jobshare:                   0
env_list:                   MANPATH=/home/software/r/r-3.0.2/share/man:/home/software/jre/jre-1.8.0-77/man:/home/sge/sge6.2u5/man:/usr/share/man:/usr/local/share/man,SHELL=/bin/bash,TERM=screen,SSH_CLIENT=18.189.17.249 64437 27687,LIBRARY_PATH=/home/software/jre/jre-1.8.0-77/lib,SGE_CELL=default,SSH_TTY=/dev/pts/11,USER=kdemuren,CPATH=/home/software/jre/jre-1.8.0-77/include/linux:/home/software/jre/jre-1.8.0-77/include,TERMCAP=SC|screen|VT 100/ANSI X3.64 virtual terminal:\,MAIL=/var/mail/kdemuren,PATH=/home/kdemuren/miniconda3/bin:/home/software/r/r-3.0.2/bin:/home/software/jre/jre-1.8.0-77/bin:/home/software/bowtie2/bowtie2-2.2.0/scripts:/home/software/bowtie2/bowtie2-2.2.0:/home/software/bedtools/bedtools-2.20.1/bin:/home/sge/sge6.2u5/bin/lx24-amd64:/usr/local/sbin:/usr/local/bin:/usr/sbin:/usr/bin:/sbin:/bin:/usr/games:/opt/454/bin:/home/kdemuren/weblogo/:/home/kdemuren/HOMER//bin/:/home/kdemuren/.bds:/home/kdemuren/bin,STY=4759.CM_CTCF.BDS,_=/home/kdemuren/.bds/bds,PWD=/home/kdemuren/TF_chipseq_pipeline,JAVA_HOME=/home/software/jre/jre-1.8.0-77,_LMFILES_=/home/software/modulefiles/bedtools/2.20.1:/home/software/modulefiles/bowtie2/2.2.0:/home/software/modulefiles/jre/1.8.0-77:/home/software/modulefiles/r/3.0.2,BOWTIE_INDEXES=/home/Genomes/bowtie_indexes,SGE_EXECD_PORT=6455,LANG=en_US.UTF-8,SGE_QMASTER_PORT=6454,MODULEPATH=/home/software/jre/jre-1.8.0-77/pkg/modulefiles:/home/software/modulefiles,SGE_ROOT=/home/sge/sge6.2u5,LOADEDMODULES=bedtools/2.20.1:bowtie2/2.2.0:jre/1.8.0-77:r/3.0.2,KRB5CCNAME=FILE:/tmp/krb5cc_103035_ej4568,SHLVL=3,HOME=/home/kdemuren,LOGNAME=kdemuren,WINDOW=0,SSH_CONNECTION=18.189.17.249 64437 18.159.2.88 27687,MODULESHOME=/usr/share/Modules,SGE_CLUSTER_NAME=koch1,NLSPATH=/usr/dt/lib/nls/msg/%L/%N.cat,XFILESEARCHPATH=/usr/dt/app-defaults/%L/Dt
script_file:                STDIN
parallel environment:  whole_nodes range: 4
usage    1:                 cpu=23:05:37, mem=75600.20397 GBs, io=223.74869, vmem=N/A, maxvmem=4.565G
scheduling info:            queue instance "bcc-bmc.q@n10" dropped because it is disabled
                            queue instance "all.q@n2" dropped because it is full
                            queue instance "all.q@n3" dropped because it is full
                            queue instance "all.q@n4" dropped because it is full

from chipseq_pipeline.

leepc12 avatar leepc12 commented on July 21, 2024

from chipseq_pipeline.

kdemuren avatar kdemuren commented on July 21, 2024

I did those things, got the same error (when I removed anaconda from the PATH, the pipeline didn't start at all)

here's the longer error message from the log file, if that helps:

05:50:25.107	Writing report file 'chipseq.bds.20180221_105925_733.report.html'
05:50:32.045	ExecutionerClusterSge 'ClusterSge[34]': Error trying to find out post-mortem info on task (PID '2952081').
	Command executed : ' qstat -f -j2952081'
	Exit code        : 1
	Stdout           : 
	Stderr           : Following jobs do not exist: 
2952081

Task failed:
	Program & line     : 'modules/postalign_bam.bds', line 812
	Task Name          : 'jsd'
	Task ID            : 'chipseq.bds.20180221_105925_733/task.postalign_bam.jsd.line_812.id_10'
	Task PID           : '2952081'
	Task hint          : 'plotFingerprint -b /home/kdemuren/ChIP_pipeline/CM_CTCF/align/rep1/D10_CTCF_rep4.nodup.bam /home/kdemuren/ChIP_pipeline/CM_CTCF/align/rep2/D10_CTCF_rep5.nodup.bam /home/kdemuren/ChIP_pipeline/CM_CTCF/align/ctl1/D10_TF_WCE.nodup.bam --JSDsample /home/kdemuren/ChIP_pipeline/CM_CTCF/align/ctl1/D10_TF_W'
	Task resources     : 'cpus: 4	mem: 7.0 GB	timeout: 21000	wall-timeout: 8640000'
	State              : 'ERROR'
	Dependency state   : 'ERROR'
	Retries available  : '1'
	Input files        : '[/home/kdemuren/ChIP_pipeline/CM_CTCF/align/rep1/D10_CTCF_rep4.nodup.bam, /home/kdemuren/ChIP_pipeline/CM_CTCF/align/rep2/D10_CTCF_rep5.nodup.bam, /home/kdemuren/ChIP_pipeline/CM_CTCF/align/ctl1/D10_TF_WCE.nodup.bam]'
	Output files       : '[/home/kdemuren/ChIP_pipeline/CM_CTCF/qc/CM_CTCF_jsd.png, /home/kdemuren/ChIP_pipeline/CM_CTCF/qc/CM_CTCF_jsd.qc]'
	Script file        : '/net/bmc-pub9/data/boyerlab/users/kdemuren/TF_chipseq_pipeline/chipseq.bds.20180221_105925_733/task.postalign_bam.jsd.line_812.id_10.sh'
	Exit status        : '1'
	Program            : 
		
		# SYS command. line 813
		
		 if [[ -f $(which conda) && $(conda env list | grep aquas_chipseq | wc -l) != "0" ]]; then source activate aquas_chipseq; sleep 5; fi;  export PATH=/net/bmc-pub9/data/boyerlab/users/kdemuren/TF_chipseq_pipeline/.:/net/bmc-pub9/data/boyerlab/users/kdemuren/TF_chipseq_pipeline/modules:/net/bmc-pub9/data/boyerlab/users/kdemuren/TF_chipseq_pipeline/utils:${PATH}:/bin:/usr/bin:/usr/local/bin:${HOME}/.bds; set -o pipefail; STARTTIME=$(date +%s)
		
		# SYS command. line 814
		
		 plotFingerprint -b /home/kdemuren/ChIP_pipeline/CM_CTCF/align/rep1/D10_CTCF_rep4.nodup.bam /home/kdemuren/ChIP_pipeline/CM_CTCF/align/rep2/D10_CTCF_rep5.nodup.bam /home/kdemuren/ChIP_pipeline/CM_CTCF/align/ctl1/D10_TF_WCE.nodup.bam --JSDsample /home/kdemuren/ChIP_pipeline/CM_CTCF/align/ctl1/D10_TF_WCE.nodup.bam \
					--labels rep1 rep2 ctl1 \
					--outQualityMetrics /home/kdemuren/ChIP_pipeline/CM_CTCF/qc/CM_CTCF_jsd.qc \
					--minMappingQuality 30 \
					-T "Fingerprints of different samples" \
					--blackListFileName /net/bmc-pub9/data/boyerlab/users/kdemuren/genome_data/mm9/mm9-blacklist.bed.gz \
					--numberOfProcessors 4 \
					--plotFile /home/kdemuren/ChIP_pipeline/CM_CTCF/qc/CM_CTCF_jsd.png
		
		# SYS command. line 823
		
		 TASKTIME=$[$(date +%s)-${STARTTIME}]; if [ ${TASKTIME} -lt 60 ]; then echo "Waiting for $[60-${TASKTIME}] seconds."; sleep $[60-${TASKTIME}]; sleep 0; fi

05:50:32.123	Wait: Task 'chipseq.bds.20180221_105925_733/task.postalign_bam.jsd.line_812.id_10' finished.
05:50:32.123	Waiting for all 'parrallel' to finish.
05:50:32.123	Waiting for parallel 'chipseq.bds.20180221_105925_733_parallel_30' to finish. RunState: FINISHED
05:50:32.123	Waiting for parallel 'chipseq.bds.20180221_105925_733_parallel_31' to finish. RunState: FINISHED
05:50:32.123	Waiting for parallel 'chipseq.bds.20180221_105925_733_parallel_32' to finish. RunState: FINISHED
05:50:32.123	Waiting for parallel 'chipseq.bds.20180221_105925_733_parallel_33' to finish. RunState: FINISHED
Fatal error: /net/bmc-pub9/data/boyerlab/users/kdemuren/TF_chipseq_pipeline/chipseq.bds, line 414, pos 3. Task/s failed.
chipseq.bds, line 91 :	main()
chipseq.bds, line 94 :	void main() { // chipseq pipeline starts here
chipseq.bds, line 101 :		jsd() // plot fingerprint and compute synthetic JS distance
chipseq.bds, line 387 :	void jsd() { // plot fingerprint
chipseq.bds, line 389 :		if ( filt_bam.hasKey("ctl1") && filt_bam.hasKey("rep1") ) {
chipseq.bds, line 414 :			wait

Creating checkpoint file: Config or command line option disabled checkpoint file creation, nothing done.
05:50:32.218	Writing report file 'chipseq.bds.20180221_105925_733.report.html'
05:50:32.226	Program 'chipseq.bds.20180221_105925_733' finished, exit value: 1, tasks executed: 1, tasks failed: 1, tasks failed names: jsd.
05:50:32.226	Finished. Exit code: 1
05:50:32.233	ExecutionerClusterSge 'ClusterSge[34]': Killed

from chipseq_pipeline.

leepc12 avatar leepc12 commented on July 21, 2024

from chipseq_pipeline.

kdemuren avatar kdemuren commented on July 21, 2024

I unloaded the python module, checked out my PATH and PYTHONPATH, and am still having the problem. There's nothing that says /mount/gensoft2 in any of those, so I'm not sure where that is coming from. Here's my .bashrc:

# .bashrc

# Source global definitions
if [ -f /etc/bashrc ]; then
        . /etc/bashrc
fi

# User specific aliases and functions
#export PYTHONPATH=$PYTHONPATH:$HOME/lib/python2.7/site-packages
#export PYTHONPATH=$PYTHONPATH:/home/software/python/python-3.5.1/lib/python3.5/site-packages
PATH=$PATH:/home/kdemuren/weblogo/
PATH=$PATH:/home/kdemuren/HOMER//bin/
export PATH=$PATH:$HOME/.bds
#module load bedtools/2.20.1
#module load bowtie2/2.2.0
module load jre/1.8.0-77
#module load r/3.0.2
#module load python/2.7.2

module unload python
unset module
`unset PYTHONPATH`
# added by Miniconda3 installer
export PATH="/home/kdemuren/miniconda3/bin:$PATH"

and here's my PATH (PYTHONPATH is empty):

/home/kdemuren/miniconda3/bin:/home/software/jre/jre-1.8.0-77/bin:/home/sge/sge6.2u5/bin/lx24-amd64:/usr/local/sbin:/usr/local/bin:/usr/sbin:/usr/bin:/sbin:/bin:/usr/games:/opt/454/bin:/home/kdemuren/weblogo/:/home/kdemuren/HOMER//bin/:/home/kdemuren/.bds:/home/kdemuren/bin

Thanks,
Kunle

from chipseq_pipeline.

leepc12 avatar leepc12 commented on July 21, 2024

Please add -no_jsd to the command line until this issue is fixed. I will take a closer look at this error later.

from chipseq_pipeline.

leepc12 avatar leepc12 commented on July 21, 2024

Also, your ~/.bashrc looks wrong unset PYTHONPATH.

from chipseq_pipeline.

kvittingseerup avatar kvittingseerup commented on July 21, 2024

I just got the same error - is there an update on the problem?

from chipseq_pipeline.

leepc12 avatar leepc12 commented on July 21, 2024

@kvittingseerup did you remove all conflicting modules and python from your bash startup scripts? and did you add unset PYTHONPATH to it?

from chipseq_pipeline.

kvittingseerup avatar kvittingseerup commented on July 21, 2024

I'll have a hard time removing python from my bashrc as I'm using it all the time. Can I specify it in the script(s) instead?

from chipseq_pipeline.

leepc12 avatar leepc12 commented on July 21, 2024

@kvittingseerup: No.
Please try unset PYTHONPATH first and then if it doesn't help go with -no_jsd.

from chipseq_pipeline.

Related Issues (20)

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. 📊📈🎉

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google ❤️ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.