|
@@ -0,0 +1,41 @@
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+for inputfile in "$@"
|
|
|
+do
|
|
|
+echo "Handle file $inputfile"
|
|
|
+hashlist=`cut -f 1 $inputfile | sort | uniq | xargs`
|
|
|
+
|
|
|
+
|
|
|
+for h in $hashlist
|
|
|
+do
|
|
|
+ echo "Handling tasks with hash = $h"
|
|
|
+ grep "^$h" $inputfile| cut -f 2 > $inputfile.$h
|
|
|
+
|
|
|
+R --no-save > /dev/null << EOF
|
|
|
+
|
|
|
+x <- scan("$inputfile.$h")
|
|
|
+hist(x[x > quantile(x,0.01) & x<quantile(x,0.99)], col="red", breaks=50, density=10)
|
|
|
+
|
|
|
+EOF
|
|
|
+mv Rplots.pdf $inputfile.$h.pdf
|
|
|
+done
|
|
|
+
|
|
|
+done
|
|
|
+echo "finished !"
|