Nothing
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)
}
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.