R/PlotXY.R

In sascha2schroeder/popEye: Analysis of eye-tracking data from reading experiments

PlotXY <- function(exp, subject, trial, sub = F, pdf = NULL, interactive = F) {
    
  # start pdf
  if (sub == F) {
    if (missing(pdf) == T) {
      par(mfrow = c(1, 1), cex = 1.25, oma = c(0, 0, 3, 0))
      if (interactive == T) par(ask = T)
    } else {
      pdf(pdf, width = 16, height = 8.5)
      par(mfrow = c(1, 1), cex = .9, oma = c(0, 0, 2, 0))
    }
  }
 
  tmp <- SelectSubjectTrial(exp, subject, trial)
 
  # smooth again
  xy2 <- SmoothData(SmoothData(tmp$xy))

  # basic plot
  # plot(xy2$x, xy2$y, type = "l", 
  #      xlim = c(min(tmp$meta$stimmat$xs) - 3*exp$setup$font$size,
  #               max(tmp$meta$stimmat$xs) + 3*exp$setup$font$size),
  #      ylim = c(max(tmp$meta$stimmat$ys) + 3*exp$setup$font$size,
  #               min(tmp$meta$stimmat$ys) - 3*exp$setup$font$size), 
  #      main = "XY Plot", xlab = "x Position (px)", ylab = "y Position (px)")
  
  # full plot
  plot(xy2$x, xy2$y, type = "l",
       xlim = c(1, exp$setup$display$resolutionX),
       ylim = c(exp$setup$display$resolutionY, 1),
       main = "XY Plot", xlab = "x Position (px)", ylab = "y Position (px)")
  
  # add saccades
  sac <- tmp$parse[tmp$parse$msg == "SAC", ]
  for (s in 1:nrow(sac)) {
    j <- sac$start[s]:sac$stop[s]
    lines(xy2$x[j], xy2$y[j], type='l', col='navyblue', lwd = 1)
  }

  # add fixations
  fix <- tmp$parse[tmp$parse$msg == "FIX", ]
  for (s in 1:nrow(fix)) {
    points(fix$xs[s], fix$ys[s], col = 'royalblue', pch = 16, cex = 1)
  }

  # add fixation number
  fix$num <- 1:nrow(fix)
  for (s in 1:nrow(fix)) {
    text(fix$xs[s], fix$ys[s] - 30, fix$num[s], col = 'royalblue', cex = .6)
  }

  # add fixation duration
  fix$dur <- fix$stop - fix$start + 1
  for (s in 1:nrow(fix)) {
    text(fix$xs[s], fix$ys[s] + 30, fix$dur[s], col = 'black', cex = .6)
  }
  
  # turn off device  
  if (sub == F) {
    if (missing(pdf) == T) {
      title(paste("Trial", tmp$meta$trialid, sep = " "), 
            outer = T, cex.main = 2)
      par(mfrow = c(1, 1), cex = 1, oma = c(0, 0, 0, 0))
      if (interactive == T) par(ask = F)
    } else {
      title(paste("Trial", tmp$meta$trialid, sep = " "), 
            outer = T, cex.main = 1.75)
      dev.off()
    }
  }

}