在R中创建我们自己的正态分布函数

在R中创建我们自己的正态分布函数,r,statistics,normal-distribution,R,Statistics,Normal Distribution,我有一个家庭作业问题,要求我创建我自己的标准正态分布函数,我已经推导出来。尽管如此,我在这样做时遇到了问题 我所做的是得到两个独立的均匀分布随机变量,并对其执行box-muller变换 尽管如此,当我绘制图表时,我没有得到我想要的图。这是我的阴谋 但我相信真正的情节应该是这样的: 我做错了吗?看来你的计划是对的。我玩了你提供的链接,它似乎工作。这可能是代码中的一些小错误 vec1 <- runif(10000) vec2 <- runif(10000) z1 <- c(NA)

我有一个家庭作业问题,要求我创建我自己的标准正态分布函数,我已经推导出来。尽管如此,我在这样做时遇到了问题

我所做的是得到两个独立的均匀分布随机变量,并对其执行box-muller变换

尽管如此,当我绘制图表时,我没有得到我想要的图。这是我的阴谋

但我相信真正的情节应该是这样的:


我做错了吗?

看来你的计划是对的。我玩了你提供的链接,它似乎工作。这可能是代码中的一些小错误

vec1 <- runif(10000)
vec2 <- runif(10000)

z1 <- c(NA)
z2 <- c(NA)

z1 <- sqrt(-2*log(vec1))*cos(2*pi*vec2)
z2 <- sqrt(-2*log(vec2))*cos(2*pi*vec1)

vec1您没有提供足够的信息来重现您的问题。但我想这最终会成为一个更合适的问题。如果您的统计数据和代码之间的链接关闭,您将在简历上获得比SO更好和/或更快的答案。请提供代码和数据。不太清楚你做了什么你能给我们看看你的代码吗?你的
z1是什么啊!谢谢你注意到这一点。他们最初有一个用途,但改变了整个设置,使他们成为多余的。