怎么编辑完 还是改不了阿？？    郁闷

#Brownian motion
T = 1                       # end of interval[0,T]
N = 300
x=0                         #initial value of the process at time t0.
t0=0                        #initial time
dt <- (T-t0)/N
t <- seq(t0,T, length=N+1)  # set up scale
X <- ts(cumsum(c(x,rnorm(N)*sqrt(dt))),start=t0, deltat=dt)
ts.plot(X)

进一步了解,请参考
package "sde"
Simulation and Inference for Stochastic Differential Equations

1# llh_xmu

你的第二种方法差一点把我的机器给弄崩溃了，在e1071包中有模拟布朗桥和维纳过程的函数：

# simulate a Wiener process on [0,1] and plot it
x <- rwiener()
plot(x,type="l")


# simulate a Brownian bridge on [0,1] and plot it
x <- rbridge()
plot(x,type="l")

画图：

是因为里面有错误 但是编辑完之后，改不了～～

Thanks～～