XX <- read.table("duration_log.txt", header=TRUE) starts <- seq(1,991,by=10) files <- XX[starts,1] digits <- matrix(nrow=100,ncol=10) durations <- matrix(nrow=100,ncol=10) for (n in 1:100){ for (m in 0:9){ digits[n,m+1] = XX[starts[n]+m,2] durations[n,m+1] = XX[starts[n]+m,3] } } counts <- 1:10; digitdurs <- 1:10; positiondurs <- 1:10 dpdurs <- matrix(nrow=10,ncol=10) dpcounts <- matrix(nrow=10,ncol=10) # Duration by digits for (digit in 0:9){ digitdurs[digit+1] <- mean(durations[digits==digit]) positiondurs[digit+1] <- mean(durations[,digit+1]) for (position in 1:10){ which <- digits[,position]==digit dpcounts[digit+1,position] <- sum(which) dpdurs[digit+1,position] <- mean(durations[which,position]) } } # # Sanity check for (n in 0:9){ counts[n+1] <- sum(digits==n) } ### ylims=range(dpdurs) par(ask=TRUE) plot(1:10,dpdurs[1,], main="3|3|4||: Mean digit duration by position", type="n", xlab="Position in 10-digit string", ylab="Mean Duration (ms)", ylim=ylims) colorlist=c("aquamarine", "black", "cyan", "darkblue", "red", "deeppink", "grey", "limegreen", "magenta", "salmon") for (digit in 0:9){ lines(1:10,dpdurs[digit+1,],type="b",pch=as.character(digit),col=colorlist[digit+1]) } posmeans = apply(durations, 2, mean) plot(1:10,posmeans, xlab="Position in string", ylab="Mean digit duration", type="b", col="red",pch="x") write.table(posmeans,file="posmeans")