Plot likelhiood ratio test statistics

Description

Plot likelhiood ratio test statistics from output of superpc.predict

Usage

superpc.plot.lrtest(object.lrtestcurv)

Arguments

Details

Value

Author(s)

Eric Bair and Robert Tibshirani

See Also

Examples

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)