gofPIOSTn tests a 2 dimensional dataset with the PIOS test for a copula. The possible copulae are "gaussian", "t", "gumbel", "clayton" and "frank". The parameter estimation is performed with pseudo maximum likelihood method. In case the estimation fails, inversion of Kendall's tau is used.gofPIOSTn(copula, x, M = 1000, param = 0.5, param.est = T, df = 4, df.est = T,
margins = "ranks", m = 1, execute.times.comp = T)"gaussian", "t", "clayton", "gumbel" and "frank".TRUE or FALSE. TRUE means that param will be estimated with a maximum likelihood estimation."t"-copula.df shall be estimated. Has to be either FALSE or TRUE, where TRUE means that it will be estimated."ranks", which is the standard approach to convert data in such a case. Alternatively can the following distributions be spM is at least 100.class gofCOP with the componentsm out of the data. The test compares then the pseudo likelihood of the data in each block with the overall parameter and with the parameter by leaving out the data in the block. By this procedure can be determined if the data in the block influence the parameter estimation significantly. The test statistic is defined as
$$T = \sum_{b=1}^B \sum_{i=1}^m [l{U_i^b;\theta_n } - l{U_i^b;\theta_n^{-b} }]$$
with the pseudo observations $U_{ij}$ for $i = 1, \dots,n$; $j = 1, \dots,d$ and
$$\theta_n = \arg \min_{\theta} \sum_{i=1}^n l(U_i; \theta)$$
and
$$\theta_n^{-b} = \arg \min_{\theta} \sum_{b^{'} \neq b}^M \sum_{i=1}^m l(U_i^{b^{'}}; \theta), b=1, \dots, M.$$
The approximate p-value is computed by the formula
$$p = \frac{1}{M} \sum_{b = 1}^M \mathbf{I}(|T_b| \geq T).$$
The applied estimation method is the two-step pseudo maximum likelihood approach, see Genest and Rivest (1995).data = cbind(rnorm(100), rnorm(100))
gofPIOSTn("gaussian", data, M = 20)Run the code above in your browser using DataLab