发布于 2014-02-26 08:49:48
我不确定您使用哪个函数将beta发行版与您的数据相匹配,但是来自fitdistr
包的MASS
函数为beta发行版的shape1
和shape2
参数提供了标准错误值:
# Obtain data to fit
set.seed(144)
data <- rbeta(500, 1, 2)
# Fit and output result
library(MASS)
fit <- fitdistr(data, "beta", start=list(shape1=0.5, shape2=0.5))
fit
# shape1 shape2
# 1.0596902 2.0406073
# (0.0602071) (0.1284133)
这里,shape1
的标准错误为0.060,shape2
为0.128。您可以使用fit$sd
获取值。
https://stackoverflow.com/questions/22045969
复制相似问题