Hello, Where is the problem in this programming. I want two way matrix, but it gives problem. The program is given below:
rep = 2
genes = 5
pred = c()
iter = array (dim = c(rep, pred))
pred =
array(dim = c(1, genes))
m = c()
l = 1
w = 1
for(m in
1:rep) {
t
= time data
c
= cens
tab
= cbind(t, c)
#
Resampling Time and Censoring
cc =
sample(c, replace = TRUE) #
Bootstrap WR Censoring Data
tt =
sample(t, replace = TRUE) #
Bootstrap WR Time Data
xx =
Surv(tt, cc)
#
Survival Times with Censored
P = predic[,1:5]
# First 5 Rows of predic
n = 310
#
Every Row has 310 values
ma =
matrix(P, n) #
P cross n, matrix
boot <- sample(seq(n), replace=TRUE) #
Bootstrap WR in seq(n)
genes5 = ma[boot, ] # 5
genes Bootstrap
Matrix
for(p in 1: genes){
i
= c()
for(i
in w:l) {
z
= genes5[,w: l]
y
= coxph(xx ~ z)
}
Pval
= summary(y) $ coeff[, 5]
if( any(Pval <=
0.01) ) sign <- 1 else sign <- 0
pred[1,
p] = sign[1]
}
iter [m, 1:5]
= pred[1: 5]
}
Thanks in advance.
Regards
Fazli Raziq
[[alternative HTML version deleted]]
______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
