setwd("/mikeData/Dropbox/rgits/baleen/") folders <- grep(pattern = "bend", x = list.files(),value = T) pdf(file = "resultsSummary.pdf") for(i in folders){ data <- read.csv(paste(i,"/Test1/Test1.steps.tracking.csv",sep="")) x=data$Displacement.E1000.Digital.Position...mm. y=data$Load.E1000.Load...N. x=x-x[1] x=-x y=-y breaking <- max(y) which.max(y) -> breakingIndex plot(x=x,y=y, type="l", xlab="Displacement (mm)", ylab="Force (N)", main=paste(i), sub=paste("max force:", breaking,"N")) points(x=x[breakingIndex], y=y[breakingIndex], col="red") } dev.off()