logLiksurface  R Documentation 
llsurface
plots the likelihood surface for distributions with two or more parameters,
llcurve
plots the likelihood curve for distributions with one or more parameters.
llsurface(data, distr, plot.arg, min.arg, max.arg, lseq = 50, fix.arg = NULL,
loglik = TRUE, back.col = TRUE, nlev = 10, pal.col = terrain.colors(100),
weights = NULL, ...)
llcurve(data, distr, plot.arg, min.arg, max.arg, lseq = 50, fix.arg = NULL,
loglik = TRUE, weights = NULL, ...)
data 
A numeric vector for non censored data or a dataframe of two columns respectively named left and right, describing each observed value as an interval for censored data. In that case the left column contains either NA for left censored observations, the left bound of the interval for interval censored observations, or the observed value for noncensored observations. The right column contains either NA for right censored observations, the right bound of the interval for interval censored observations, or the observed value for noncensored observations. 
distr 
A character string "name" naming a distribution for which the corresponding density function dname and the corresponding distribution function pname must be classically defined. 
plot.arg 
a twoelement vector with the names of the two parameters that will
vary for 
min.arg 
a twoelement vector with lower plotting bounds for

max.arg 
a twoelement vector with upper plotting bounds for

lseq 
length of sequences of parameters. 
fix.arg 
a named list with fixed value of other parameters. 
loglik 
a logical to plot loglikelihood or likelihood function. 
back.col 
logical (for llsurface only). Contours are plotted with a background gradient of colors if TRUE. 
nlev 
number of contour levels to plot (for llsurface only). 
pal.col 
Palette of colors. Colors to be used as back (for llsurface only). 
weights 
an optional vector of weights to be used in the fitting process.
Should be 
... 
Further graphical arguments passed to graphical functions. 
These two function are not intended to be called directly but is internally called in
llplot
.
llsurface
plots the likelihood surface for distributions with two
varying parameters and other parameters fixed.
When back.col
, image
(2Dplot) is used.
When nlev > 0
, contour
(2Dplot) is used to add
nlev
contours.
llcurve
plots the likelihood curve for distributions with one
varying parameter and other parameters fixed.
MarieLaure DelignetteMuller and Christophe Dutang.
DelignetteMuller ML and Dutang C (2015), fitdistrplus: An R Package for Fitting Distributions. Journal of Statistical Software, 64(4), 134.
See llplot
for an automatic (log)likelihood plots (surface ou curve)
of an object of class "fitdist"
or "fitdistcens"
and plot
, contour
,
image
for classic plotting functions.
# (1) loglikelihood or likelihood curve
#
n < 100
set.seed(1234)
x < rexp(n)
llcurve(data = x, distr = "exp", plot.arg = "rate", min.arg = 0, max.arg = 4)
llcurve(data = x, distr = "exp", plot.arg = "rate", min.arg = 0, max.arg = 4,
loglik = FALSE)
llcurve(data = x, distr = "exp", plot.arg = "rate", min.arg = 0, max.arg = 4,
main = "loglikelihood for exponential distribution", col = "red")
abline(v = 1, lty = 2)
# (2) loglikelihood surface
#
x < rnorm(n, 0, 1)
llsurface(data =x, distr="norm", plot.arg=c("mean", "sd"),
min.arg=c(1, 0.5), max.arg=c(1, 3/2), back.col = FALSE,
main="loglikelihood for normal distribution")
llsurface(data =x, distr="norm", plot.arg=c("mean", "sd"),
min.arg=c(1, 0.5), max.arg=c(1, 3/2),
main="loglikelihood for normal distribution",
nlev = 20, pal.col = heat.colors(100),)
points(0, 1, pch="+", col="red")
llsurface(data =x, distr="norm", plot.arg=c("mean", "sd"),
min.arg=c(1, 0.5), max.arg=c(1, 3/2),
main="loglikelihood for normal distribution",
nlev = 0, back.col = TRUE, pal.col = rainbow(100, s = 0.5, end = 0.8))
points(0, 1, pch="+", col="black")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.