¿Alguien sabe de algún código bien escrito (en Matlab o R) para MCMC de salto reversible? Preferiblemente, una aplicación de demostración simple para complementar los documentos sobre el tema, que sería útil para comprender el
Se puede usar una base de datos de (población, área, forma) para mapear la densidad de población asignando un valor constante de población / área a cada forma (que es un polígono como un bloque del Censo, tramo, condado, estado, lo que sea). Sin embargo, las poblaciones generalmente no están...
He estado tratando de discernir qué significa exactamente la salida "coef" y "(exp) coef" de coxph. Parece que el "(exp) coef" son comparaciones de la primera variable en el modelo de acuerdo con el grupo asignado en el comando. ¿Cómo llega la función coxph a los valores de "coef" y "(exp)...
En R, el drop1comando genera algo ordenado. Estos dos comandos deberían obtener algún resultado: example(step)#-> swiss drop1(lm1, test="F") El mío se ve así: > drop1(lm1, test="F") Single term deletions Model: Fertility ~ Agriculture + Examination + Education + Catholic +...
Necesito agrupar unidades en grupos para minimizar la suma de cuadrados dentro del grupo (WSS), pero necesito asegurarme de que cada grupo contenga al menos unidades. ¿Alguna idea de si alguna de las funciones de agrupación de R permite la agrupación en agrupaciones sujetas a una restricción de...
Espero que esta no sea una pregunta tonta. Digamos que tengo una distribución continua arbitraria. También tengo una estadística, y me gustaría usar esta distribución arbitraria para obtener un valor p para esta estadística. Me doy cuenta de que en R es fácil hacer esto siempre que su distribución...
Estoy luchando por encontrar un método para reducir el número de categorías en datos nominales u ordinales. Por ejemplo, supongamos que quiero construir un modelo de regresión en un conjunto de datos que tenga varios factores nominales y ordinales. Si bien no tengo problemas con este paso, a...
Cerrado. Esta pregunta está fuera de tema . Actualmente no está aceptando respuestas. ¿Quieres mejorar esta pregunta? Actualice la pregunta para que sea sobre el tema de Cross Validated. Cerrado hace 2 años . A veces necesito obtener solo la primera fila...
Estoy a punto de sumergirme en el aprendizaje de R y mi proyecto de aprendizaje implicará aplicar una regresión de efectos mixtos o aleatorios a un conjunto de datos para desarrollar una ecuación predictiva. Comparto la preocupación del escritor en esta publicación ¿Cómo elegir la biblioteca nlme...
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 7 meses . Es fácil encontrar un área de cálculo de...
Estoy tratando de descubrir cómo controlar los parámetros de suavizado en un modelo mgcv: gam. Tengo una variable binomial que intento modelar principalmente como una función de las coordenadas xey en una cuadrícula fija, además de algunas otras variables con influencias menores. En el pasado, he...
Pasé tres días incursionando tmdespués de leer un borrador de un amigo donde exploró un corpus de texto con UCINET, mostrando nubes de texto, gráficos de red de dos modos y descomposición de valor único (con gráficos, usando Stata). Me encontré con una gran cantidad de problemas: en Mac OS X, hay...
Soy un novato en estadística y R y tengo problemas para usar la función Levene (me gustaría verificar la igualdad de varianza de dos muestras). La documentación dice que debería ejecutar: levene.test (y, grupo) ¿Pero no tengo idea de lo que debo poner como y y grupo? Tengo dos muestras...
¿Puede alguien decirme cómo hacer que R calcule el punto de ruptura en un modelo lineal por partes (como un parámetro fijo o aleatorio), cuando también necesito estimar otros efectos aleatorios? He incluido un ejemplo de juguete a continuación que se ajusta a una regresión de palo de hockey /...
¿Alguna recomendación para elegir una biblioteca de optimización restringida adecuada para mi función de optimización? Estoy minimizando ai) la función no lineal con restricciones lineales de igualdad y desigualdad, y ii) tengo disponible el gradiente y el hessian de la función. Si ayuda, la...
Estoy tratando de comprender la medida dinámica de deformación del tiempo para comparar series de tiempo juntas. Tengo tres series de datos de series de tiempo como esta: T1 <- structure(c(0.000213652387565, 0.000535045478866, 0, 0, 0.000219346347883, 0.000359669104424, 0.000269469145783,...
Soy un estudiante graduado en ciencias de la computación. He estado haciendo un análisis factorial exploratorio para un proyecto de investigación. Mis colegas (que lideran el proyecto) usan SPSS, mientras que yo prefiero usar R. Esto no importó hasta que descubrimos una gran discrepancia entre los...
Voy a ver los videos del curso gratuito de aprendizaje automático en línea de Andrew Ng en Stanford. Analiza el Descenso de degradado como un algoritmo para resolver regresión lineal y funciones de escritura en Octave para realizarlo. Presumiblemente podría reescribir esas funciones en R, pero mi...
Basado en una muestra de tiempos de supervivencia, me gustaría estimar la probabilidad de sobrevivir al tiempo , para algunos específicos , utilizando el estimador de Kaplan-Meier. ¿Es posible hacer esto en ? Tenga en cuenta que no es necesariamente un tiempo de
Considere el siguiente biplot de PCA: library(mvtnorm) set.seed(1) x <- rmvnorm(2000, rep(0, 6), diag(c(5, rep(1,5)))) x <- scale(x, center=T, scale=F) pc <- princomp(x) biplot(pc) Hay un montón de flechas rojas trazadas, ¿qué significan? Sabía que la primera flecha etiquetada con...