In gkeele/Diploffect.INLA: INLA Implementation of the Diploffect Model

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

library(Diploffect.INLA)

data(exampleCC)
data(locusmatrix)

inla.diploffect <- run.diploffect.inla(formula = y~1+(1|strain)+(1|dose.date), 
                                       add.on = FALSE,                     
                                       data = exampleCC, 
                                       K = NULL, 
                                       prob.matrix = locusmatrix,
                                       num.draws = 10, 
                                       use.dip.lincomb = TRUE, 
                                       seed = 1, 
                                       gamma.rate = 1, 
                                       impute.on = "CCline"
                                       )

inla.diploffect.summary <- run.diploffect.inla.summary.stats(inla.diploffect)

plot.straineff.ci(inla.diploffect.summary, flip = FALSE)

plot.deviation.ci(inla.diploffect.summary, flip = FALSE)

plot.diplotype.ci(inla.diploffect.summary, flip = FALSE)

plot.varexp.ci(inla.diploffect.summary, add.numbers = TRUE)

gkeele/Diploffect.INLA documentation built on May 17, 2023, 8:37 a.m.