Skip to content

Instantly share code, notes, and snippets.

@wbazant
Created September 20, 2017 10:42
Show Gist options
  • Select an option

  • Save wbazant/f936724e0d07b0ce79cd49567bb625dc to your computer and use it in GitHub Desktop.

Select an option

Save wbazant/f936724e0d07b0ce79cd49567bb625dc to your computer and use it in GitHub Desktop.
grep GO:0003713 /nfs/production3/ma/home/atlas3-production/bioentity_properties/archive/ensembl_89_36/mus_musculus.ensgene.go.tsv | cut -f 1 | sort -u > $ATLAS_PROD/go-unstable/pvalue-by-hand/genes-for-GO_0003713.txt
head -n 1 $ATLAS_EXPS/E-GEOD-59831/E-GEOD-59831-analytics.tsv > $ATLAS_PROD/go-unstable/pvalue-by-hand/E-GEOD-59831-analytics.tsv
join -1 1 -2 1 <(sort -k1 $ATLAS_EXPS/E-GEOD-59831/E-GEOD-59831-analytics.tsv) $ATLAS_PROD/go-unstable/pvalue-by-hand/genes-for-GO_0003713.txt >> $ATLAS_PROD/go-unstable/pvalue-by-hand/E-GEOD-59831-analytics.tsv
# (with headers)
wc -l $ATLAS_EXPS/E-GEOD-59831/E-GEOD-59831-analytics.tsv $ATLAS_PROD/go-unstable/pvalue-by-hand/E-GEOD-59831-analytics.tsv
# 45514 /nfs/public/ro/fg/atlas/experiments/E-GEOD-59831/E-GEOD-59831-analytics.tsv
# 177 /nfs/production3/ma/home/atlas3-production/go-unstable/pvalue-by-hand/E-GEOD-59831-analytics.tsv
tail -n +1 /nfs/public/ro/fg/atlas/experiments/E-GEOD-59831/E-GEOD-59831-analytics.tsv | cut -f 3 | grep -v NA | wc -l
# 20483
tail -n +1 $ATLAS_PROD/go-unstable/pvalue-by-hand/E-GEOD-59831-analytics.tsv | cut -f 3 | grep -v NA | wc -l
# 158
cut -f 1 /nfs/production3/ma/home/atlas3-production/bioentity_properties/archive/ensembl_89_36/mus_musculus.ensgene.go.tsv | sort -u | wc -l
# 51826
# (s,k,M,N) = (20483, 158, 177, 51826)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment