anm.loglik {asbio} | R Documentation |

## Animated plots of log-likelihood functions

### Description

Plots the normal, exponential, Poisson, binomial, and "custom" log-likelihood functions. By definition, likelihoods for parameter estimates are calculated by holding data constant and varying estimates. For the normal distribution a fixed value for the parameter which is not being estimated (`\mu`

or `\sigma^2`

) is established using MLEs.

### Usage

```
anm.loglik(X, dist = c("norm", "poi", "bin", "exp", "custom"),
plot.likfunc = TRUE, parameter = NULL, func = NULL, poss = NULL,
plot.density = TRUE, plot.calc = FALSE, xlab = NULL, ylab = NULL,
conv = diff(range(X))/70, anim = TRUE, est.col = 2, density.leg = TRUE,
cex.leg = 0.9, interval = 0.01, ...)
loglik.norm.plot(X, parameter = c("mu", "sigma.sq"), poss = NULL,
plot.likfunc = TRUE, plot.density = TRUE, plot.calc = FALSE,
xlab = NULL, ylab = NULL, conv = 0.01, anim = TRUE, est.col = 2,
density.leg = TRUE, cex.leg = 0.9, interval = 0.01, ...)
loglik.pois.plot(X, poss = NULL, plot.likfunc = TRUE,
plot.density = TRUE, plot.calc = FALSE, xlab = NULL, ylab = NULL,
conv = 0.01, anim = TRUE, interval = 0.01, ...)
loglik.binom.plot(X, poss = NULL, xlab = NULL, ylab = NULL,
plot.likfunc = TRUE, plot.density = TRUE, conv = 0.01, anim = TRUE,
interval = 0.01, ...)
loglik.exp.plot(X, poss = NULL, plot.likfunc = TRUE,
plot.density = TRUE, plot.calc = FALSE, xlab = NULL, ylab = NULL,
conv = 0.01, anim = TRUE, est.col = 2, density.leg = TRUE,
cex.leg = 0.9, interval = 0.01, ...)
loglik.custom.plot(X, func, poss, anim = TRUE, interval = 0.01,
xlab, ylab, ...)
anm.loglik.tck()
```

### Arguments

`X` |
A vector of quantitative data. The function does not currently handle extremely large datasets, |

`parameter` |
The parameter for which ML estimation is desired in |

`poss` |
An optional vector containing a sequence of possible parameter estimates. Elements in the vector must be distinct. If |

`dist` |
The type of assumed distribution there are currently five possibilities: |

`plot.likfunc` |
A logical command for indicating whether a graph of the log-likelihood function should be created. |

`plot.density` |
A logical command for indicating whether a second graph, in which densities are plotted on the pdf, should be created. |

`plot.calc` |
A logical command for indicating whether a third graph, in which log-densities are added to one another, should be created. |

`xlab` |
Optional |

`ylab` |
Optional |

`conv` |
Precision of likelihood function. Decreasing |

`anim` |
A logical command indicating whether animation should be used in plots. |

`est.col` |
Color used in depicting estimation. |

`density.leg` |
Logical. Should the legend for density be shown? |

`cex.leg` |
Character expansion for legend. |

`interval` |
Speed of animation, in seconds per frame. May not work in all systems; see |

`func` |
Custom likelihood function to be specified when using |

`...` |
Additional arguments from |

### Details

These plots are helpful in explaining the workings of ML estimation for parameters. Animation is included as an option to further clarify processes.
When specifying `poss`

be sure to include the estimate that you "want" the log-likelihood function to maximize in the vector of possibilities, e.g. `mean(X)`

for estimation of `\mu`

.

### Value

Three animated plots can be created simultaneously. The first plot shows the normal, Poisson, exponential, binomial, or custom log-likelihood functions. The second plot shows the pdf with ML estimates for parameters.
On this graph densities of observations are plotted as pdf parameters are varied. By default these two graphs will be created simultaneously on a single graphics device.
By specifying `plot.calc = TRUE`

a third plot can also be created which shows that log-likelihood is the sum of the log-densities.
Animation in this third plot will be automatically sped up, using a primitive routine, for large datasets, and slowed for small datasets.
The third plot will not be created for the binomial pdf because there will only be a single outcome from the perspective of likelihood (e.g. 10 successes out of 22 trials).
The second and third plots will not be created for custom likelihood functions.
Loading package tcltk allows use of the function `anm.loglik.tck`

which provides an interactive GUI to run `anm.loglik`

.

### Author(s)

Ken Aho

### See Also

### Examples

```
## Not run:
##Normal log likelihood estimation of mu.
X<-c(11.2,10.8,9.0,12.4,12.1,10.3,10.4,10.6,9.3,11.8)
anm.loglik(X,dist="norm",parameter="mu")
##Add a plot describing log-likelihood calculation.
anm.loglik(X,dist="norm",parameter="mu",plot.calc=TRUE)
##Normal log likelihood estimation of sigma squared.
X<-c(11.2,10.8,9.0,12.4,12.1,10.3,10.4,10.6,9.3,11.8)
anm.loglik(X,dist="norm",parameter="sigma.sq")
##Exponential log likelihood estimation of theta
X<-c(0.82,0.32,0.14,0.41,0.09,0.32,0.74,4.17,0.36,1.80,0.74,0.07,0.45,2.33,0.21,
0.79,0.29,0.75,3.45)
anm.loglik(X,dist="exp")
##Poisson log likelihood estimation of lambda.
X<-c(1,3,4,0,2,3,4,3,5)
anm.loglik(X,dist="poi")
##Binomial log likelihood estimation of p.
X<-c(1,1,0,0,0,1,0,0,0,0)#where 1 = a success
anm.loglik(X,dist="bin",interval=.2)
##Custom log-likelihood function
func<-function(X=NULL,theta)theta^5*(1-theta)^10
anm.loglik(X=NULL,func=func,dist="custom",poss=seq(0,1,0.01),
xlab="Possibilities",ylab="Log-likelihood")
##Interactive GUI, requires package 'tcltk'
anm.loglik.tck()
## End(Not run)
```

*asbio*version 1.9-7 Index]