如何使用黎曼和计算 J=P(0<=Z<=1),其中 Z~N(0,1)?
所以 J 是函数 (1/sqrt(2*pi))*exp^((-x^2)/2) 从 0 到 1 的积分
这是我在 R 中实现此功能的方法
m<-5000
a<-0
b<-1
w<-(b-a)/m
x<-seq(a+(w/2),b-(w/2),w)
h<-(1/sqrt(2*pi))*exp^((-x^2)/2)
# Error in exp^((-x^2)/2) : non-numeric argument to binary operator
sum(h*w)
#Error: object 'h' not found
我不知道为什么标记这样的错误,我输入 is.numeric(x)
并返回TRUE
那么如果我只组合数字,问题到底出在哪里?
最佳答案
删除 exp
函数后面的 ^
。尝试:
m<-5000
a<-0
b<-1
w<-(b-a)/m
x<-seq(a+(w/2),b-(w/2),w)
h<-(1/sqrt(2*pi))*exp((-x^2)/2)
sum(h*w)
[1] 0.3413447
关于r - 使用黎曼和计算积分,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/56517919/