inst/md/example_sets.md

In BGLR: Bayesian Generalized Linear Regression

Fitting models with window-specific variances of effects

 library(devtools)
 install_git('https://github.com/gdlc/BGLR-R/')
 library(BGLR)

 data(mice)
 X=scale(mice.X,center=T,scale=F)
 QTLs=seq(from=100,to=10000,length=10)
 signal=rowSums(X[,QTLs])
 signal=signal/sd(signal)
 y=signal+rnorm(nrow(X))

 mrkGroups=rep(1:2000,each=10)[1:ncol(X)]
 fm=BGLR(y=y,ETA=list(list(X=X,model='BRR_sets',sets=mrkGroups))) # Note: method and sets

 plot(fm$ETA[[1]]$varB,cex=.1,col=4,type='o')
 abline(v=QTLs,col=2,lty=2,lwd=.5)

BGLR documentation built on May 12, 2022, 1:06 a.m.