superpc.plot.lrtest {superpc} | R Documentation |
Plot likelhiood ratio test statistics from output of superpc.predict
superpc.plot.lrtest(object.lrtestcurv, call.win.metafile = FALSE)
object.lrtestcurv |
Output from superpc.lrtest.curv |
call.win.metafile |
For use by PAM Excel interface |
Eric Bair and Robert Tibshirani
set.seed(332) #generate some data x<-matrix(rnorm(1000*40),ncol=40) y<-10+svd(x[1:60,])$v[,1]+ .1*rnorm(40) ytest<-10+svd(x[1:60,])$v[,1]+ .1*rnorm(40) censoring.status<- sample(c(rep(1,30),rep(0,10))) censoring.status.test<- sample(c(rep(1,30),rep(0,10))) featurenames <- paste("feature",as.character(1:1000),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<- superpc.predict(a, data, data.test, threshold=1.0, n.components=1, prediction.type="continuous") bb<-superpc.lrtest.curv(a,data,data.test) superpc.plot.lrtest(bb)