Comments (9)
I've checked the class of expTPM. After the line "expTPM<-expDatList[['TPM']]", class(expTPM) is matrix. But after the line "expTPM<-expTPM[sameProbes,]", class(expTPM) is numeric. Then when it gets to "ncol=ncol(expTPM)", the "ncol" is NULL.
from cellnet.
How many fastqs are you trying to process?
from cellnet.
@pcahan1 Sorry, I realized I was trying to test with only one fastq file.
But when I tried to test with three fastq files, I got another error at cn_salmon:
./salmonRes_FKO_mESC/quant.sf
./salmonRes_WT/quant.sf
./salmonRes_FKO_diff/quant.sf
Error in geneIndexList[[i]] : subscript out of bounds
Actually I didn't got any error when I test on my Mac. But I couldn't figure out what's going on when I tried to test it on server. I used the latest version of salmon and build the index myself.
from cellnet.
I wonder where I can download the fasta file you used to build index: Mus_musculus.GRCm38.cdna.ncrna.exo.fa
from cellnet.
Hi, We are working on this and hope to have a fix very soon. Will post here as soon as we can. Thanks.
from cellnet.
The latter issue
./salmonRes_FKO_mESC/quant.sf
./salmonRes_WT/quant.sf
./salmonRes_FKO_diff/quant.sf
Error in geneIndexList[[i]] : subscript out of bounds
is likely due to using a different index that what we have provided.
We have fixed the earlier error in gene_expr_sum that pops up when the pipeline is applied to one query sample.
Please let us know if you have any more issues. Thanks!
from cellnet.
@pcahan1 Thanks!
I then used cdna and ncrna files from Ensembl to rebuild the index with the latest salmon software and it works throughout the cn_salmon function.
I noticed that the index file that you provided (that is "salmon.index.mouse.122116") doesn't work well with the latest version of salmon(0.8.2), which would raise a error so that we need to rebuild the index with latest version of salmon.
However, I got another problem at the "Reconstruction of GRNs" step when using the "cn_make_grn" function:
esc : 5
liver : 3
neuron : 5
VE : 1
Number of samples per CT: 1
Error in expDat[, rownames(stGRN)] : incorrect number of dimensions
Calls: cn_make_grn
In addition: Warning message:
In if (is.na(tfs)) { :
the condition has length > 1 and only the first element will be used
I was trying to use my own data to train a new CellNet, btw. The training metadata table is in ".rda" format in your paper. I just read my csv file and save as ".rda" in R and then use it for cn_make_grn function together with the output of cn_salmon function. Is there anything wrong with this conversion ? I hope that I describe my situation clearly.
from cellnet.
The output file of cn_salmon is a list, while the parameter for cn_make_grn need to be a matrix ? Step 13-16 in your Nature Protocol paper doesn't describe how to build a new CellNet in details. So far, I have just finished collecting fast files, creating a training metadata table in "csv" format and run through the cn_salmon function...
from cellnet.
Hi, I am going to move this last error
Error in expDat[, rownames(stGRN)] : incorrect number of dimensions
to another issue and close this one as it is not related to the original bug. I cannot recommend using so few samples per cell type or tissue. In the original CellNet paper, the minimum was around 60 samples per cell type, and in the rna-seq we have gone as low as 30-40.
from cellnet.
Related Issues (20)
- cn_salmon Error: cannot open compressed file 'geneToTrans_Homo_sapiens.GRCh38.80.exo_Jul_04_2015.R' HOT 5
- Problem running cn_make_processor HOT 4
- Error in "cn_make_grn" HOT 1
- Example on running locally HOT 1
- cn_salmon HOT 1
- Question: NIS scores for single cell data
- Origin of TF list for Homo sapiens HOT 1
- cn_apply() Error in predict.randomForest HOT 4
- Error at step 8 of the PROCEDURE HOT 1
- Error while running cn_barplot_grnSing
- Beg you to help us reconstruct GRNs HOT 1
- interpreting NIS plot HOT 2
- Can't download the reference example HOT 5
- cn_apply: does not match
- STEP 6: Error in geneIndexList[[i]] : subscript out of bounds
- Unused argument in cn_salmon HOT 2
- Human training data HOT 2
- cn_salmon error HOT 1
- cn_salmon Error:arguments implying differing number of rows HOT 8
- Using CellNet downstream of salmon HOT 4
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from cellnet.