How to simulate a null distribution of paralog ages without ancient WGDs to use in K-S goodness-of-fit test
Entering edit mode
4 months ago

Hi there,

I would like to determine whether my gene age distribution plot contains potential whole-genome duplications. I would like to do with a K-S goodness-of-fit test by comparing the paralog ages in my plot to a simulated null distribution (without ancient WGDs).

My question is as follows: Is there a way for an R novice like myself to model this null distribution in order to compare it with my gene age distribution plot? I found the following paper describing how to model the age distribution of paralogs under a constant birth–death model but I'm finding it very difficult to understand and implement: Widespread genome duplications throughout the history of flowering plants

Some advice or useful resources would be extremely appreciated.

Thank you for your time

WGD distributions plots age KS gene • 129 views

Login before adding your answer.

Traffic: 1428 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6