H2 = 0
H3 - H999 = =SQRT(EXP(0,05)*ABS($I2-50))*NORMINV(RAND();0;1)
I2 - I999 = =$E$2 $D$2*EXP($A$2*($G2-$B$2) $C$2)/(1 EXP($A$2*($G2-$B$2) $C$2)) H2
嗨，大家好，
我在excel中創建了上面的公式。我想將這個實作到 R。任何人都可以幫我解決這個問題嗎？
我已經嘗試使用公式和 data.frame 來實作。
問題是 sigma 取決于 t 和 X(t)-1。
有人可以幫我嗎？
謝謝，馬克斯
uj5u.com熱心網友回復：
這應該讓你開始：
A <- 2
B <- 5
C <- 3
D <- 30
E <- 100
delta <- 0.05
t <- seq(0.05, by = 0.05, length.out = 998)
sigma <- c(0, sqrt(exp(0.05)*abs(t[-998] - 50))*rnorm(997))
Xt <- E D*exp(A*(t - B) C)/(1 exp(A*(t - B) C)) sigma
轉載請註明出處，本文鏈接：https://www.uj5u.com/qukuanlian/524524.html
標籤：r公式微分方程
上一篇：動態物件的數字日歷向量