вычислить стандартную ошибку параметра бета-функции

Question

я изучаю р используя книгу под названием «Открытие статистики с помощью R».Это здорово, но, кажется, некоторые места пропускаются.

Итак, у меня есть функция в р который вычисляет параметры a, b из следующего PDF-файла бета-функция:

Моя функция возвращает следующие параметры, найденные в выборке из 500:

  [1] 1.028316 2.095143 #a b

Я пытаюсь вычислить стандартную ошибку параметров.

Мне интересно, как это можно реализовать в р?

Насколько я могу найти в Интернете, стандартные ошибки рассчитываются на основе выборки, а не параметров.Итак, все, что я реализовал, это здесь:

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

Заранее спасибо.

Answer 1

Я не уверен, какую функцию вы используете для подбора бета-дистрибутива к вашим данным, но 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 для shape1 и 0,128 для shape2.Вы можете получить значения с помощью fit$sd.