`R/prml_analysis.R`

`prml_filter.Rd`

PRML filter: Bayes factor of Poisson versus Poisson mixtures.

prml_filter(xs_bn, n_gq = 20, n_per = 100, alpha = 0.5)

xs_bn | A vector. Spike counts of repeated single-stimulus trial data. |
---|---|

n_gq | A number. 20 by default. Number of grids in Gaussion quadrature. |

n_per | A number. 100 by default. Permutation of likihood estimation to obtain the order-invariant estimator. |

alpha | 0.5 by default. The range of the spike counts estimator \( [Y_{0.25}-\alpha IQR,Y_{0.75}+\alpha IQR] \) |

A number. Bayes factor of Poisson versus Poisson mixtures by PRML algorithm.