clonal diversity d50 calculation error

[rbonifer]

clonalDiversity(tcr_G1, cloneCall = "aa", group.by = "donor_condition")

This initial code properly computes the default diversity plots. Although I was hoping to also see d50, so when I include the metric argument for d50, I get:

clonalDiversity(tcr_G1, cloneCall = "aa", group.by = "donor_condition", metrics = "d50")
Error in combine_vars():
! At least one layer must contain all faceting variables: variable
✖ Plot is missing variable
✖ Layer 1 is missing variable
✖ Layer 2 is missing variable

Is there an additional argument that I am missing? I was initially confused because this d50 computation just seems derivative of the clonalProportion function, which functions as advertised. Apologies if this was asked previously, but I could not find any issues revolving around d50. Let me know if there is any further information I can provide.

Thanks,
Riley

[ncborcherding]

Hey @rbonifer,

Thanks for reaching out - the d50 calculation is a recent addition from a user request. I am not sure it makes sense generally, especially if the question is diversity of the repertoire.

In terms of the actual error - can you give me the output of sessionInfo() and some more information on your data?

A similar error was seen in #415 and it may be a result of not having the correct number of variables.

Nick