如何提高R中参数值不能用向量表示的分布函数的计算速度

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

比如对于pnorm（x,a,b）函数，x,a,b可以同时表示成向量带入计算，但是对于pstable（x,a,b,c,d）函数a,b,c,d就不能表示成向量，
那这个时候如果我要计算10^6组(x,a,b,c,d)所对应的pstable值，如果用for语句把每一组的pstable值循环求一遍速度太慢，请问有没有提高运算速度的其他方法？谢谢！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：分布函数 如何提高 不能用 Stable Table 分布函数 如何

1. pstabley<-function(y) pstable(y[1],y[2],y[3],y[4])

复制代码

这回用pstabley（7）就可以了

nuomin 发表于 2015-3-4 20:30
这回用pstabley（7）就可以了

不是呀，假如y是n*5的矩阵，每一行对应的是（x,a,b,c,d）的一组值，一共有n=10^6组，那我就要对应的求出n个pstable的值，而pstable（x,a,b,c,d）函数规定啊a,b,c,d四个参数length==1，也就是不能为向量，不能直接带入y[,2],y[,3],y[,4]和y[,5]，那我怎样以快速的方法求出这n个pstable的值？

过眼云烟83 发表于 2015-3-4 21:08
不是呀，假如y是n*5的矩阵，每一行对应的是（x,a,b,c,d）的一组值，一共有n=10^6组，那我就要对应的求出n ...

用mapply（）函数，这是R的区别于其他语言的本质特征

nuomin 发表于 2015-3-5 13:48
用mapply（）函数，这是R的区别于其他语言的本质特征

谢谢，虽然速度上并不能有所提高，但是的确简化了编程语言