2 files changed, 102 insertions, 5 deletions
diff --git a/code/rtt/rttests.py b/code/rtt/rttests.py
index 6356397..fdc99bb 100644
--- a/code/rtt/rttests.py
+++ b/code/rtt/rttests.py
@@ -46,7 +46,7 @@ def writeestRTTs(datapath, samples, estimates):
 def calcEST(samples):
     ests = []
     K = 4
-    G = 0.1 # =100 msec
+    G = 0.001 # =100 msec
     alpha = 1/8.0
     beta = 1/4.0
 
diff --git a/code/rtt/rttplot.R b/code/rtt/rttplot.R
index 1d9b871..b4f3379 100644
--- a/code/rtt/rttplot.R
+++ b/code/rtt/rttplot.R
@@ -1,12 +1,109 @@
+options("scipen" = 10)
 data = read.csv(file="../../data/rtt/numpak/r1est.csv", header=FALSE, sep=",")
 
-frame_number = log(data$V1)
+frame_number = data$V1
 rtt_sample = data$V2
 rtt_est = data$V3
 
-plot(frame_number, rtt_sample, type="l", xlab="log_10 Frame numbers", ylab="time (s)", main="Num. Packets - Rank 1 - RTT Sample and RTT Estimate")
-lines(frame_number, rtt_est, col="blue")
-legend(13,0.15, legend=c("Sample", "Estimate"), col=c("black", "blue"), lty=1)
+plot(frame_number, rtt_sample, type="l", xlab="Frame numbers", ylab="time (s)", main="Num. Packets - Rank 1 - RTT Sample and RTT Estimate")
+lines(frame_number, rtt_est, col="blue", type="l")
+legend('topright', legend=c("Sample", "Estimate"), col=c("black", "blue"), lty=1)
 
 dev.print(png, '../../plots/rtt/numpak/r1.png', width=500)
+dev.off()
 
+data = read.csv(file="../../data/rtt/numpak/r2est.csv", header=FALSE, sep=",")
+frame_number = data$V1
+rtt_sample = data$V2
+rtt_est = data$V3
+
+plot(frame_number, rtt_sample, type="l", xlab="Frame numbers", ylab="time (s)", main="Num. Packets - Rank 2 - RTT Sample and RTT Estimate")
+lines(frame_number, rtt_est, col="blue", type="l")
+legend('topright', legend=c("Sample", "Estimate"), col=c("black", "blue"), lty=1)
+
+dev.print(png, '../../plots/rtt/numpak/r2.png', width=500)
+dev.off()
+
+data = read.csv(file="../../data/rtt/numpak/r3est.csv", header=FALSE, sep=",")
+frame_number = data$V1
+rtt_sample = data$V2
+rtt_est = data$V3
+
+plot(frame_number, rtt_sample, type="l", xlab="Frame numbers", ylab="time (s)", main="Num. Packets - Rank 3 - RTT Sample and RTT Estimate")
+lines(frame_number, rtt_est, col="blue", type="l")
+legend('topright', legend=c("Sample", "Estimate"), col=c("black", "blue"), lty=1)
+
+dev.print(png, '../../plots/rtt/numpak/r3.png', width=500)
+dev.off()
+
+data = read.csv(file="../../data/rtt/bytes/r1est.csv", header=FALSE, sep=",")
+frame_number = data$V1
+rtt_sample = data$V2
+rtt_est = data$V3
+
+plot(frame_number, rtt_sample, type="l", xlab="Frame numbers", ylab="time (s)", main="Total Bytes - Rank 1 - RTT Sample and RTT Estimate")
+lines(frame_number, rtt_est, col="blue", type="l")
+legend('topright', legend=c("Sample", "Estimate"), col=c("black", "blue"), lty=1)
+
+dev.print(png, '../../plots/rtt/bytes/r1.png', width=500)
+dev.off()
+
+data = read.csv(file="../../data/rtt/bytes/r2est.csv", header=FALSE, sep=",")
+frame_number = data$V1
+rtt_sample = data$V2
+rtt_est = data$V3
+
+plot(frame_number, rtt_sample, type="l", xlab="Frame numbers", ylab="time (s)", main="Total Bytes - Rank 2 - RTT Sample and RTT Estimate")
+lines(frame_number, rtt_est, col="blue", type="l")
+legend('topright', legend=c("Sample", "Estimate"), col=c("black", "blue"), lty=1)
+
+dev.print(png, '../../plots/rtt/bytes/r2.png', width=500)
+dev.off()
+
+data = read.csv(file="../../data/rtt/bytes/r3est.csv", header=FALSE, sep=",")
+frame_number = data$V1
+rtt_sample = data$V2
+rtt_est = data$V3
+
+plot(frame_number, rtt_sample, type="l", xlab="Frame numbers", ylab="time (s)", main="Total Bytes - Rank 3 - RTT Sample and RTT Estimate")
+lines(frame_number, rtt_est, col="blue", type="l")
+legend('topright', legend=c("Sample", "Estimate"), col=c("black", "blue"), lty=1)
+
+dev.print(png, '../../plots/rtt/bytes/r3.png', width=500)
+dev.off()
+
+data = read.csv(file="../../data/rtt/dur/r1est.csv", header=FALSE, sep=",")
+frame_number = data$V1
+rtt_sample = data$V2
+rtt_est = data$V3
+
+plot(frame_number, rtt_sample, type="l", xlab="Frame numbers", ylab="time (s)", main="Duration - Rank 1 - RTT Sample and RTT Estimate")
+lines(frame_number, rtt_est, col="blue", type="l")
+legend('topright', legend=c("Sample", "Estimate"), col=c("black", "blue"), lty=1)
+
+dev.print(png, '../../plots/rtt/dur/r1.png', width=500)
+dev.off()
+
+data = read.csv(file="../../data/rtt/dur/r2est.csv", header=FALSE, sep=",")
+frame_number = data$V1
+rtt_sample = data$V2
+rtt_est = data$V3
+
+plot(frame_number, rtt_sample, type="l", xlab="Frame numbers", ylab="time (s)", main="Duration - Rank 2 - RTT Sample and RTT Estimate")
+lines(frame_number, rtt_est, col="blue", type="l")
+legend('topright', legend=c("Sample", "Estimate"), col=c("black", "blue"), lty=1)
+
+dev.print(png, '../../plots/rtt/dur/r2.png', width=500)
+dev.off()
+
+data = read.csv(file="../../data/rtt/dur/r3est.csv", header=FALSE, sep=",")
+frame_number = data$V1
+rtt_sample = data$V2
+rtt_est = data$V3
+
+plot(frame_number, rtt_sample, type="l", xlab="Frame numbers", ylab="time (s)", main="Duration - Rank 3 - RTT Sample and RTT Estimate")
+lines(frame_number, rtt_est, col="blue", type="l")
+legend('topright', legend=c("Sample", "Estimate"), col=c("black", "blue"), lty=1)
+
+dev.print(png, '../../plots/rtt/dur/r3.png', width=500)
+dev.off()