# R/VuongClarke.r In SemiParSampleSel: Semi-Parametric Sample Selection Modelling with Continuous or Discrete Response

#### Documented in VuongClarke

```VuongClarke <- function(obj1,obj2,sig.lev=0.05){

l1 <- obj1\$logLik
l2 <- obj2\$logLik

if(l1==l2) stop("The two competing models have identical log-likelihoods!")

p1 <- obj1\$t.edf
p2 <- obj2\$t.edf
n  <- obj1\$n
li12 <- obj1\$fit\$l.par - obj2\$fit\$l.par
w  <- sqrt(var(li12)*(n-1))
vt <- (l1-l2-(p1-p2)/2*log(n))/w
li12b <- li12-(p1-p2)/(2*n)*log(n)
b     <- sum(li12b>0)

if(abs(vt) <= qnorm(1 - sig.lev/2)) dvt <- " Vuong's test: it is not possible to discriminate between the two models."
if(vt >  qnorm(1 - sig.lev/2))      dvt <- " Vuong's test: Model 1 is preferred over Model 2."
if(vt < -qnorm(1 - sig.lev/2))      dvt <- " Vuong's test: Model 2 is preferred over Model 1."

db <- "Clarke's test: it is not possible to discriminate between the two models."

if(b >= n/2){ pvalue <- 2 * (1 - pbinom(b - 1, n, 0.5))
if(pvalue <= sig.lev) db <- "Clarke's test: Model 1 is preferred over Model 2."
}
if(b < n/2) { pvalue <- 2 * (pbinom(b, n, 0.5))
if(pvalue <= sig.lev) db <- "Clarke's test: Model 2 is preferred over Model 1."
}

cat("\n",dvt,"\n",sep="")
cat(db,"\n\n",sep="")

}
```

## Try the SemiParSampleSel package in your browser

Any scripts or data that you put into this service are public.

SemiParSampleSel documentation built on May 29, 2017, 7:54 p.m.