Double loop in R with a matrix -


I know the error I make, but I can not find any solution for the monte carlo method. Programming a double loop

  set.seed (-1256, generic .kind = "box-muller") a & lt; - For matrix (nsimul, 85) (1 in kernel: Nsimul) {r = c () [1] = r0_CIR S = c () S [1] = I0 A [, 1] = for r0_CIR (in Jammu 1: NumPassi) {epsilon = rnorm (2,0,1) r [j + 1] = r [ja] + alphaastar * (gamaster-r [ja]) * deltat + ro * sqrt (r [j]) * Epicelone [1] * SQLT (Delatat) if (R [J + 1] <0) r [j + 1] = abs (r [j + 1]) epsilon_S = epsilon [1] + sqrt (1-core = 2) * Epicelon [2] s [j + 1] = s [j] * XP (r [j] -sigma ^ 2/2-div) * delatat + sigma * epicelonic * scert (deltat)) [a , J + 1] = R [J + 1]}}   

I get this error when I try to run the code

  error [[& lt; -` (`* tmp *`, j + 1, value = 0.0102279735166489):   

Comments

Popular posts from this blog

php - PDO bindParam() fatal error -

logging - How can I log both the Request.InputStream and Response.OutputStream traffic in my ASP.NET MVC3 Application for specific Actions? -

java - Why my included JSP file won't get processed correctly? -