I see what you mean that this method isn’t a perfect way to compare segments but I’ve had a very hard time finding other similarity metrics to use. Any idea what might be a better way to quantify or visualize how similar our data sets are?
My lab has some RNA secondary structure data on a number of virus RNA segments from a method called SHAPE-MaP, which gives each nucleotide in a sequence a reactivity value. We would like to make a 2D plot that clusters segments with similar reactivity value profiles together. However, each segment has a different number of nucleotides, and so we can't figure out how to use standard tools like PCA or t-SNE because each dataset has a different number of dimensions. Is there a way to perform dimensionality reduction on a data set where each point has a different number of dimensions?
1 answer
Not sure whether what you want to do will be informative. That aside, UMAP works with sparse data. Simply insert missing values for shorter sequences to make the length identical to the longest sequence.
You can try DNA/RNA language models and extract their embeddings based on your sequences. Those should be of the same length and presumably can be used for dimensionality reduction.
Log in to answer this question.