Spike times

A CSV file is used to store at least two columns: Channel containing the channel name, and Time containing the time (in seconds) at which a spike was detected on that channel. The rows of the CSV do not need to be ordered by channel, but the times for each channel should be ordered, earliest spike first.

Channel positions

A second CSV contains the (x,y) spatial location each channel. If the recording is from a multi-well array, then an extra Well column determines the name of the well. (If no Well column is provided, the package assumes a plate with one well.)

Plate information

The package by default has information about just two plates, the Axion 12 and 48 well systems. As shown in the examples below, in addition we need to provide information regarding the plate. If any information is missing, the system will revert to defaults.

Example 1: a single-well MEA recording

Here there is just a single MEA recording, from a hexagonal MEA (Wong et al., 1993). As there is only one well, the Well information is absent from the position file. The data files are provided within the package, and the top of each file looks as follows:

times = system.file("extdata/textreader/wong1993_p0.times", package="meaRtools")
pos = system.file("extdata/textreader/wong1993_p0.pos", package="meaRtools")
show_top_file(times)

Channel,Time
c1,14.51755
c1,14.56795
c1,78.75835
c1,78.7723
c1,78.7951
c1,78.80975
c1,78.83945
c1,78.86245
c1,78.90175

show_top_file(pos)

"Channel","x","y"
"c1",70,-242.48
"c2",0,-242.48
"c3",0,-242.48
"c4",-140,-242.48
"c5",175,-181.86
"c6",175,-181.86
"c7",105,-181.86
"c8",-35,-181.86
"c9",-35,-181.86

hex_platelayout = list(n_well = 1, #number of wells 
                        wells = c("w1"), #names of those wells.
                        n_well_r = 1, # number of wells / row
                        n_well_c = 1, # number of wells / col
                        layout = c(1, 1), # layout when plotting
                        n_elec_r = 8,
                        n_elec_c = 8,
                        xlim = c(-400, 400), # xlimits for plotting
                        ylim = c(-400, 400), # ylimits for plotting
                        spacing = 50,  # distance (um) separating electrodes
                        corr_breaks = 0 # vector of correlation distances
                        )

add_plateinfo("hex-1well", hex_platelayout)

## [1] "Axion 48 well" "Axion 12 well" "hex-1well"

s = read_spikelist_text(times, pos, array="hex-1well")
meaRtools:::.plot_mealayout(s$layout, use_names = TRUE, cex=0.3)

meaRtools:::.plot_meanfiringrate(s, main = "Mean Firing Rate by Plate (Hz)")

Example 2: a synthetic multi-well recording

This second dataset is a composite of six recordings from P9 and P11 mouse retina (Demas et al, 2013), synthesised to make a 6-well plate.

demas_platelayout = list(n_well = 6,
                        wells = paste0("w", 1:6),
                        n_well_r = 2,
                        n_well_c = 3,
                        layout = c(3, 2),
                        n_elec_r = 8,
                        n_elec_c = 8,
                        xlim = c(-100, 7200),
                        ylim = c(0, 6000),
                        spacing = 200,
                        corr_breaks = 0
                        )
add_plateinfo("demas-6well", demas_platelayout)

## [1] "Axion 48 well" "Axion 12 well" "hex-1well"    
## [4] "demas-6well"

times = system.file("extdata/textreader/demas.times", package="meaRtools")
pos = system.file("extdata/textreader/demas.pos", package="meaRtools")
s = read_spikelist_text(times, pos, array="demas-6well")
meaRtools:::.plot_mealayout(s$layout, use_names = TRUE, cex=0.3)