Comments (4)
Hi Simon,
This isn't something I've really done before, since most projects I've worked on have just had a single input library.
I think the median approach you suggest should work? If you wanted to get really fancy, you could probably put errors/deviations on these medians (median absolute deviations across the replicates is one idea?) and use these in some kind of significance calculation. Whether or not this is overkill probably depends on how certain you need to be in your results downstream - from my experience looking at essential genes from TraDIS experiments in different bacteria, I suspect this might take a lot of work to get to more or less the same result in the end, though I'd be quite happy to be proved wrong.
Of course the easiest thing if you just want a quick result would just be to run the script three times on your three libraries, check that they agree, and maybe exclude things that are borderline or switch categories in different experiments.
Feel free to email me directly if you want any advice on specific data, or just a second set of eyes to sanity-check.
-Lars
from bio-tradis.
Hi Lars,
if you have only worked with a single input library, I am wondering from where you got the statistical power for
tradis_comparison.R
? I would expect edgeR to require at least 2 replicates of a group, right?
I am having 3 independent biological replicates per group.
Well, the median approach was the best one I could think of.
So instead of taking the median I take the MAD, did I understand that right? And find a significance how?
If you could elaborate on that (or email me personally) I would be very happy.
I am not quite sure if running the thing 3 times in a row would be statistically correct. Of course we could perform as you suggest, but we would not able to take any technical or human bias into account? I would expect that each replicate might behave differently.
Lets say I have 0 | 1 | 25 as insertion counts for one replicate and gene. The median would be 1, but the question is, would that be close to the biological truth? If your experience tells you that, then I would be happy to follow that way. However, how would I justify that in a paper? You think that would fly?
I have to ask, if I can forward you any data, then I will come back to you! Thanks for the offer.
Best,
Simon
from bio-tradis.
Hi Simon,
Generally the way we (and I think most labs, at least for in vitro work where population bottlenecks are not an issue) do these sorts of experiments is to generate a single dense base transposon mutant library, then grow this up in independent replicates and compare these. Examples of comparisons might be something like the input culture and bacteria recovered from an organ for an infection model, or growth in rich media over some time period with and without an antibiotic. So the replication is done at the level of the experiment, and not at the level of generating the mutant library.
As you've found, tradis_essentiality.R only deals with the single base library. If you've created multiple independent dense transposon mutant libraries, it's probably possible to create something more robust than tradis_essentiality.R for essentiality calling. This is not something I've done, and I'm not aware of anyone who has. The question is really whether it makes sense to spend the time to do this, and this depends on the scope of the statement you want to make based on this data. If all you want is to be able to say "we found X genes were essential, and Y% are shared with E. coli" (or similar), just running the analysis 3 times and verifying it is consistent is probably sufficient. If you're planning to write an entire paper just about the gene essentiality, then it might be worth spending the time thinking about different ways of approaching this data.
Anyway, I think this is getting off-topic for an issue report -- please close the issue & email me directly if you want to discuss this further.
from bio-tradis.
Thanks!
Sure, I will come back to you.
from bio-tradis.
Related Issues (20)
- installing biotradis on MacOS HOT 4
- test failures with samtools v1.10 HOT 1
- Essentiality analysis code HOT 4
- ftp download link in biotradistutorial.pdf no longer exists
- Possible issue with calculation of read start and ends in module Cigar.pm HOT 2
- Problem with insertion plots (files longer than they should be) HOT 9
- tradis_essentiality.R output a file containing change points?
- tradis_comparison issue HOT 3
- -mm paramater not working as expected HOT 1
- Unable to install biotradis
- Essentiality function not working HOT 1
- Issues on running BioTraDIS on multiple contigs HOT 1
- EdgeR not working with new R version?
- how is the insertion index calculated? HOT 6
- `tradis_merge_plots` cannot find `tradisfind` HOT 3
- bacteria_tradis tsv extension bug
- Possible bug handling reads with clipping when converting mapped reads to insertion plot HOT 1
- Problems with bacteria_tradis mapping HOT 1
- conda install提示需要perl-moose HOT 1
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 bio-tradis.