zzh
/
gatk-3.8
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Added instantaneous rate plotting routine to performance log plotter 

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@4926 348d0f76-0448-11de-a6fe-93d51630548a
This commit is contained in:
depristo 2011-01-03 13:33:35 +00:00
parent 8d2c14b29c
commit 54adbd2581
1 changed files with 26 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

27
R/plot_GATK_performance_log.R
View File

@ -13,10 +13,35 @@ if ( onCMDLine ) {
if ( onCMDLine ) pdf(OUTPUT_PDF)
vec.margin <- function(x) {
l = length(x)
d = x[-1] - x[1:(l-1)]
c(x[1], d[1:(l-1)])
}
l = length(d$units.processed)
d$units.processed.margin = vec.margin(d$units.processed)
#prev = 0
#for ( i in 1:l ) {
# cur = d$units.processed[i]
# d[i,]$units.processed.margin = cur - prev
# prev = cur
#}
generateOneReport <- function(d) {
qs = quantile(d$processing.speed, probs = c(0.05, 0.5, 0.95))
qs = quantile(d$processing.speed, probs = c(0.01, 0.5, 0.99))
# unit processing time
plot(d$elapsed.time, d$processing.speed, main=DESCRIPTION, xlab="Elapsed time (sec)", ylab="Processing speed (seconds per 1M units)", ylim=c(qs[1], qs[3]), type="b", col="cornflowerblue", lwd=2)
abline(h=qs[2], lty=2)
# instantaneous processing speed
running_median_window = 101
rm = runmed(d$units.processed.margin, running_median_window)
POINT_COL = "#0000AA33"
plot(d$elapsed.time, d$units.processed.margin, main=DESCRIPTION, xlab="Elapsed time (sec)", ylab="Units processed in last timing interval", type="p", cex = 0.5, col=POINT_COL)
lines(d$elapsed.time, rm, lwd=3, col="red")
legend("topleft", c("Observations", "101-elt running median"), fill=c(POINT_COL, "red"))
}
if ( RUNME ) {