Preguntas etiquetadas con r

Estoy leyendo A. Agresti (2007), Introducción al análisis de datos categóricos , 2do. edición, y no estoy seguro si entiendo este párrafo (p.106, 4.2.1) correctamente (aunque debería ser fácil): En la Tabla 3.1 sobre ronquidos y enfermedades cardíacas en el capítulo anterior, 254 sujetos informaron ronquidos todas las noches, de …

9 r  generalized-linear-model  logistic 

El modelo Bradley – Terry – Luce (BTL) establece que , donde es la probabilidad de que el objeto se considere "mejor", más pesado, etc., que el objeto , y , y son parámetros.pji=logit−1(δj−δi)pji=logit−1(δj−δi)p_{ji} = logit^{-1}(\delta_j - \delta_i)pijpijp_{ij}jjjiiiδiδi\delta_iδjδj\delta_j Esto parece ser un candidato para la función glm, con family = …

9 r  logistic  categorical-data  paired-data  bradley-terry-model 

Mis datos son una serie temporal de población ocupada, L, y el lapso de tiempo, año. n.auto=auto.arima(log(L),xreg=year) summary(n.auto) Series: log(L) ARIMA(2,0,2) with non-zero mean Coefficients: ar1 ar2 ma1 ma2 intercept year 1.9122 -0.9567 -0.3082 0.0254 -3.5904 0.0074 s.e. NaN NaN NaN NaN 1.6058 0.0008 sigma^2 estimated as 1.503e-06: log likelihood=107.55 …

9 r  regression  arima 

Sé que esta es probablemente una pregunta básica ... Pero parece que no encuentro la respuesta. Estoy adaptando un GLM con una familia de Poisson, y luego intenté echar un vistazo a las predicciones, sin embargo, el desplazamiento parece tenerse en cuenta: model_glm=glm(cases~rhs(data$year,2003)+lhs(data$year,2003), offset=(log(population)), data=data, subset=28:36, family=poisson()) predict (model_glm, type="response") …

9 r  generalized-linear-model  prediction  offset 

@whuber ha demostrado cómo simular resultados multivariados ( y1y1y_1 , y2y2y_2 e ) para un punto de tiempo.y3y3y_3 Como sabemos, los datos longitudinales a menudo ocurren en estudios médicos. Mi pregunta es cómo simular resultados multivariados de medidas repetidas en R? Por ejemplo, medimos repetidamente , e en 5 puntos …

9 r  repeated-measures  simulation 

Necesito un paquete que me pueda dar la ecuación para un modelo SVM lineal. Actualmente estoy usando e1071 así: library(e1071) m = svm(data, labels, type='C', kernel='linear', cost=cost, probability=FALSE, scale=scale) w = t(m$coefs) %*% data[m$index,] #Weight vector b = -model$rho #Offset Sin embargo, no estoy seguro de cómo e1071::svm()selecciona las clases …

9 r  svm  e1071 

Tengo problemas al usar las funciones cor()y cor.test(). Solo tengo dos matrices (solo valores numéricos y el mismo número de filas y columnas) y quiero tener el número de correlación y el valor p correspondiente. Cuando uso cor(matrix1, matrix2)obtengo los coeficientes de correlación para todas las celdas. Solo quiero un …

9 r  correlation 

Estoy usando el paquete kernlab en R para construir un SVM para clasificar algunos datos. El SVM está funcionando bien ya que proporciona 'predicciones' de una precisión decente, sin embargo, mi lista de variables de entrada es más grande de lo que quisiera y no estoy seguro de la importancia …

9 r  machine-learning  svm  genetic-algorithms 

Digamos que tiene datos de supervivencia como este: obs <- data.frame( time = c(floor(runif(100) * 30), floor((runif(100)^2) * 30)), status = c(rbinom(100, 1, 0.2), rbinom(100, 1, 0.7)), group = gl(2,100) ) Para realizar una prueba de rango de registro estándar, se puede usar survdiff(Surv(time, status) ~ group, data = obs, …

9 r  survival  wilcoxon-signed-rank 

Como pasatiempo secundario, he estado explorando las series temporales de pronósticos (en particular, usando R). Para mis datos, tengo el número de visitas por día, por cada día que se remonta a casi 4 años. En estos datos hay algunos patrones distintos: De lunes a viernes tiene muchas visitas (mayor …

9 r  time-series  seasonality  forecasting 

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 hace 3 años . Estoy muy interesado en aprender a crear informes recurrentes a partir de mi código R y …

9 r 

Estoy buscando un ejemplo de código simple de cómo ejecutar un filtro de partículas en R. El paquete pomp parece admitir el bit matemático del espacio de estado, pero los ejemplos son un poco difíciles de seguir programáticamente para un desarrollador OO simple como yo, particularmente cómo cargar los datos …

9 r 

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 hace 8 meses . El problema mencionado en esta pregunta se solucionó en la versión 1.7.3 del paquete R glmnet. …

9 r  multinomial  glmnet 

He creado mi propia versión ligeramente mejorada del termplot que uso en este ejemplo, puedes encontrarla aquí . He publicado anteriormente en SO, pero cuanto más lo pienso, creo que esto probablemente esté más relacionado con la interpretación del modelo de riesgos proporcionales de Cox que con la codificación real. …

9 r  survival  cox-model 

Tengo algunos datos de series temporales en los que la variable medida es enteros positivos discretos (recuentos). Quiero probar si hay una tendencia al alza con el tiempo (o no). La variable independiente (x) está en el rango 0-500 y la variable dependiente (y) está en el rango 0-8. Pensé …

9 r  regression  python