search for: comparar

Estimados miembros de la lista, Estoy haciendo una análisis de supervivencia con R. Adjunto mis datos. Quiero analizar la supervivencia de 5 grupos diferentes y compararla. Para ello estoy utilizando el paquete survival. > s = Surv(c$tiempo, c$estado) > f = survfit(s ~ tratamiento, data = c) > d = survdiff(s ~ tratamiento, data = c) > d Call: survdiff(formula = s ~ tratamiento, data = c) N Observed Expected (O-E)^2/E (O-E)^2/V tra...

Hola, Manuel, No entiendo tu pregunta (la repito aqui para que sea mas explicito): hay alguna forma de comparar la probabilidad de supervivencias (en este caso anual) entre grupos sin utilizar un chi-cuadrado y un valor de P. Entiendo que lo que hace survdiff es comparar las curvas de supervivencia, pero yo quiero comparar la probabilidad de supervivencia entre grupos al final del estudio. Con survdiff (log...

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 realmente obs<-table(...

Tengo dos dataframe que corresponden a los valores horarios de los contaminantes y el otro a los valores horarios de las variables meteorológicas para el mismo periodo de tiempo. La variable temporal date en ambos dataframe tiene el formato dmy_hms. Los dataframe tienen diferente número de filas. dataframe contaminantes contiene 175247 filas y el dataframe meteorología 175321. Quiero saber

Quiero comparar varias dosis letales 50% (LD50) usando análisis probit. He seguido un ejemplo que viene en paquete DRC, pero no obtengo el resultado esperado. Lo que quiero es saber si las LD50s, son diferentes y si la diferencias son estadísticamente significativas. Gracias de antemano. José Arturo e-mail. j...

Buenas, ¿alguien sabe qué puede ser este error? data es una df y las columnas 23 y 24 son las que quiero comparar; están bien. testK<-kappa2(data[,c(23,24)], "equal") Error in sort.list(y) : 'x' must be atomic for 'sort.list' Have you called 'sort' on a list? -- Dr Manuel Mendoza Department of Biogeography and Global Change National Museum of Natural History (MNCN) Spani...

Esta vez no es un problema, pero quizás podáis aconsejarme sobre qué gráfico elegir, pues hice algunas pruebas y no queda muy bien. Quiero comparar 9 rectas (o lo que sea) hechas de tan solo dos puntos. Para 9 grupos de especies distintos tengo un dato para 2050 y otro para 2070, y pretendo que se aprecie bien la evolución (de tan solo dos fechas) de cada uno de los 9 grupos, y las diferencias entre grupos. Podría ser en un mismo gráfico o un...

Hola a todos: quiero consultarles para estar seguro de que estoy entendiendo bien el funcionamiento de la función smoothSpline() del paquete 'timeSeries'. Tengo una serie temporal con datos mensuales a la cual quiero suavizar usando splines para, por ejemplo, comparar con otras series temporales. Por lo que estuve viendo, me conviene usar la función smoothSpline() que se basa en smooth.spline() del paquete 'stats'. Pero me quedan algunas dudas respecto de los argumentos que se utilizan en ambas funciones. Principalmente, el que más me interesa es 'nk...

....e. 10) Se observa la respuesta de las especies a cada tratamiento. Para cada especie tendremos varios valores (20) para C y varios valores (20) para H. Hacemos un modelo lineal en el que se prueba el efecto del tratamiento teniendo en cuenta las especies como un factor aleatorio: NO me interesa comparar entre especies, solo avisar al modelo que debido a que son especies diferentes habrán diferencias. Después, independientemente del test anterior, hacemos para cada especie el test para evaluar el efecto del tratamiento. p.e. Para la especie 1 comparamos los 20 valores de C con los 20 valores de...

...n ggplot, con 29 categorías, por lo que tengo que utilizar library(RColorBrewer) para disponer de suficientes colores. El problema es que al hacerlo acorde a dos variables distintas (color=var1 y color=var2), cuyas 29 categorías son obviamente las mismas, les da distintos colores y no puedo comparar los mapas. Muchas gracias, Manuel . -- Dr Manuel Mendoza Department of Biogeography and Global Change National Museum of Natural History (MNCN) Spanish Scientific Council (CSIC) C/ Serrano 115bis, 28006 MADRID Spain

...sv original de los autores no hace el análisis porque algo falta y no nos percatamos de que es. Adjunto tabla original Saludos cordiales #ESTO ES LO QUE VIENE EN LA AYUDA Y TRABAJA BIEN rm(list = ls()) library(FactoMineR) data(mortality) library(epicalc) des(mortality) #Agregue esto para comparar las dos datas res<-MFA(mortality,group=c(9,9),type=c("f","f"), # REALIZA EL ANÁLISIS BIEN name.group=c("1979","2006")) plot(res,choix="freq",invisible="ind",habillage="group") #CON LA TABLA ORGINAL en csv NO N...

Hola. Tengo que representar en un mapa, valores de proporciones de extranjeros con tarjetas de residencia en vigor, según nacionalidad y país de procedencia. Los datos son de este estilo Rumania Andalucía 0,200365823 Aragón 0,377192073 Asturias 0,205353012 I. Baleares 0,06592893 Canarias 0,03050312 Cantabria 0,190397393 Castilla-León 0,211670421 Castilla-La Mancha

...maneras, creo que coges días al azar para calcular to "proxy". Aunque yo personalmente cogería días consecutivos porque probablemente el consumo en muchos productos no sea independiente temporalmente. Por otro lado, de cara a la implementación no sé si sería mejor coger 6 o siete días o comparar ambas distribuciones. Más que nada porque si se te acaba el stock en la +1 después de que puedas hacer tu pedido entonces el tiempo real serían 7 días. Espero que te parezca bien el feedback. Estaría encantado de discutir esto con los compañeros. Un abrazo Javier -- [[alternative HTML vers...

...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(zeroinfl, type="prob") [,...

...a a tod en s, estoy trabajando con modelos de regresión mixta con el paquete "nlme", concretamente con la función "lme". 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

Hola Lorena: Lo prometido es deuda. Dije que miraría si iba este asunto: ggplot(data = dLSaa, aes(x = factor(AA), y = AD, colour = factor(AA))) + geom_boxplot() + xlab( "AA" ) + ylab( "AD" ) + labs(colour = "Leyenda") + ggtitle(label = "Comparación de AD frente a AA") Y este es el resultado: [image:

...rdiales > > > > > > > > #ESTO ES LO QUE VIENE EN LA AYUDA Y TRABAJA BIEN > > > > rm(list = ls()) > > > > library(FactoMineR) > > > > data(mortality) > > > > library(epicalc) > > > > des(mortality) #Agregue esto para comparar las dos datas > > > > res<-MFA(mortality,group=c(9,9),type=c("f","f"), # REALIZA EL ANÁLISIS > > BIEN > > > > name.group=c("1979","2006")) > > > > plot(res,choix="freq",invisible="ind"...

Saludos: Quisiera saber si se pueden comparar dos documentos con R para detectar si uno está plagiado del otro. Los dos documentos los tengo en pdf. Gracias Javier Gómez González [[alternative HTML version deleted]]

Bom dia eu chamo-me Ana, estou a tentar fazer matriz de correla??o e covariancia para comparar 4 variaveis e saber quais s?o as mehores... fiz o segundo codigo mas devo tar a fazer alguma confus?o para nao me aparecer o output correcto: cor(dados[,2:5],method=c("pearson")) cor_with_p_test(dados[2:5]) ct <- cor.test(dados[,2:5]) ct$p ct$v #matriz de covariancias cov(dados[,2:5...

...1 5 0 0 0 2 0 2 4 3 1 2 6 0 4 0 2 0 0 0 0 12 0 3 2 10 1 0 0 0 0 0 0 0 0 11 0 3 3 0 0 0 7 5 3 0 0 4 0 0 0 3 4 0 0 0 1 5 1 NA 0 5 0 0 0 3 5 2 0 0 5 5 0 0 0 0 3 0 2 3 6 4 2 0 4 2 0 0 0 4 6 0 0 necesito crear un boxplot donde x siempre será Ep (episodios), e Y considere por ejemplo GS y GI, para poder comparar dichos comportamientos. De preferencia, me gustaría hacerlo simplemente con boxplot(), pero estoy abierta a intentarlo con ggplot de igual forma. Muchas gracias, saludos -- *Lorena Saavedra A.**Ing. Recursos Naturales Renovables* *+56 9 9880 2972* [[alternative HTML version deleted]]