La mayoría de las distribuciones estándar en R tienen una familia de comandos: pdf / pmf, cdf / cmf, cuantil, desviaciones aleatorias (por ejemplo, dnorm, pnorm, qnorm, rnorm). Sé que es bastante fácil utilizar algunos comandos estándar para reproducir estas funciones para las distribuciones...
Es útil estudiar el código de análisis de datos de expertos. Recientemente he estado examinando github y hay varias personas que comparten código de análisis de datos allí. Esto incluye algunos paquetes R (que por supuesto están disponibles directamente de CRAN), pero también varios ejemplos de...
Utilizo el paquete "boot" para calcular un valor p de arranque de 2 lados aproximado , pero el resultado está demasiado lejos del valor p de usar t.test. No puedo entender qué hice mal en mi código R. ¿Puede alguien darme una pista para esto? time = c(14,18,11,13,18,17,21,9,16,17,14,15,...
¿Hay alguna medida de similitud o distancia entre dos matrices de covarianza simétrica (ambas tienen las mismas dimensiones)? Estoy pensando aquí en análogos a la divergencia KL de dos distribuciones de probabilidad o la distancia euclidiana entre vectores, excepto que se aplica a las matrices. Me...
[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...
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...
Mira esta imagen: Si extraemos una muestra de la densidad roja, se espera que algunos valores sean inferiores a 0.25, mientras que es imposible generar dicha muestra a partir de la distribución azul. Como consecuencia, la distancia Kullback-Leibler de la densidad roja a la densidad azul es...
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...
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 conseguir...
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...
A menudo leo sobre una función que es 'altamente no lineal'. En mi opinión, hay "lineal" y "no lineal", entonces, ¿de qué se trata esto 'altamente'? ¿Hay una diferencia formal de no lineal? ¿Cómo se
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 . Tengo un diagrama de dispersión. ¿Cómo puedo...
Quizás esta pregunta sea ingenua, pero: Si la regresión lineal está estrechamente relacionada con el coeficiente de correlación de Pearson, ¿hay alguna técnica de regresión estrechamente relacionada con los coeficientes de correlación de Kendall y
Bloqueado . Esta pregunta y sus respuestas están bloqueadas porque la pregunta está fuera de tema pero tiene un significado histórico. Actualmente no acepta nuevas respuestas o interacciones. Me preguntaba si es posible hacer cálculos simbólicos en R? Por ejemplo,...
Me sorprende que esto no se haya preguntado antes, pero no puedo encontrar la pregunta en stats.stackexchange. Esta es la fórmula para calcular la varianza de una muestra distribuida normalmente: ∑(X−X¯)2n−1∑(X−X¯)2n−1\frac{\sum(X - \bar{X}) ^2}{n-1} Esta es la fórmula para calcular el error...
Me pregunto si hay una manera simple de producir una lista de variables usando un bucle for y dar su valor. for(i in 1:3) { noquote(paste("a",i,sep=""))=i } En el código anterior, trato de crear a1, a2, a3, que asignan a los valores de 1, 2, 3. Sin embargo, R da un mensaje de error. Gracias por...
Esta es una pregunta sobre terminología. ¿Es un "vago previo" lo mismo que un previo no informativo, o hay alguna diferencia entre los dos? Mi impresión es que son iguales (al buscar juntos vagos y no informativos), pero no puedo estar
Considere un Jeffreys anterior donde , donde es la información de Fisher.yop ( θ ) ∝ | i ( θ ) |----√p(θ)∝|i(θ)|p(\theta) \propto \sqrt{|i(\theta)|}yoii Sigo viendo esto antes mencionado como un previo no informativo, pero nunca vi un argumento por qué no es informativo. Después de todo, no es un...
Reflexiono sobre la discusión en torno a esta pregunta y, en particular, sobre el comentario de Frank Harrell de que la estimación de la varianza en un modelo reducido (es decir, uno a partir del cual se probaron y rechazaron varias variables explicativas) debería usar los grados de libertad...
Necesito usar variables binarias (valores 0 y 1) en k-means. Pero k-means solo funciona con variables continuas. Sé que algunas personas todavía usan estas variables binarias en k-means ignorando el hecho de que k-means solo está diseñado para variables continuas. Esto es inaceptable para...