... | ... |
@@ -47,10 +47,10 @@ Z <- prepareTensorfromList(Multi,10) |
47 | 47 |
Z <- aperm(Z,c(2,1,3)) |
48 | 48 |
require(RTCGA.clinical) |
49 | 49 |
Clinical <- list(BLCA.clinical,BRCA.clinical,CESC.clinical,COAD.clinical) |
50 |
-Multi_sample <- list(BLCA.rnaseq[seq_len(100),1,drop=F], |
|
51 |
- BRCA.rnaseq[seq_len(100),1,drop=F], |
|
52 |
- CESC.rnaseq[seq_len(100),1,drop=F], |
|
53 |
- COAD.rnaseq[seq_len(100),1,drop=F]) |
|
50 |
+Multi_sample <- list(BLCA.rnaseq[seq_len(100),1,drop=FALSE], |
|
51 |
+ BRCA.rnaseq[seq_len(100),1,drop=FALSE], |
|
52 |
+ CESC.rnaseq[seq_len(100),1,drop=FALSE], |
|
53 |
+ COAD.rnaseq[seq_len(100),1,drop=FALSE]) |
|
54 | 54 |
#patient.stage_event.tnm_categories.pathologic_categories.pathologic_m |
55 | 55 |
k <- c(770,1482,773,791) |
56 | 56 |
#patient.bcr_patient_barcode |
... | ... |
@@ -63,7 +63,8 @@ HOSVD <- computeHosvd(Z) |
63 | 63 |
cond<- attr(Z,"sampleData") |
64 | 64 |
index <- selectFeatureProj(HOSVD,Multi,cond,de=1e-3,input_all=3) #Batch mode |
65 | 65 |
head(tableFeatures(Z,index)) |
66 |
-genes <-unlist(lapply(strsplit(tableFeatures(Z,index)[,1],"|",fixed=T),"[",1)) |
|
66 |
+genes <-unlist(lapply(strsplit(tableFeatures(Z,index)[,1],"|", |
|
67 |
+ fixed=TRUE),"[",1)) |
|
67 | 68 |
``` |
68 | 69 |
|
69 | 70 |
# Enrichr |
... | ... |
@@ -86,7 +86,7 @@ input_all <- selectSingularValueVectorLarge(HOSVD,cond,input_all=c(12,1)) |
86 | 86 |
Finally, we perform the following function to select features in individual |
87 | 87 |
omics profiles in an interative mode |
88 | 88 |
``` |
89 |
-HOSVD$U[[1]] <- HOSVD$U[[2]] #selectFeatureSquareのHOSVD$U[[1]]を[[2]]にすればこれは不要 |
|
89 |
+ |
|
90 | 90 |
index_all <- selectFeatureSquare(HOSVD,input_all,CLL_data, |
91 | 91 |
de=c(0.5,0.1,0.1,1)) |
92 | 92 |
``` |