get_ncp: Compute ncp for BF null distribution

In william-denault/ffw: A package to perform Fast functional association analysis based on wavelets

Description

Compute non central parameter of the Bayes factors distribution.

Usage

get_ncp(Y, Counfounder, WCs, sigma_b)

Arguments

Y	a vector of the variable of interest.
Counfounder	the counfounding matrix with same a number of line equal to the length of Y. The intercept should not be included, if missing will generate a intercept matrix.
WCs	the vector of a specific wavelet coefficient for each inidivual.
sigma_b	value of the prior used in the Wavelet screaming.

Details

Compute non central parameter of the Bayes factors distribution using its closed form which is provided in the paper of Quan Zhou and Yongtao Guan.

References

Quan Zhou and Yongtao Guan, On the Null Distribution of Bayes Factors in linear Regression, Journal of the American Statistical Association, 518, 2017.

See Also

get_lambda1

Examples

## Not run: 
Y <- rnorm(n=1000)
WCs <- rnorm(n=1000)
sigma <-0.2
get_ncp(Y=Y,WCs=WCs,sigma_b=sigma)

## End(Not run)

william-denault/ffw documentation built on Jan. 20, 2021, 8:28 p.m.