R/plotP300evg.R

In Dardare/prjct:

plotP300evg <- function(target, nont){
  target <- tmp[[1]]
  nontarget <- tmp[[2]]
  samples1 <- target
  samples2 <- nontarget


  pl <- function(A, title){
    A1 <- abind(A, along=3)
    AM <- rowMeans(A1, dims=2)
    rownames(AM) <- -249:250/500


    print(
      ggplot(melt(AM), aes(x=Var1, y=value, colour=as.factor(Var2))) + geom_line() + labs(title=title)
    )
  }

  pl(samples1[[3]], 'P300')
  pl(samples2[[3]], 'NON P300')
}