Combining distance matricies between multiple aligned genes into final distance matrix
0
0
Entering edit mode
2.9 years ago

I have alignments of Busco genes present across all my assemblies. I want to get a final distance matrix based on these individual gene alignments so I can plot an MDS.

I plan to generate either a similarity matrix with Rs Seqinr package or a distance matrix with RaxML for each gene. Firstly are similarity and distance here just the same with reversed given that the distances are euclidean I believe?

Secondly, how can I combines all these matrices into a single one?

phylogeny distance alignment • 796 views
ADD COMMENT
0
Entering edit mode

Wouldn't it make more sense to concatenate all alignments first and then proceed with your phylogenetic analysis based on a single MSA? Gene-wise calculation of distances gets rid of length information of different gene-based alignments unless you somehow bring them back into the calculus.

ADD REPLY
0
Entering edit mode

I am not sure I understand what you mean. I am taking all (for example) BUSCO 929at5338 genes across all 46 samples and performing an MSA. I am unsure how I could concatenate those alignments with all others into one and and use it for distance analysis.

ADD REPLY
0
Entering edit mode

Could I not just multiple the distance value by gene length to weight every distance by thee length then sum all the dataframes together?

ADD REPLY

Login before adding your answer.

Traffic: 1487 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6