Confidence bands for rqss quantile smoothed effects
Theo Mouton
Hi all,
I am using the “rqss” method in geom_quantile to plot additive quantile regression curves on a variety of plots. In the quantreg package an option is available to plot confidence bands for the smoothed effects. Is there any way I could plot these same confidence bands on a ggplot?
I have more than 40 graphs to plot and would really like to be able to use ggplots instead of normal R plots.
Thanks,
TM
| Theo Mouton |
| Freshwater Research Assistant |
| +64-7-859-1889 | Gate 10 Silverdale Road, Hillcrest, Hamilton | www.niwa.co.nz |
 |
| To ensure compliance with legal requirements and to maintain cyber security standards, NIWA's IT systems are subject to ongoing monitoring, activity logging and auditing. This monitoring and auditing service may be provided by third parties. Such third parties can access information transmitted to, processed by and stored on NIWA's IT systems. |
Ben Bolker
--
--
You received this message because you are subscribed to the ggplot2 mailing list.
Please provide a reproducible example: https://github.com/hadley/devtools/wiki/Reproducibility
To post: email ggp...@googlegroups.com
To unsubscribe: email ggplot2+unsubscribe@googlegroups.com
More options: http://groups.google.com/group/ggplot2
---
You received this message because you are subscribed to the Google Groups "ggplot2" group.
To unsubscribe from this group and stop receiving emails from it, send an email to ggplot2+unsubscribe@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
Theo Mouton
Thanks Ben for your answer.
Attached is the summary of my rqss model. I am not familiar with using geom_ribbon and when using this code bellow, geom_ribbon applies to my points, not the geom_quantile line. Any help would be greatly appreciated.
pe<-ggplot(data=H1a, aes(x = shade,y = V1))
pe<- pe + geom_point()
pe<- pe + geom_quantile(method="rqss", quantiles = 0.99, lambda=15, linejoin = "round", lineend = "round")
pe <- pe + geom_ribbon(aes(ymin = V1-2.47, ymax = V1+2.47), alpha=0.2)
| Theo Mouton |
| Freshwater Research Assistant |
| +64-7-859-1889 | Gate 10 Silverdale Road, Hillcrest, Hamilton | www.niwa.co.nz |
 |
| To ensure compliance with legal requirements and to maintain cyber security standards, NIWA's IT systems are subject to ongoing monitoring, activity logging and auditing. This monitoring and auditing service may be provided by third parties. Such third parties can access information transmitted to, processed by and stored on NIWA's IT systems. |
From: Ben Bolker [mailto:bbo...@gmail.com]
Sent: Tuesday, 30 May 2017 11:51 a.m.
To: Theo Mouton <Theo....@niwa.co.nz>
Cc: ggp...@googlegroups.com
Subject: Re: Confidence bands for rqss quantile smoothed effects
figure out how to extract the y-values of the lower and upper confidence intervals for every predicted point, and take a look at ?geom_ribbon ...
On Mon, May 29, 2017 at 7:47 PM, Theo Mouton <Theo....@niwa.co.nz> wrote:
Hi all,
I am using the “rqss” method in geom_quantile to plot additive quantile regression curves on a variety of plots. In the quantreg package an option is available to plot confidence bands for the smoothed effects. Is there any way I could plot these same confidence bands on a ggplot?
I have more than 40 graphs to plot and would really like to be able to use ggplots instead of normal R plots.
Thanks,
TM
Theo Mouton
Freshwater Research Assistant
+64-7-859-1889 | Gate 10 Silverdale Road, Hillcrest, Hamilton | www.niwa.co.nz
To ensure compliance with legal requirements and to maintain cyber security standards, NIWA's IT systems are subject to ongoing monitoring, activity logging and auditing. This monitoring and auditing service may be provided by third parties. Such third parties can access information transmitted to, processed by and stored on NIWA's IT systems.
--
--
You received this message because you are subscribed to the ggplot2 mailing list.
Please provide a reproducible example: https://github.com/hadley/devtools/wiki/Reproducibility
To post: email ggp...@googlegroups.com
To unsubscribe: email ggplot2+u...@googlegroups.com
More options: http://groups.google.com/group/ggplot2
---
You received this message because you are subscribed to the Google Groups "ggplot2" group.
To unsubscribe from this group and stop receiving emails from it, send an email to ggplot2+u...@googlegroups.com.