作者lofu (lofu)
看板R_Language
标题关於自订函数
时间Sat Apr 25 21:23:52 2015
各位先进大家好小弟有个写R上碰到的问题想请教各位大大?!
--R入门阶段
> 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