Estoy implementando un clasificador SVM no lineal con kernel RBF. Me dijeron que la única diferencia con respecto a un SVM normal era que simplemente tenía que reemplazar el producto de puntos con una función del núcleo: Sé cómo funciona un SVM lineal normal, es decir, después de resolver el...
Tengo dos grupos de 10 participantes que fueron evaluados tres veces durante un experimento. Para probar las diferencias entre los grupos y entre las tres evaluaciones, ejecuté un ANOVA de diseño mixto 2x3 con group(control, experimental), time(primero, segundo, tres) y group x time. Ambos timey...
Logré crear una red neuronal de mis datos. Pero no estoy tan seguro sobre la interpretación de la salida R. Usé el siguiente comando para crear una red neuronal: > net=nnet(formula = category~iplen+date_time, size=0,skip=T,lineout=T) # weights: 3 initial value 136242.000000 final value...
He medido dos parámetros (carbono orgánico disuelto DOC = y, y descarga = x). Cuando estas dos variables se trazan una contra la otra, obtenemos un ciclo de histéresis (vea el ejemplo de código y la imagen). Ahora, para un análisis más detallado, me gustaría determinar el área de este circuito...
Quiero calcular el parámetro de la distribución exponencial partir de una población de muestra sacada de esta distribución en condiciones sesgadas. Hasta donde sé, para una muestra de n valores, el estimador habitual es . Sin embargo, mi muestra está sesgada de la siguiente manera:λλ\lambdami- λ...
Mi conjunto de datos tiene las siguientes variables: Tratamiento (4 tipos- fijo) Ubicación (8 ubicaciones- fijo) Posición en la ubicación (3 posiciones por ubicación, fija) Se toman muestras en cada posición (3 muestras por posición aleatoria) Tiempo (dos tiempos de muestreo - fijo) Tasa de...
He estado investigando el paquete de ratones , y aún no he descubierto una forma de usar las imputaciones múltiples para hacer un modelo de Cox, luego validar ese modelo con la función del paquete rmsvalidate() . Aquí hay un código de muestra de lo que tengo hasta ahora, usando el conjunto de datos...
Como seguimiento a la sugerencia del Modelo para una regresión de Cox con covariables dependientes del tiempo, aquí está la gráfica de Kaplan Meier que explica la naturaleza dependiente del tiempo de los embarazos. En otras palabras, el conjunto de datos ahora se divide en un conjunto de datos...
¿Alguien tiene un breve ejemplo educativo breve sobre cómo usar las redes neuronales ( nneten R por ejemplo) con el propósito de predicción? Aquí hay un ejemplo, en R, de una serie temporal T <- seq(0,20,length=200) Y <- 1 + 3*cos(4*T+2) +.2*T^2 + rnorm(200) plot(T,Y,type="l") Esto es solo...
Actualmente estoy usando R para predecir una serie temporal con estas instrucciones: X <- ts(datas, frequency=24) X.arima <- Arima(X, order=c(2,1,0), seasonal=c(1,1,1)) pred <- predict(X.arima, n.ahead=24) plot.ts(pred$pred) Como puede ver, tengo datos cada hora, y elegí el período...
Estoy tratando de realizar un análisis de crecimiento de clase latente (LCGA) y / o modelos de mezcla de crecimiento (GMM) en R. Los datos que estoy usando es un número cada vez mayor de tenedores de repositorios git (variable discreta, no categórica), como puedes ver en este conjunto de datos...
Los datos consisten en espectros ópticos (intensidad de luz contra frecuencia) registrados en diferentes tiempos. Los puntos se adquirieron en una cuadrícula regular en x (tiempo), y (frecuencia). Para analizar la evolución del tiempo a frecuencias específicas (un aumento rápido, seguido de una...
Tengo problemas para encontrar una solución con respecto a cómo ejecutar una prueba post-hoc (Tukey HSD) después de un ANOVA de medidas repetidas de 2 factores (ambos dentro de los sujetos) en R. Para el ANOVA, he usado la función aov: summary(aov(dv ~ x1 * x2 + Error(subject/(x1*x2)),...
Cuando se trata de datos con factores, R puede usarse para calcular las medias para cada grupo con la función lm (). Esto también proporciona los errores estándar para las medias estimadas. Pero este error estándar difiere de lo que obtengo de un cálculo a mano. Aquí hay un ejemplo (tomado de aquí...
Estoy usando el lme4paquete en R para hacer algunos modelos logísticos de efectos mixtos. Entiendo que la suma de cada efecto aleatorio debe ser cero. Cuando hago modelos mixtos lineales de juguete usando lmer, los efectos aleatorios generalmente son < confirmando mi creencia de que el But en...
Tengo datos que se parecen un poco a esto: ID Status 01 A 02 G 03 E ... ... 100 G Entiendes la idea, creo. Tengo estos datos de dos poblaciones separadas (cohortes) y quiero comparar la distribución de la variable de estado en una población con la distribución en otra. La pregunta que estoy...
Tengo una pregunta sobre cómo un estadístico interpretaría normalmente una salida anova. Digamos que tengo salida anova de R. > summary(fitted_data) Call: lm(formula = V1 ~ V2) Residuals: Min 1Q Median 3Q Max -2.74004 -0.33827 0.04062 0.44064 1.22737 Coefficients: Estimate Std. Error t...
Intentando replicar los resultados del artículo publicado recientemente, Aghion, Philippe, John Van Reenen y Luigi Zingales. 2013. "Innovación y propiedad institucional". American Economic Review, 103 (1): 277-304. (Los datos y el código de estado están disponibles en
Sabemos que si un estimador es un estimador imparcial de theta y si su varianza tiende a 0 cuando n tiende al infinito, entonces es un estimador consistente para theta. Pero esta es una condición suficiente y no necesaria. Estoy buscando un ejemplo de un estimador que sea consistente pero cuya...
Estoy tratando de ajustar un modelo logarítmico lineal a una gran cantidad de variables de los datos de la encuesta. Hay algunas razones por las que podría ser preferible ajustar las regresiones logísticas a esos datos. Varias autoridades sugieren que estos son equivalentes. Sin embargo, tengo...