저는 r에서 화합물 포아송 과정을 시뮬레이션하려고합니다. 이 과정은 $ \ sum_ {j = 1}^{N_t} Y_j $에 의해 정의됩니다. $ Y_n $은 i.i.d와 독립적입니다. $ N (0,1) $ values와 $ N_t $는 $ 1 $ 매개 변수가있는 포아송 프로세스입니다. 난 행운을 빌어 r에서 이것을 시뮬레이트하려하고있다. Simutale 0부터 T까지 CPP :r에서 화합물 포아송 과정을 시뮬레이트합니다
초급자 : I는 다음과 같이 계산하는 알고리즘이 $ K = 0 $
반복하면서 $ \ sum_ {I = 1}^K T_i < T $으로
세트 $에서 K = K + 1 $
시뮬레이션 $ T_k \ 심 EXP (\ 람다) $ (내 케이스 $의 \ 람다 = 1 $)
시뮬레이션 $ Y_k \ SIM N (0 , 1) $ (이것은 특별한 경우 일 뿐이므로 이것을 배포본으로 변경할 수 있습니다.)
궤적은 다음과 같이 주어진다. 궤적은 다음과 같이 주어진다. 여기서, $ N (t) = sup (k : \ sum_ {i = 1}^k T_i \ leq t) $
프로세스를 계획 할 수 있도록 누군가를 r에서 시뮬레이션하도록 도와 줄 수 있습니까? 나는 시도했으나 끝내지 못했습니다.
'rnorm','rexp '및'while'을 사용할 수 있습니까? 속도는 느리지 만 다른 프로그래밍 언어와는 전혀 다른 것입니다. 너 뭐 해봤 니? – AdamO