[R-es] Leyenda gráfico combinado

Carlos Ortega co| @end|ng |rom qu@||tyexce||ence@e@
Lun Oct 5 00:28:06 CEST 2020


Hola,

Usando el paquete que comentaba y cogiendo el data.frame con el paquete
"datapasta"...
He aplicado un ajuste polinómico de grado 3.

#--------------------------------------------
library(data.table)
library(ggplot2)

datIn <- data.table(
  Individuo = c(1L,2L,3L,4L,5L,6L,7L,8L,9L,10L,
                11L,12L,13L,14L,15L,16L,17L,18L,19L,20L,21L,22L,23L,
                24L,25L,26L,27L,28L,29L,30L),
  Observada = c(200.21,154.817,514.919,234.068,232.191,
                379.53,297.466,450.94,317.84,308.16,178.317,530.919,
                368.224,251.352,154.866,309.39,291.637,513.398,322.25,
                451.932,612.163,574.745,92.629,451.836,449.408,689.179,
                885.454,484.617,876.734,506.156),
  Estimada = c(180.75719,110.00147,455.28532,226.17628,
               218.58544,346.95982,309.05514,421.53012,276.81604,
               305.13638,182.79552,492.87962,347.28844,255.40259,143.53278,
               307.57602,288.5641,454.39712,323.53048,442.87195,555.30366,
               518.25317,56.79211,413.1442,445.26074,741.37221,735.90541,
               450.47909,742.06702,521.66028)
)

datIn_lg <- melt(datIn, id=1)
names(datIn_lg) <- c('Individuo', 'Tipo', 'Valor')

library(ggpmisc)
formula <- y ~ poly(x, 3, raw = TRUE)
ggplot(datIn_lg, aes( Individuo, Valor, group = Tipo, color = Tipo)) +
  geom_point() +
  geom_smooth(method = "lm", formula = formula) +
  stat_poly_eq(aes(label =  paste(stat(eq.label), stat(adj.rr.label),
                                  sep = "*\", \"*")),
               formula = formula, parse = TRUE) +
  labs(x = expression(italic(x)), y = expression(italic(y))) +
  facet_wrap(~Tipo) +
  theme_bw()
#----------------------------

Saludos,
Carlos Ortega
www.qualityexcellence.es

El dom., 4 oct. 2020 a las 23:16, Jimmy Erney Reyes Velasco (<
jimmyreyesvelasco using gmail.com>) escribió:

> Lo siento creo que olvide poner los datos.
> mis datos son:
>
> Individuo Observada Estimada
> 1 200.210 180.75719
> 2 154.817 110.00147
> 3 514.919 455.28532
> 4 234.068 226.17628
> 5 232.191 218.58544
> 6 379.530 346.95982
> 7 297.466 309.05514
> 8 450.940 421.53012
> 9 317.840 276.81604
> 10 308.160 305.13638
> 11 178.317 182.79552
> 12 530.919 492.87962
> 13 368.224 347.28844
> 14 251.352 255.40259
> 15 154.866 143.53278
> 16 309.390 307.57602
> 17 291.637 288.5641
> 18 513.398 454.39712
> 19 322.250 323.53048
> 20 451.932 442.87195
> 21 612.163 555.30366
> 22 574.745 518.25317
> 23 92.629 56.79211
> 24 451.836 413.1442
> 25 449.408 445.26074
> 26 689.179 741.37221
> 27 885.454 735.90541
> 28 484.617 450.47909
> 29 876.734 742.06702
> 30 506.156 521.66028 por otro lado intenté poner la ecuación para un
> modelo potencial pero no pude, no pude saber cómo escribir el argumento
> para la fórmula.
>
> El dom., 4 de oct. de 2020 a la(s) 14:09, Carlos Ortega (
> cof using qualityexcellence.es) escribió:
>
>> Hola,
>>
>> Además de la alternativa de Emilio, como estás buscando incluir una
>> fórmula y el valor de ajuste sobre el gráfico, quizás te interese usar las
>> funcionalidades de este paquete:
>>
>>    - *ggpmisc*:
>>    https://cran.r-project.org/web/packages/ggpmisc/vignettes/user-guide.html
>>    - Y en particular la función *stat_poly_eq()*.
>>    - Puedes ver ejemplos prácticos de esta función sobre a mitad de la
>>       viñeta.
>>       - Y que por cierto, permite además ajuste y representación del
>>       tipo *nls()* sobre tus datos .
>>
>> Saludos,
>> Carlos Ortega
>> www.qualityexcellence.es
>>
>> El dom., 4 oct. 2020 a las 16:07, Jimmy Erney Reyes Velasco (<
>> jimmyreyesvelasco using gmail.com>) escribió:
>>
>>> Hola buenos días
>>> hice un gráfico combinado de líneas, puntos y barras en ggplot2, pero no
>>>>>> 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="white")+
>>>   labs(y = expression(paste("Biomasa ", (g/cm^{2}))))+
>>>   ggtitle("Espeletia standleyana")+ theme_minimal()+
>>>   annotate(geom = 'text',
>>>            x = 0,
>>>            y = 750,
>>>            hjust = 0,
>>>            label = "Biomasa=220,774(IAF)+39,759(Aprom)-163,883")+
>>>   annotate(geom = 'text',
>>>            x = 0,
>>>            y = 700,
>>>            hjust = 0,
>>>            label = "R² = 0,982      P=3,8187E-17")
>>>
>>> ¿Alguien sabe cómo hacerlo? Agradezco mucho cualquier información.
>>> que tengan buen día
>>>
>>> <
>>> https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail
>>> >
>>> Libre
>>> de virus. www.avast.com
>>> <
>>> https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail
>>> >
>>> <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
>>>
>>>         [[alternative HTML version deleted]]
>>>
>>> _______________________________________________
>>> R-help-es mailing list
>>> R-help-es using r-project.org
>>> https://stat.ethz.ch/mailman/listinfo/r-help-es
>>>
>>
>>
>> --
>> Saludos,
>> Carlos Ortega
>> www.qualityexcellence.es
>>
>

-- 
Saludos,
Carlos Ortega
www.qualityexcellence.es

	[[alternative HTML version deleted]]



Más información sobre la lista de distribución R-help-es