Contour plot for two parameters likelihood
Create a contour plot (superimposed with a heat map)
lik.contour(x, y, z, levels = NULL, nlevels = 11, heat = TRUE, col.heat = NULL, ...)
x, y, z |
As in |
levels |
As in |
nlevels |
As in |
heat |
If |
col.heat |
Vector of heat colors |
... |
Additional arguments to |
This function is a wrapper for contour, with a different method to compute
a default value for levels. If heat = TRUE, a heatmap produced by image is added to the plot.
See contour for details on parameters.
Hervé Perdry and Claire Dandine-Roulland
data(AGT) x <- as.bed.matrix(AGT.gen, AGT.fam, AGT.bim) # Compute Genetic Relationship Matrix K <- GRM(x) # eigen decomposition of K eiK <- eigen(K) # simulate a phenotype set.seed(1) y <- 1 + lmm.simu(tau = 1, sigma2 = 2, eigenK = eiK)$y # Likelihood TAU <- seq(0.5,2.5,length=30) S2 <- seq(1,3,length=30) lik1 <- lmm.diago.likelihood(tau = TAU, s2 = S2, Y = y, eigenK = eiK) lik.contour(TAU, S2, lik1, heat = TRUE, xlab = "tau", ylab = "sigma^2")
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.