Файл: Correlation examples2.svg

Резюме

#Title: Пример корреляции x и y для различных распределений пар (x, y)# Теги: математика; Статистика; Корреляция# Автор: Дени Бойджело# Необходимые пакеты: mvtnorm (rmvnorm), RSVGTipsDevice (devSVGTips)# Как использовать: output ()## Это переведенная на R версия кода Matematica 6 от Imagecreator.библиотека (mvtnorm)библиотека (RSVGTipsDevice)MyPlot <- function (xy, xlim = c (-4, 4), ylim = c (-4, 4), eps = 1e-15) { title = round (cor (xy [, 1], xy [, 2]), 1) if (sd (xy [, 2]) <eps) title = "" # корр. коэфф. не определено plot (xy, main = title, xlab = "", ylab = "", col = "темно-синий", pch = 16, cex = 0,2, xaxt = "n", yaxt = "n", bty = "n", xlim = xlim, ylim = ylim)}MvNormal <- функция (n = 1000, cor = 0,8) { for (i in cor) { sd = матрица (c (1, i, i, 1), ncol = 2) х = rmvnorm (n, c (0, 0), sd) MyPlot (x) }}вращение <- function (t, X) return (X% *% matrix (c (cos (t), sin (t), -sin (t), cos (t)), ncol = 2))RotNormal <- функция (n = 1000, t = pi / 2) { sd = матрица (c (1, 1, 1, 1), ncol = 2) х = rmvnorm (n, c (0, 0), sd) для (я в т) MyPlot (вращение (i, x))}Другие <- функция (n = 1000) { х = runif (n, -1, 1) y = 4 * (x ^ 2 - 1/2) ^ 2 + runif (n, -1, 1) / 3 MyPlot (cbind (x, y), xlim = c (-1, 1), ylim = c (-1/3, 1 + 1/3)) y = runif (n, -1, 1) xy = вращение (-pi / 8, cbind (x, y)) lim = sqrt (2 + sqrt (2)) / sqrt (2) MyPlot (xy, xlim = c (-lim, lim), ylim = c (-lim, lim)) xy = вращение (-pi / 8, xy) MyPlot (xy, xlim = c (-sqrt (2), sqrt (2)), ylim = c (-sqrt (2), sqrt (2)))  y = 2 * x ^ 2 + runif (n, -1, 1) MyPlot (cbind (x, y), xlim = c (-1, 1), ylim = c (-1, 3)) y = (x ^ 2 + runif (n, 0, 1/2)) * образец (seq (-1, 1, 2), n, replace = TRUE) MyPlot (cbind (x, y), xlim = c (-1,5, 1,5), ylim = c (-1,5, 1,5)) y = cos (x * pi) + rnorm (n, 0, 1/8) х = грех (х * пи) + rnorm (п, 0, 1/8) MyPlot (cbind (x, y), xlim = c (-1,5, 1,5), ylim = c (-1,5, 1,5)) xy1 = rmvnorm (п / 4, с (3, 3)) xy2 = rmvnorm (n / 4, c (-3, 3)) xy3 = rmvnorm (n / 4, c (-3, -3)) xy4 = rmvnorm (n / 4; c (3; -3)) MyPlot (rbind (xy1, xy2, xy3, xy4), xlim = c (-3-4, 3 + 4), ylim = c (-3-4, 3 + 4))}output <- function () { devSVGTips (width = 7, height = 3.2) # удалить первую и последнюю строку, чтобы не экспортировать svg par (mfrow = c (3, 7), oma = c (0,0,0,0), mar = c (2,2,2,0)) MvNormal (800, c (1,0, 0,8, 0,4, 0,0, -0,4, -0,8, -1,0)); RotNormal (200, c (0, пи / 12, пи / 6, пи / 4, пи / 2-пи / 6, пи / 2-пи / 12, пи / 2)); Другое (800) dev.off () # удалить первую и последнюю строку, чтобы не экспортировать svg}

Лицензирование

Я, владелец авторских прав на это произведение, публикую его под следующей лицензией:

Этот файл доступен по лицензии Creative Commons CC0 1.0 Universal Public Domain Dedication .

Лицо, связавшее произведение с этим документом, посвятило произведение общественному достоянию , отказавшись от всех своих прав на произведение во всем мире в соответствии с законом об авторском праве, включая все смежные и смежные права, в той степени, в которой это разрешено законом. Вы можете копировать, изменять, распространять и выполнять работу даже в коммерческих целях, не спрашивая разрешения. http://creativecommons.org/publicdomain/zero/1.0/deed.enCC0Creative Commons Zero, посвящение в общественное достояниеложныйложный