He llevado a cabo una agrupación de puntos de coordenadas (longitud, latitud) y he encontrado resultados sorprendentes y adversos de los criterios de agrupación para el número óptimo de agrupaciones. Los criterios se toman del clusterCrit()paquete. Los puntos que estoy tratando de agrupar en una...
Tengo un conjunto de datos meteorológicos diarios, que tiene, como era de esperar, un efecto estacional muy fuerte. Adapté un modelo ARIMA a este conjunto de datos usando la función auto.arima del paquete de pronóstico. Para mi sorpresa, la función no aplica ninguna operación estacional:...
Mi diseño es el siguiente. yyy es la respuesta de Bernoulli X1X1x_1 es una variable continua X2X2x_2 es una variable categórica (factor) con dos niveles El experimento está completamente dentro de los sujetos. Es decir, cada sujeto recibe cada combinación deX1X1x_1 y X2X2x_2. Esta es una...
Quiero calcular los errores estándar de una distribución hiperbólica ajustada. En mi notación, la densidad viene dada por Estoy usando el paquete HyperbolicDistr en R. Calculo los parámetros mediante el siguiente
Tengo algunas preguntas básicas sobre PCA (análisis de componentes principales) y LDA (análisis discriminante lineal): En PCA hay una manera de calcular la proporción de varianza explicada. ¿También es posible para LDA? ¿Si es así, cómo? ¿Es la salida "Proporción de traza" de la ldafunción (en...
http://cran.r-project.org/web/packages/quadprog/quadprog.pdf El paquete R quadprogparece ser capaz de resolver el problema de programación cuadrática solo cuando la matrizreDD Es positivo definitivo. Sin embargo, hay un caso cuando la matriz reDDNo es positivo definido. como min (X2+y2- 6 x...
¿Es posible realizar una selección aproximada de hiperparámetros completamente bayesianos (1) (por ejemplo, escala de covarianza) con el código GPML, en lugar de maximizar la probabilidad marginal (2)? Creo que usar métodos MCMC para resolver las integrales que involucran hiperparámetros antes...
Entonces, me parece que la función de pesos en lm da más peso a las observaciones cuanto mayor sea el valor de "peso" de la observación asociada, mientras que la función lme en lme hace exactamente lo contrario. Esto se puede verificar con una simple simulación. #make 3 vectors- c is used as an...
Tengo dos preguntas, Pregunta 1: ¿Cómo puedo demostrar que la distribución posterior es una distribución beta si la probabilidad es binomial y la anterior es beta? Pregunta 2: ¿Cómo afectan las elecciones los parámetros anteriores a los posteriores? ¿No deberían ser todos iguales? ¿Es posible...
Sé que es posible cuando realiza una regresión lineal múltiple calcular el intervalo de confianza en el R cuadrado y en el R cuadrado ajustado. ¿Alguien sabe cómo hacerlo con
Tengo un conjunto de datos que contiene unos cientos de transacciones de tres proveedores que operan en más de 100 países durante un período de tres años. Hemos descubierto que el país de ventas no es un factor significativo en los precios alcanzados (los productos son más o menos productos...
Esta pregunta ha surgido en otro hilo que comencé, así que pensé que obtendría la opinión de más personas al respecto. Mi pregunta es ¿Es el residual, e, un estimador del error, ?ϵϵ\epsilon La razón por la que pregunto es la siguiente. En OLS, la varianza de los residuos, , se conoce como la...
Esta es una derivación de esta pregunta: ¿Cómo comparar dos grupos con múltiples mediciones para cada individuo con R? En las respuestas allí (si entendí correctamente) aprendí que la varianza dentro del sujeto no afecta las inferencias hechas sobre las medias grupales y que está bien simplemente...
Tengo el siguiente tipo de datos (codificados en R): v.a = c('cat', 'dog', 'dog', 'goat', 'cat', 'goat', 'dog', 'dog') v.b = c(1, 2, 1, 2, 1, 2, 1, 2) v.c = c('blue', 'red', 'blue', 'red', 'red', 'blue', 'yellow', 'yellow') set.seed(12) v.d = rnorm(8) aov(v.a ~ v.b + v.c + v.d) # Error Me...
Estoy tratando de hacer SVM de una clase en R. He estado tratando de usar el paquete kernlab e1071 / ksvm. Pero no estoy seguro si lo estoy haciendo correctamente. ¿Hay algún ejemplo de trabajo para SVM de una clase en R? También, Estoy dando una gran matriz de predictores como X. Dado que se...
Estoy ejecutando un CFA y obtengo buenos índices de ajuste (CFI = .99, RMSEA = .01) para una escala unidimensional. Sin embargo, cuando pruebo la consistencia interna, obtengo los s ( ) de Cronbach pobres . He intentado todo, desde eliminar valores atípicos hasta soltar elementos y aún así terminar...
Estoy realizando regresiones cuantiles en R usando el paquete quantreg. Mi conjunto de datos incluye 12,328 observaciones que van desde 0,12 a 330. Los puntos de tiempo para mis datos no son exactamente continuos; Todos los datos caen en una de las pocas docenas de contenedores que van desde 73 a...
He calculado los valores de ICC para dos grupos y ahora me gustaría comparar los valores de ICC para determinar si los grupos difieren en su repetibilidad. En la literatura, la gente simplemente ha utilizado las pruebas t para comparar la repetibilidad, pero no me queda claro cómo hacerlo. Por...
¿Hay alguna forma de que pueda realizar LASSO con regresión binomial negativa en R? Estoy realizando una regresión binomial negativa en mi conjunto de datos porque los datos están demasiado dispersos para imponer la regresión de Poisson. Mientras tanto, también me enfrento a un problema de...
Tengo la siguiente tabla en R df <- structure(list(x = structure(c(12458, 12633, 12692, 12830, 13369, 13455, 13458, 13515), class = "Date"), y = c(6080, 6949, 7076, 7818, 0, 0, 10765, 11153)), .Names = c("x", "y"), row.names = c("1", "2", "3", "4", "5", "6", "8", "9"), class =...