Preguntas etiquetadas con r

Aquí está mi contexto para esta pregunta: por lo que puedo decir, no podemos ejecutar una regresión de mínimos cuadrados ordinarios en R cuando usamos datos ponderados y el surveypaquete. Aquí, tenemos que usar svyglm(), que en su lugar ejecuta un modelo lineal generalizado (¿qué puede ser lo mismo? Estoy …

13 r  generalized-linear-model  least-squares  r-squared  deviance 

Mis predicciones provenientes de un modelo de regresión logística (glm en R) no están delimitadas entre 0 y 1 como esperaba. Comprendo la regresión logística es que los parámetros de entrada y modelo se combinan linealmente y la respuesta se transforma en una probabilidad utilizando la función de enlace logit. …

13 r  regression  logistic  generalized-linear-model 

Si construyo una matriz 2D compuesta completamente de datos aleatorios, esperaría que los componentes PCA y SVD esencialmente no explicaran nada. En cambio, parece que la primera columna SVD parece explicar el 75% de los datos. ¿Cómo puede ser esto posible? ¿Qué estoy haciendo mal? Aquí está la trama: Aquí …

13 r  pca  svd 

El bildpaquete parece ser un excelente paquete para respuestas binarias en serie. Pero es por tiempo discreto. Me gustaría especificar una función de tiempo suave para la conexión de odds ratio de la respuesta actual Y con respuestas binarias medidas en momentos anteriores, o al menos una versión de Markov …

13 r  repeated-measures  binary-data  panel-data 

Estoy buscando un paquete R general, limpio y rápido (es decir, usando rutinas C ++) para simular rutas desde una difusión no lineal no homogénea como (1) usando el esquema de Euler-Maruyama, el esquema de Milstein (o cualquier otro). Está destinado a integrarse en un código de estimación más grande …

13 r  simulation  stochastic-processes  markov-process 

Quiero predecir Tree Heights en un área determinada usando algunas variables obtenidas a través de la teledetección. Como la biomasa aproximada, etc. Quiero usar primero una regresión lineal (sé que no es la mejor idea, pero es un paso obligado para mi proyecto). Quería saber qué tan mal puede afectarlo …

13 r  multiple-regression  spatial  autocorrelation 

"Machine Learning For Hackers" de O'Reilly dice que cada componente principal representa un porcentaje de la variación. He citado la parte relevante de la página a continuación (capítulo 8, p. 207). Hablando con otro experto, acordaron que es el porcentaje. Sin embargo, los 24 componentes suman 133.2095%. ¿Como puede ser? …

13 r  pca 

Estoy tratando de ejecutar un logit bayesiano en los datos aquí . Estoy usando bayesglm()el armpaquete en R. La codificación es bastante sencilla: df = read.csv("http://dl.dropbox.com/u/1791181/bayesglm.csv", header=T) library(arm) model = bayesglm(PASS ~ SEX + HIGH, family=binomial(link="logit"), data=df) summary(model) da el siguiente resultado: Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) …

13 r  bayesian  p-value 

Estas dos funciones existen en R pero no sé sus diferencias. Parece que solo devuelven los mismos valores p cuando se llama wilcox.testcon correct=FALSE, y wilcox_test(en el paquete de monedas) con distribution="aymptotic". Para otros valores, devuelven valores p diferentes. También wilcox.testsiempre devuelve W = 0 para mi conjunto de datos, …

13 r  hypothesis-testing  p-value  wilcoxon-mann-whitney 

Me encantaría realizar una prueba post-hoc de TukeyHSD después de mi Anova bidireccional con R, obteniendo una tabla que contenga los pares ordenados agrupados por diferencia significativa. (Perdón por la redacción, todavía soy nuevo con las estadísticas). Me gustaría tener algo como esto: Entonces, agrupados con estrellas o letras. ¿Alguna …

13 r  anova  multiple-comparisons  post-hoc  tukey-hsd 

Estoy tratando de modelar algunos datos usando el glmnetpaquete en R. Digamos que tengo los siguientes datos training_x <- data.frame(variable1 = c(1, 2, 3, 2, 3), variable2 = c(1, 2, 3, 4, 5)) y <- c(1, 2, 3, 4, 5) (Esto es una simplificación; mis datos son mucho más complicados). …

13 r  glmnet 

Estoy tratando de hacer predicciones usando un modelo de bosque aleatorio en R. Sin embargo, obtengo errores ya que algunos factores tienen valores diferentes en el conjunto de prueba que en el conjunto de entrenamiento. Por ejemplo, un factor Cat_2tiene valores 34, 68, 76, etc., en el conjunto de prueba …

13 r  categorical-data  random-forest 

En respuesta a una pregunta sobre la selección de modelos en presencia de multicolinealidad , Frank Harrell sugirió : Ponga todas las variables en el modelo pero no pruebe el efecto de una variable ajustada por los efectos de las variables competidoras ... Las pruebas fragmentarias de las variables competidoras …

13 r  model-selection  multicollinearity 

Me preguntaba cómo generaría datos a partir de una ecuación de regresión de Poisson en R? Estoy un poco confundido sobre cómo abordar el problema. Entonces, si supongo que tenemos dos predictores y X 2 que se distribuyen N ( 0 , 1 ) . Y la intersección es 0 …

13 r  regression  poisson-distribution  simulation 

Esta también es una pregunta que se relaciona mucho con Python como un banco de trabajo de estadísticas y sobresale como un banco de trabajo de estadísticas . Sé que hay una gran discusión sobre Ruby versus Python, pero este no es el punto en esta pregunta. Pensé que Ruby …

13 r  python  software  ruby