BGLR (version 1.1.3)

getVariances: getVariances

Description

Computes the sample-variance (var()) for sets of markers as well as the total variance.

Usage

getVariances(X, B, sets, verbose=TRUE)

Value

A matrix with variances for markers as well as the total.

Arguments

X

(numeric, \(n \times p\)) incidence matrix for \(\boldsymbol \beta\).

B

(numeric), object returned by the function readBinMat().

sets

(numeric).

verbose

(logical), if TRUE it shows progress information in the console.

Author

Gustavo de los Campos.

Examples

if (FALSE) {
#Demos

library(BGLR)
data(wheat)
y=wheat.Y[,1] ; X=scale(wheat.X)
dir.create('test_saveEffects')
setwd('test_saveEffects')
fm=BGLR(y=y,ETA=list(list(X=X,model='BayesB',saveEffects=TRUE)),nIter=12000,thin=2,burnIn=2000)
B=readBinMat('ETA_1_b.bin')
plot(B[,1],type='o',col=4)
VAR=getVariances(B=B,X=X,sets=sample(1:20,size=1279,replace=T))
head(VAR)
plot(VAR[,"total"])

}