library(ggplot2) library(tcltk) speed <- read.csv("speed.csv", header=FALSE) means <- apply(speed,1,mean) speed_norm <- speed/means std_speed_norm <- apply(speed_norm,1,sd) # ggplot format legend_theme <- theme(plot.title = element_text(size=20, hjust = 0.5, face="bold"), axis.title.x = element_text(size=14, face="bold"), axis.title.y = element_text(size=14, face="bold"), legend.title=element_text(size=14), legend.text=element_text(size=12)) x11() # windows() gg <- ggplot() + geom_point(aes(x=means*3.6, y=std_speed_norm), color="blue", size = 0.1) + geom_smooth(aes(x=means*3.6, y=std_speed_norm), color="red", method='lm') + ggtitle("Standard deviation of normalized Speed") + xlab("Average speed (km/h)") + ylab("Standard deviation") + xlim(min(means*3.6), max(means*3.6)) + ylim(0, 0.4) + legend_theme print(gg) capture <- tk_messageBox(message = "hit spacebar to close plots", detail = "close now") # speed[1:2,1:2] # apply(speed[1:2,1:2],1,mean) # speed[1:2,1:2]/c(1,2)