mirror of
https://github.com/agdamsbo/daDoctoR.git
synced 2024-11-22 03:40:23 +01:00
27 lines
718 B
R
27 lines
718 B
R
#' Summarising function of HWE calculation
|
|
#'
|
|
#' For easy printing.
|
|
#' @param a1 Allele 1.
|
|
#' @param a2 Allele 2.
|
|
#' @param f factor for grouping.
|
|
#' @keywords hardy-weinberg-equllibrium
|
|
#' @export
|
|
#' @examples
|
|
#' hwe_sum()
|
|
|
|
hwe_sum<-function (a1, a2, f) {
|
|
require(daDoctoR)
|
|
lst <- list()
|
|
df <- data.frame(a1, a2)
|
|
|
|
for (i in 1:length(ls <- split(df, f))) {
|
|
grp <- names(ls)[i]
|
|
obs <- data.frame(hwe_allele(ls[[i]][, 1], ls[[i]][,2])[[c("observed.dist")]])
|
|
pval <- round(hwe_allele(ls[[i]][, 1], ls[[i]][, 2])[[c("p.value")]],3)
|
|
prnt <- paste0(obs[1, ], " (", obs[2, ], ")")
|
|
names(prnt) <- names(obs)
|
|
lst <- list(lst, grp, obs.dist = obs, print = prnt,hwe.pv = pval)
|
|
}
|
|
return(lst)
|
|
}
|