R语言布朗运动模拟股市、物种进化树状图、二项分布可视化

本文模拟了在连续和离散时间布朗演化一些简单的方法。

由Kaizong Ye,Weilong Zhang撰写

布朗运动的数学模型(也称为随机游动)也可以用来描述许多现象以及微小颗粒的随机运动, 如股市的波动和在化石中的物理特性的演变。


布朗运动是随机模式,即改变了从一次到下一个是随机从正态分布绘制均值为0.0,方差为σ2×ΔT。换句话说,根据布朗运动的预期方差通过时间与瞬时差σ2线性增加。


股市模拟

首先,模拟股市一个实例为100的离散时间布朗运动,其中,扩散过程的方差为σ2=0.01。


## 首先模拟随机数
x <- rnorm(n = length(t) - 1, sd = sqrt(sig2)) ## 计算累加和 x <- c(0, cumsum(x))
image.png

视频

马尔可夫链蒙特卡罗方法MCMC原理与R语言实现

探索见解

去bilibili观看

探索更多视频

画图

我们得出各t的时间间隔的随机正偏离改变;然后在每个时间间隔,我们计算累积总和。从而可以看出布朗运动的变化的分布是不变的,并且不依赖于时间的状态。

1)  
X <- cbind(rep(0, nsim), t(apply(X, 1, cumsum)))  
plot(t, X[1, ], xlab = "time", ylab = "phenotype", ylim = c(-2, 2), typ
image.png

 

为了看到这结果如何取决于σ2,我们比较除以10 的SIG2的结果:

X <- matrix(rnorm(n = nsim * (length(t) - 1), sd = sqrt(sig2/10)), nsim,

如上所述,根据布朗运动的预期方差只是σ2。然后,我将使用模拟10000个相同的条件下的结果,以“理顺”我们的结果是:

v <- apply(X, 2, var)  
plot(t, v, type = "l", xlab = "time", ylab = "variance among simulation

图片

PYTHON 用几何布朗运动模型和蒙特卡罗MONTE CARLO随机过程模拟股票价格可视化分析耐克NKE股价时间序列数据

阅读文章


image.png

物种进化

然后,我们尝试用布朗运动模拟物种进化树状图。查看数据的变化情况:

t <- 100  # 总时间  
n <- 30  # 总分支  
b <- (log(n) - log(2))/t


随时关注您喜欢的主题


image.png

现在,来模拟树,我们只需要分别模拟在每个分支的所有分支,然后由最终状态“转向”每个子分支它的父节点。 因为在每个时间步布朗进化的结果是独立于其它所有时间步长。

沿着每条边模拟进化

yy <- sapply(yy, function(x, y) y[[x]][length(y[[x]])], y = X)  
text(x = max(H), y = yy, tree$tip.label)
image.png

在现实中,布朗运动的大部分模拟使用连续的而不是离散的时间进行。

这是因为布朗运动意味着不同物种之间的协方差之间的预期差异。

关于布朗进化的一些其他特点:

在某些情况下,在树的不同部分的布朗进化有可能存在不同的速率。

因此可以简单的模仿不同部门的不同的速率布朗运动。


tree <- sim.history(tree, Q, anc = "1")
image.png

els = TRUE,  
    spread.cost = c(1, 0))

1111.png 下面模拟不同的树从而通过散点图证实相同的父节点产生的树拥有相似的协方差。

plot(tree, edge.width = 2, direction = "downwards")
image.png
rplotMatrix(t(X))
image.png

布朗运动不假定在其下个体谱系移动的过程是高斯过程。其结果将服从高斯分布 – 和中心极限定理。

t <- 0:100    
sig2 <- 0.01  
nsim <- 1000

二项分布的布朗运动

我们模拟二项分布的布朗运动 并查看方差是否和之前一样等于1

apply(X[2:nsim, ], 1, function(x, t) lines(t, x), t = t)
image.png
image.png
image.png

布朗运动,一般认为是没有趋势;然而它(在某些情况下)可以模拟一个模型的趋势。这里是一个模拟(使用如上述相同的一般方法)趋势的一个例子。

X <- matrix(rnorm(mean = 0.02, n = nsim * (length(t) - 1), sd = sqrt(sig2/4)),  
    nsim, length(t) - 1)
image.png

对于布朗运动的数学模型的简单形式的形式: S_T= eS_t-1  其中e是从概率分布绘制。因此,后续还有更多的应用值得进一步研究。


可下载资源

关于作者

Kaizong Ye拓端研究室(TRL)的研究员。

本文借鉴了作者最近为《R语言数据分析挖掘必知必会 》课堂做的准备。

​非常感谢您阅读本文,如需帮助请联系我们!

 
QQ在线咨询
售前咨询热线
15121130882
售后咨询热线
0571-63341498

关注有关新文章的微信公众号


永远不要错过任何见解。当新文章发表时,我们会通过微信公众号向您推送。

技术干货

最新洞察

This will close in 0 seconds