Zacząłem się uczyć
Rkę i jestem zachwycony tym językiem (platformą).
Naprawdę fajne rzeczy można zrobić. Np. do symulowania
procesu Wienera
D2_Wiener <- function() {
dev.new(width = 10, height = 4)
par(mfrow = c(1, 3), oma = c(0, 0, 2, 0))
for(i in 1:3) {
W1 <- cumsum(rnorm(100000))
W2 <- cumsum(rnorm(100000))
plot(W1,W2, type= "l", ylab = "", xlab = "")
}
mtext("2-dimensional Wiener-processes with no correlation",
outer = TRUE, cex = 1.5, line = -1)
}
Kiedy będziemy chcieli dodać skorelowane zmienne to wynik symulacji jest poniżej:
Correlated_Wiener <- function(cor) {
dev.new(width = 10, height = 4)
par(mfrow = c(1, 3), oma = c(0, 0, 2, 0))
for(i in 1:3) {
W1 <- cumsum(rnorm(100000))
W2 <- cumsum(rnorm(100000))
W3 <- cor * W1 + sqrt(1 - cor^2) * W2
plot(W1, W3, type= "l", ylab = "", xlab = "")
}
mtext(paste("2-dimensional Wiener-processes (",cor," correlation)",
sep = ""), outer = TRUE, cex = 1.5, line = -1)
}
Poniżej zamieszkam ściągi do Rki: