Una regla de puntuación adecuada es una regla que se maximiza mediante un modelo 'verdadero' y no permite 'cubrir' o jugar con el sistema (informar deliberadamente diferentes resultados, como es la verdadera creencia del modelo para mejorar la puntuación). El puntaje de Brier es correcto, la...
Digamos que tenemos un modelo mod <- Y ~ X*Condition + (X*Condition|subject) # Y = logit variable # X = continuous variable # Condition = values A and B, dummy coded; the design is repeated # so all participants go through both Conditions # subject = random effects for different subjects...
Según tengo entendido, las escuelas del Reino Unido enseñan que la desviación estándar se encuentra usando: mientras que las escuelas de los Estados Unidos enseñan: (en un nivel básico de todos modos). Esto ha causado problemas a varios de mis alumnos en el pasado, ya que buscaron en...
Me toma la idea de la contracción de James-Stein (es decir, que una función no lineal de una sola observación de un vector de normales posiblemente independientes puede ser un mejor estimador de las medias de las variables aleatorias, donde 'mejor' se mide por error al cuadrado ) Sin embargo, nunca...
Estoy trabajando en un proyecto que involucra 14 variables y 345,000 observaciones para datos de vivienda (cosas como año de construcción, pies cuadrados, precio de venta, condado de residencia, etc.). Me preocupa tratar de encontrar buenas técnicas gráficas y bibliotecas R que contengan buenas...
Tengo un experimento que intentaré resumir aquí. Imagina que arrojo tres piedras blancas frente a ti y te pido que juzgues su posición. Grabo una variedad de propiedades de las piedras y su respuesta. Hago esto en varias materias. Genero dos modelos. Una es que la piedra más cercana a ti predice tu...
Al programar en R, he usado el paquete multinúcleo varias veces. Sin embargo, nunca he visto una declaración sobre cómo maneja sus números aleatorios. Cuando uso openMP con C, tengo cuidado de usar un RNG paralelo apropiado, pero con R supongo que sucede algo sensato. ¿Alguien puede confirmar que...
Cerrado. Esta pregunta está fuera de tema . Actualmente no está aceptando respuestas. Cerrado hace 4 años . Bloqueado . Esta pregunta y sus respuestas están bloqueadas porque la pregunta está fuera de tema pero tiene un significado histórico. Actualmente no...
Dadas dos matrices x e y, ambas de longitud n, calzo un modelo y = a + b * x y quiero calcular un intervalo de confianza del 95% para la pendiente. Esto es (b - delta, b + delta) donde b se encuentra de la manera habitual y delta = qt(0.975,df=n-2)*se.slope y se.slope es el error estándar en la...
Hice un modelo de regresión logística usando glm en R. Tengo dos variables independientes. ¿Cómo puedo trazar el límite de decisión de mi modelo en el diagrama de dispersión de las dos variables? Por ejemplo, ¿cómo puedo trazar una figura como:
Actualmente estoy recopilando datos para un experimento sobre las características psicosociales asociadas con la experiencia del dolor. Como parte de esto, estoy recopilando mediciones de GSR y BP electrónicamente de mis participantes, junto con varias medidas implícitas y de autoinforme. Tengo...
Estoy usando GNU R en una PC Ubuntu-Lucid que tiene 4 CPU. Para utilizar las 4 CPU, instalé el paquete "r-cran-multicore". Como el manual del paquete carece de ejemplos prácticos que entiendo, necesito consejos sobre cómo optimizar mi script para utilizar las 4 CPU. Mi conjunto de datos es un...
Quiero usar en la distribución para modelar retornos de activos de intervalo corto en un modelo bayesiano. Me gustaría estimar tanto los grados de libertad (junto con otros parámetros en mi modelo) para la distribución. Sé que los rendimientos de los activos son bastante no normales, pero no sé...
¿Cómo obtendrías los medios por hora para múltiples columnas de datos, durante un período diario, y mostrarías los resultados de doce "Hosts" en el mismo gráfico? Es decir, me gustaría graficar cómo se ve un período de 24 horas, para una semana de datos. El objetivo final sería comparar dos...
Jugué con algunas pruebas de raíz unitaria en R y no estoy completamente seguro de qué hacer con el parámetro k lag. Utilicé la prueba aumentada de Dickey Fuller y la prueba de Philipps Perron del paquete tseries . Obviamente, el parámetro predeterminado (para ) depende solo de la longitud de 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 2 años . ¿Podría alguien idear el código R para trazar...
Me interesaría encontrar formas en R para actualizar eficientemente un modelo lineal cuando se agrega una observación o un predictor. biglm tiene una capacidad de actualización al agregar observaciones, pero mis datos son lo suficientemente pequeños como para residir en la memoria (aunque tengo una...
Digamos que tengo un simple experimento factorial 2x2 en el que quiero hacer ANOVA. Así, por ejemplo: d <- data.frame(a=factor(sample(c('a1','a2'), 100, rep=T)), b=factor(sample(c('b1','b2'), 100, rep=T))); d$y <- as.numeric(d$a)*rnorm(100, mean=.75, sd=1) + as.numeric(d$b)*rnorm(100,...
Cuando se realiza un GLM y se obtiene el error "no definido debido a singularidades" en la salida de anova, ¿cómo se contrarresta este error? Algunos han sugerido que se debe a la colinealidad entre covariables o que uno de los niveles no está presente en el conjunto de datos (ver: interpretación...
Estoy tratando de usar la función ' densidad ' en R para hacer estimaciones de densidad del núcleo. Tengo algunas dificultades para interpretar los resultados y comparar varios conjuntos de datos, ya que parece que el área bajo la curva no es necesariamente 1. Para cualquier función de densidad de...