求助，R怎么做积分

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我想用R做一个积分，可是程序老是说出现意外的符号，怎么都找不到问题，请高人帮我看看吧。



积分.txt (216 Bytes)

2013-2-8 14:58:31 上传

这是我写的程序

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：怎么做 积分

你这个语句写的错误很多，是自己写的吗？

你的代码错误很多喔。
1.0.0358X当中应该有乘号
2.正态密度函数应该为dnorm()
3.为什么要写一个pt的函数？这个造成了myfun()中出现非数值的错误。把这个直接写入就可以了。
4.integrate(myfun2(x,2),0,29.47)，这里就是引用函数名myfun2，你的原函数设置没有第二个参数，怎么会多个2出来呢？
4.integrate()完成后是个结构体，需要引用其中的value，才能进行加减计算
可行结果：

> myfun1=function(x){(x/14.04)*(0.0358*x+2949.38)^(-2)*dnorm(x,8.84,0.91)}
> myfun2=function(x){(x/14.04)*(0.0229*x+2949.76)^(-2)*dnorm(x,8.84,0.91)}
> integrate(myfun1,29.47,Inf)$value+integrate(myfun2,0,29.47)$value
[1] 7.235227e-08

parazhu 发表于 2013-2-9 08:50
你的代码错误很多喔。
1.0.0358X当中应该有乘号
2.正态密度函数应该为dnorm()

我是写的对数正太，就是的dlnorm。我试试您写的啊，刚学R，不是特别熟。谢了啊

parazhu 发表于 2013-2-9 08:50
你的代码错误很多喔。
1.0.0358X当中应该有乘号
2.正态密度函数应该为dnorm()

您好，我还想再问个问题。R能不能做带参数的积分啊？就是说积分出来的结果是另一个参数的函数。就是比如说一个函数是关于x和y两个变量的，能不能只对x做定积分啊。

就我所知，一元积分和使用cubature软件包的多元积分，由于使用的是数值方法，符号积分是无法完成的。
可能有人会开发符号积分的软件包，如果哪位知道可以跟帖分享。但个人认为，如果你想使用符号运算，那么使用matlab和mathmatica，会十分方便，毕竟软件各有所长。我在做几个G的大容量数据处理时候，就用SAS，做符号运算就用mathmatica，调自己写的大程序就用matlab，使用统计和计量分析用SPLUS和R，几个混合使用，依靠各个软件的长处。

parazhu 发表于 2013-2-12 12:17
就我所知，一元积分和使用cubature软件包的多元积分，由于使用的是数值方法，符号积分是无法完成的。
可能 ...

噢噢，那我试试mathmatica好了，谢啦~~非常感谢^^

F=function(x){plnorm(x,mu,sigma)}应该可以