Compute PSP copula (Nelsen, 2006, p. 23) is named by the author (Asquith) for the copBasic package and is
$$\mathbf{PSP}(u,v) = \frac{\mathbf{\Pi}}{\mathbf{\Sigma} - \mathbf{\Pi}} = \frac{uv}{u + v - uv}\mbox{,}$$
where \(\mathbf{\Pi}\) is the indpendence or product copula (P
) and \(\mathbf{\Sigma}\) is the sum \(\mathbf{\Sigma} = u + v\). The \(\mathbf{PSP}(u,v)\) copula is a special case of the \(\mathbf{N4212}(u,v)\) copula (N4212cop
). The \(\mathbf{PSP}\) is included in copBasic because of its simplicity and for pedagogical purposes. The name “PSP” comes from Product, Summation, Product to loosely reflect the mathematical formula shown. Nelsen (2006, p. 114) notes that the PSP copula shows up in several families and designates it as “\(\mathbf{\Pi}/(\mathbf{\Sigma}-\mathbf{\Pi})\).” The PSP is undefined for \(u = v = 0\) but no internal trapping is made; calling functions will have to intercept the NaN
so produced for \(\{0, 0\}\). The \(\mathbf{PSP}\) is left internally untrapping NaN
so as to be useful in stressing other copula utility functions herein.
PSP(u, v, ...)
Nonexceedance probability \(u\) in the \(X\) direction;
Nonexceedance probability \(v\) in the \(Y\) direction; and
Additional arguments to pass, which for this copula are not needed, but given here to support flexible implementation.
Value(s) for the copula are returned.
Nelsen, R.B., 2006, An introduction to copulas: New York, Springer, 269 p.
# NOT RUN {
PSP(0.4,0.6)
PSP(0,0)
PSP(1,1)
# }
Run the code above in your browser using DataLab