Estoy tratando de obtener el previo de Jeffreys para una distribución binomial negativa. No puedo ver dónde me equivoco, así que si alguien pudiera ayudar a señalar eso, sería apreciado. Bien, entonces la situación es la siguiente: debo comparar las distribuciones anteriores obtenidas usando un...
Tengo dos series temporales de datos diarios. Uno es sign-upsy el otro terminationsde suscripciones. Me gustaría predecir esto último utilizando la información contenida en ambas variables. Al observar el gráfico de estas series, es obvio que las terminaciones están correlacionadas con múltiplos...
Versión corta: ¿hay una introducción a los escritos (documentos y libros) de Ronald Fisher sobre estadísticas que está dirigida a aquellos con poca o ninguna experiencia en estadísticas? Estoy pensando en algo así como un "lector Fisher anotado" dirigido a los no estadísticos. Expongo la...
He estado leyendo sobre la estimación de máxima verosimilitud y la estimación máxima a posteriori y hasta ahora he encontrado ejemplos concretos solo con la estimación de máxima verosimilitud. He encontrado algunos ejemplos abstractos de estimación máxima a posteriori, pero nada concreto aún con...
Estoy tratando de usar el neuralnetpaquete de R (documentación aquí ) para la predicción. Aquí lo que estoy tratando de hacer: library(neuralnet) x <- cbind(runif(50, min=1, max=500), runif(50, min=1, max=500)) y <- x[, 1] * x[, 2] train <- data.frame(x, y) n <- names(train) f <-...
No sé cómo se llaman tales tramas y, por lo tanto, le di a esta pregunta un título estúpido. Digamos que tengo un conjunto de datos ordenado de la siguiente manera 4253 4262 4270 4383 4394 4476 4635 ... Cada número corresponde a la cantidad de publicaciones que un determinado usuario contribuyó...
¿Existe una estrategia para elegir la cantidad de árboles en un GBM? Específicamente, el ntreesargumento en Rla gbmfunción de. No veo por qué no debe establecer ntreesel valor razonable más alto. Me di cuenta de que una mayor cantidad de árboles reduce claramente la variabilidad de los resultados...
Usando Amelia en R, obtuve múltiples conjuntos de datos imputados. Después de eso, realicé una prueba de medidas repetidas en SPSS. Ahora, quiero agrupar los resultados de las pruebas. Sé que puedo usar las reglas de Rubin (implementadas a través de cualquier paquete de imputación múltiple en R)...
Usé el siguiente código R para ajustar un modelo probit: p1 <- glm(natijeh ~ ., family=binomial(probit), data=data1) stepwise(p1, direction='backward/forward', criterion='BIC') ¿Quiero saber qué hace stepwisey qué hace backward/forwardexactamente y cómo seleccionar las
Tengo una pregunta sobre la interpretación de la llamada tsboot en R. Revisé la documentación de Kendall y el paquete de arranque, pero no soy más inteligente que antes. Cuando ejecuto un bootstrap usando, por ejemplo, el ejemplo en el paquete Kendall, donde la estadística de prueba es la tau de...
En un modelo jerárquico de datos donde \ lambda \ sim \ textrm {Gamma} (\ alpha, \ beta) parece ser típico en la práctica elegir valores ( \ alpha, \ beta) de modo que la media y la varianza de la distribución gamma coincidan aproximadamente con la media y la varianza de los datos y (por...
Estoy ejecutando una simulación rápida para comparar diferentes métodos de agrupación, y actualmente me encuentro con un problema tratando de evaluar las soluciones de agrupación. Conozco varias métricas de validación (muchas de ellas se encuentran en cluster.stats () en R), pero supongo que es...
Apliqué algunos datos para encontrar la mejor solución de variables del modelo de regresión usando la regresión de cresta en R. He usado lm.ridgey glmnet(cuándo alpha=0), pero los resultados son muy diferentes, especialmente cuando lambda=0. Supongamos que ambos estimadores de parámetros tienen los...
Para los no estadísticos como yo, es muy difícil captar la idea de la VImétrica (variación de la información) incluso después de leer el documento relevante de Marina Melia " Comparación de agrupamientos: una distancia basada en la información " (Journal of Multivariate Analysis, 2007). De hecho,...
Me gustaría comprobar Rsi mis datos se ajustan a las distribuciones log-normal o Pareto. ¿Cómo podría hacer eso? Quizás ks.testpodría ayudarme a hacer eso, pero ¿cómo podría obtener los parámetros y para la distribución de Pareto para mis
He escrito un programa para simular un volado de reproducción aleatoria tarjeta. Cada carta está numerada, con el palo que va CLUBS, DIAMONDS, HEARTS, SPADESy el rango de Dos a Diez y luego Jack, Reina, Rey y As. Por lo tanto, el Two of Clubs tiene un Número de 1, el Three of Clubs un 2 ... As of...
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...
Estoy tratando de calcular una media y una desviación estándar de 2 percentiles para una distribución lognormal. Tuve éxito en realizar el cálculo para una distribución normal usando X = mean + sd * Zy resolviendo para media y sd. Creo que me falta una ecuación cuando intento hacer lo mismo para...
Estoy reproduciendo un ejemplo de modelos generalizados, lineales y mixtos . Mi MWE está abajo: Dilution <- c(1/128, 1/64, 1/32, 1/16, 1/8, 1/4, 1/2, 1, 2, 4) NoofPlates <- rep(x=5, times=10) NoPositive <- c(0, 0, 2, 2, 3, 4, 5, 5, 5, 5) Data <- data.frame(Dilution, NoofPlates,...
Me gustaría corregir manualmente un cierto coeficiente, digamos , luego ajustar los coeficientes a todos los demás predictores, mientras en el modelo.β 1 = 1.0β1=1.0β1=1.0\beta_1=1.0β1=1.0β1=1.0\beta_1=1.0 ¿Cómo puedo lograr esto usando R? En particular, me gustaría trabajar con LASSO ( glmnet) si...