dtsoucas / DWLS

Gene expression deconvolution using dampened weighted least squares

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

Error when using buildSignatureMatrixMAST()

RK1912 opened this issue · comments

commented

Hi
I am using the buildSignatureMatrixMAST() function and I am getting the following error :

This is the traceback :

Error in dimnames(x) <- dn : 
  length of 'dimnames' [2] not equal to array extent 
3.
`colnames<-`(`*tmp*`, value = unique(id)) 
2.
`colnames<-`(`*tmp*`, value = unique(id)) at Deconvolution_functions.R#384
1.
buildSignatureMatrixMAST(scdata = data, id = as.character(phenodata$cellType), 
    path = ".", diff.cutoff = 0.5, pval.cutoff = 0.01) 

The scdata is a gene x cell name matrix .
id is the cell type that each of the cell belongs to .

Please help me fix this error !

Thanks ,
RK

commented

Hi
After investigating the script , I realized that there are no differentially expressed genes identified out of 18 cell types. All the "_MIST.Rdata" files show 0 observations. My data are Log2CPM values.
I have a few questions :

  1. what is the best data transformation to use for dwls ? ( raw counts , CPM, etc)
  2. I did not get any genes even after increasing the p-val cutoffs. So how can I get this data ?

Thanks,
RK

2. I did not get any genes even after increasing the p-val cutoffs. So how can I get this data ?

Following. I am facing similar issues too and I have the same doubts.