¿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...
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...
Varias fuentes sugieren que hay muchas consecuencias negativas de la discretización (categorización) de variables continuas antes del análisis estadístico (muestra de referencias [1] - [4] a continuación). Por el contrario [5] sugiere que se sabe que algunas técnicas de aprendizaje automático...
El proceso de generación de datos es:y= sin ( x + I( d= 0 ) ) + sin ( x + 4 ∗ I( d= 1 ) ) + I( d= 0 )z2+ 3 I( d= 1 )z2+ N ( 0 , 1 )y=sin(x+I(d=0))+sin(x+4∗I(d=1))+I(d=0)z2+3I(d=1)z2+N(0,1)y = \text{sin}\Big(x+I(d=0)\Big) + \text{sin}\Big(x+4*I(d=1)\Big) + I(d=0)z^2 + 3I(d=1)z^2 +...
Primero simulemos algunos datos para una regresión logística con partes fijas y aleatorias: set.seed(1) n <- 100 x <- runif(n) z <- sample(c(0,1), n, replace=TRUE) b <- rnorm(2) beta <- c(0.4, 0.8) X <- model.matrix(~x) Z <- cbind(z, 1-z) eta <- X%*%beta + Z%*%b pr <-...
Algunos de ustedes podrían haber leído este bonito artículo: O'Hara RB, Kotze DJ (2010) No log-transforma los datos de recuento. Methods in Ecology and Evolution 1: 118–122. Klick . Actualmente estoy comparando modelos binomiales negativos con modelos gaussianos en datos transformados. A...
Estoy buscando una medida adecuada de la "varianza explicada" de un Poisson GLM (usando una función de enlace de registro). He encontrado una cantidad de recursos diferentes (tanto en este sitio como en otros lugares) que analizan una serie de diferentes pseudo-R2R2R^2 medidas, pero casi todos los...
Me gustaría especificar un modelo de regresión logística donde tenga la siguiente relación: mi[YyoEl |Xyo] = f( βXyo 1+β2Xyo 2)E[Yi|Xi]=f(βxi1+β2xi2)E[Y_i|X_i] = f(\beta x_{i1} + \beta^2x_{i2}) donde es la función de logit inversa.Fff ¿Existe una forma "rápida" de hacer esto con funciones R...
He estado usando la función glm.fit en R para ajustar parámetros a un modelo de regresión logística. Por defecto, glm.fit utiliza mínimos cuadrados repesados de forma iterativa para ajustar los parámetros. ¿Cuáles son algunas razones por las cuales este algoritmo no podría converger, cuando se...
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,...
Quiero ajustar un modelo lineal por R con family=binomial(link="identity"), sin embargo, la familia binomial no tiene un enlace de identidad. ¿Qué tengo que
Así que estoy jugando con la idea de escribir un algoritmo que crezca y elimine un árbol de regresión de los datos y luego, en los nodos terminales del árbol, se ajuste a un GLM. He estado tratando de leer sobre la idea, pero parece que no puedo encontrar ningún nombre consistente para la técnica....
Así que aquí estoy estudiando modelos lineales generalizados. Sé que esta pregunta es bastante ingenua y simple, pero no sé exactamente por qué la función canónica de enlace es tan útil. ¿Podría alguien proporcionarme una intuición sobre este
Enseño física a estudiantes de secundaria y me gustaría que mis estudiantes realicen una comparación rudimentaria del modelo bayesiano para obtener datos de sus experimentos. Descubrí una forma para que lo hagan (ver más abajo), pero no estoy seguro de que sea correcto. Agradecería mucho cualquier...
Cerrado. Esta pregunta está fuera de tema . Actualmente no está aceptando respuestas. ¿Quieres mejorar esta pregunta? Actualice la pregunta para que esté en el tema de Cross Validated. Cerrado el año pasado . Estoy tratando de ajustar el modelo lineal...
Primero déjame darte algunos antecedentes; Resumiré mis preguntas al final. La distribución Beta, parametrizada por su media yμμ\muϕϕ\phi, tiene Var( Y) = V( μ ) / ( ϕ + 1 )Var(Y)=V(μ)/ /(ϕ+1)\operatorname{Var}(Y) = \operatorname{V}(\mu)/(\phi+1), donde es la función de variación.V( μ ) = μ ( 1...