d <- read.table("data_v00.txt",header=T)
load ("Normal_v00_R.gz")

winm <- c(0:999)*.05
wins <- c(0:999)*.02

pdf(file="plots.pdf")

for (i in (1:50)){
    plot(winm, -s[,1,i],
      type="l", 
      main="http:/www.neurotraces.com/edf",
      sub=sprintf("Median nerve, recording APB, stimulus wrist, case %d",i), 
      xlab="ms",
      ylab="mV",ylim=c(-20,20))
    l1 <- d[7*(i-1)+1,]$lat1
    l2 <- d[7*(i-1)+1,]$lat2
    l3 <- d[7*(i-1)+1,]$lat3
    l4 <- d[7*(i-1)+1,]$lat4
    lines (c(l1,l1),c(-20,20),col="blue")
    lines (c(l2,l2),c(-20,20),col="pink")
    lines (c(l3,l3),c(-20,20),col="green")
    lines (c(l4,l4),c(-20,20),col="light blue")
}

for (i in (1:50)){
    plot(winm, -s[,2,i],
      type="l", 
      main="http:/www.neurotraces.com/edf",
      sub=sprintf("Median nerve, recording APB, stimulus elbow, case %d",i), 
      xlab="ms",
      ylab="mV",
      ylim=c(-20,20))
    l1 <- d[7*(i-1)+2,]$lat1
    l2 <- d[7*(i-1)+2,]$lat2
    l3 <- d[7*(i-1)+2,]$lat3
    l4 <- d[7*(i-1)+2,]$lat4
    lines (c(l1,l1),c(-20,20),col="blue")
    lines (c(l2,l2),c(-20,20),col="pink")
    lines (c(l3,l3),c(-20,20),col="green")
    lines (c(l4,l4),c(-20,20),col="light blue")
}


for (i in (1:50)){
    plot(wins, -s[,3,i],type="l", 
      main="http:/www.neurotraces.com/edf",
      sub=sprintf("Median nerve. recording 2d finger, stimulus wrist, case %d",i), 
      xlab="ms", 
      ylab="uV",
      ylim=c(-100,100))
    l1 <- d[7*(i-1)+3,]$lat1
    l2 <- d[7*(i-1)+3,]$lat2
    l3 <- d[7*(i-1)+3,]$lat3
    l4 <- d[7*(i-1)+3,]$lat4
    lines (c(l1,l1),c(-100,100),col="blue") 
    lines (c(l2,l2),c(-100,100),col="pink")
    lines (c(l3,l3),c(-100,100),col="green")
    lines (c(l4,l4),c(-100,100),col="light blue")
} 

for (i in (1:50)){
    plot(wins, -s[,4,i],type="l", 
      main="http:/www.neurotraces.com/edf",      
      sub=sprintf("Median nerve, recording 2d finger, stimulus elbow, case %d",i), 
      xlab="ms", 
      ylab="uV",
      ylim=c(-100,100))
    l1 <- d[7*(i-1)+4,]$lat1
    l2 <- d[7*(i-1)+4,]$lat2
    l3 <- d[7*(i-1)+4,]$lat3
    l4 <- d[7*(i-1)+4,]$lat4
    lines (c(l1,l1),c(-100,100),col="blue")
    lines (c(l2,l2),c(-100,100),col="pink")
    lines (c(l3,l3),c(-100,100),col="green")
    lines (c(l4,l4),c(-100,100),col="light blue")
}

for (i in (1:50)){
    plot(wins, -s[,5,i],type="l", 
      main="http:/www.neurotraces.com/edf",
      sub=sprintf("Ulnar nerve, recording 5th finger, stimulus wrist, case %d",i), 
      xlab="ms",
      ylab="uV",
      ylim=c(-100,100))
    l1 <- d[7*(i-1)+5,]$lat1
    l2 <- d[7*(i-1)+5,]$lat2
    l3 <- d[7*(i-1)+5,]$lat3
    l4 <- d[7*(i-1)+5,]$lat4
    lines (c(l1,l1),c(-100,100),col="blue")
    lines (c(l2,l2),c(-100,100),col="pink")
    lines (c(l3,l3),c(-100,100),col="green")
    lines (c(l4,l4),c(-100,100),col="light blue")
}

for (i in (1:50)){
    plot(wins, -s[,6,i],type="l", 
      main="http:/www.neurotraces.com/edf",
      sub=sprintf("Radial nerve, recording 1st finger, stimulus wrist, case %d",i), 
      xlab="ms",
      ylab="uV",
      ylim=c(-100,100))
    l1 <- d[7*(i-1)+6,]$lat1
    l2 <- d[7*(i-1)+6,]$lat2
    l3 <- d[7*(i-1)+6,]$lat3
    l4 <- d[7*(i-1)+6,]$lat4
    lines (c(l1,l1),c(-100,100),col="blue")
    lines (c(l2,l2),c(-100,100),col="pink")
    lines (c(l3,l3),c(-100,100),col="green")
    lines (c(l4,l4),c(-100,100),col="light blue")
}

for (i in (1:50)){
    plot(wins, -s[,7,i],type="l", 
      main="http:/www.neurotraces.com/edf",
      sub=sprintf("Median nerve, recording 1st finger, stimulus wrist, case %d",i), 
      xlab="ms",
      ylab="uV",
      ylim=c(-100,100))
    l1 <- d[7*(i-1)+7,]$lat1
    l2 <- d[7*(i-1)+7,]$lat2
    l3 <- d[7*(i-1)+7,]$lat3
    l4 <- d[7*(i-1)+7,]$lat4
    lines (c(l1,l1),c(-100,100),col="blue")
    lines (c(l2,l2),c(-100,100),col="pink")
    lines (c(l3,l3),c(-100,100),col="green")
    if (!is.na(l4)){
      lines (c(l4,l4),c(-100,100),col="light blue")
    }
}

dev.off()