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

Re: Using Mathematica for Electrophysiology Data Analysis

• To: mathgroup at smc.vnet.net
• Subject: [mg110421] Re: Using Mathematica for Electrophysiology Data Analysis
• From: "Kevin J. McCann" <Kevin.McCann at umbc.edu>
• Date: Thu, 17 Jun 2010 02:04:43 -0400 (EDT)
• References: <hva62h$qf2$1@smc.vnet.net>

Francois,

My usual approach to histograms is to first let Mathematica do it, then, once I
have a sense of the data structure, I use BinCounts to pick the bins
that I want and get the counts in them.

As far as the multivariate Gaussian fit to a bimodal distribution. Hmmm.
I am not sure why I would think that that is a good idea. Perhaps a
relevance vector machine approach would provide a better fit.

Kevin

Fran=E7ois Blumenfeld wrote:
> 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
>