问利用R生成右删失Copula数据

EN

    library(methods)
    library(copula)
    library(muhaz)
    library(foreach)
    library(survival)
    library (MASS)

tau = 0.2
nsample = 100
err_dist = "Normal"
generate.dat <- function(tau,nsample,err_dist)
  {
  C <- rexp(nsample, rate = 0.8) # Censoring

  cc <- claytonCopula(copClayton@iTau(tau))
  u <- rCopula(nsample, cc)

  e <- qnorm(u,mean=0,sd=1) # Error term

  b1 <- 0.5 ; b2 <- 0.2  # using the old betas 
  beta <- matrix(c(b1, 0, 0, b2), 2) # as Beta got to be 2x2

  RHO <- 0  # Indpendence
  MU1 <- 0; S1 <- 1 # Std. Normal
  MU2 <- 0; S2 <- 1 # Std. Normal
  MU <- c(MU1, MU2) # mean vector
  SIGMA <- matrix (c(S1^2, S1*S2*RHO, S1*S2*RHO, S2^2), 2) # Var-Covar matrix
  Z <- mvrnorm(nsample, mu=MU, Sigma=SIGMA)

  T <- exp(-Z%*%beta + e)
  T1 <- T[,1]
  T2 <- T[,2]

  x1 <- pmin(T1,C)
  x2 <- pmin(T2,C); 
  d1 <- ifelse(T1<=C, 1, 0)
  d2 <- ifelse(T2<=C, 1, 0)

  dat <- data.frame(x1, x2, d1, d2)  
  return(dat)
  print(dat)
}