All functions

ClusterLhnData()

A function for fitting Shahar's LHN data to the linear-nonlinear-poisson model.

DetectVariableBaselineUsingBayesianModelSelection()

A function that uses model selection to detect variability in the baseline firing rate of a cell. The algorithm assumes Poisson spiking in the baseline period and compares the posterior on a single rate explaining the data to having two rates.

EstimateResponseDelayByMomentMatching()

Estimates the response delay by cross correlating the observed response with the odor profile convolved with an exponential filter.

addnacols()

Reorder odour response matrix adding nas as necessary

anyid2shortid() anyid2stack() anyid2longid()

Convert any identifier / file path to Shahar's short id or stack id

baseline_subtract_allfreqs()

Subtract baseline spike rate from list of smoothed psth data

basesubtract_heatmap_cor_dist()

An attaempt to create a basesubtraction option first step at calculating lifetime sparsness

common_odours_for_cells()

Find the common set of odours for which all cells have trials

createSummarySpikesMat() createSummarySpikesArray()

Convert raw spike summary array into matrix or array without cells/odours missing data

create_raw_summary_array() create_simple_summary_arrays()

Create the raw summary array for all spikes

heatmap_anatomy()

heatmap for set of cells on nblast anatomy distance

heatmap_cor_dist() spike_cor_dist()

heatmap for set of cells and odours based on correlation distance

jet.colors()

Return a colour palette function

physplit.analysis

Packaged and versioned analysis functions for Shahar's cells

plotcellsf()

Plot spike/Vm trace for a cell using defaults appropriate for Frechter et al

poissonTestOdoursSF()

Carry out Poisson test on absolute number of spikes in odour response

prop.ci()

Approximate (1-alpha)100% confidence interval for proportion of a population

required.sample.size()

Estimate sample size to find population proportion with given tolerance

sample_finite_population()

Sample from finite population with known number of true positives

truepos_given_sample() summary(<truepos>)

Estimate distribution of true positives given sampling resuts