View source: R/kernel_methods.R
Phi_fn | R Documentation |
This function calculates the Φ matrix, which is a second order Stein operator applied to a polynomial. See South et al (2020) for further details. This function is not required for estimation but may be useful when evaluation samples are not initially available since estimators using heldout samples are of the form mean(f - f_hat) + b[1] where f_hat = K0*a + Phi*b for heldout K0 and Phi matrices.
Phi_fn(samples, derivatives, polyorder = NULL, apriori = NULL)
samples |
An N by d matrix of samples from the target |
derivatives |
An N by d matrix of derivatives of the log target with respect to the parameters |
polyorder |
(optional) The order of the polynomial to be used in the parametric component, with a default of 1. We recommend keeping this value low (e.g. only 1-2). |
apriori |
(optional) A vector containing the subset of parameter indices to use in the polynomial. Typically this argument would only be used if the dimension of the problem is very large or if prior information about parameter dependencies is known. The default is to use all parameters 1:d where d is the dimension of the target. |
An N by Q matrix (where Q is determined by the polynomial order and the apriori).
Leah F. South
South, L. F., Karvonen, T., Nemeth, C., Girolami, M. and Oates, C. J. (2020). Semi-Exact Control Functionals From Sard's Method. https://arxiv.org/abs/2002.00033
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.