Phi_fn: Phi matrix calculation

Description Usage Arguments Value Author(s) References

View source: R/kernel_methods.R


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)



An N by d matrix of samples from the target


An N by d matrix of derivatives of the log target with respect to the parameters


(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).


(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.

ZVCV documentation built on July 2, 2020, 2:38 a.m.

Related to Phi_fn in ZVCV...