top_sites: top_sites

View source: R/functions.R

top_sitesR Documentation

top_sites

Description

Using mitochondrial variants frequency matrix as input, the 20/50 highest variable sites were called by identifying the rows with highest standard deviation (sd) across columns(cells).

Usage

top_sites(d)

Arguments

d

A dataframe containing mitochondrial variants frequency matrix where a column represented a single cell and a row represented variants frequency of a specific mitochondrial genotype. Required.

Details

Highly variable sites identification

Value

A dataframe containing 2 modified mitochondrial genotype matrix with 20/50 highest variable sites.

Examples

data("TF1_clones")
data=TF1_clones$data
d=data_prepare(data)
tops=data.frame()
tops2=data.frame()
for(i in 1:12){
    di=d[[i]]
    tops=rbind(tops,top_sites(di)[[1]])
    tops2=rbind(tops2,top_sites(di)[[2]])
}


songjiajia2018/LINEAGE documentation built on Oct. 17, 2022, 6:17 a.m.