Мне было интересно, возможно ли генерировать коррелированные случайные биномиальные переменные, следуя подходу линейного преобразования? Ниже я попробовал что-то простое в R, и оно дает некоторую корреляцию. Но мне было интересно, есть ли принципиальный способ сделать это? X1 = rbinom(1e4, 6, .5) ;...
Я пробовал обычный поиск в Google и т. Д., Но большинство ответов, которые я нахожу, либо несколько неоднозначны, либо специфичны для языка / библиотеки, например Python или C ++ stdlib.hи т. Д. Я ищу независимый от языка математический ответ, а не специфику библиотеки. В качестве примера многие...
Как работает метод инверсии? Скажем , у меня случайную выборку X1,X2,...,XnX1,X2,...,XnX_1,X_2,...,X_n с плотностью f(x;θ)=1θx(1−θ)θf(x;θ)=1θx(1−θ)θf(x;\theta)={1\over \theta} x^{(1-\theta)\over \theta} более...
Например, в R, MASS::mvrnorm()функция полезна для генерации данных, чтобы продемонстрировать различные вещи в статистике. Он принимает обязательный Sigmaаргумент, который представляет собой симметричную матрицу, определяющую ковариационную матрицу переменных. Как бы я создал симметричную матрицу с...
Как я могу сделать выборку из распределения смеси, и в частности из смеси нормальных распределений в R? Например, если я хотел сделать выборку из: 0,3× N( 0 , 1 )+0,5× N( 10 , 1 )+0.2× N( 3 , .1 )0.3×N(0,1)+0.5×N(10,1)+0.2×N(3,.1) 0.3\!\times\mathcal{N}(0,1)\; + \;0.5\!\times\mathcal{N}(10,1)\; +...
Я пытаюсь написать программу на R, которая имитирует псевдослучайные числа из распределения с помощью кумулятивной функции распределения: F( х ) = 1 - опыт( - а х - бр + 1Икср + 1) ,х ≥ 0F(Икс)знак равно1-ехр(-aИкс-бп+1Иксп+1),Икс≥0F(x)= 1-\exp \left(-ax-\frac{b}{p+1}x^{p+1}\right), \quad x \geq 0...
Я читаю об адаптивном MCMC (см., Например, главу 4 « Справочника цепи Маркова Монте-Карло» , изд. Brooks et al., 2011; а также Andrieu & Thoms, 2008 ). Основной результат Roberts and Rosenthal (2007) состоит в том, что если схема адаптации удовлетворяет исчезающему условию адаптации (плюс...
Я пытаюсь сгенерировать коррелированную случайную последовательность со средним значением = , дисперсией = , коэффициентом корреляции = . В приведенном ниже коде я использую & как стандартные отклонения, а & как средство.0001110.80.80.8s1s2m1m2 p = 0.8 u = randn(1, n) v = randn(1, n) x = s1...
Я использую разложение Холецкого для моделирования коррелированных случайных величин с учетом матрицы корреляции. Дело в том, что результат никогда не воспроизводит структуру корреляции так, как он задан. Вот небольшой пример на Python, чтобы проиллюстрировать ситуацию. import numpy as np n_obs =...
Каков наилучший способ попробовать дистрибуцию Кантора ? У него есть только cdf, и мы не можем его
Здесь «усечение» подразумевает снижение точности случайных чисел, а не усечение последовательности случайных чисел. Например, если у меня есть действительно случайных чисел (взятых из любого распределения, например, нормальных, равномерных и т. Д.) С произвольной точностью, и я обрезаю все числа,...
Я хотел бы нарисовать целые числа от 1 до некоторого определенного , бросая некоторое количество справедливых шестигранных кубика (d6). Хороший ответ объяснит, почему его метод производит однородные и независимые целые числа.NNN В качестве иллюстративного примера было бы полезно объяснить, как...
Я хотел бы создать образцы из синей области, определенной здесь: Наивным решением является использование выборки отбраковки на единицу площади, но это обеспечивает эффективность только (~ 21,4%).1−π/41-π/41-\pi/4 Есть ли какой-нибудь способ, которым я могу сделать выборку более эффективно?...
Я написал программу, которая генерирует случайные данные. Если программа работает правильно, эти данные должны соответствовать определенному, известному распределению вероятности. Я хотел бы запустить программу, сделать некоторые расчеты по результату и получить значение p. Прежде чем кто-либо еще...
Мне нужно генерировать случайные числа после нормального распределения в пределах интервала . (Я работаю в Р.)(a,b)(a,b)(a,b) Я знаю, что функция rnorm(n,mean,sd)будет генерировать случайные числа после нормального распределения, но как установить пределы интервала в этом? Для этого есть какие-то...
Я хотел бы нарисовать образец x∼N(0,Σ)x∼N(0,Σ)\mathbf{x} \sim N\left(\mathbf{0}, \mathbf{\Sigma} \right) . Википедия предлагает использовать либо разложение Холецкого, либо Собственное , то есть или Σ=D1DT1Σ=D1D1T \mathbf{\Sigma} = \mathbf{D}_1\mathbf{D}_1^T Σ = Q Λ QTΣзнак равноQΛQT...
В Libre Office Calc rand()доступна функция, которая выбирает случайное значение от 0 до 1 из равномерного распределения. Я немного заржавел в своей вероятности, поэтому, когда я увидел следующее поведение, я был озадачен: A = 200x1 столбец rand()^2 B = 200x1 столбец rand()*rand() mean(A) знак равно...
Как и почему генераторы случайных чисел (ГСЧ) важны в вычислительной статистике? Я понимаю, что случайность важна при выборе выборок для многих статистических тестов, чтобы избежать смещения в отношении любой гипотезы, но есть ли другие области вычислительной статистики, где важны генераторы...
Из статистической случайности Википедии : Глобальная случайность и локальная случайность различны. Большинство философских концепций случайности являются глобальными, потому что они основаны на идее, что «в долгосрочной перспективе» последовательность выглядит действительно случайной, даже если...
Предположим, у нас есть Икс1∼ униф ( n , 0 , 1 ) ,Икс1~UNIF(N,0,1),X_1 \sim \textrm{unif}(n,0,1), Икс2~ unif (n,0,1),Икс2~UNIF(N,0,1),X_2 \sim \textrm{unif}(n,0,1), где - равномерная случайная выборка размером n, иunif (n,0,1 )UNIF(N,0,1)\textrm{unif}(n,0,1) Y= Х1,Yзнак равноИкс1,Y=X_1, Z= 0,4 х1+...