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.