Re: [R-es] [Grupo Usuarios de R de Madrid]: Siguiente reunión 15-octubre....(*Material Disponible*)
Hola, Ya está disponible todo el material que presentamos durante la reunión del pasado miércoles. Los detalles los podéis encontrar aquí: http://madrid.r-es.org/miercoles-15-de-octubre-2014/ Saludos, Carlos Ortega www.qualityexcellence.es El 15 de octubre de 2014, 10:06, Carlos Ortega c...@qualityexcellence.es escribió: Hola, Un a breve nota para recordar que hoy tendremos la reunión del Grupo de Usuarios de R de Madrid. - La agenda es la misma que anunciamos. - La reunión será en la Facultad de Educación (UNED) en la sala Fernández Huertas, planta 1ª. Los detalles de la convocatoria los podéis encontrar en el portal del grupo: http://madrid.r-es.org/ Gracias, Carlos Ortega www.qualityexcellence.es El 8 de octubre de 2014, 10:47, Carlos Ortega c...@qualityexcellence.es escribió: Buenas a todos, Ya tenemos programada la agenda de la siguiente reunión del grupo. De forma escueta: - *Presentaciones:* - Pedro Valverde https://www.linkedin.com/profile/view?id=232087352authType=NAME_SEARCHauthToken=qJG6locale=es_ESsrchid=51802411410559294801srchindex=1srchtotal=139trk=vsrp_people_res_nametrkInfo=VSRPsearchId%3A51802411410559294801%2CVSRPtargetId%3A232087352%2CVSRPcmpt%3Aprimary : “Modelización de los salarios en España”. - Pedro Concejero: https://www.linkedin.com/profile/view?id=8293681authType=NAME_SEARCHauthToken=VdNxlocale=en_USsrchid=51802411412718027835srchindex=1srchtotal=3trk=vsrp_people_res_nametrkInfo=VSRPsearchId%3A51802411412718027835%2CVSRPtargetId%3A8293681%2CVSRPcmpt%3Aprimary“decidiR – teoría del prospecto – paquete pt http://cran.r-project.org/web/packages/pt/index.html “. - Carlos Ortega: https://www.linkedin.com/in/carlosortegafernandez “Otros diveRtimentos sobre Money-is-a-gas”. - *Píldoras:* - Jorge Ayuso https://www.linkedin.com/profile/view?id=229197095authType=NAME_SEARCHauthToken=GvWulocale=es_ESsrchid=51802411412718420247srchindex=1srchtotal=33trk=vsrp_people_res_nametrkInfo=VSRPsearchId%3A51802411412718420247%2CVSRPtargetId%3A229197095%2CVSRPcmpt%3Aprimary : “OpenCpu http://cran.r-project.org/web/packages/opencpu/index.html y Docker https://www.docker.com/“. Más detalles aquí: http://madrid.r-es.org/miercoles-15-de-octubre-2014/ -- Saludos, Carlos Ortega www.qualityexcellence.es -- Saludos, Carlos Ortega www.qualityexcellence.es -- Saludos, Carlos Ortega www.qualityexcellence.es [[alternative HTML version deleted]] ___ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es
[R-es] Fwd: Question about GEE and interaction plot
Hello all! I am trying to analyze a dataset using the GEE. This is what I have so far C=read.table(Cnarinofinal4.txt,header=T) attach(C) names(C) library(geepack) Form - formula(Conteo~Grupo+Tiempo+Comportamiento+Comportamiento*Grupo+Comportamiento*Hora) M4 - geeglm(Form,family = poisson, id =Individuo, corstr = exchangeable) anova(M4) 1) I want to know if someboy knows how to perform post hoc comparisons for significant terms, in this case I have a significance result for the interaction Comportamiento*Hora and the factor Comportamiento. 2) I want to perform and interaction plot, I tried with the following code: interaction.plot(C$Comportamiento,C$Tiempo,C$Conteo) Nevertheless, I would like to plot the graph in a decreasing order by the mean and also add the error bars. For the first approach I do not have available info, and for the second one I found that data must be sorted before using the function summarySE from the ggplot2 library, nevertheless it does not seem to work anymore. Any help you can provide will be welcome! HoraIndividuo Comportamiento Tiempo Grupo Conteo 0 1 Exploración NOCHE MACHO 1 0 1 interacciones NOCHE MACHO 0 0 1 Mantenimiento NOCHE MACHO 0 0 1 AlimentaciónNOCHE MACHO 5 0 1 ReproducciónNOCHE MACHO 0 0 1 Reposo NOCHE MACHO 2 0 1 Exacavación NOCHE MACHO 0 0 2 Exploración NOCHE MACHO 0 0 2 interacciones NOCHE MACHO 0 0 2 Mantenimiento NOCHE MACHO 0 0 2 AlimentaciónNOCHE MACHO 0 0 2 ReproducciónNOCHE MACHO 0 0 2 Reposo NOCHE MACHO 7 0 2 Exacavación NOCHE MACHO 0 0 3 Exploración NOCHE MACHO 1 0 3 interacciones NOCHE MACHO 0 0 3 Mantenimiento NOCHE MACHO 0 0 3 AlimentaciónNOCHE MACHO 0 0 3 ReproducciónNOCHE MACHO 0 0 3 Reposo NOCHE MACHO 7 0 3 Exacavación NOCHE MACHO 0 0 4 Exploración NOCHE MACHO 3 0 4 interacciones NOCHE MACHO 0 0 4 Mantenimiento NOCHE MACHO 0 0 4 AlimentaciónNOCHE MACHO 1 0 4 ReproducciónNOCHE MACHO 0 0 4 Reposo NOCHE MACHO 5 0 4 Exacavación NOCHE MACHO 0 0 5 Exploración NOCHE MACHO 1 0 5 interacciones NOCHE MACHO 0 0 5 Mantenimiento NOCHE MACHO 0 0 5 AlimentaciónNOCHE MACHO 1 0 5 ReproducciónNOCHE MACHO 0 0 5 Reposo NOCHE MACHO 8 0 5 Exacavación NOCHE MACHO 0 0 6 Exploración NOCHE MACHO 3 0 6 interacciones NOCHE MACHO 0 0 6 Mantenimiento NOCHE MACHO 0 0 6 AlimentaciónNOCHE MACHO 0 0 6 ReproducciónNOCHE MACHO 0 0 6 Reposo NOCHE MACHO 2 0 6 Exacavación NOCHE MACHO 0 0 7 Exploración NOCHE HEMBRA 3 0 7 interacciones NOCHE HEMBRA 0 0 7 Mantenimiento NOCHE HEMBRA 0 0 7 AlimentaciónNOCHE HEMBRA 0 0 7 ReproducciónNOCHE HEMBRA 0 0 7 Reposo NOCHE HEMBRA 5 0 7 Exacavación NOCHE HEMBRA 0 0 8 Exploración NOCHE HEMBRA 2 0 8 interacciones NOCHE HEMBRA 0 0 8 Mantenimiento NOCHE HEMBRA 0 0 8 AlimentaciónNOCHE HEMBRA 0 0 8 ReproducciónNOCHE HEMBRA 0 0 8 Reposo NOCHE HEMBRA 2 0 8 Exacavación NOCHE HEMBRA 0 0 9 Exploración NOCHE HEMBRA 2 0 9 interacciones NOCHE HEMBRA 0 0 9 Mantenimiento NOCHE HEMBRA 0 0 9 AlimentaciónNOCHE HEMBRA 0 0 9 ReproducciónNOCHE HEMBRA 0 0 9 Reposo NOCHE HEMBRA 5 0 9 Exacavación NOCHE HEMBRA 0 0 10 Exploración NOCHE HEMBRA 6 0 10 interacciones NOCHE HEMBRA 0 0 10 Mantenimiento NOCHE HEMBRA 0 0 10 AlimentaciónNOCHE HEMBRA 0 0 10 ReproducciónNOCHE HEMBRA 0 0 10 Reposo NOCHE HEMBRA 2 0 10 Exacavación NOCHE HEMBRA 0 0 11 Exploración NOCHE HEMBRA 1 0 11 interacciones NOCHE HEMBRA 0 0 11 Mantenimiento NOCHE HEMBRA 0 0 11 AlimentaciónNOCHE HEMBRA 2 0 11 ReproducciónNOCHE HEMBRA 0 0 11 Reposo NOCHE HEMBRA 4 0 11 Exacavación NOCHE HEMBRA 0 0 12 Exploración NOCHE HEMBRA 2 0 12 interacciones NOCHE HEMBRA 0 0 12
Re: [R-es] Ayuda prueba de significancia funcion nls()
Hola, ¿qué tal? Hu... no conozco ese modelo que tratas de ajustar demasiado bien, pero creo que el parámetro A significa demasiadas cosas (y depende de demasiadas cosas) a la vez. Por ejemplo, seguro, está muy correlacionado con alfa porque opera como una especie de término independiente de lo que hay dentro de exp y su valor depende mucho, por ejemplo, de si normalizas o no la edad. Yo buscaría una formulación o reparametrización alternativa. Pero en este caso en particular, parecía claro que el problema estaba en la falta de ajuste en esa cola poco relevante. Un saludo, Carlos J. Gil Bellosta http://www.datanalytics.com El día 16 de octubre de 2014, 22:12, Boris Polanco borisc...@gmail.com escribió: Hola, he vuelto a realizar el ajuste con la sugerencia que me diste y ha funcionado, pero mi pregunta es si para con otro modelo similar no me sale un parámetro no significativo, cómo interpreto el no rechazo de la hipótesis nula? saludos El 15 de octubre de 2014, 17:44, Carlos J. Gil Bellosta c...@datanalytics.com escribió: Hola, ¿qué tal? Yo eliminaría las observaciones por encima de los noventaitantos (¿o cien?) años. Son solo ruido, como puede observarse. Estarán basadas en un porcentaje nimio de la muestra. Igual entonces cambian las cosas. Un saludo, Carlos J. Gil Bellosta http://www.datanalytics.com El día 15 de octubre de 2014, 23:18, Boris Polanco borisc...@gmail.com escribió: Hola, estoy haciendo un estudio de la desaceleración de la mortalidad y tengo dudas con la función nls(), adjunto un archivo donde he redactado mi inquietud. saludos cordiales Boris Polanco ___ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es ___ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es