Tutorial Contents

Phase-plane

Event-specific analysis

See also

Contents

Phase-plane analysis

Most data visualization in neuroscience involves plotting a time-dependent variable such as the membrane potential against time itself. However, in phase-plane analysis, one time-dependent variable is plotted against another time-dependent variable. This provides an alternative visualizationIn dynamical systems analysis the state variables of the underlying differential equations are plotted against each other. This provides more than just an alterative visualization - it allows mathematical investigation of features such as system stability. However, this is more used in computational than experimental neuroscience, and is beyond the scope of this tutorial (and this author). of the shapes of features such as spikes or PSPs, typically by plotting the membrane potential (V) against its time derivative (dV/dt) at each moment in time over the duration of the feature. The phase-plane plot can reveal subtle changes in shape over time that are difficult to pick up in an extended record.

 

a
Phase plot of sine wave chirp
b
3D phase-plane of chirp
Phase plot of a sine wave with gradually increasing frequency (a "chirp"). The colours blue → red indicate increasing time in the source data. a. 2D plot. b. 3d plot.

The phase plot shows a line that spirals clockwise and is colour-coded for time, with blue being earlier, red being later. At the earliest time the main view shows that the membrane potential is 0, but the rate-of-rise of the sine wave is maximal (for that frequency), so the bluest phase plot point is in the middle of the X axis, representing 0 voltage, but above the middle on the Y axis, representing a positive rate of change (voltage derivative: dV/dt). As the sine wave climbs to its maximum, the phase plot moves rightwards indicating increasing voltage, but downwards indicating a decreasing derivative. At the peak of the sine wave the derivative is zero. Thus the maximum and minimum on the X axis (mV) of the phase-plane represent the peaks and troughs of the sine wave, and at these voltages the value on the Y axis (dV/dt) is zero. The sine wave has constant amplitude, so the X maxima and minima stay the same. The the maximum and minimum on the Y axis represent the most rapidly changing voltages, which occur when the sine wave itself crosses the zero value. However, the sine wave is a chirp of increasing frequency, so the Y maxima and minima increase with each cycle. Hence the spiral.

The phase progression can be seen more easiy in 3D:

What do real neurophysiological data look like?

Follow the line clockwise from the blue to the red end. The line first moves to the right, indicating depolarization, and then back to the left, indicating hyperpolarization. There is an initial small hump which is generated by the EPSP. The EPSP almost flattens out at its peak, which is why the hump is hump-shaped! - the derivative almost returns to zero at the flattish top of the EPSP. Then the spike takes off. The maximum rate of rise is about 51 mV/ms, which occurs when the spike has a membrane potential of about -8 mV (you can see the values in the X and Y readouts when you hover the mouse over the plot). The spike peaks at about +15 mV, when the derivative is 0. are two loops visible in the plot. The maximum fall rate is about -26.5 mV/ms, which occurs at a membrane potential of -3 mV.

a
spikes
b
3D phase-plane plot
Phase plot of an intracellular recording. a. The recording. b. A 3-D plot of the recording.

The phase plot is now generated from all the data in the file. The 7 large loops represent the 7 spikes displayed in the record. The small loops at the left represent subthreshold events. The EPSP underlying the single spike failure is the largest of these small loops, drawn in green-yellow (rotating the graph will help visualize this event). The main display shows that the 1st and 4th spikes take off rather late in the PSP, and the blue and yellow loops in the phase plane initially follow the failure trajectory, but just manage to make threshold. A small orange loop is visible in the 3D display, and this is caused by the small EPSP following immediately after the 5th spike in the record.

You can visualize these subthreshold events more clearly by zooming the 3D display.

Event-specific analysis

Dataview allows phase-plane analysis of discontiguous sections of data if they are defined by events. This can be useful if the regions of interest are widely dispersed in a recording - it avoids the computational expense of dealing with unimportant data.

The file file spike shape already has suitable events in channel a that were defined using simple threshold-crossing.

The phase plot is similar to the one from the whole recording, but it only includes data within the events.

See also

Spike shape: Tool for extracting numerical parameters describing spike shape.