4.11

Выборка из распределения

[18/100%]
Показать
LaTeX
Задача 4.11.1

Если UU равномерно распределена на [0,1][0,1], каково распределение X=⌊nU⌋+1X = \lfloor n U\rfloor +1?

?
Задача 4.11.2

Дана последовательность S0S_{0} первых nn целых чисел в произвольном порядке. Поступайте следующим образом:

(a) выберите случайно любую позицию P0P_{0} из {1,2,…,n}\left\{ 1,2, \ldots , n\right\} и поменяйте местами число в этой позиции S0S_{0} с числом на nn-м месте S0S_{0}, получив S1S_{1}.

(b) выберите случайно любую позицию P1P_{1} из {1,2,…,n−1}\left\{ 1,2, \ldots , n-1\right\} и поменяйте местами число в этой позиции S1S_{1} с числом на (n−1)(n-1)-м месте S1S_{1}, получив S2S_{2},

(c) на (r−1)(r-1)-м шаге число в позиции Pr−1P_{r-1}, выбранной случайно из {1,2,…,n−r+1}\left\{ 1,2, \ldots , n-r+1\right\}, меняется местами с числом на (n−r+1)(n-r+1)-м месте последовательности Sr−1S_{r-1}. Покажите, что Sn−1S_{n-1} равновероятно может быть любой из nn ! перестановок {1,2,…,n}\left\{ 1,2, \ldots , n\right\}.

?
Задача 4.11.3

Используйте метод отбраковки для моделирования выборки из гамма-плотности Γ(λ,t)\Gamma (\lambda , t), где t(≥1)t( \geq 1) не обязательно целое. [Указание: возможно, стоит начать с экспоненциальной случайной величины с параметром 1/t1 / t.]

?
Задача 4.11.4

Покажите, как получить выборку из бета-плотности β(α,β)\beta (\alpha , \beta ), где α,β≥1\alpha , \beta \geq 1. [Указание: используйте Упражнение (4.11.3).]

?
Задача 4.11.5

Опишите три различных метода получения выборки из плотности распределения f(x)=6x(1−x),0≤x≤1f(x) = 6 x(1-x), 0 \leq x \leq 1.

?
Задача 4.11.6

Конечный вещественный вектор называется вероятностным вектором, если все его элементы неотрицательны и в сумме дают 1. Покажите, что вероятностный вектор p\mathbf{p} длины nn может быть записан в виде

p=1n−1∑r=1nvr \mathbf{p} = \frac{1}{n-1} \sum _{r = 1}^{n} \mathbf{v}_{r}

где каждый vr\mathbf{v}_{r} — вероятностный вектор, имеющий не более двух ненулевых элементов. Опишите метод, основанный на этом наблюдении, для генерации выборки из p\mathbf{p}, рассматриваемого как функция вероятности.

?
Задача 4.11.7

Пусть U1U_{1} и U2U_{2} независимы и равномерно распределены на [0,1][0,1], и пусть Ti=2Ui−1T_{i} = 2 U_{i}-1. Покажите, что при условии события R=T12+T22≤1R = \sqrt{T_{1}^{2}+T_{2}^{2}} \leq 1,

X=T1R−2log⁡R2,Y=T2R−2log⁡R2 X = \frac{T_{1}}{R} \sqrt{-2 \log R^{2}}, \quad Y = \frac{T_{2}}{R} \sqrt{-2 \log R^{2}}

являются независимыми стандартными нормальными случайными величинами.

?
Задача 4.11.8

Пусть UU равномерно распределена на [0,1][0,1] и 0<q<10 < q < 1. Покажите, что X=1+⌊log⁡U/log⁡q⌋X = 1+\lfloor \log U / \log q\rfloor имеет геометрическое распределение.

?
Задача 4.11.9

Точка (X,Y)(X, Y) выбирается равномерно случайным образом в полукруге x2+y2≤1,x≥0x^{2}+y^{2} \leq 1, x \geq 0. Каково распределение Z=Y/XZ = Y / X?

?
Задача 4.11.10

Пусть XX — неотрицательная целочисленная случайная величина с h(r)=P(X=r∣X≥r)h(r) = \mathbb {P}\left(X = r \mid X \geq r\right). Если {Ui:i≥0}\left\{ U_{i}: i \geq 0\right\} независимы и равномерно распределены на [0,1][0,1], покажите, что Z=min⁡{n:Un≤h(n)}Z = \min \left\{ n : U_{n} \leq h(n)\right\} имеет то же распределение, что и XX.

?
Задача 4.11.11

Пусть g(x1,x2,…,xn)g\left(x_{1}, x_{2}, \ldots , x_{n}\right) — функция, возрастающая по каждой из своих переменных, и пусть {Ur:r≥1}\left\{ U_{r}: r \geq 1\right\} — независимые одинаково распределённые случайные величины, имеющие равномерное распределение на [0,1][0,1]. Покажите, что

cov⁡{g(U1,U2,…,Un),g(1−U1,1−U2,…,1−Un)}≤0 \operatorname {cov}\left\{ g\left(U_{1}, U_{2}, \ldots , U_{n}\right), g\left(1-U_{1}, 1-U_{2}, \ldots , 1-U_{n}\right)\right\} \leq 0

[Указание: используйте неравенство ФКГ из задачи (3.11.18).] Объясните, как это может помочь в эффективном оценивании I=∫01g(x)dxI = \int_{0}^{1} g(\mathbf{x}) d \mathbf{x}.

?
Задача 4.11.12

Мы хотим оценить I=∫g(x)fX(x)dx=E[g(X])I = \int g(x) f_{X}(x) d x = \mathbb {E}\left[g(X\right]), где либо трудно производить выборку из плотности fXf_{X}, либо g(X)g(X) имеет очень большую дисперсию. Пусть fYf_{Y} эквивалентна fXf_{X}, то есть для всех x,fX(x)=0x, f_{X}(x) = 0 тогда и только тогда, когда fY(x)=0f_{Y}(x) = 0. Пусть {Yi:0≤i≤n}\left\{ Y_{i}: 0 \leq i \leq n\right\} — независимые случайные величины с плотностью fYf_{Y}, и определим

J=1n∑r=1ng(Yr)fX(Yr)fY(Yr) J = \frac{1}{n} \sum _{r = 1}^{n} \frac{g\left(Y_{r}\right) f_{X}\left(Y_{r}\right)}{f_{Y}\left(Y_{r}\right)}

Метод, изобретённый Дж. М. Хаммерсли и К. У. Мортоном в 1956 году.

Покажите, что:

?
(a)

E[J]=I=E[g(Y)fX(Y)fY(Y)]\mathbb {E}\left[J\right] = I = \mathbb {E}\left[\frac{g(Y) f_{X}(Y)}{f_{Y}(Y)}\right],

(b)

Var⁡(J)=1n[E[g(Y)2fX(Y)2fY(Y)2]−I2]\operatorname {Var}\left(J\right) = \frac{1}{n}\left[\mathbb {E}\left[\frac{g(Y)^{2} f_{X}(Y)^{2}}{f_{Y}(Y)^{2}}\right]-I^{2}\right],

(c)

J→ п.н. IJ \xrightarrow {\text{ п.н. }} I при n→∞n \rightarrow \infty. (См. главу 7 для описания сходимости.)

Идея состоит в том, что из fYf_{Y} должно быть легко производить выборку, и она должна быть выбрана, по возможности, так, чтобы Var⁡(J)\operatorname {Var}\left(J\right) было намного меньше, чем n−1[E[g(X)2]−I2]n^{-1}\left[\mathbb {E}\left[g(X)^{2}\right]-I^{2}\right]. Функция fYf_{Y} называется плотностью важности.

Задача 4.11.13

Постройте два различных метода генерации выборки из арксинусной плотности распределения

f(x)=2π1−x2,0≤x≤1 f(x) = \frac{2}{\pi \sqrt{1-x^{2}}}, \quad 0 \leq x \leq 1
?
Задача 4.11.14

Пусть (X,Y)(X, Y) — случайная точка, выбранная равномерно на единичном диске. Покажите, что

(U,V,W)=(2X1−X2−Y2,2Y1−X2−Y2,1−2(X2+Y2)) (U, V, W) = \left(2 X \sqrt{1-X^{2}-Y^{2}}, 2 Y \sqrt{1-X^{2}-Y^{2}}, 1-2\left(X^{2}+Y^{2}\right)\right)

равномерно распределена на единичной сфере.

?
Задача 4.11.15

Требуется генерировать выборку из плотности fXf_{X}, удовлетворяющей fX=f1+f2f_{X} = f_{1}+f_{2}, где f1,f2:R→[0,∞)f_{1}, f_{2}: \mathbb {R} \rightarrow [0, \infty ), и f1≤fYf_{1} \leq f_{Y} для некоторой случайной величины YY. Пусть UU равномерно распределена на [0,1][0,1] и независима от YY. Из YY и UU мы генерируем другую случайную величину ZZ следующим образом. Если U>f1(Y)/fY(Y)U > f_{1}(Y) / f_{Y}(Y), пусть ZZ имеет плотность f2(z)/∫−∞∞f2(u)duf_{2}(z) / \int_{-\infty }^{\infty } f_{2}(u) d u; в противном случае положим Z=YZ = Y. Покажите, что ZZ имеет функцию плотности fXf_{X}.

?
Задача 4.11.16

Функцию F−1F^{-1}, определённую в теореме (4.11.1), можно назвать квантильной функцией рассматриваемой случайной величины.

?
(a)

Для случайной величины UU с равномерным распределением на (0,1)(0,1) и непрерывной возрастающей функции G:(0,1)→RG:(0,1) \rightarrow \mathbb {R} покажите, что квантильная функция G(U)G(U) есть GG.

(b)

Для случайной величины XX с функцией распределения FF и плотностью ff определим HX(v)=E[F−1(U)∣U>v]H_{X}(v) = \mathbb {E}\left[F^{-1}(U) \mid U > v\right], где UU такая же, как в пункте (a). Преобразованием Харди-Литлвуда величины XX называется HX(V)H_{X}(V), где VV равномерно распределена на (0,1)(0,1) и независима от UU. Покажите, что квантильная функция HX(V)H_{X}(V) есть (1−v)−1∫v∞F−1(w)dw(1-v)^{-1} \int_{v}^{\infty } F^{-1}(w) d w для v∈(0,1)v \in (0,1).

Задача 4.11.17

Для произвольной (не обязательно непрерывной) функции распределения F=FXF = F_{X} определим

G(x,u)=P(X<x)+uP(X=x)=F(x−)+u(F(x)−F(x−)),u∈[0,1] \begin{aligned} G(x, u) & = \mathbb {P}\left(X < x\right)+u \mathbb {P}\left(X = x\right) \\ & = F(x-)+u(F(x)-F(x-)), \quad u \in [0,1] \end{aligned}

Пусть UU равномерно распределена на (0,1)(0,1) и независима от XX. Назовём G(X,U)G(X, U) обобщённым распределительным преобразованием величины XX. Покажите, что W=G(X,U)W = G(X, U) равномерно распределена на (0,1)(0,1), и что случайная величина F−1(W)F^{-1}(W) имеет функцию распределения FF.

?
Задача 4.11.18

Пусть (X1,X2,…,Xn)\left(X_{1}, X_{2}, \ldots , X_{n}\right) — случайный вектор с совместной функцией распределения FF и непрерывными маргинальными функциями распределения F1,F2,…,FnF_{1}, F_{2}, \ldots , F_{n}. Покажите, что существует функция распределения «копульного вида» G(u1,u2,…,un)G\left(u_{1}, u_{2}, \ldots , u_{n}\right), такая что F(x1,x2,…,xn)=G(F1(x1),F2(x2),…,Fn(xn))F\left(x_{1}, x_{2}, \ldots , x_{n}\right) = G\left(F_{1}\left(x_{1}\right), F_{2}\left(x_{2}\right), \ldots , F_{n}\left(x_{n}\right)\right). [Функция распределения «копульного вида» — это такая функция, все маргинальные распределения которой равномерны на [0,1][0,1]].

?