Si hay 0 en la tabla de contingencia y estamos ajustando modelos poisson / loglineales anidados (usando la glmfunción de R ) para una prueba de razón de probabilidad, ¿necesitamos ajustar los datos antes de ajustar los modelos glm (por ejemplo, agregar 1/2 a todos los recuentos)? Obviamente,...
¿Cómo interpreta una curva de supervivencia del modelo de riesgo proporcional de Cox? En este ejemplo de juguete, supongamos que tenemos un modelo de riesgo proporcional de Cox ageen kidneydatos variables y generamos la curva de supervivencia. library(survival) fit <- coxph(Surv(time,...
Los modelos penalizados se pueden usar para estimar modelos donde el número de parámetros es igual o incluso mayor que el tamaño de la muestra. Esta situación puede surgir en modelos log-lineales de grandes tablas dispersas de datos categóricos o de conteo. En estos entornos, a menudo también es...
La situación Tengo un conjunto de datos con una dependiente una variable independiente . Quiero ajustar una regresión lineal continua por partes con puntos de interrupción conocidos / fijos que ocurren en . Los breakpoins se conocen sin incertidumbre, por lo que no quiero estimarlos. Luego ajusto...
Tengo datos del siguiente diseño experimental: mis observaciones son recuentos del número de éxitos ( K) del número correspondiente de ensayos ( N), medidos para dos grupos, cada uno compuesto de Iindividuos, de Ttratamientos, donde en cada combinación de factores hay Rréplicas . Por lo tanto, en...
Estoy tratando de ajustar un modelo logarítmico lineal a una gran cantidad de variables de los datos de la encuesta. Hay algunas razones por las que podría ser preferible ajustar las regresiones logísticas a esos datos. Varias autoridades sugieren que estos son equivalentes. Sin embargo, tengo...
Podemos simular una regresión lineal sin aleatoriedad, lo que significa que hacemos lugar de . Entonces, si ajustamos un modelo lineal, los coeficientes serán idénticos a la "verdad fundamental". Aquí hay un ejemplo.y= Xβy=Xβy=X\betay= Xβ+ ϵy=Xβ+ϵy=X\beta+\epsilon set.seed(0) n <- 1e5 p <-...
Tengo dos grupos de 10 participantes que fueron evaluados tres veces durante un experimento. Para probar las diferencias entre los grupos y entre las tres evaluaciones, ejecuté un ANOVA de diseño mixto 2x3 con group(control, experimental), time(primero, segundo, tres) y group x time. Ambos timey...
Hemos ajustado una función no lineal a los datos observados. El siguiente paso debe ser la evaluación de la bondad de ajuste de esta función (como para modelos lineales).R2R2R^2 ¿Cuáles son las formas habituales de medir esto? Editar 1: El ajuste se realizó de la siguiente manera: Realizar...
Estoy implementando un clasificador SVM no lineal con kernel RBF. Me dijeron que la única diferencia con respecto a un SVM normal era que simplemente tenía que reemplazar el producto de puntos con una función del núcleo: Sé cómo funciona un SVM lineal normal, es decir, después de resolver el...
Estoy leyendo un texto, "Probabilidad y estadística" de Devore. Estoy viendo 2 ítems en la página 740: el valor esperado y la varianza de la estimación de β1β1\beta_1 , que es el parámetro de pendiente en la regresión lineal Yi=β0+β1Xi+ϵiYi=β0+β1Xi+ϵiY_i = \beta_0 + \beta_1 X_i + \epsilon_i ....
Tengo datos que describen con qué frecuencia tiene lugar un evento durante una hora ("número por hora", nph) y cuánto duran los eventos ("duración en segundos por hora", dph). Estos son los datos originales: nph <- c(2.50000000003638, 3.78947368414551, 1.51456310682008, 5.84686774940732,...