Skip to content
Snippets Groups Projects
Commit 3b0613c8 authored by Aaron's avatar Aaron
Browse files

added preliminary tests for map

parent e620be4c
No related branches found
No related tags found
No related merge requests found
chr1 0 100
chr1 100 200
chr2 0 100
chr2 100 200
chr3 0 100
chr3 100 200
BT=../../bin/bedtools
check()
{
if diff $1 $2; then
echo ok
else
echo fail
fi
}
###########################################################
# Test defaults
############################################################
echo " map.t01...\c"
echo \
"chr1 0 100 30
chr1 100 200 1
chr2 0 100 .
chr2 100 200 .
chr3 0 100 6
chr3 100 200 4" > exp
$BT map -a ivls.bed -b values.bed > obs
check obs exp
rm obs exp
###########################################################
# Test sum
############################################################
echo " map.t02...\c"
echo \
"chr1 0 100 30
chr1 100 200 1
chr2 0 100 .
chr2 100 200 .
chr3 0 100 6
chr3 100 200 4" > exp
$BT map -a ivls.bed -b values.bed -o sum > obs
check obs exp
rm obs exp
###########################################################
# Test count
############################################################
echo " map.t03...\c"
echo \
"chr1 0 100 3
chr1 100 200 1
chr2 0 100 0
chr2 100 200 0
chr3 0 100 3
chr3 100 200 1" > exp
$BT map -a ivls.bed -b values.bed -o count > obs
check obs exp
rm obs exp
###########################################################
# Test mean
############################################################
echo " map.t04...\c"
echo \
"chr1 0 100 10
chr1 100 200 1
chr2 0 100 .
chr2 100 200 .
chr3 0 100 2
chr3 100 200 4" > exp
$BT map -a ivls.bed -b values.bed -o mean > obs
check obs exp
rm obs exp
###########################################################
# Test max
############################################################
echo " map.t05...\c"
echo \
"chr1 0 100 15
chr1 100 200 1
chr2 0 100 .
chr2 100 200 .
chr3 0 100 3
chr3 100 200 4" > exp
$BT map -a ivls.bed -b values.bed -o max > obs
check obs exp
rm obs exp
###########################################################
# Test min
############################################################
echo " map.t06...\c"
echo \
"chr1 0 100 5
chr1 100 200 1
chr2 0 100 .
chr2 100 200 .
chr3 0 100 1
chr3 100 200 4" > exp
$BT map -a ivls.bed -b values.bed -o min > obs
check obs exp
rm obs exp
###########################################################
# Test mode
############################################################
echo " map.t07...\c"
echo \
"chr1 0 100 5
chr1 100 200 1
chr2 0 100 .
chr2 100 200 .
chr3 0 100 1
chr3 100 200 4" > exp
$BT map -a ivls.bed -b values2.bed -o mode > obs
check obs exp
rm obs exp
###########################################################
# Test anti-mode
############################################################
echo " map.t08...\c"
echo \
"chr1 0 100 10
chr1 100 200 1
chr2 0 100 .
chr2 100 200 .
chr3 0 100 1
chr3 100 200 4" > exp
$BT map -a ivls.bed -b values2.bed -o antimode > obs
check obs exp
rm obs exp
chr1 0 10 a1 10 +
chr1 10 20 a2 5 +
chr1 20 30 a3 15 +
chr1 120 130 a4 1 +
chr3 0 10 a5 1 +
chr3 10 20 a6 2 +
chr3 20 30 a7 3 +
chr3 120 130 a8 4 +
\ No newline at end of file
chr1 0 10 a1 10 +
chr1 10 20 a2 5 +
chr1 11 21 a3 5 +
chr1 20 30 a4 15 +
chr1 120 130 a5 1 +
chr3 0 10 a6 1 +
chr3 10 20 a7 2 +
chr3 20 30 a8 3 +
chr3 120 130 a9 4 +
chr3 120 130 a9 4 +
chr1 0 10 a1 10 +
chr1 10 20 a2 5 +
chr1 11 21 a3 5 +
chr1 20 30 a4 15 +
chr1 20 30 a5 15 +
chr1 20 30 a6 15 +
chr1 120 130 a7 1 +
chr3 0 10 a8 1 +
chr3 10 20 a9 2 +
chr3 20 30 a10 3 +
chr3 120 130 a11 4 +
chr3 120 130 a12 4 +
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment