summary.SemiParBIVProbit: SemiParBIVProbit summary

Description

It takes a fitted SemiParBIVProbit object produced by SemiParBIVProbit() and produces some summaries from it.

Usage

## S3 method for class 'SemiParBIVProbit':
summary(object, n.sim=1000, s.meth="svd", sig.lev=0.05, thrs1=0.5, thrs2=0.5,...)

Arguments

object

A fitted SemiParBIVProbit object as produced by SemiParBIVProbit().

n.sim

The number of simulated coefficient vectors from the posterior distribution of the estimated model parameters. This is used to calculate `confidence' intervals for $\rho$.

s.meth

Matrix decomposition used to determine the matrix root of the covariance matrix. See the documentation of mvtnorm for further details.

sig.lev

Significance level used for intervals obtained via posterior simulation.

thrs1,thrs2

Threshold to calculate the correct prediction ratios for the two binary responses. Default value is 0.5. It only works when selection=FALSE.

...

Other arguments.

Value

tableP1Table containing parametric estimates, their standard errors, z-values and p-values for equation 1.
tableP2As above but for equation 2.
tableNP1Table of nonparametric summaries for each smooth component including estimated degrees of freedom, estimated rank, approximate Wald statistic for testing the null hypothesis that the smooth term is zero and corresponding p-value, for equation 1.
tableNP2As above but for equation 2.
nSample size.
rhoEstimated correlation parameter between the two equations.
formula1Formula used for equation1.
formula2Formula used for equation2.
l.sp1Number of smooth components in equation 1.
l.sp2Number of smooth components in equation 2.
t.edfTotal degrees of freedom of the estimated bivariate probit model.
CIrs`Confidence' intervals for $\rho$.
selThis is used for internal calculations.
n.selNumber of selected observations in the sample selection case.
RE,RE.typeThese are used for internal calculations.
massesIf RE=TRUE and RE.type="NP", then it returns the probabilities associated with the biavariate mass points.
table.REIf RE=TRUE then it returns a table representing the estimated parameters for the employed random effect distribution.
table.RJoint absolute frequency table of observed binary value combinations. This as well as all remaining quantities below are not provided when selection=TRUE.
table.PJoint absolute frequency table of predicted binary value combinations.
table.FIt is given by table.P/table.R.
MRBased on the comparison between table.P and table.R, it reports the matching rate as percentage.
P1,P2Marginal probabilities for the two binary outcomes y1=1 and y2=1.
QPS1,QPS2Quadratic probability scores for the two binary responses as suggested by Diebold and Rudebusch (1989). Values are on the interval [0,2], with 0 indicating a perfect fit.
CR1,CR2Correct prediction ratios for the two binary responses as percentage.

Details

As in the package mgcv, based on the results of Marra and Wood (2012), `Bayesian p-values' are returned for the smooth terms. These have better frequentist performance than their frequentist counterpart. Let $\hat{\bf f}$ and ${\bf V}_f$ denote the vector of values of a smooth term evaluated at the original covariate values and the corresponding Bayesian covariance matrix, and let ${\bf V}_f^{r-}$ denote the rank $r$ pseudoinverse of ${\bf V}_f$. The statistic used is $T=\hat{\bf f}^\prime {\bf V}_f^{r-} \hat{\bf f}$. This is compared to a chi-squared distribution with degrees of freedom given by $r$, which is obtained by biased rounding of the estimated degrees of freedom. See Marra (2013) for further details. Note that covariate selection can also be achieved using a single penalty shrinkage approach as shown in Marra and Wood (2011).

References

Diebold F.X. and Rudebusch G.D. (1989), Scoring the Leading Indicators. Journal of Business, 62(3), 369-391. Marra G. (2013), On P-values for Semiparametric Bivariate Probit Models. Statistical Methodology, 10(1), 23-28. Marra G. and Wood S.N. (2011), Practical Variable Selection for Generalized Additive Models. Computational Statistics and Data Analysis, 55(7), 2372-2387. Marra G. and Wood S.N. (2012), Coverage Properties of Confidence Intervals for Generalized Additive Model Components. Scandinavian Journal of Statistics, 39(1), 53-74.

Examples

Run this code

## see examples for SemiParBIVProbit

Run the code above in your browser using DataLab