Read in a directory of spike times from Marla Feller and create a "spikes" data structure.

Share:

Description

Marla Feller has provided MEA data from an hexagonal array. The data from each electrode is stored in a separate file. This function reads those text files and creates a "spikes" data structure.

Usage

1
2
feller.read.spikes(filename, ids=NULL, time.interval=1,
                   beg=NULL, end=NULL)

Arguments

filename

Name of the text file to be read in.

ids

Which electrodes should be kept.

time.interval

How often to estimate the firing rate.

beg

Start time of the recording.

end

End time of the recording.

Details

The data from each electrode is not in seconds; divide by 20000 to get the time in seconds.

Value

Return the data structure 's'.

Author(s)

Stephen Eglen

See Also

jay.read.spikes sanger.read.spikes

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
data.file <- system.file('examples', '2002-11-26-1', package='sjemea')
s <- feller.read.spikes(data.file)
ejcmealayout.check()
fourplot(s)

## This file contains the picture of the array for comparison
system.file("examples/ejc_ElectrodeNumbering.jpg",package="sjemea")

## next two examples show electrodes being rejected/accepted.
s <- feller.read.spikes(data.file, ids=c('-', 'a7a', 'b1a', 'h2a'))
s <- feller.read.spikes(data.file, ids=c(     'a7a', 'b1a', 'h2a'))