Hi, I am new to the package. I noticed that "def...
# airrflow
j
Hi, I am new to the package. I noticed that "define_clones" folder results of clonal_diversity.tsv files from all_reps_clone_report folder (containing all samples) and an individual sample report folder are different ,such as "d" diversity score with q=0 as an example. Are there any transformation done on "d" from individual sample to be included in all_reps_clone_report? Which 'd' score, should a user to choose ? Thanks so much for any suggestion!
g
Hi Jin, could you share an example of this happening? I think it is due to the bootstrapping sample used in each case to compute diversity, as it is taken as the number of sequences in the sample with least sequences, so when considering all repertoires together, this number will be different.
j
I noticed it even in airrflow demo example here
Here is the detail for sample SRR1383451 at q=0, d is 742.125 for that sample in the individual folder. But the same sample SRR1383451 at q=0 from all_reps_clone_report folder is 61.31. Always is about 1 order difference in my own samples.
g
Yes this is due to the boostrapping sample: when computing diversity, all samples used for diversity computation are bootsrapped without replacement to the number of sequences of the sample with the fewest amount of sequences. So if one sample has less sequences than another, when including them together for the diversity calculation the results can differ.
It would be recommended to exclude from the diversity analysis samples with too few sequences, or many less sequences than others
j
👍
Is there a source or publication describing this algorithm since the two 'd 'scores are quite different? Many thanks!!
g
We use the Hill Diversity index for diversity calculation
Screenshot 2024-09-23 at 20.33.48.png
In a bootstrapped fashion as mentioned above