I want to calculate if the sample has whole genome duplication event or not?

Question

[Hyper_Odin]

I have the copy number data, and i am taking into consideration the mean count of the major allele from each chromosome

the method is published here.

and following is the script, I want to confirm is taking the mean of each chromosome correct?

major_cn <- diffa %>%
  filter(chromosome %in% 1:22) %>%
  group_by(chromosome) %>%
  summarize(major_cn = mean(nMajor))

prop_major_cn_two_or_more <- sum(major_cn$major_cn >= 2) / nrow(major_cn)

if (prop_major_cn_two_or_more >= 0.5) {
  cat("whole genome duplication.\n")
} else {
  cat("no evidence\n")
}

Thanks