Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Summarizer: re-order results #16

Open
swvanderlaan opened this issue Aug 22, 2018 · 0 comments
Open

Summarizer: re-order results #16

swvanderlaan opened this issue Aug 22, 2018 · 0 comments
Assignees
Labels

Comments

@swvanderlaan
Copy link
Owner

There is an issue with the re-ordering of the nominal results: these are empty after attempting to re-order using the following statements:

if [[ ${QTL_TYPE} == "CIS" ]]; then
	
		if [[ ${CLUMP} == "Y" ]]; then
			echo "* Clumped ${QTL_TYPE}-QTL results."
			gzip -fv ${SUMMARY}/${STUDYNAME}_QC_qtlnom_clumped_summary.txt
			gzip -fv ${SUMMARY}/${STUDYNAME}_QC_qtlperm_clumped_summary.txt
			
			echo "* And re-ordering based on p-value."
			zcat ${SUMMARY}/${STUDYNAME}_QC_qtlnom_clumped_summary.txt.gz | (head -n 1 && tail -n +3  | sort -t , -k 24) > ${SUMMARY}/${STUDYNAME}_QC_qtlnom_clumped_summary.txt
			gzip -fv ${SUMMARY}/${STUDYNAME}_QC_qtlnom_clumped_summary.txt
		
		elif [[ ${ANALYSIS_TYPE} == "MQTL" ]]; then
			echo "* Non-clumped ${QTL_TYPE}-QTL results."
			gzip -fv ${SUMMARY}/${STUDYNAME}_QC_qtlperm_summary.txt
			
			echo "* And re-ordering based on p-value."
			cat ${SUMMARY}/${STUDYNAME}_QC_qtlnom_summary.txt | (head -n 1 && tail -n +3  | sort -t , -k 32) > ${SUMMARY}/${STUDYNAME}_QC_qtlnom_summary.txt
			gzip -fv ${SUMMARY}/${STUDYNAME}_QC_qtlnom_summary.txt
			
		else
			echo "* Non-clumped ${QTL_TYPE}-QTL results."
			gzip -fv ${SUMMARY}/${STUDYNAME}_QC_qtlperm_summary.txt
			
			echo "* And re-ordering based on p-value."
			cat ${SUMMARY}/${STUDYNAME}_QC_qtlnom_summary.txt | (head -n 1 && tail -n +3  | sort -t , -k 24) > ${SUMMARY}/${STUDYNAME}_QC_qtlnom_summary.txt
			gzip -fv ${SUMMARY}/${STUDYNAME}_QC_qtlnom_summary.txt
			
		fi

	elif [[ ${QTL_TYPE} == "TRANS" ]]; then
		gzip -fv ${SUMMARY}/${STUDYNAME}_QC_qtlnom_summary.txt

	fi
@swvanderlaan swvanderlaan self-assigned this Aug 22, 2018
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

1 participant