FieldTrip

Note that this reference documentation is identical to the help that is displayed in Matlab when you type “help ft_spike_jpsth”.

  FT_spike_JPSTH computes the joint peristimulus histograms for spiketrains
  and a shift predictor (for example see Aertsen et al. 1989).
 
  The shift predictor is computed in consecutive trials in a symmetric way.
  For example, we compute the jpsth for chan 1 in trial 1 versus chan 2 in
  trial 2, but also for chan 1 in trial 2 versus chan 2 in trial 1. This
  gives (nTrials-1)*2 jpsth matrices for individual trials. Picking
  consecutive trials and computing the shift predictor in a symmetric way
  ensures that slow changes in the temporal structure do not affect the
  shift predictor.
 
  Use as
    [jpsth] = ft_spike_jpsth(cfg,psth)
 
  The input PSTH should be organised as the input from FT_SPIKE_PSTH,
  FT_SPIKE_DENSITY or TIMELOCKANALYSIS containing a field PSTH.trial and
  PSTH.time. In any case, one is expected to use cfg.keeptrials = 'yes' in
  these functions.
 
  Configurations:
    cfg.normalization    = 'no' (default), or 'yes'.  If requested (see cfg.normalization), the joint
                            psth is normalized as in van Aertsen et al. (1989), by
                            D(u,v)/sqrt(D(u,u)*D(v,v) w
                            here D(u,v) is the difference of the average jpsth with the predicted jpsth
                            (see ref. for details), giving a quantity between -1 and 1.
                            Since this method normalizes by the mean across all trials, it can be
                            confounded by latency drifs over trials.
    cfg.shiftpredictor   = 'no' (default) or 'yes'. If 'yes', then JPSTH.AVG, JPSTH.VAR
                            AND JPSTH.DOF will apply to the shiftpredictor, not the jpsth
                            proper.
    cfg.channelcmb       =  Mx2 cell-array with selection of channel pairs (default = {'all' 'all'}),
                            see FT_CHANNELCOMBINATION for details
    cfg.trials           = 'all' (default) or numerical or logical array of to be selected trials.
    cfg.latency          = [begin end] in seconds, 'maxperiod' (default), 'prestim'(t<=0), or
                           'poststim' (t>=0).
    cfg.keeptrials       = 'yes' or 'no' (default).