Estadísticas y Big Data

¿Existe un conjugado previo para la distribución de Laplace ? Si no, ¿existe una expresión de forma cerrada conocida que se aproxime a la posterior para los parámetros de la distribución de Laplace? He buscado mucho en Google sin éxito, así que mi suposición actual es "no" en las preguntas …

13 bayesian  conjugate-prior  laplace-distribution 

He escrito un código que puede hacer el filtrado de Kalman (utilizando varios filtros de tipo Kalman diferentes [Filtro de información et al.]) Para el Análisis de espacio de estado lineal gaussiano para un vector de estado n-dimensional. Los filtros funcionan muy bien y estoy obteniendo una buena salida. Sin …

13 maximum-likelihood  algorithms  kalman-filter  likelihood 

¿Se supone que las bandas de confianza y predicción alrededor de una regresión no lineal son simétricas alrededor de la línea de regresión? Lo que significa que no toman la forma del reloj de arena como en el caso de las bandas para la regresión lineal. ¿Porqué es eso? Aquí …

13 confidence-interval  nonlinear-regression  prediction-interval 

¿Cuál es la distribución máxima de entropía para una variable continua positiva, dados sus primeros y segundos momentos? Por ejemplo, una distribución gaussiana es la distribución de entropía máxima para una variable ilimitada, dada su media y desviación estándar, y una distribución Gamma es la distribución de entropía máxima para …

13 distributions  standard-deviation  mean  maximum-entropy 

Estoy usando la biblioteca R 'multcomp' ( http://cran.r-project.org/web/packages/multcomp/ ) para calcular la prueba de Dunnett. Estoy usando el script a continuación: Group <- factor(c("A","A","B","B","B","C","C","C","D","D","D","E","E","F","F","F")) Value <- c(5,5.09901951359278,4.69041575982343,4.58257569495584,4.79583152331272,5,5.09901951359278,4.24264068711928,5.09901951359278,5.19615242270663,4.58257569495584,6.16441400296898,6.85565460040104,7.68114574786861,7.07106781186548,6.48074069840786) data <- data.frame(Group, Value) aov <- aov(Value ~ Group, data) summary(glht(aov, linfct=mcp(Group="Dunnett"))) Ahora, si ejecuto este script a través de la Consola R …

13 r  multiple-comparisons 

Después de pasar por algunas matemáticas ligeramente concisas, creo que tengo una ligera intuición de la estimación de la densidad del núcleo. Pero también soy consciente de que estimar la densidad multivariada para más de tres variables podría no ser una buena idea, en términos de las propiedades estadísticas de …

13 nonparametric  pdf  kernel-smoothing  bivariate  density-estimation 

Suponga que son variables aleatorias iid que siguen la distribución de Poisson con media . ¿Cómo puedo probar que no hay un estimador imparcial de la cantidad ?X0,X1, ... ,XnorteX0,X1,...,Xnorte X_{0},X_{1},\ldots,X_{n} λλ \lambda 1λ1λ \dfrac{1}{\lambda}

13 probability  poisson-distribution  unbiased-estimator  estimators  iid 

Introducción: Habiendo notado la atención recibida hoy por esta pregunta, " ¿Puede ANOVA ser significativo cuando ninguna de las pruebas t por pares lo es? ", Pensé que podría volver a enmarcarlo de una manera interesante que mereciera su propio conjunto de respuestas. . Una variedad de resultados incongruentes (al …

13 hypothesis-testing  anova  statistical-significance  t-test  multiple-comparisons 

He ejecutado una regresión múltiple en la que el modelo en su conjunto es significativo y explica aproximadamente el 13% de la varianza. Sin embargo, necesito encontrar la cantidad de varianza explicada por cada predictor significativo. ¿Cómo puedo hacer esto usando R? Aquí hay algunos datos de muestra y código: …

13 r  regression  variance 

En R (2.15.2) instalé una vez un ARIMA (3,1,3) en una serie de tiempo y una vez un ARMA (3,3) en la serie de tiempo una vez diferenciada. Los parámetros ajustados difieren, lo que atribuí al método de ajuste en ARIMA. Además, ajustar un ARIMA (3,0,3) en los mismos datos …

13 r  time-series  arima  fitting  arma 

He estado trabajando en un modelo logístico y tengo algunas dificultades para evaluar los resultados. Mi modelo es un logit binomial. Mis variables explicativas son: una variable categórica con 15 niveles, una variable dicotómica y 2 variables continuas. Mi N es grande> 8000. Estoy tratando de modelar la decisión de …

13 r  logistic  goodness-of-fit  roc 

Estoy trabajando en algo como el siguiente problema. Tengo un montón de usuarios y N libros. Cada usuario crea una clasificación ordenada de todos los libros que ha leído (que probablemente sea un subconjunto de los N libros), por ejemplo, Libro 1> Libro 40> Libro 25. Ahora quiero convertir estas …

13 multiple-comparisons  ranking  recommender-system  bradley-terry-model 

Aparentemente, es el caso de que si , entoncesXi∼N(0,1)Xi∼N(0,1)X_i \sim N(0,1) X1X2+X3X4∼Laplace(0,1)X1X2+X3X4∼Laplace(0,1)X_1 X_2 + X_3 X_4 \sim \mathrm{Laplace(0,1)} He visto artículos sobre formas cuadráticas arbitrarias, que siempre resultan en horribles expresiones de chi-cuadrado no central. La simple relación anterior no me parece nada obvia, así que (si es cierto) ¿alguien …

13 normal-distribution  proof  quadratic-form  laplace-distribution 

Frecuentemente solía diagnosticar mis datos multivariados usando PCA (datos ómnicos con cientos de miles de variables y docenas o cientos de muestras). Los datos a menudo provienen de experimentos con varias variables categóricas independientes que definen algunos grupos, y a menudo tengo que pasar por algunos componentes antes de poder …

13 pca  linear-model  r-squared  discriminant-analysis 

¿Es posible realizar un análisis de potencia para la prueba U de Kruskal-Wallis y Mann-Whitney? En caso afirmativo, ¿hay algún paquete / función de R que lo realice?

13 r  nonparametric  power-analysis  kruskal-wallis