The issue is with the calculate_custom_traits() function, which uses pivot_wider() to create columns for each glycan. This does not work when there are analytes without a glycan.
Will have to rewrite the function. Could probably be written more efficiently as well.
For now this can be solved by creating separate block files in LaCyTools for non-glycosylated peptides. That way they are not assigned to the same cluster as glycopeptides.
The issue is with the
calculate_custom_traits()
function, which usespivot_wider()
to create columns for each glycan. This does not work when there are analytes without a glycan.Will have to rewrite the function. Could probably be written more efficiently as well.