Comments (3)
Hi, could you share what you put to sgrna_stat
? Thank you
from cb2.
`sgrna_stat <- measure_sgrna_stats(J82_count, df_design,group_a = "T0",group_b = "G",delim = "_",ge_id = "Gene",sg_id = "sgRNA")
head(sgrna_stat)
Gene sgRNA n_a n_b phat_a vhat_a phat_b vhat_b cpm_a cpm_b
1 COL15A1 sgCOL15A1_2 1 3 1.831101e-05 0 2.158151e-05 3.582704e-13 18.31101 22.17798
2 LACC1 sgLACC1_1 1 3 1.675262e-05 0 1.862168e-05 3.459494e-12 16.75262 18.57624
3 BOD1L1 sgBOD1L1_3 1 3 1.465855e-05 0 1.632593e-05 2.691174e-11 14.65855 16.32042
4 ITFG1 sgITFG1_1 1 3 2.873270e-06 0 1.092631e-05 5.610771e-12 2.87327 10.88184
5 KCTD16 sgKCTD16_2 1 3 2.459324e-05 0 2.277680e-05 3.781127e-13 24.59324 24.60295
6 ADD3 sgADD3_1 1 3 1.090869e-05 0 1.185389e-05 2.556015e-11 10.90869 11.84388
logFC t_value df p_ts p_pa p_pb fdr_ts fdr_pa fdr_pb
1 0.2633311870 5.4639798 2 0.03190101 0.98404950 0.01595050 0.1001397 0.9997659 0.1503497
2 0.1410711114 1.0048866 2 0.42077345 0.78961327 0.21038673 0.5329416 0.9997659 0.5300289
3 0.1455238811 0.3214133 2 0.77837816 0.61081092 0.38918908 0.8370835 0.9997659 0.8109563
4 1.6171339532 3.3997608 2 0.07669632 0.96165184 0.03834816 0.1780190 0.9997659 0.2326820
5 0.0005473835 -2.9540115 2 0.09803526 0.04901763 0.95098237 0.2077280 0.1735382 0.9999973
6 0.1090668677 0.1869572 2 0.86894153 0.56552923 0.43447077 0.9079486 0.9997659 0.8791015`
from cb2.
have you resolved the problem? otherwise, could you send the r object to [email protected]?
Thank you,
from cb2.
Related Issues (20)
- run_sgrna_quant fails HOT 3
- inconsistent between .fq.gz and .fastq? HOT 7
- calc_mappability HOT 1
- export join_count_and_design HOT 2
- plot_count_distribution add export option HOT 1
- gene_stat cpms HOT 2
- Clarity around logFC for gene_stats HOT 3
- Problem with gene-level statistic HOT 1
- Error in arising in run_sgrna_quant HOT 4
- How to handle paired-end fastqs HOT 5
- Cluster setup failed. 31 of 31 workers failed to connect. HOT 2
- run_sgrna_quant report the wrong sequences associated to sgRNA names HOT 2
- maximum guide length HOT 2
- Error in cb2_count() HOT 4
- Is it possible to analyze by inserting one mismatch guide RNA barcode? HOT 3
- Interaction Terms and Complex Designs HOT 1
- Fasta file issue HOT 16
- CB2 vhat estimation with variable total read count
- Choice of count normalisation
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 cb2.