inst/md/mixedModel.md

In gdlc/BGLR-R: Bayesian Generalized Linear Regression

Using BGLR to fitmodels with multiple sets of effects

  library(BGLR)
  data(mice)

  X=scale(mice.X) # markers
  A=mice.A        # pedigree-relationships
  pheno=mice.pheno

  fm=BGLR(y=pheno$Obesity.BMI,
          ETA=list(
            fixed=list(~factor(GENDER)+factor(Litter),data=pheno,model='FIXED'),
                        cage=list(~factor(cage),data=pheno,model='BRR'),
                        ped=list(K=A,model='RKHS'),
                        mrk=list(X=X,model='BayesB')
        )
    )

    fm$ETA$fixed$b
    fm$ETA$cage$b
    fm$ETA$ped$u
    fm$ETA$mrk$b

Back to examples

gdlc/BGLR-R documentation built on April 23, 2024, 11:01 p.m.