similar to: Simulation of Stochastic processes

Dear R Helpers, I have searched for any R package or code for simulating multivariate fractional Brownian motion (mFBM) or multivariate fractional Gaussian noise (mFGN) when a covariance matrix are given. Unfortunately, I could not find such a package or code. Can you suggest any solution for multivariate FBM and FGN simulation? Thank you for your help. Best Regards, Ryan ----- Wonsang You

I would like to plot multiple random walks onto the same graph. My p variable dictates how may random walks there will be. par(mfrow=c(1,1)) p <- 100 N <- 1000 S0 <- 10 mu <- 0.03 sigma <- 0.2 nu <- mu-sigma^2/2 x <- matrix(rep(0,(N+1)*p),nrow=(N+1)) y <- matrix(rep(0,(N+1)*p),nrow=(N+1)) t<- (c(0:N))/N for (j in 1:p) { z <- rnorm(N,0,1) x[1,j] <- 0 y[1,j]

Hi, I have been using maxLik to do some MLE of Geometric Brownian Motion Process and everything has been going fine, but know I have tried to do it with jumps. I have create a vector of jumps and then added this into my log-likelihood equation, know I am getting a message: NA in the initial gradient My codes is hear # n<-length(combinedlr) j<-c(1,2,3,4,5,6,7,8,9,10)

Hello all, could anyone explain the exact meaning of parameter nproc? Why different values of nproc give so different critical values, i.e. meanL2BB$computeCritval(0.05,nproc=3) [1] 0.9984853 meanL2BB$computeCritval(0.05,nproc=1) [1] 0.4594827 The strucchange-package description gives "integer specifying for which number of processes Brownian motions should be simulated" - do I need

Is there a package that will allow me to fit Brownian motion and Ornstein?Uhlenbeck models of evolution for discrete traits? I know that geiger and ouch have commands for fitting these models for continuous traits, but these aren't suitable for discrete trait evolution, correct? -- View this message in context:

Computational Aspects of Simulation and Inference for Stochastic Processes and the YUIMA Project This two-day workshop is aimed at presenting the latest results on simulation and inference for stochastic processes and their current and prospect implementation within the YUIMA R package. Dates: March 27 (Tue) 9:15-12:15, 14:45-17:45 March 28 (Wed) 9:15-12:15, 14:45-17:45 Location:

Hi list, I'm trying to generate fractional brownian surfaces in R. [A detailed description with respect to various techniques to generating these neutral landscapes has been described by Timothy Keitt (spectral representation of neutral landscapes, Landscape Ecology, 2000) and probably various other authours.] Are there any packages that deal with this or related items? I've been

Hi, I'm writing a short course tutorial to Browninan Motion/ Langevin Equation. At the end of the theory section I wanted to add a short GNU R example, so the students can play a little around. I already looked in the MASS book (by Venables and Ripley) but I couldn't find any Brownian Motion/ Langevin Equation package. Are there any good packages or tutorials available which cover R and

I''d like to be able to use a has_many :through association and treat the associations as a list but I''m getting this error when I try an use an acts_as_list method: NoMethodError: undefined method `move_to_bottom'' I''m using edge rails r6786. Here are my domain rules: Activities are things students can do. Units consists of a sequenced list of

Hello, I would appreciate if someone could help me with this query. I would like to plot a line chart of all of the points in a "for" loop. I would also like to plot the final point with a symbol (to show where the random walk ends). Here is the code I am using: Brownian.fn <- function(mu, sigma, T, N){ dt <- T/N t <- c(rep(NA, N)) B1 <- c(rep(NA, N)) B2 <- c(rep(NA,

To R programming experts, I am a undergraduate student, and now doing research personally. I apply diagonal bivariate poisson (R package "bivpois") with stochatics weighted function (refer to dixoncoles97 section 4.5 to 4.7). However I dont know how to fit this stochatical weighted function to the completed bivariate poisson model. I know that some other references for dynamic soccer

Hi, i need help building a program for the evaluation of a cash or nothing option. The option is written on a stock that today has a price of X. Nine months before i will have this situation: If a<X<b the option pays 3 dollars If X<a or X>b the option pays nothing The price of the title is described by a geometric brownian motion, so i need to start a Montecarlo simulation to find

hi all, I have a simple question that does power spectral analysis related to capacity dimension, information dimension, lyapunov exponent, hurst exponent. If yes then please show me the way. I am newbie in the world of chaos. Sayonara With Smile & With Warm Regards :-) G a u r a v Y a d a v Senior Executive Officer, Economic Research & Surveillance Department, Clearing

Dear all, I'm trying to get a Brownian bridge kernel (kernelbb) for each combination of two consecutive animal locations (see commands below) and put them, with a loop, inside a list. It works well at the beginning but after 42 runs, it appears the following warning : >Error in seq.default(yli[1], yli[2], by = diff(xg[1:2])) :  >  invalid (to - from)/by in seq(.) I looked at the

Hi, this is my first time using R. I want to simulate the following process: "in a population of size N, there are i individuals bearing genotype A, the number of those bearing A is j in the next generation, which following a binominal distribution (choose j from 2*N, the p is i/2*N), to plot the probability of the next generations, my script is as follows. It cannot run successfully,

I am trying to learn how to use the graphics from the lattice package ( and am very new to R). I am trying to replicate the example plot referenced below, by using the lattice xyplot & lpolygon to create panels. I get what appears to be the correct shape of the filled region, but cannot get the position to overlay properly. I have attempted with various settings of position. ( i.e.

Hello all, I wanted to share with you this announcement for an upcoming workshop on the analysis of neural data to be held this summer at woods hole, MA. The workshop has been extremely productive in the past years and we anticipate another good working group this year. This workshop is applicable for any working in time series analysis whether it involves point processes or sampled continuous

Hello everybody, I need your help about the package SN and the skew student distribution. Il will be very grateful if I have the solution. I construct a stochastic model with a white noise not gaussian but following a skew student distribution. I fit the noise on monthly data to obtain the four parameters. The question is : how to annualize the parameters to use my model for simulate daily data

Hi all. A friend of mine asked me to assist with an R presentation to the Penn State LUG in a few weeks. The audience will be predominately computer scientists and/or engineers (I suspect) and they are just looking for presentations of good open source software and my friend volunteered R. We are both relatively new to R, so I thought it would be best to ask the list for your opinion....what do

Dear all, I know that my question is somewhat special but I tried several times to solve the problems on my own but I am unfortunately not able to compute the following test statistic using the quantreg package. Well, here we go, I appreciate every little comment or help as I really do not know how to tell R what I want it to do^^ My situation is as follows: I have a data set containing a