Annalyzing data using R

Analyze.R

#Author : Muhammad Ehsan ul Haq

#Read the data from the text file
dataAlpha = read.table("dataAlpha.txt", header = TRUE) 
dataBravo = read.table("dataBravo.txt", header = TRUE)
#number of no-responses
noResponseAlpha = sum(dataAlpha == -1)
noResponseBravo = sum(dataBravo == -1)

print(sprintf("Number of no Responses in Alpha = %d", noResponseAlpha))
print(sprintf("Number of no Responses in Bravo = %d", noResponseBravo))

#finding standard deviation
sdAlpha = sd(unlist(dataAlpha["RTT"]))
sdBravo = sd(unlist(dataBravo["RTT"]))

print(sprintf("Standard Deviation of RTT of Alpha = %f", sdAlpha))
print(sprintf("Standard Deviation of RTT of Bravo = %f", sdBravo))

#finding mean of the data
mAlpha = mean(unlist(dataAlpha["RTT"]))
mBravo = mean(unlist(dataBravo["RTT"]))

print(sprintf("Mean of RTT of Alpha = %f", mAlpha))
print(sprintf("Mean of RTT of Bravo = %f", mBravo))

#plotting the data
plot(data.frame(dataAlpha["StartTime"], dataAlpha["RTT"]), "Time", type = "p", main = "Alpha")
plot(data.frame(dataBravo["StartTime"], dataBravo["RTT"]), "Time", type = "p", main = "Bravo", col = "red")

plot(data.frame(dataAlpha["StartTime"], dataAlpha["RTT"]), "Time", type = "p", main = "Alpha and Bravo")
points(data.frame(dataBravo["StartTime"], dataBravo["RTT"]), pch = 4, col = "red")
legend("topright", legend = c("Alpha", "Bravo"), col = c("black", "red"), pch = c(1, 4), cex = 0.6)