Re: [R-es] Leyenda gráfico combinado

[Carlos Ortega]

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 (<
jimmyreyesvela...@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 (
> c...@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 (<
>> jimmyreyesvela...@gmail.com>) escribió:
>>
>>> 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="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_source=link_campaign=sig-email_content=webmail
>>> >
>>> Libre
>>> de virus. www.avast.com
>>> <
>>> https://www.avast.com/sig-email?utm_medium=email_source=link_campaign=sig-email_content=webmail
>>> >
>>> <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
>>>
>>>   

Re: [R-es] Leyenda gráfico combinado

[Jimmy Erney Reyes Velasco]

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 (
c...@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 (<
> jimmyreyesvela...@gmail.com>) escribió:
>
>> 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="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_source=link_campaign=sig-email_content=webmail
>> >
>> Libre
>> de virus. www.avast.com
>> <
>> https://www.avast.com/sig-email?utm_medium=email_source=link_campaign=sig-email_content=webmail
>> >
>> <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
>>
>> [[alternative HTML version deleted]]
>>
>> ___
>> R-help-es mailing list
>> R-help-es@r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-help-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

Re: [R-es] Leyenda gráfico combinado

[Carlos Ortega]

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 (<
jimmyreyesvela...@gmail.com>) escribió:

> 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="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_source=link_campaign=sig-email_content=webmail
> >
> Libre
> de virus. www.avast.com
> <
> https://www.avast.com/sig-email?utm_medium=email_source=link_campaign=sig-email_content=webmail
> >
> <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-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

Re: [R-es] Leyenda gráfico combinado

[Emilio L. Cano]

Hola,

Tal vez tu problema esté en el rango del eje X. No puedo comprobarlo porque tu 
código no es reproducible, pero lo siguiente funciona:

ggplot(data.frame(x=101:200, y = rnorm(10, 100)), 
   aes(x, y)) + geom_bar(stat="identity") +  
  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")

Si x=individuo, entonces el annotate puede quedar fuera del rango (aunque 
debería ajustarse automáticamente como en el código que pego)

Con un ejemplo reproducible será más fácil ayudarte, bastará con una muestra de 
los datos (fácil con dput, mira aquí):
https://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example
 


Un saludo,
Emilio




> El 4 oct 2020, a las 16:07, Jimmy Erney Reyes Velasco 
>  escribió:
> 
> 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="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
> 
> 
> Libre
> de virus. www.avast.com
> 
> <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
> 
>   [[alternative HTML version deleted]]
> 
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-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] Leyenda gráfico combinado

[Jimmy Erney Reyes Velasco]

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="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


Libre
de virus. www.avast.com

<#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es