R code to find the probability of presence of a gene using its presence absence data from the current day species