# \donttest{
n <- 100
rho <- .6
cens <- 0.25
nu <- 4
set.seed(20200527)
w <- cbind(1,runif(n,-1,1),rnorm(n))
x <- cbind(w[,1:2])
family <- "T"
c <- qt(cens, df=nu)
sigma2 <- 1
beta <- c(1,0.5)
gamma<- c(1,0.3,-.5)
gamma[1] <- -c*sqrt(sigma2)
data <- rHeckman(x,w,beta,gamma,sigma2,rho,nu,family=family)
# }
Run the code above in your browser using DataLab