作者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/m.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