Вопросы с тегом «random-generation»

11
Как создать матрицу случайной корреляции, которая имеет приблизительно нормально распределенные недиагональные записи с заданным стандартным отклонением?

Я хотел бы создать матрицу случайной корреляции так, чтобы распределение ее недиагональных элементов выглядело примерно как нормальное. Как я могу это сделать? Мотивация такая. Для набора данных из временных рядов распределение корреляции часто выглядит достаточно близким к нормальному. Я хотел бы...

11
Какие числа наименее вероятно будут выбраны людьми в лотерее?

Mega Millions сегодня составляет более 500 миллионов долларов. Я помню, как читал газету JSTOR о некоторых числах, которые вряд ли будут выбраны. Например, многие люди выбирают 7, потому что это их счастливое число, и я хочу противоположность этому. Однако мое членство в JSTOR закончилось. Какие...

11
Генерация случайных чисел по Коши

Мне нужно нарисовать случайные числа из логарифмического распределения, которое имеет плотность: Может кто-нибудь помочь мне или указать мне книгу / бумагу, которая может показать мне, как?е( х ; μ , σ) = 1х πσ[ 1 + ( l n ( x ) - μσ)2],е(Икс;μ,σ)знак...

11
Как симулировать цензурированные данные

Мне интересно, как я могу смоделировать выборку из n времен жизни распределения Вейбулла, которые включают наблюдения типа I с правой цензурой. Например, пусть n = 3, shape = 3, scale = 1 и уровень цензуры = .15, а время цензуры = .88. Я знаю, как сгенерировать выборку Вейбулла, но я не знаю, как...

11
Ссылки и лучшие практики для установки семян в генерации псевдослучайных чисел

В этом документе , касающемся команды «set seed», сотрудники Stata обсуждают вопросы, связанные с установкой seed при генерации псевдослучайных чисел. Примечательным «не» является «не использовать последовательно последовательность натуральных чисел в качестве начальных чисел, потому что это имеет...

11
Почему runif не генерирует один и тот же результат каждый раз?

Почему генераторы случайных чисел, такие как runif()в R, не генерируют один и тот же результат каждый раз? Например: X <- runif(100) X генерирует разные выходы каждый раз. В чем причина создания разных выходов каждый раз? Какие функции у него есть в фоновом режиме, чтобы сделать...

10
Модель истории дискретного времени (выживания) в R

Я пытаюсь вписать модель с дискретным временем в R, но я не уверен, как это сделать. Я читал, что вы можете организовать зависимую переменную в разных строках, по одной для каждого временного наблюдения, и использовать glmфункцию со ссылкой logit или cloglog. В этом смысле, у меня есть три колонки:...

10
Лучший способ заполнить N независимых генераторов случайных чисел от 1 значения

В моей программе мне нужно запустить N отдельных потоков, каждый с собственным RNG, который используется для выборки большого набора данных. Мне нужно иметь возможность заполнить весь этот процесс одним значением, чтобы я мог воспроизвести результаты. Достаточно ли просто последовательно...

10
Учитывая монету с неизвестным уклоном, эффективно генерировать отличия от справедливой монеты

Учитывая монету с неизвестным смещением , как я могу генерировать переменные - настолько эффективно, насколько это возможно - которые распределены Бернулли с вероятностью 0,5? То есть, используя минимальное количество сальто на сгенерированный...

10
Почему Anova () и drop1 () предоставили разные ответы для GLMM?

У меня есть GLMM формы: lmer(present? ~ factor1 + factor2 + continuous + factor1*continuous + (1 | factor3), family=binomial) Когда я использую drop1(model, test="Chi"), я получаю другие результаты, чем если бы я использовал Anova(model, type="III")из пакета автомобиля или summary(model). Последние...

10
С адреса электронной почты до квази-случайного числа [закрыто]

Закрыто. Этот вопрос не по теме . В настоящее время он не принимает ответы. Хотите улучшить этот вопрос? Обновите вопрос, чтобы он соответствовал теме перекрестной проверки. Закрыто 2 года назад . Моя цель: Я хотел бы иметь функцию, которая берет адрес электронной почты и выводит квазислучайное...

10
RNG, R, mclapply и кластер компьютеров

Я запускаю симуляцию на R и кластере компьютеров и имею следующую проблему. На каждом из компьютеров X я запускаю: fxT2 <- function(i) runif(10) nessay <- 100 c(mclapply(1:nessay, fxT2), recursive=TRUE) Есть 32 компьютера, каждый с 16 ядрами. Тем не менее, около 2% случайных чисел идентичны....

10
Эффективная выборка порогового бета-распределения

Как мне выбрать образец из следующего дистрибутива? x ∼ B ( α , β) , х > к x∼B(α,β), x>k x \sim B(\alpha, \beta),\space x > k Если не слишком велико, то выборка отклонения может быть лучшим подходом, но я не уверен, как действовать, когда k велико. Возможно, есть какое-то асимптотическое...

9
Монте-Карло == применяется случайный процесс?

У меня никогда не было курса официальной статистики, но из-за моего направления исследований я постоянно сталкиваюсь со статьями, которые применяют несколько статистических концепций. Часто я вижу описание процесса Монте-Карло, применяемого к данной ситуации, и то, что я могу собрать 9 из 10 раз,...

9
Что подразумевается под «шумом Лапласа»?

В настоящее время я пишу алгоритм для дифференциальной конфиденциальности с использованием механизма Лапласа. К сожалению, у меня нет опыта в статистике, поэтому многие термины мне неизвестны. Так что теперь я спотыкаюсь о термине: шум Лапласа . Чтобы сделать дифференциальный набор данных закрытым,...

9
Реальные примеры неэффективного генератора случайных чисел

Все мы знаем, что генераторы случайных чисел в компьютерах не генерируют истинные случайные числа, а вместо этого генерируют псевдослучайные числа. Кроме того, некоторые RNG лучше, чем другие, а некоторые реализованы лучше, чем другие. Каковы некоторые примеры того, когда использовался плохой ГСЧ...

9
Имитация данных в соответствии с моделью посредничества

Я заинтересован в поиске процедуры для моделирования данных, которые соответствуют указанной модели посредничества. В соответствии с общей структурой модели линейных структурных уравнений для тестирования моделей посредничества, впервые описанной Barron и Kenny (1986) и описанной в других местах,...

9
Генерация случайных величин, удовлетворяющих ограничениям

Мне нужно сгенерировать список случайных величин учетом ограничений, которые можно выразить в виде где - это матрица если имеет записей. Во всех случаях, с которыми я имею дело, , например будет около 14 000, а будет 50. Я не уверен, какой метод я буду использовать для случайной выборки, будь то...

9
Как отобрать усеченное многочленное распределение?

Мне нужен алгоритм для выборки усеченного полиномиального распределения. Это, Икс⃗ ∼ 1ZпИкс11… РИксККИкс1! … ХК!Икс→~1Zп1Икс1...пКИксКИкс1!...ИксК!\vec x \sim \frac{1}{Z} \frac{p_1^{x_1} \dots p_k^{x_k}}{x_1!\dots x_k!} где - нормализационная константа, имеет положительных компонент и . Я...