birth: RJMCMC: Birth of a change point.

Description Usage Arguments Value

View source: R/PWC.R

Description

RJMCMC: Birth of a change point.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
birth(
  partition,
  heights,
  d.spikes,
  hyper.param,
  cvalue,
  kappa,
  mu,
  kappa0,
  k.max,
  sum.of.prob,
  lambda,
  T.min = NULL,
  max.T.min = NA,
  ISI.type = "Gamma",
  do.log,
  which.heights = "independent"
)

Arguments

partition

Partition of time experiment time. A vector of K values where the first entry is 0 and the last is the end of experiment time.

heights

Vector of (K-1) heights corresponding to the partition.

d.spikes

Data frame containing spike sequences.

hyper.param

ISI parameter.

cvalue

parameter c in the RJMCMC, to choose which event to perform in MCMC iteration.

kappa

parameter of priors height

mu

parameter of priors height

kappa0

parameter of priors height

k.max

Maximum number of change points in RJMCMC.

sum.of.prob

sum^k.max_i=0 Poi(k | lambda).

lambda

Parameter of the Poisson distribution.

T.min

Refractory period.

max.T.min

Maximal allowed T.min values from d.spikes.

ISI.type

The ISI distribution.

do.log

Flag, where if do.log is true the calculations are computed on the log scale.

which.heights

Method for prior distribution of heights, either 'independent' or 'martingale'.

Value

new partition and heights


JPNotts/Package documentation built on Oct. 5, 2021, 2:04 p.m.