|
0
|
1 #!/bin/bash
|
|
|
2
|
|
|
3 tool_path=$(dirname $0)
|
|
|
4 ped=$1
|
|
|
5 map=$2
|
|
4
|
6 fileout_label=$(date "+%Y%m%d%H%M%S")
|
|
|
7 fileout_matrix=$3
|
|
|
8 fileout_plot=$4
|
|
|
9 fileout_log=$5
|
|
8
|
10 groups=$6
|
|
|
11
|
|
0
|
12
|
|
|
13 rsync -a $ped input.ped
|
|
|
14 rsync -a $map input.map
|
|
8
|
15 if [ -f $groups ]
|
|
|
16 then
|
|
|
17 cp -rf $groups input.individual_info.txt
|
|
|
18 fi
|
|
0
|
19
|
|
|
20 perl $tool_path/MDSbasedOnIBSmatrix.pl --in input --out $fileout_label
|
|
|
21
|
|
|
22 rm -f input.ped input.map
|
|
|
23
|
|
|
24 cp $fileout_label.ibs_matrix.txt $fileout_matrix
|
|
|
25 cp $fileout_label.mds_plot.txt $fileout_plot
|
|
|
26 cp input.plink.log $fileout_log
|
|
|
27
|
|
|
28
|
|
|
29 rm -f $fileout_label.ibs_matrix.txt $fileout_label.mds_plot.txt input.plink.log
|
