¿Cómo puedo calcular la sensibilidad y la especificidad (o medidas análogas) de una prueba de diagnóstico continuo para predecir un resultado continuo (por ejemplo, presión arterial) sin dicotomizar el resultado? ¿Algunas ideas? Parece que los investigadores han hecho esto usando modelos de...
Tengo problemas para encontrar una solución con respecto a cómo ejecutar una prueba post-hoc (Tukey HSD) después de un ANOVA de medidas repetidas de 2 factores (ambos dentro de los sujetos) en R. Para el ANOVA, he usado la función aov: summary(aov(dv ~ x1 * x2 + Error(subject/(x1*x2)),...
Estoy usando el lme4paquete en R para hacer algunos modelos logísticos de efectos mixtos. Entiendo que la suma de cada efecto aleatorio debe ser cero. Cuando hago modelos mixtos lineales de juguete usando lmer, los efectos aleatorios generalmente son < confirmando mi creencia de que el But en...
Estoy tratando de entender cuál es la diferencia entre un HMM estándar y un HMM bayesiano. Wikipedia solo menciona brevemente cómo se ve el modelo, pero necesito un tutorial más detallado. ¿Alguien sabe de un documento o una implementación que pueda ver? También tengo problemas con la terminología...
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...
De Wikipedia El modelo lineal general (GLM) es un modelo lineal estadístico. Puede escribirse como 1 Y=XB+U,Y=XB+U, \mathbf{Y} = \mathbf{X}\mathbf{B} + \mathbf{U}, dónde YYY es una matriz con series de mediciones multivariadas, XXX es una matriz que podría ser una matriz de diseño, BBB es...
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...
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...
Mi pregunta general es: ¿por qué usar en bayesglmlugar de otros métodos de clasificación? Nota: Solo me interesa la predicción. Tengo una cantidad decente de datos (~ 100,000 obs). Siento que el tamaño de la muestra es lo suficientemente grande como para que los parámetros de una regresión...
La subitización es la enumeración rápida y precisa de las pantallas de baja numerosidad, que se distingue del conteo por una no linealidad aguda en la gráfica de tiempos de respuesta. A continuación se muestra una trama representativa, de Watson, DG, Maylor, EA y Bruce, LAM (2007). Observe que los...
Estoy tratando de llegar a un modelo utilizando la regresión binomial negativa (GLM binomial negativa). Tengo un tamaño de muestra relativamente pequeño (mayor que 300), y los datos no están escalados. Noté que hay dos formas de medir la bondad de ajuste: una es la desviación y la otra es la...
¿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 =...
Tengo el siguiente modelo lineal generalizado. El objeto glmDVse modela como una proporción de éxitos sobre las pruebas totales. Los objetos x_ison variables continuas. ¿Cómo se ve esto en notación matemática? winp.glm = glm(glmDV ~ x1 + x2 + x3 + x4 + x5 + x6 + x7, data=myData,...
Tengo un conjunto de datos en el que estoy tratando de hacer una regresión, y estoy fallando. La situación: Miles de operadores de robots de batalla están librando batallas entre ellos utilizando robots de batalla. Algunos robots de batalla son fuertes y poderosos, y otros son débiles; los...
Mi proyecto actual puede requerir que construya un modelo para predecir el comportamiento de un determinado grupo de personas. el conjunto de datos de entrenamiento contiene solo 6 variables (la identificación es solo para fines de identificación): id, age, income, gender, job category, monthly...
Estoy trabajando en modelos de datos de recuento cero inflados usando el psclpaquete. ¡Me pregunto por qué no hay desarrollo de modelos para modelos de datos de conteo inflados! Además, ¿por qué no hay desarrollo de modelos de datos de conteo bimodales, digamos cero y 2 inflados! Una vez que generé...
Dejar E(h)E(h)\mathcal{E(h)} una función que dada alguna hipótesis hhh devuelve el error de generalización para ese fijo hhh. Estaba leyendo algunas notas sobre la selección del modelo y el error de generalización y decía: "Si tuviéramos acceso a E(h)E(h)\mathcal{E(h)}, tampoco habría problema...
Entiendo el modelo AR (p): su entrada es la serie temporal que se está modelando. Estoy completamente atascado cuando leo sobre el modelo MA (q): su aportación es innovación o choque aleatorio, ya que a menudo se formula. El problema es que no puedo imaginar cómo obtener un componente de...
Esta es una pregunta bastante genérica: suponga que quiero construir un modelo para predecir la próxima observación basada en la anterior. norteNN observaciones (norteNNpuede ser un parámetro para optimizar experimentalmente). Básicamente, tenemos una ventana deslizante de características de...