search for: observados

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 observados...

...no siempre se pueden combinar > una librería con la otra. > > Javier Marcuzzi > > El 06/11/14 a las 06:48, Lorena Tudela Marco escibió: > >> 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(z...

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

Dear R users, I have been plotting several graphics in only one page. I intend to plot only one legend in each page. Since the graphics are inside 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")

...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 ad...

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=&q...

Estimados todos: Me es grato escribir a esta lista de ayuda para R, ya que comparto 100% la filosofía de Software libre en especial software en 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

...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

Estimados: Estoy entrampado con el siguiente problema. Poseo un dataframe (f), con datos x, y, z, los cuales no son coordenadas ni alturas, sino niveles de tres variables cualquiera. Lo que busco es presentar en un gráfico es: a) el scatterplot de las variables ''x'' e ''y'', y además b) las isolineas de ''z'' para los distintos niveles de

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

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). P se puede estimar mediante las

Porque con "less" estás hacienddo un test de una cola. En tu ejemplo estás haciendo un test de dos colas: > 2*pt(-5.4349,df=21.982) [1] 1.855473e-05 El 21/12/2013 21:33, HERNANDEZ CORONADO JORGE escribió: > t.test(1:10,y=c(7:20)) > > Welch Two Sample t-test > > data: 1:10 and c(7:20) > t = -5.4349, df = 21.982, p-value = 1.855e-05 > alternative hypothesis:

Si miras el código de t.test() verás que para calcular el p-valor llama a pt(). En la página de ayuda de pt() se cuenta que: "For the central case of pt, a normal approximation in the tails, otherwise via pbeta. For the non-central case of pt based on a C translation of Lenth, R. V. (1989). Algorithm AS 243 ? Cumulative distribution function of the non-central t distribution, Applied

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