1

我有一个程序可以模拟一个公平的骰子 100 次。使用这个程序，我需要在直方图上有一个钟形曲线，显示正态/高斯分布。在我的代码中，代表实验结果的变量是 z。钟形曲线的定义包含两个参数：ave 确定曲线在 z 轴上必须具有最大值的位置，sig 确定峰的宽度。

 set.seed(123)

 x <- sample(1:6, size=100, replace = TRUE)

 hist(x,
 main="100 fair rolls",
 xlab = "Dice Result",
 ylab = "Probability",
 xlim=c(0.5,6.5), 
 breaks=-1:100+.5,
 prob=TRUE )

我到目前为止的高斯分布方程：

    bell<-function(z,ave,sig){
exp(-(z-ave)^2/(2*sig))/sqrt(2*pi*sig)
    }

    points<-0:600

    lines(points,bell(points,350,292))
4

0 回答 0