[Date Index] [Thread Index] [Author Index]

Using Mathematica for Electrophysiology Data Analysis

• To: mathgroup at smc.vnet.net
• Subject: [mg110386] Using Mathematica for Electrophysiology Data Analysis
• From: François Blumenfeld <hypotyposis at me.com>
• Date: Wed, 16 Jun 2010 05:38:25 -0400 (EDT)

Hi,

I'm playing around with Mathematica to analyze some electrophysiology  
data; this is on my recreation time, so apologies if the following  
questions are very basic.

I've imported data (pairs of time points and amplitudes) from pClamp  
into Mathematica as lists. The lists are pretty huge (around 150,000  
elements per list-I've separated times and amplitudes in two lists).
I've created a histogram of the amplitudes list, letting Mathematica  
decide which binning method to use and the width of the bin. Without  
user input, it actually did a great job at picking bin width and  
binning method; but I don't know how to retrieve that information:  
i.e. what bin width was picked and what binning method used for a  
particular histogram?

My second question is about fitting a multivariate Gaussian  
distribution to a multimodal (usually bimodal) histogram. I've payed  
around with NonLinearModelFit, but for some reason I don't think I'm  
doing it correctly. My third question is related to this: how do I  
retrieve the processed data from the histogram to construct a  
frequency table? I'm guessing that this might in fact be the  
intermediate step I need to complete before being able to fit a curve  
to the histogram data.

Any suggestions welcome!

Francois Blumenfeld-Kouchner