Me pregunto cómo puedo simular una muestra de n vidas de distribución de Weibull que incluya observaciones censuradas a la derecha de Tipo I. Por ejemplo, tengamos n = 3, forma = 3, escala = 1 y la tasa de censura = .15, y el tiempo de censura = .88. Sé cómo generar una muestra de Weibull pero no...
Estoy buscando generar 450 puntos de datos en R. Hay tres conjuntos distintos 150 de cada uno distribuidos en una banda circular con diferentes radios (en 1, 2.8 y 5). En particular, estoy buscando reproducir el primer gráfico en p546 de Los elementos del aprendizaje estadístico. Estaría muy...
Sé cómo generar una secuencia con media . Por ejemplo, en Matlab, si quiero generar una secuencia de longitud , es:0 ± 1 10000±1±1\pm 1000±1±1\pm 1100001000010000 2*(rand(1, 10000, 1)<=.5)-1 Sin embargo, ¿cómo generar una secuencia con una media de , es decir, con siendo ligeramente...
En este documento , que se refiere al comando "establecer semilla", las personas de Stata discuten cuestiones relacionadas con la configuración de semillas al generar números pseudoaleatorios. Un notable "no" es "no use en serie la secuencia de números naturales como semillas, porque tiene un...
¿Por qué los generadores de números aleatorios como runif()en R no generan el mismo resultado cada vez? Por ejemplo: X <- runif(100) X está generando diferentes salidas cada vez. ¿Cuál es la razón para generar diferentes salidas cada vez? ¿Qué funcionalidades tiene en segundo plano para...
Suponga un archivo de datos con más de 80 millones de unos y ceros, generados aleatoriamente. A partir de este archivo, queremos crear una lista de enteros decimales aleatorios. Este es el plan para hacer esta conversión. Divida los 80 millones de dígitos en grupos de 4 dígitos binarios....
Cerrado. Esta pregunta está fuera de tema . Actualmente no está aceptando respuestas. ¿Quieres mejorar esta pregunta? Actualice la pregunta para que esté en el tema de Cross Validated. Cerrado hace 2 años . Mi objetivo: Me gustaría tener una función que...
¿Cómo debo muestrear eficientemente de la siguiente distribución? x ∼ B ( α , β) , x > k x∼B(α,β), x>k x \sim B(\alpha, \beta),\space x > k Si no es demasiado grande, entonces el muestreo de rechazo puede ser el mejor enfoque, pero no estoy seguro de cómo proceder cuando es grande. ¿Quizás...
Estoy ejecutando una simulación en R y un grupo de computadoras y tengo el siguiente problema. En cada una de las X computadoras que ejecuto: fxT2 <- function(i) runif(10) nessay <- 100 c(mclapply(1:nessay, fxT2), recursive=TRUE) Hay 32 computadoras, cada una con 16 núcleos. Sin embargo,...
En mi programa, necesito ejecutar N subprocesos separados, cada uno con su propio RNG, que se utiliza para muestrear un gran conjunto de datos. Necesito poder sembrar todo este proceso con un solo valor para poder reproducir los resultados. ¿Es suficiente simplemente aumentar secuencialmente la...
Necesito crear vectores aleatorios de números reales a_i que satisfagan las siguientes restricciones: abs(a_i) < c_i; sum(a_i)< A; # sum of elements smaller than A sum(b_i * a_i) < B; # weighted sum is smaller than B aT*A*a < D # quadratic multiplication with A smaller than D where...
Necesito dibujar números aleatorios de una distribución log-cauchy que tiene densidad: ¿Alguien puede ayudarme o señalarme un libro / documento que pueda mostrarme cómo?F( x ; μ , σ) = 1x πσ[ 1 + ( l n ( x ) -
Tengo un GLMM de la forma: lmer(present? ~ factor1 + factor2 + continuous + factor1*continuous + (1 | factor3), family=binomial) Cuando lo uso drop1(model, test="Chi"), obtengo resultados diferentes a los que uso Anova(model, type="III")del paquete del automóvil o summary(model). Estos dos...
Dada una moneda con sesgo desconocido , ¿cómo puedo generar variables, de la manera más eficiente posible, que estén distribuidas por Bernoulli con probabilidad 0.5? Es decir, usar el número mínimo de vueltas por variante
Estoy tratando de ajustar un modelo de tiempo discreto en R, pero no estoy seguro de cómo hacerlo. He leído que puede organizar la variable dependiente en diferentes filas, una para cada observación de tiempo, y usar la glmfunción con un enlace logit o cloglog. En este sentido, tengo tres...
He estado usando secuencias de baja discrepancia durante un tiempo para Distribuciones uniformes, ya que he encontrado que sus propiedades son útiles (principalmente en gráficos de computadora por su apariencia aleatoria y su capacidad para cubrir densamente [0,1] de manera incremental). Por...
Necesito generar una lista de variables aleatorias sujetas a restricciones que se pueden expresar en la forma donde es una matriz si tiene entradas. En todos los casos con los que estoy tratando, , por ejemplo, será de alrededor de 14,000 será de 50. No estoy seguro de qué método utilizaré para el...
Actualmente estoy escribiendo un algoritmo para privacidad diferencial usando el mecanismo de Laplace. Lamentablemente no tengo antecedentes en estadísticas, por lo tanto, muchos términos son desconocidos para mí. Así que ahora estoy tropezando con el término: ruido de Laplace . Para que un...
Nunca tuve un curso formal de estadística, pero debido a mi línea de investigación constantemente encuentro artículos que aplican varios conceptos estadísticos. A menudo veré una descripción de un proceso de Monte Carlo aplicado a una situación dada, y por lo que puedo reunir 9 de cada 10 veces se...
Supongamos que tiene una moneda justa que puede lanzar tantas veces como desee (posiblemente infinitamente contable). ¿Es posible generar la distribución uniforme discreta en , donde(1,2,...,k)(1,2,...,k)(1,2,...,k)kkk NO es una potencia de 2? ¿Como lo harias? Si esto es demasiado general,...