redirect stdout to the output file, fix generation of the R plot
depend on libreadline, again
add initial callpeak test