mirror of
https://github.com/agdamsbo/daDoctoR.git
synced 2024-11-22 20:00:22 +01:00
35 lines
859 B
R
35 lines
859 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){
|
||
|
|
## HWE summarising function, for several groups defined by factor f. Alleles are provided as vectors, a1 and a2.
|
||
|
|
|
||
|
|
source("https://raw.githubusercontent.com/agdamsbo/research/master/hwe_allele.R")
|
||
|
|
lst<-list()
|
||
|
|
df<-data.frame(cbind(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)
|
||
|
|
}
|