extract_tree: Extract sub result from the result of the wavelet_screening

In william-denault/WaveletScreaming: A package to perform Wavelet Screening on GWAS data

Description

Function to performed zoomed analysis of the wavelet screening output

Usage

extract_tree(res, lev_res, thresh)

Arguments

res	Output of Wavelet_screening.
lev_res	the maximum level of resolution needed, has to be less or equal to the request level of resolution in the Wavelet_screening.
thresh	Minimal value of the Bayes Factor to defined a sub-region if missing set as 1.

Value

A vector corresponding to the subtree for the zoomed analysis.

Examples

## Not run: 

#using res for the Wavelet_screening exemple


sub_analysis <- function(res, lev_res )
{
 sub <- extract_tree(res,lev_res=lev_res)
 my_pi <- adaptative_EM_Lambda(sub)
 out <-  adaptative_Lambda (my_pi, sub)
 return(out)
}


sub_analysis(res, 6)


## End(Not run)

william-denault/WaveletScreaming documentation built on Jan. 23, 2021, 12:34 p.m.