Description Usage Arguments Details Examples
These are methods for the maxLik related objects. See also the documentation for the corresponding generic functions
1 2 3 4 5 6 7 8 |
object |
a ‘maxLik’ object ( |
k |
numeric, the penalty per parameter to be used; the default ‘k = 2’ is the classical AIC. |
x |
a ‘maxLik’ object |
eigentol |
The standard errors are only calculated if the ratio of the smallest and largest eigenvalue of the Hessian matrix is less than “eigentol”. Otherwise the Hessian is treated as singular. |
... |
other arguments for methods |
calculates Akaike's Information Criterion (and other information criteria).
extracts the estimated parameters (model's coefficients).
extracts standard errors (using the Hessian matrix).
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | ## estimate mean and variance of normal random vector
set.seed(123)
x <- rnorm(50, 1, 2)
## log likelihood function.
## Note: 'param' is a vector
llf <- function( param ) {
mu <- param[ 1 ]
sigma <- param[ 2 ]
return(sum(dnorm(x, mean=mu, sd=sigma, log=TRUE)))
}
## Estimate it. Take standard normal as start values
ml <- maxLik(llf, start = c(mu=0, sigma=1) )
coef(ml)
stdEr(ml)
AIC(ml)
|
Loading required package: miscTools
Please cite the 'maxLik' package as:
Henningsen, Arne and Toomet, Ott (2011). maxLik: A package for maximum likelihood estimation in R. Computational Statistics 26(3), 443-458. DOI 10.1007/s00180-010-0217-1.
If you have questions, suggestions, or comments regarding the 'maxLik' package, please use a forum or 'tracker' at maxLik's R-Forge site:
https://r-forge.r-project.org/projects/maxlik/
mu sigma
1.068807 1.833129
mu sigma
0.2592487 0.1833165
[1] 206.4963
attr(,"df")
[1] 2
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.