作者lofu (lofu)
看板R_Language
标题关於自订函数
时间Sat Apr 25 21:23:52 2015
各位先进大家好小弟有个写Ｒ上碰到的问题想请教各位大大？！
--Ｒ入门阶段
> cum = function(v){
for(i in 1:v){
n = which(cumsum(runif(10))>1)[1]
print(n)
}
}
>cum(3)
[1] 3
[1] 2
[1] 3
产生出来的结果是我想要的但是我想把结果as.dataframe
却无法，一直想出不解答想请各位大大帮帮忙～
谢大家
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.119.157.48
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1429968235.A.4B7.html
1F：→ andrew43: v 的内容是什麽？ 04/25 21:42
2F：→ andrew43: 只是一个整数吗？ 04/25 21:43
3F：→ celestialgod: cum = data.frame(V = sapply(1:v, function(x){ 04/25 21:54
4F：→ celestialgod: sum(cumsum(runif(10)) < 1 ) + 1})) 04/25 21:55
5F：→ celestialgod: 因为你用print，不是把值回传 04/25 21:55
6F：→ celestialgod: 这是一个模拟，算inf_k (unif(0,1) k个和>1) 04/25 21:58
7F：→ celestialgod: 根据蒙提卡罗法应该可以得到是exp(1) 04/25 21:58
9F：→ celestialgod: 以前硕班机率论有证过大於1的case 04/25 22:01
10F：→ celestialgod: 应该是在随机过程的部分，N(t)~sum of uniform dist 04/25 22:03
11F：→ celestialgod: 那个证明是满难的 for general case 04/25 22:05
12F：→ celestialgod: 如果是> 1的话 我记得积分没几行 04/25 22:05
13F：→ lofu: 感恩，问题搞定了 ！ 04/26 11:50