# NOT RUN {
set.seed(332)
#generate some data
x <- matrix(rnorm(50*20), ncol=20)
y <- 10 + svd(x[1:10,])$v[,1] + .1*rnorm(20)
ytest <- 10 + svd(x[1:10,])$v[,1] + .1*rnorm(20)
censoring.status <- sample(c(rep(1,15), rep(0,5)))
censoring.status.test <- sample(c(rep(1,15), rep(0,5)))
featurenames <- paste("feature", as.character(1:50), sep="")
data <- list(x=x,
y=y,
censoring.status=censoring.status,
featurenames=featurenames)
data.test <- list(x=x,
y=ytest,
censoring.status=censoring.status.test,
featurenames=featurenames)
a <- superpc.train(data, type="survival")
aa <- superpc.cv(a, data)
fit.red <- superpc.predict.red(a,
data,
data.test,
threshold=.6)
fit.redcv <- superpc.predict.red.cv(fit.red,
aa,
data,
threshold=.6)
# }
Run the code above in your browser using DataLab