Create a gist now

Instantly share code, notes, and snippets.

What would you like to do?
This script shows how to load and plot fnirs data (.nirs format)
# install the package if not available with:
# install.packages("R.matlab")
# load the nirs file
data <- readMat("Simple_Probe.nirs")
# get the available names
# format description:
# find the dimensions of data: samples and channels
# number of channels
# create a custom plotting function
myPlot <- function(y){
plot(data$t, y,
type = "l",
xlab = "Time", ylab = "Intensity",
ylim = c(range(data$d)),
main = paste("Channel ", toString(y)))
# add the lines when triggers occur
lines(data$t, data$s*1500, # multiply the trigger value to make it visible
type = "l", xlab = "Time", col = "red")
# apply the plotting function to all channels
apply(data$d[, 1:ncol(data$d)], 2, myPlot)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment