R/plot.proLik.R

In nadiv: Functions to construct non-additive genetic relatedness matrices

plot.proLik <- function(x, CL = TRUE, alpha = 0.05, type = "l", ...)
{
  if(CL == TRUE & is.null(alpha)) stop("'alpha' must be non-null to include CL in the plot")

  plot(x$lambdas ~ x$var.estimates, 
     xlab = x$component, ylab = "LRT statistic", 
     type = "l", lwd = 2, ...)
     if(CL){  
        chi <- (-0.5 * qchisq(alpha, df = 1, lower.tail = FALSE))
        abline(h = chi, lty = "dotted", col = "red", lwd = 2)
    }  
}