Description Usage Arguments Details Value Examples
View source: R/riskscore_fct.R
Main input function for SurvRank.
1 2  | 
cv.ob | 
 output of the   | 
data | 
 same list used as input in   | 
th | 
 Defaults to 0.5. Threshold of used features. th=0.5 majority vote approach  | 
surv.tab | 
 Defaults to c(0.5). Calculates for selected features survival curves.   | 
f | 
 Defaults to NA. ranking approach function. One of   | 
fix.var | 
 Defauts to NA. not NA, fixed number of features is calculated  | 
list.t | 
 Defauls to "weighted". Which toplist should be chosen? Possible choices are "weighted", "unweighted", "rank", "top1se","cluster" or "final"  | 
ncl | 
 Defaults to 1. Number of clusters for parallel execution.  | 
plt | 
 Default=F. Should plot of survival curves be generated?  | 
... | 
 arguments that can be passed to underlying functions, not used now  | 
details to follow
Output of the riskscore_fct, basically a list containing the following elements
 | 
 toplist of features that have been chosen  | 
 | 
 Matrix of survival AUCs with fixed number of features, but not fixed features!! (could also be calculated before)  | 
 | 
 cox model output for selected features, according to   | 
 | 
 AIC criterion of cox model  | 
 | 
 summary object of the fitted cox model  | 
 | 
 concordance measure of fitted cox model  | 
 | 
 survfit object of the cox model)  | 
 | 
 predictions of the cox model (fitted values)  | 
 | 
 survfit object according to   | 
 | 
 Cox model on the groups generated by   | 
 | 
 surfdiff: Tests if there is a difference between two or more survival curves using the G-rho family of tests, or for a single curve against a known alternative  | 
Additionally two plots are generated: if f is not NA, a boxplot of the survival AUCs, averaged for cross-validation iterations. The second plot shows the resulting survival curves according to surv.tab.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17  | ## Simulating a survival data set
N=100; p=10; n=3
x=data.frame(matrix(rnorm(N*p),nrow=N,p))
beta=rnorm(n)
mx=matrix(rnorm(N*n),N,n)
fx=mx[,seq(n)]%*%beta/3
hx=exp(fx)
ty=rexp(N,hx)
tcens=1-rbinom(n=N,prob=.3,size=1)
y=Surv(ty,tcens)
data=list()
data$x<-x; data$y<-y
out<-CVrankSurv_fct(data,2,3,3,fs.method="cox.rank")
## Using the weighted toplist
risk<-riskscore_fct(out,data,list.t="weighted")
## Selected names
risk$selnames
 | 
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.