حساب الخطأ القياسي لمعلمة دالة بيتا

Question

أنا أتعلم ر باستخدام كتاب اسمه "اكتشاف الإحصائيات باستخدام R".إنه أمر رائع ولكن يبدو أنه يتخطى بعض المناطق.

لذا، لدي وظيفة في ر الذي يحسب المعلمات a, b من قوات الدفاع الشعبي التالية وظيفة بيتا:

تقوم وظيفتي بإرجاع المعلمات التالية التي تم العثور عليها من عينة مكونة من 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.