superpc.plot.lrtest {superpc} | R Documentation |
Plot likelhiood ratio test statistics from output of superpc.predict
superpc.plot.lrtest(object.lrtestcurv)
object.lrtestcurv |
Output from superpc.lrtest.curv |
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) status<- sample(c(rep(1,30),rep(0,10))) status.test<- sample(c(rep(1,30),rep(0,10))) featurenames <- paste("feature",as.character(1:1000),sep="") data<-list(x=x,y=y, status=status, featurenames=featurenames) data.test<-list(x=x,y=ytest, status=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)