plot_heatmapsOnSel didn't work for dataset already annotated
0
0
Entering edit mode
xubeisi • 0
@xubeisi-11767
Last seen 2.5 years ago
United States
 

 

library("MineICA")

dat <- data.frame(matrix(rnorm(10000),ncol=10,nrow=1000))
rownames(dat) <- paste("g", 1:1000, sep="")
colnames(dat) <- paste("s", 1:10, sep="")
## build a data.frame containing sample annotations
annot <- data.frame(type=c(rep("a",5),rep("b",5)))
rownames(annot) <- colnames(dat)
## run ICA
resJade <- runICA(X=dat, nbComp=3, method = "JADE")
## build params
params <- buildMineICAParams(resPath="toy/")
## build IcaSet object
icaSettoy <- buildIcaSet(params=params, A=data.frame(resJade$A), S=data.frame(resJade$S),
dat=dat, pData=annot, alreadyAnnot=TRUE)
params <- icaSettoy$params
icaSettoy <- icaSettoy$icaSet

keepvv = c("type")
keepcc <- 1

plot_heatmapsOnSel(icaSet = icaSettoy, selCutoff = 2, level = "genes", keepVar = keepvv, keepComp=keepcc, doSamplesDendro = TRUE, doGenesDendro = TRUE, heatmapCol = maPalette(low = "blue",high = "red", mid = "yellow", k=44), file = "heatmapWithDendro_zval3.pdf")

 

will report

 

Error in { : task 1 failed - "non-numeric matrix extent"
In addition: Warning message:
executing %dopar% sequentially: no parallel backend registered

 

if I change 

keepvv = c("type1")

the phenotype data won't draw but the heatmap seems fine. 

I tried several debug but still can't figure out how to make it work like using 

data(icaSetCarbayo)

 

Thanks!

 

 

MineICA • 934 views
ADD COMMENT

Login before adding your answer.

Traffic: 765 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6