decodebiology

echo `date`;
dir=(`dirname $1`)

base=(`basename $1`)



echo "1/7 Preparing files";
cat $1 | sed 's/ "/\t/g' | sed 's/"; /\t/g' | cut -f1,4,5,7,10,12,16,18,22 | sed 's/";//' | awk '!x[$6]++' | awk '{print $6"\t"$0;}' > $dir/enst_annotation.tmp;

decodebiology

# Bland-Altman plot R function.
# Author: jmmateos@mce.hggm.es
baplot <- function(m1, m2, ...) {
    # m1 and m2 are the measurements
    means <- (m1 + m2) / 2
    diffs <- m1 - m2
    mdiff <- mean(diffs)
    sddiff <- sd(diffs)
    
    # Compute the figure limits