对于计算ml估计,我应用了包mle2 {bbmle}。我的输出为S4模式。我无法选择“ $”和“ @”参数估计。 我要求选择例如土地估价= 0.31148。 感谢您的回应。 输出: “最大似然估计。
Call:
mle2(minuslogl = nsdlo, start = list(landa = 0.01, teta = 1.001,
alpha = 1.001, beta = 1.002), data = list(x1 = dataset1,
x2 = dataset2))
Coefficients:
Estimate Std. Error z value Pr(z)
landa 0.31148 9.83330 0.0317 0.9747
teta 0.95531 28.32446 0.0337 0.9731
alpha 0.93491 12.68511 0.0737 0.9412
beta 0.99756 NA NA NA
-2 log L: -0.00330078``
答案 0 :(得分:0)
要使用coef
函数访问系数,有关以下示例,请参见?bbmle::mle2
:
x <- 0:10
y <- c(26, 17, 13, 12, 20, 5, 9, 8, 5, 4, 8)
d <- data.frame(x,y)
LL <- function(ymax=15, xhalf=6)
-sum(stats::dpois(y, lambda=ymax/(1+x/xhalf), log=TRUE))
fit1F <- mle2(LL, fixed=list(xhalf=6))
获取第一个系数:
coef(fit1F)[1]
或使用@
:
fit1F@fullcoef[1]