[R-es] Dibujar rectas sobre varios un grafico multipanel

[Oscar Perpiñan]

Hola,

Prueba esto:

xyplot(mean ~ dia|sol, data = dats,
       type=c('p', 'r'), layout=c(4, 1),
       ylim=extendrange(c(min(dats$ll), max(dats$ul))),
       ll = dats$ll, ul = dats$ul,
       panel=function(x, y, subscripts, ll, ul, ...) {
           panel.segments(as.numeric(x), ll[subscripts],
                          as.numeric(x), ul[subscripts], ...)
           panel.xyplot(x, y, ...)
       })

Otra opción es latticeExtra::segplots pero no es fácil añadir las
líneas de regresión.

Saludos.

Oscar.
Oscar Perpiñán Lamigueiro
Grupo de Sistemas Fotovoltaicos (IES-UPM)
Dpto. Ingeniería Eléctrica (EUITI-UPM)
URL: http://procomun.wordpress.com
Twitter: @oscarperpinan


El día 24 de junio de 2013 21:10, neo <ericconchamunoz en gmail.com> escribió:
> -----BEGIN PGP SIGNED MESSAGE-----
> Hash: SHA1
>
> Estimada comunidad, tengo que hacer un grafico con 4 paneles (4 graficos
> pegados juntos) y sobre este tengo que dibujar 4 rectas ajustadas a los
> datos de cada panel. Intente hacerlo con lattice usando el siguiente
> codigo :
>
> library(lattice)
> library(latticeExtra)
>
>
> xYplot(Cbind(mf[,3],mf[,5],mf[,6])~dia|factor(sol),data=mf,xlab="Time
> (days)",ylab="Unsaturation index", method="bars",
> keys="lines",layout=c(4,1),abline=list(list(a=4.532e-04,b=-2.138e-05),list(a=7.137e-05,b=2.699e-05),list(a=3.627e-04,b=2.552e-05),list(a=6.251e-04,
> b=4.938e-05)))
>
> y obtengo el grafico que adjunto, el que no responde a mi necesidad pues
> debo dibujar una recta en cada panel.
>
> Alguna idea de como hacerlo ?
>
> Adjunto los datos para que prueben si lo desean.
>
> Muchas gracias por vuestro tiempo,
>
> Eric.
> -----BEGIN PGP SIGNATURE-----
> Version: GnuPG v1.4.12 (GNU/Linux)
> Comment: Using GnuPG with Thunderbird - http://www.enigmail.net/
>
> iQEcBAEBAgAGBQJRyJm0AAoJEDoUMoeY1RuM9OIH/0u0H4iZ1aKlWb09TPKFEMNz
> DZ6fyI6w0Bl5mwavAP/DLlPf4IIikhOR7d0nn/fWGRXOJdPDw34LDWJa/BxPf1CY
> 1tP5TfDAIfaNX9xIaaOuqJgztL1uA7+FDb0LGZY4MkCv/nhIW9NPbSz8M+u+VZw+
> KI+J/Gcymm/uHEH9T/dvb4ihQQ17Lg8xmDIRUYiVtxG9eaCqZaBnosxsH5oLC7oF
> /9aZBC09VnhtpSMtMk3jGMWtCKVHwt7rYrRzgNNVXb2xtl+Kq5ti7mgQ9WQfaGE9
> ATEgKxJzBZFf66FEdXuXra99lXfCN1/qBqqMvD6cw0B02C0MR+3mhzgbG3cDhys=
> =S8jM
> -----END PGP SIGNATURE-----
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>