n.CT <- 100
n.CC <- 50
nevent.C <- 100
n.ECp <- 200
nevent.ECp <- 180
accrual <- 16
out.mevent.CT <- 6
out.mevent.CC <- 6
driftHR <- 1
cov.C <- list(list(dist="norm",mean=0,sd=1,lab="cov1"),
list(dist="binom",prob=0.4,lab="cov2"))
cov.cor.C <- rbind(c( 1,0.1),
c(0.1, 1))
cov.EC <- list(list(dist="norm",mean=0,sd=1,lab="cov1"),
list(dist="binom",prob=0.4,lab="cov2"))
cov.cor.EC <- rbind(c( 1,0.1),
c(0.1, 1))
cov.effect <- c(0.9,0.9)
indata <- trial.simulation.t2e(
n.CT=n.CT, n.CC=n.CC, nevent.C=nevent.C,
n.ECp=n.ECp, nevent.ECp=nevent.ECp, accrual=accrual,
out.mevent.CT, out.mevent.CC, driftHR,
cov.C=cov.C, cov.cor.C=cov.cor.C,
cov.EC=cov.EC, cov.cor.EC=cov.cor.EC, cov.effect=cov.effect)
n.EC <- 50
method.whomatch <- "conc.treat"
method.matching <- "optimal"
method.psorder <- NULL
out.psmatch <- psmatch(
study~cov1+cov2, data=indata, n.EC=n.EC,
method.whomatch=method.whomatch, method.matching=method.matching,
method.psorder=method.psorder)
indata.match <- rbind(indata[indata$study==1,],indata[out.psmatch$subjid.EC,])
method.borrow <- list(list(prior="cauchy",scale=2.0),
list(prior="normal",scale=0.5))
commensurate.t2e(
survival::Surv(time,status)~cov1+cov2,data=indata.match,
method.borrow=method.borrow,chains=1,iter=100)
Run the code above in your browser using DataLab