calc.burst.distributions: calculate and plot burst featues distributions
In IGM.MEA: IGM MEA Analysis

Description Usage Arguments Details Note Author(s) Examples

View source: R/burstDistributions_functions.R

The function calculates normalized distributions of selected bursting features and plots distribution graphs of all treatments in a recording. The function also prints csv output in the /Analysis directory for downstream stats such as permutation test of treatment labels

1
2
3

calc.burst.distributions(s, minVals = 1, xlimit = 25, binsInSec = 5, 
feature = "non", filterValuesByMin = 0, minValues = 0, perWell = 0, outputdir = getwd(),
min.electrodes=4, timeStamp="DATE_TIME")

`s`	MEA data structure
`minVals`	minimum values number per electrode, electrodes with a smaller number of values than that are discarded
`xlimit`	max limit of values, for example: xlimit = 25 for IBI analysis means that IBIs longer than 25 seconds will not be part of distribution calculations
`binsInSec`	how many bins to cut each of the segments. For example: IBI analysis has 25 seconds as xlimit, to analyse in a 0.1 sec resolution binsInSec should be set to 10, for 1 sec resolution set binsInsec to 1
`feature`	what feature to analyze, options are "IBI", "ISI, "nspikesInBurst", "duration", "spikesDensityInBurst"
`filterValuesByMin`	should analysis disregard values with lower then filterValuesByMin number of values ? (0/1, default is 0). For example, if set to 1 for duration analysis, should analysis consider also bursts shorter than filterValuesByMin ?
`minValues`	disregards values with lower then filterValuesByMin , only if filterValuesByMin set to 1
`perWell`	should distribution analysis be performed by testing treatment differences on well level means (1) or electrode level means(0)
`outputdir`	output directory
`min.electrodes`	minimum electrodes for an active well
`timeStamp`	time stamp for the output files

Plot distributions calculates normalized distributions of bursting features. 'Normalized distribution' are a way to handle biases caused by noisy electrodes/wells. The function will calculate a normalized histogram (values 0-1) of each feature for each electrode. Next, it will average histogram values either per well and then average all wells per treatment, or directly per treatment. All comparisons between treatments will be then made by plotting the normalized histograms of each treatment and running a K-S test between them.

Output is a made of: 1) Plots of all selected burst features distributions. 2) CSV files ending with _distributions.csv that harbor all electrodes per treatment for all the recordings loaded in the IGM.MEA pipeline for a specific MEA plate

Sahar Gelfman

# Load exapmple of recording Robject (MEA data structure)
data("S")  
feature="IBI"; 
#calc.burst.distributions(S, minVals = 15, xlimit = 20, binsInSec = 5, 
#feature = feature, perWell = 0, outputdir = "/Analysis")