在R上创建柯西分布

Question

我正在执行的c）部分，这需要柯西分布。这是我的尝试：

counter <- 0
slopes <- c()

while(counter < 1000) {
  xvals <- runif(n = 100, min = -1, max = 1)
  evals <- rcauchy(100)
  y <- 5 + 3 * xvals + evals
  fit <- lm(y ~ xvals)
  slopes <- c(slopes, fit$coefficients[[2]])
  counter <- sum(counter, 1)
}
mean(slopes)
print(evals)
plot(evals)

但是，我得到了一些非常奇怪的结果，我怀疑它们是否正确：

我也尝试设置evals

Answer 1

如注释中所述，您的代码正在做“正确的事情”，但是您可能需要调整直觉，以了解尾部发生的事情的重要性

至于功能之间的区别：

• [rcauchy从柯西分布中提取随机值
• [dcauchy是柯西分布的“ probability density function”

有许多与此相关的代码位，rnorm和dnorm用于正态分布，但是您可以在https://stats.stackexchange.com/q/157662/17060中找到有关此差异的更多注释