全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
7314 7
2013-02-08
我想用R做一个积分,可是程序老是说出现意外的符号,怎么都找不到问题,请高人帮我看看吧。
电脑老是说有问题


积分.txt
大小:(216 Bytes)

 马上下载



二维码

扫码加我 拉你入群

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

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

全部回复
2013-2-8 17:52:32
你这个语句写的错误很多,是自己写的吗?
二维码

扫码加我 拉你入群

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

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

2013-2-9 08:50:02
你的代码错误很多喔。
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


二维码

扫码加我 拉你入群

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

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

2013-2-10 02:28:45
parazhu 发表于 2013-2-9 08:50
你的代码错误很多喔。
1.0.0358X当中应该有乘号
2.正态密度函数应该为dnorm()
我是写的对数正太,就是的dlnorm。我试试您写的啊,刚学R,不是特别熟。谢了啊
二维码

扫码加我 拉你入群

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

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

2013-2-12 09:45:57
parazhu 发表于 2013-2-9 08:50
你的代码错误很多喔。
1.0.0358X当中应该有乘号
2.正态密度函数应该为dnorm()
您好,我还想再问个问题。R能不能做带参数的积分啊?就是说积分出来的结果是另一个参数的函数。就是比如说一个函数是关于x和y两个变量的,能不能只对x做定积分啊。
二维码

扫码加我 拉你入群

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

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

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

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群