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) } ############################################ myname <- "myl" # Your PennName here durationfile <- paste(myname,"durations",sep="_") digitfile <- paste(myname,"digits", sep="_") # write.table(round(durations),file=durationfile) write.table(digits,file=digitfile)