Preguntas etiquetadas con r

¿Hay alguna manera en R (una función incorporada) para calcular la matriz de transición para una cadena de Markov a partir de un conjunto de observaciones? Por ejemplo, ¿tomar un conjunto de datos como el siguiente y calcular la matriz de transición de primer orden? dat<-data.frame(replicate(20,sample(c("A", "B", "C","D"), size = …

29 r  markov-process 

Utilizo la función auto.arima () en el paquete de pronóstico para ajustar los modelos ARMAX con una variedad de covariables. Sin embargo, a menudo tengo una gran cantidad de variables para seleccionar y generalmente termino con un modelo final que funciona con un subconjunto de ellas. No me gustan las …

29 r  time-series  lasso  regularization  elastic-net 

Soy bastante nuevo en estadísticas y necesito tu ayuda. Tengo una pequeña muestra, como sigue: H4U 0.269 0.357 0.2 0.221 0.275 0.277 0.253 0.127 0.246 Ejecuté la prueba de Shapiro-Wilk usando R: shapiro.test(precisionH4U$H4U) y obtuve el siguiente resultado: W = 0.9502, p-value = 0.6921 Ahora, si supongo que el nivel …

29 r  distributions  interpretation  goodness-of-fit  normality-assumption 

Siempre uso lm()en R para realizar una regresión lineal de en . Esa función devuelve un coeficiente tal queyyyXxxββ\betay= βx .y=βx.y = \beta x. Hoy aprendí sobre mínimos cuadrados totales y esa princomp()función (análisis de componentes principales, PCA) se puede utilizar para realizarla. Debería ser bueno para mí (más preciso). …

29 r  pca  least-squares  deming-regression  total-least-squares 

Podemos usar lm()para predecir un valor, pero aún necesitamos la ecuación de la fórmula del resultado en algunos casos. Por ejemplo, agregue la ecuación a las gráficas.

29 r  regression  lm 

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 2 años . Estoy usando caret para ejecutar un bosque aleatorio validado cruzado sobre un conjunto de datos. La …

29 r  random-forest  caret  regression  prediction  fitting  social-science  poisson-distribution  distributions  characteristic-function  bayesian  prior  regression  normal-distribution  interaction  nonparametric  skewness  svm  standard-deviation  standard-error  regression-coefficients  igraph  natural-language  word2vec  word-embeddings  regression  machine-learning  sampling  r  regression  machine-learning  random-forest  ensemble  sampling  unbiased-estimator  proof  estimators  mse  probability  conditional-probability  bayes  anova  missing-data  neural-networks  recommender-system  r  confidence-interval  sample  multiple-imputation  r  time-series  forecasting  mase 

Estoy usando el paquete robustbase para ejecutar una estimación de glm. Sin embargo, cuando lo hago, aparece el siguiente error: Error in solve.default(crossprod(X, DiagB * X)/nobs, EEq) : system is computationally singular: reciprocal condition number = 1.66807e-16 ¿Qué significa / indica esto? ¿Y cómo puedo depurarlo? PD. Si necesita algo …

29 r  generalized-linear-model  robust 

Estoy usando R (3.1.1) y modelos ARIMA para pronosticar. Me gustaría saber cuál debería ser el parámetro de "frecuencia", que se asigna en la ts()función , si estoy usando datos de series temporales que son: separado por minutos y se extiende por 180 días (1440 minutos / día) separado por …

28 r  time-series  arima  multiple-seasonalities  mape 

Estoy usando glmnet para calcular las estimaciones de regresión de crestas. Obtuve algunos resultados que me hicieron sospechar que glmnet realmente está haciendo lo que creo que hace. Para verificar esto, escribí un script R simple donde comparo el resultado de la regresión de cresta realizada por resolver y el …

28 r  ridge-regression  glmnet 

Me gustaría obtener un valor p y un tamaño de efecto de una variable categórica independiente (con varios niveles), que es "general" y no para cada nivel por separado, como es la salida normal de lme4R. Es igual que lo que la gente informa cuando ejecuta un ANOVA. ¿Cómo puedo …

28 r  hypothesis-testing  anova  mixed-model  lme4-nlme 

[El título inicial "Medición de similitud para árboles de agrupamiento jerárquico" fue cambiado más tarde por @ttnphns para reflejar mejor el tema] Estoy realizando una serie de análisis de conglomerados jerárquicos en un marco de datos de registros de pacientes (por ejemplo, similar a http://www.biomedcentral.com/1471-2105/5/126/figure/F1?highres=y ) Estoy experimentando con diferentes …

28 r  clustering  distance-functions  similarities  dendrogram 

Entiendo que esta pregunta es bastante amplia, pero me pregunto cuáles deberían ser los puntos decisivos para decidir crear (o no) un nuevo paquete para R. Para ser más específico, agregaría que la pregunta no se trata de las razones para use R en sí mismo, más sobre la decisión …

28 r  software 

Espero que no les importe esta pregunta, pero necesito ayuda para interpretar el resultado de un resultado de modelo de efectos mixtos lineales que he tratado de aprender a hacer en R. Soy nuevo en el análisis de datos longitudinales y la regresión de efectos mixtos lineales. Tengo un modelo …

28 r  mixed-model  interpretation  panel-data 

Como ejemplo, considere el ChickWeightconjunto de datos en R. La varianza obviamente crece con el tiempo, así que si uso una regresión lineal simple como: m <- lm(weight ~ Time*Diet, data=ChickWeight) Mis preguntas: ¿Qué aspectos del modelo serán cuestionables? ¿Se limitan los problemas a extrapolar fuera del Timerango? ¿Qué tan …

28 r  regression  heteroscedasticity  assumptions 

Encontré una fórmula para pseudo en el libro Extending the Linear Model with R, Julian J. Faraway (p. 59).R2R2R^2 1 - Desviación residualDesviación nula1-Desviación residualDesviación nula1-\frac{\text{ResidualDeviance}}{\text{NullDeviance}} . ¿Es esta una fórmula común para pseudo para GLM?R2R2R^2

28 r  regression  generalized-linear-model  r-squared