search for: observado

Buenos días a todxs, Estoy comparando la predicción de los valores (0, 1, 2, 3,.....hasta 13) frente a los observados. Con la idea de comparar el modelo Zero inflated y el Binomial negativo y ver cual presenta mas distancia frente a las predicciones observadas. Para ello introduzco los códigos en la consola: #Modelo ZIM pred<-round(colSums(predict(zeroinfl, type="prob") [,1:14])) #Valores observado...

Hola Javier, Si, cuando hablo de valor observado me refiero al valor real en campo y el predicho al que estiman los modelos. Disculpa, que no lo detallase así desde el principio. En mi caso trabajo con dos diferentes: Zero inflated y Binomial Negativo y me gustaría comprobar que diferencia (distancia) existe entre cada uno de ellos y la realidad...

Buenos días; Estoy intentando aprender algo de sobre los modelos lineales mixtos. Estoy siguiendo el libro "Multilevlel analysis: techniques and applications" de Joop Hox. Estoy intentado reproducir las tablas y gráficos que trae el libro con el paquete lmer. En este libro se indica que en la comprobación de los supuestos del modelo hay que estudiar los gráficos de residuales tanto

...de a loop I have to do it automatically. The main problem is that inside the loop I work with differente types of variables, I mean , the range is different for each one. I have already tried : par(xpd=TRUE) pos<-par("usr") legend(pos[1],pos[3]-(pos[3]*.05),c("Modelo","Observado"),col = c("blue","orange"),cex=1.0, pch=15) I would like to be able to plot the legend outer the graphics region. Something like I've done with title. You can see one figure an attached file. I really appreciate any help. Nilza Barros -- Abra?o, Nilza Barros

...leta, cex.lab=0.7, cex.axis=0.7, cex.main=0.7, xlab="Coord X", ylab="Coord Y", main=paste("Espécie: ",splist[i]), xlim=c(210000,240000) ) col_lab <- c(max(quad_N_sp$x),min(quad_N_sp$x)) #! Vector com os limites min e max do N de indivíduos observados color.legend(248000,120000,250000,128000,col_lab,sort(unique(paleta)),gradie nt="y",cex=0.6)#! Legenda text(245300,130500,"Nº Indivíduos",cex=0.6) plot(blocos,ol="grey40",fg=NA,add=T) I'd like to replace the grey shade by other colors. Thanks in a...

Hola buenos días hice un gráfico combinado de líneas, puntos y barras en ggplot2, pero no sé cómo puedo poner la leyenda de eso gráfico para que me represente para las líneas con puntos los valores estimados por un modelo y observados. este es mi código: ggplot(MLM,aes(x=Individuo)) + geom_bar(aes(y=Observada), stat = "identity", color= "gray") + geom_line(aes(y=Estimada), stat = "identity", color="blue", size=1.5) + geom_point(aes(y=Estimada), shape=21, fill="blue", color=&...

...n Estadística ya que es la carrera que estoy siguiendo. El motivo de este mensaje es por un par de dudas que no pude resolver: 1. He utilizado las funciones para realizar pruebas de normalidad (kolmogorov-smirnov, cuando n>50) y (Shapiro, cuando n<50); y he obtenido los resultados ("D observado" en el caso de kolmogorov) pero no he podido hallar una función que me genere los valores de "D tabulado" para hacer el contraste y tomar la decisión de rechazar/no rechazar la hipótesis de normalidad. ¿Alguien sabe una función que me pueda ayudar?. 2. Además alguien me puede ayuda...

...t;. Una vez que hago mi modelo y el summary del mismo obtengo por ejemplo la significación (valor p) de la variable independiente, así como el valor del AIC para comparar mi modelo con otro similar, pero no sé cómo obtener un valor del grado de ajuste de la recta de regresión resultante a los datos observados, algo así como el R al cuadrado ajustado que me da el modelo "lm". ¿Alguien sabe cómo obtenerlo? Muchas gracias de antemano, Salu2, Javier

R-Amigos Mis mejores deseos en estas fiestas y que el año nuevo premie la generosidad con la que los miembros de este grupo comparten su sabiduría y sus trucos. Felicidades -- Sent from my IPAD [[alternative HTML version deleted]]

Estimados compañeros: tengo una duda sobre R; cuando ejecuto la prueba t, R calcula una t experimental y le calcula su probabilidad asociada y la compara con P=0.05 (valor por defecto). Mi duda es si para calcular este area, tipifica primera a z y luego entra en la tabla de z y la calcula y la presenta en pantalla. De otra manera no veo como asociar un valor de t a una probabilidad (la

...################################### # Definición del Data Frame f <- c(x, y, z) # Modelo lineal definido entre las variables f.lm <- lm(z ~ x + y, data=f) # generación de la grilla de las combinaciones para los distintos niveles de ''x'' e ''y'' en el rango observado f.fit <- expand.grid(x=seq(min(f$x),max(f$x),0.5), y=seq(min(f$y),max(f$y),0.5)) # Prediccion de z para la grilla generada usando el modelo lineal obtenido anteriormente z <- predict(f.lm, newdata=f.fit) f.fit$z <- as.numeric(z) f.fit$z <- round(f.fit$z, 2) # Gráfico ggplot para las...

Si te refieres a t.test, lo que hace R es calcular el estadístico t y después compara el valor observado con la distribución de probabilidad de la t para los grados de libertad de tu prueba. Ej: > t.test(1:10, y = c(7:20), alternative="less") Welch Two Sample t-test data: 1:10 and c(7:20) t = -5.4349, df = 21.982, p-value = 9.276e-06 alternative hypothesis: true difference...

Hola Jesús, Si no entiendo mal lo que planteas, has observado unas frecuencias N = c(12, 0, 8, 6, 4, 2) de pedidos acaecidos en 1:6 días, en total n = sum(N) pedidos. Llamemos P = (p1, p2, ?, p6) a las probabilidades ?teóricas? de pedido en cada día. En principio el vector N se puede considerar una realización de una distribución multinomial M(n; p1, ?, p6)....

...Why? > > > > > El 21 de diciembre de 2013, 21:07, Marcelino de la Cruz > <marcelino.delacruz en upm.es <mailto:marcelino.delacruz en upm.es>> escribió: > > Si te refieres a t.test, lo que hace R es calcular el estadístico t > y después compara el valor observado con la distribución de > probabilidad de la t para los grados de libertad de tu prueba. Ej: > > > t.test(1:10, y = c(7:20), alternative="less") > > Welch Two Sample t-test > > data: 1:10 and c(7:20) > t = -5.4349, df = 21.982, p-va...

...o, es decir, consultado las tablas. Gracias El 21 de diciembre de 2013, 21:07, Marcelino de la Cruz <marcelino.delacruz en upm.es <mailto:marcelino.delacruz en upm.es>> escribió: Si te refieres a t.test, lo que hace R es calcular el estadístico t y después compara el valor observado con la distribución de probabilidad de la t para los grados de libertad de tu prueba. Ej: > t.test(1:10, y = c(7:20), alternative="less") Welch Two Sample t-test data: 1:10 and c(7:20) t = -5.4349, df = 21.982, p-value = 9.276e-06 alternative...

Saludo estimados compañeros y compañeras Tengo una matriz de datos de 51 filas por 160 columnas, Cada columna es una serie de tiempo. Existe alguna libreria q me calcule la posicion de la fila donde se hallan estos maximos y minimos o los eventos de estas curvas Agradezco la atencion CARLOS ANDRES [[alternative HTML version deleted]]

Hola: tengo un pequeño problemita que es más bien estadístico que específico de R, pero de todos modos espero me puedan dar una mano. Quiero encontrar una forma de reportar la Bondad de ajuste de un modelo No Lineal. Estuve buscando y leyendo y parece que tal cosa no es posible o al menos no es fiable, desde el punto de vista estadístico. El tema es que quiero mandar un trabajo a una

Hola a todos, -Estoy estudiando el efecto de dos genotipos (~tratamientos) en la aparición de síndrome metabólico (MetS) con datos longitudinales recogidos a tiempo 0,7,10,15,20 y 25 años. -He hecho un dataframe con las siguientes variables MetS: Síndrome Metabólico (Si=1,No=0) bmi: Indice de masa corporal (IMC) cuando se produce la conversión a MetS+ . Para los que permancen MetS-, esta variable

Hola: Desde hace varias versiones no puedo utilizar la ayuda HTML de R. Cuando, por ejemplo, lanzo ?nlm me abre el navegador con ERROR The requested URL could not be retrieved While trying to retrieve the URL: http://127.0.0.1:28683/library/stats/html/nlm.html The following error was encountered: Connection to 127.0.0.1 Failed The system returned: (111) Connection refused The remote host