如何在R中使用带有glm()的t分布?

时间:2017-12-04 21:31:39

标签: r glm

我已经确定地确定了我用标准lm()回归生成的残差可以通过具有6-ish自由度的t分布合理地建模。我希望glm()使用该错误模型,但我没有看到t适合其中一个系列。关于glm()的替代方案与t一起发挥作用的任何建议,或者一个能够替代(或超集)t的家庭的任何建议?

1 个答案:

答案 0 :(得分:2)

library(heavy) data(ereturns) fit <- heavyLm(m.marietta ~ CRSP, data = ereturns, family = Student(df = 6)) summary(fit) # Linear model under heavy-tailed distributions # Data: ereturns; Family: Student(df = 2.83727) # # Residuals: # Min 1Q Median 3Q Max # -0.142237 -0.036156 0.003433 0.041310 0.546533 # # Coefficients: # Estimate Std.Error Z value p-value # (Intercept) -0.0072 0.0082 -0.8876 0.3748 # CRSP 1.2637 0.1902 6.6459 0.0000 # # Degrees of freedom: 60 total; 58 residual # Scale estimate: 0.002520795 # Log-likelihood: 71.81294 on 3 degrees of freedom 可以执行t-student回归模型。以下是文档中的示例:

HTTPRedirectHandler