Sé que R no es particularmente útil para analizar grandes conjuntos de datos dado que R carga todos los datos en la memoria, mientras que algo como SAS hace análisis secuenciales. Dicho esto, hay paquetes como bigmemory que permiten a los usuarios realizar análisis de datos grandes (análisis...
Utilizando el excelente paquete de pronósticos de Rob Hyndman, me encontré con la necesidad no solo de tener intervalos de predicción, sino también de simular una serie de caminos futuros, dadas las observaciones pasadas de una serie temporal con estacionalidades complejas. Hay algo para series de...
Tengo un conjunto variable de respuestas que se expresan como un intervalo, como la muestra a continuación. > head(left) [1] 860 516 430 1118 860 602 > head(right) [1] 946 602 516 1204 946 688 donde izquierda es el límite inferior y derecha es el límite superior de la respuesta. Quiero...
Estoy tratando de calcular esta distribución
En la predicción de series de tiempo utilizando varios modelos como AR, MA, ARMA, etc., generalmente nos enfocamos en el modelado de los datos en el cambio de tiempo. Pero cuando tenemos 2 series de tiempo en las que el coeficiente de correlación de Pearson muestra que están altamente...
Necesito calcular el Dirichlet CDF , pero solo puedo encontrar implementaciones del PDF . ¿Saben ustedes de alguna biblioteca (preferiblemente en R)
Tengo algunos conjuntos de datos de EEG que estoy probando en dos clases. Puedo obtener una tasa de error decente de LDA (las distribuciones condicionales de clase no son gaussianas, pero tienen colas similares y una separación suficientemente buena), por lo que quiero trazar el ROC del predictor...
Esta es una pregunta para principiantes sobre un ejercicio en "Computación bayesiana con R" de Jim Albert. Tenga en cuenta que si bien esto podría ser tarea, en mi caso no lo es, ya que estoy aprendiendo métodos bayesianos en R porque creo que podría usarlo en mis futuros análisis. De todos modos,...
Escribo esta pregunta con referencia a un ejemplo en p138-142 del siguiente documento: ftp://ftp.software.ibm.com/software/analytics/spss/documentation/amos/20.0/en/Manuals/IBM_SPSS_Amos_User_Guide.pdf . Aquí hay figuras ilustrativas y una tabla: Entiendo que la variable latente no tiene una...
Cuando ejecuto este código: require(nlme) a <- matrix(c(1,3,5,7,4,5,6,4,7,8,9)) b <- matrix(c(3,5,6,2,4,6,7,8,7,8,9)) res <- lm(a ~ b) print(summary(res)) res_gls <- gls(a ~ b) print(summary(res_gls)) Obtengo los mismos coeficientes y la misma significación estadística en los...
Me gustaría hacer análisis de poder para pruebas de hipótesis de (no) igualdad de proporciones en las cuales las proporciones son muy pequeñas. Me gustaría hacerlo sin usar aproximaciones normales (o Poisson) de la distribución binomial. Hay varios tipos generales de preguntas de poder que me...
Permítame hacer una pregunta básica. Entiendo la mecánica de Naive Bayes para variables discretas, y puedo rehacer los cálculos "a mano". (código de HouseVotes84 hasta el final por debajo). Sin embargo, me cuesta ver cómo funcionan los mecanismos para las variables continuas (código de ejemplo a...
Tengo un parámetro que se encuentra entre . Digamos que puedo ejecutar un experimento y obtener , donde es un gaussiano estándar. Lo que necesito es una estimación de que es 1) imparcial 2) casi seguramente acotada. El requisito (2) es crucial para
Estoy tratando de ejecutar una regresión de Cox en un conjunto de datos de 2,000,000 de filas de la siguiente manera usando solo R. Esta es una traducción directa de un PHREG en SAS. La muestra es representativa de la estructura del conjunto de datos original. ## library(survival) ### Replace...
Estoy trabajando en un modelo que depende de una función parametrizada fea que actúa como una función de calibración en una parte del modelo. Utilizando una configuración bayesiana, necesito obtener antecedentes no informativos para los parámetros que describen mi función. Sé que, idealmente,...
Tengo 2 estudios que analizan la respuesta del paciente al mismo medicamento. El estudio 1 encontró 10,000 genes expresados sobre el fondo y 500 de ellos se expresan de manera diferencial y se conocen como la firma de respuesta al medicamento. El estudio 2 encontró 1,000 genes que representan la...
Necesito un único conjunto de datos imputados (por ejemplo, para crear un grupo ficticio del grupo imputado a partir de los datos de ingreso per cápita del país imputado). R ofrece paquetes de paquetes para crear múltiples datos imputados (por ejemplo, Amelia) y combinar resultados de múltiples...
Quiero incluir pesos de muestra en mi modelo de regresión cuantil, pero no estoy seguro de cómo hacerlo. Ya he definido mi peso, que son los pesos replicados que ya figuran en el conjunto de datos de la encuesta (calculado en el paquete de la encuesta):
Yo mismo siempre usaría la media geométrica para estimar una mediana lognormal. Sin embargo, en el mundo de la industria, a veces usar la mediana de la muestra da mejores resultados. Por lo tanto, la pregunta es, ¿existe un rango / punto de corte a partir del cual la mediana de la muestra se pueda...
¿Alguien sabe un método para comparar dos variables con una prueba de chi cuadrado si las variables son de diferentes encuestas con diferentes svydesign()declaraciones? Estoy buscando probar una diferencia en una distribución variable en dos ondas de una encuesta, pero la svychisq()declaración se...