Soy muy nuevo en mínimos cuadrados parciales (PLS) y trato de entender la salida de la función R plsr()en el plspaquete. Simulemos datos y ejecutemos el PLS: library(pls) n <- 50 x1 <- rnorm(n); xx1 <- scale(x1) x2 <- rnorm(n); xx2 <- scale(x2) y <- x1 + x2 + rnorm(n,0,0.1); yy...
Estoy ejecutando la siguiente prueba de raíz unitaria (Dickey-Fuller) en una serie de tiempo usando la ur.df()función en el urcapaquete. El comando es: summary(ur.df(d.Aus, type = "drift", 6)) El resultado es: ############################################### # Augmented Dickey-Fuller Test Unit...
En R, si configuro.seed () y luego uso la función de muestra para aleatorizar una lista, ¿puedo garantizar que no generaré la misma permutación? es decir... set.seed(25) limit <- 3 myindex <- seq(0,limit) for (x in seq(1,factorial(limit))) { permutations <- sample(myindex)...
Tengo una muestra de datos generados Rpor rnorm(50,0,1), por lo que los datos obviamente adquieren una distribución normal. Sin embargo, Rno "conoce" esta información de distribución sobre los datos. ¿Existe algún método Rque pueda estimar de qué tipo de distribución proviene mi muestra? Si no,...
Me han hecho la siguiente pregunta como parte de la tarea: Diseñe e implemente un estudio de simulación para examinar el rendimiento del bootstrap para obtener intervalos de confianza del 95% en la media de una muestra univariada de datos. Su implementación puede ser en R o SAS. Los aspectos...
Tengo un montón de conjuntos de datos relacionados. Las correlaciones de Pearson entre pares de ellas suelen ser definitivamente mayores que las correlaciones de Spearman. Eso sugiere que cualquier correlación es lineal, pero uno podría esperar que incluso si Pearson y Spearman fueran iguales. ¿Qué...
Estoy tratando de analizar datos de medidas repetidas y estoy luchando para que funcione R. Mis datos son esencialmente los siguientes, tengo dos grupos de tratamiento. Todos los sujetos de cada grupo se evalúan todos los días y se les asigna una puntuación (el porcentaje correcto en una prueba)....
¿Por qué incluso tener antecedentes no informativos? No proporcionan información sobre . Entonces, ¿por qué usarlos? ¿Por qué no solo usar previos informativos? Por ejemplo, supongamos que θ ∈ [ 0 , 1 ] . Entonces es θ ~ T ( 0 , 1 ) una previa no informativa para θ ?θθ\thetaθ∈[0,1]θ∈[0,1] \theta...
En una publicación anterior, me preguntaba cómo lidiar con los puntajes EQ-5D . Recientemente me topé con la regresión logística de cuantiles sugerida por Bottai y McKeown que introduce una forma elegante de lidiar con resultados acotados. La fórmula es simple: l o gi t ( y) = l o g( y- ym i nym a...
Nota: esta pregunta es una nueva publicación, ya que mi pregunta anterior tuvo que ser eliminada por razones legales. Al comparar PROC MIXED de SAS con la función lmedel nlmepaquete en R, me topé con algunas diferencias bastante confusas. Más específicamente, los grados de libertad en las...
He usado rjags para ejecutar MCMC en un modelo, especificado en el lenguaje JAGS. ¿Hay una buena manera de extraer ese modelo y realizar predicciones con él (usando las distribuciones posteriores de mis parámetros)? Puedo volver a especificar el modelo en R y conectar los modos de mis parámetros...
Dado un conjunto de datos: x <-
Estoy tratando de probar varios enfoques de análisis de datos funcionales. Idealmente, me gustaría probar el panel de enfoques que tengo sobre datos funcionales simulados. Intenté generar FD simulada utilizando un enfoque basado en un ruido gaussiano sumador (código a continuación), pero las curvas...
Tengo problemas de rendimiento al usar el MCMCglmmpaquete en R para ejecutar un modelo de efectos mixtos. El código se ve así: MC1<-MCMCglmm(bull~1,random=~school,data=dt,family="categorical" , prior=list(R=list(V=1,fix=1), G=list(G1=list(V=1, nu=0))) , slice=T, nitt=iter, ,burnin=burn,...
Tengo la siguiente función de probabilidad: Prob=11+e−zProb=11+e−z\text{Prob} = \frac{1}{1 + e^{-z}} dónde z=B0+B1X1+⋯+BnXn.z=B0+B1X1+⋯+BnXn.z = B_0 + B_1X_1 + \dots + B_nX_n. Mi modelo se parece Pr(Y=1)=11+exp(−[−3.92+0.014×(bid)])Pr(Y=1)=11+exp(−[−3.92+0.014×(bid)])\Pr(Y=1) = \frac{1}{1 +...
¿Es mejor diferenciar una serie (suponiendo que la necesite) antes de usar un Arima O mejor usar el parámetro d dentro de Arima? Me sorprendió cuán diferentes son los valores ajustados dependiendo de qué ruta se tome con el mismo modelo y datos. ¿O estoy haciendo algo
Los investigadores a menudo usan dos medidas que tienen elementos muy similares y argumentan que miden cosas diferentes (por ejemplo, "Siempre me preocupo cuando estoy cerca de autos"; "Tengo miedo de los autos"). Llamemos a las medidas hipotéticas Medida y ansiedad del miedo a los automóviles de...
Digamos que tengo los siguientes datos y estoy ejecutando un modelo de regresión: df=data.frame(income=c(5,3,47,8,6,5), won=c(0,0,1,1,1,0), age=c(18,18,23,50,19,39), home=c(0,0,1,0,0,1)) Por un lado, ejecuto un modelo lineal para predecir los ingresos: md1 = lm(income ~ age + home + home,...
Estoy usando el paquete randomForest en R y usando los datos del iris, el bosque aleatorio generado es una clasificación, pero cuando uso un conjunto de datos con alrededor de 700 características (las características son cada píxel en una imagen de 28x28 píxeles) y la columna de etiqueta se llama...
Estoy tratando de descomponer una matriz de covarianza basada en un conjunto de datos dispersos / breves. Me doy cuenta de que la suma de lambda (varianza explicada), según se calcula con svd, se amplifica con datos cada vez más vacíos. Sin huecos, svdy eigenarrojar los mismos resultados. Esto no...