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...
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...
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...
¿Existe un principio general sobre si uno debe calcular la correlación de Pearson para dos variables aleatorias X e Y antes de tomar su transformación logarítmica o después? ¿Existe un procedimiento para probar cuál es más apropiado? Producen valores similares pero diferentes, ya que la...
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...
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...
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...
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...
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...
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...
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
Tenemos un conjunto de datos con dos covariables y una variable de agrupación categórica y queremos saber si existen diferencias significativas entre la pendiente o la intercepción entre las covariables asociadas con las diferentes variables de agrupación. Hemos utilizado anova () y lm () para...
@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...
He estado usando el Root Mean Squared Error(RMSE) para medir la precisión de los valores pronosticados usando un modelo. Entiendo que el valor devuelto está usando las unidades de mis medidas (en lugar de un porcentaje). Sin embargo, me gustaría citar mis valores como un porcentaje. El enfoque que...
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...
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...
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...
Me preguntaba si alguien tenía experiencia en el uso de la función de ratones, como se describe en ratones: ¿Imputación multivariada por ecuaciones encadenadas en R (JSS 2011 45 (3))? Tengo un conjunto de datos con varias variables, cada una con diferentes grados de datos faltantes. Mi pregunta...
Estoy usando el paquete randomForest en R (R versión 2.13.1, randomForest versión 4.6-2) para la regresión y noté un sesgo significativo en mis resultados: el error de predicción depende del valor de la variable de respuesta. Los valores altos se predicen por debajo y los valores bajos se predicen...
Los siguientes injertos se toman de este artículo . Soy novato en bootstrap e intento implementar el bootstrapping paramétrico, semiparamétrico y no paramétrico para el modelo mixto lineal con R bootpaquete. Código R Aquí está mi