代码之家  ›  专栏  ›  技术社区  ›  Reanimation

计算β函数参数的标准误差

  •  0
  • Reanimation  · 技术社区  · 10 年前

    我在学习 R 使用名为“使用R发现统计数据”的书。这很好,但似乎跳过了一些领域。

    所以,我在 R 计算参数的 a, b 以下pdf格式 beta function :

    func

    我的函数返回从500个样本中找到的以下参数:

      [1] 1.028316 2.095143 #a b
    

    我试图计算参数的标准误差。

    我想知道如何在 R ?

    据我所知, standard errors 根据样本而不是参数计算。所以,我所实现的就是:

        stderr <- function(x) sqrt(var(x)/length(x))
    

    提前感谢。

    1 回复  |  直到 10 年前
        1
  •  2
  •   josliber Martin Ballet    10 年前

    我不确定您使用的是哪一个函数来将测试版分布适合您的数据,但 fitdistr 函数 MASS 包为 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)
    

    这里,标准误差为0.060 形状1 和0.128 形状2 。您可以使用 fit$sd .