r - 使用黎曼和计算积分

标签 r math statistics integration simulation

如何使用黎曼和计算 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/

相关文章:

r - 按日期子集箱线图,按月排列 x 轴

r - 用最少的步数链接 r 中的记录

java - 计算整数的幂

字符串长度与字符串长度列表的按行比较

xml - 如何将XML中具有相同名称的多个节点转换为R中的df/list?

java - Java 中的数学问题;计算得出的不可思议的 NaN

VBA .TopPadding 和点到英寸转换数学问题

c++ - 如何从 C++ 中的双变量正态分布和学生 T 分布生成随机样本?

r - 尝试比较两个分布

performance - 有没有办法加快 R 中的库加载?