R/utils.R

In variani/finemapr: R wrapper to fine-mappers

#------------------
# data processing
#------------------

prepare_zscore_writing <- function(tab)
{
  filter(tab, finemap) %>% select(snp, zscore)
}

merge_tab_snp <- function(tab, snp)
{
    select(tab, rank_z, snp) %>%
    left_join(snp, by = "snp") %>%
    select(snp, rank_z, rank_pp, everything()) %>%
    arrange(rank_pp)
}

#----------------
# LD
#----------------

simulate_ld_diag <- function(snps)
{
  ld <- diag(length(snps))
  rownames(ld) <- snps
  colnames(ld) <- snps
  
  return(ld)
}