我使用最大化似然法来拟合一些分布。 R“ maxLogL”中实现的方法是一种了不起的工具,效果很好。文档:
AIC和BIC带有摘要功能,但我想将这些值保留在变量中。您可以轻松复制以下代码:
library(EstimationTools)
set.seed(10)
z <- rnorm(n = 1000, mean = 0.1, sd = 1)
fit1 <- maxlogL(x = z, dist = 'dnorm', start=c(0, 2), lower = 0, upper = 2)
a <-summary(fit1)
将打印以下内容:
Optimization routine: nlminb
Standard Error calculation: Hessian from optim
AIC BIC
2824.494 2820.494
Estimate Std. Error
mean 0.011375 0.0313
sd 0.991346 0.0222
我的问题是:
谢谢
PD:我特别想使用maxLogL函数。
我对问题1提出的更接近的解决方案是:
results <- capture.output(summary(fit1))
将控制台输出保留为变量。但这是一个烦人的解决方案,因为那时我需要处理文本。应该有一些简单更好的东西。