PreMean <- matrix(ncol=3, nrow=2, byrow=TRUE)<br>
PreMean[1, ] <- c(1, -1, 0)<br>
PreMean[2, ] <- c(-0.8, 0, 0.8)<br>
PreZeta <- c(0.4, 1.8, 1)<br>

GenExpression <- function(x, cancertype, sdv) {<br>
    rnorm(1, PreMean[x,cancertype], sdv)<br>
}<br>

GenSurvival <- function(funtype, cancertype, x1, x2, cancerprob) {<br>
#    rm(.Random.seed)<br>
	 S <- runif(1, 0, 1)<br>
	 
    if(cancertype==1) {  h <- 0.3*x1 + 0.4*x2 - 0.1  }<br>
    if(cancertype==2) {  h <- 0.5*(x1-x2)^2 + 1      }<br>
	 if(cancertype==3) {  h <- 0.006*x1^3 - 0.2*x2    }<br>
	 
	 if(funtype=="1") {  -(log(S)/exp(-h))  }<br>
	 
	 #survival function is Weibull<br>
	 exp((-log(S)/exp(-h))^(1/PreZeta[cancertype]))<br>
}<br>