[R] lm() with spearman corr option ?

Cade, Brian cadeb at usgs.gov
Fri Apr 29 16:38:38 CEST 2016


I think you would just need to replace the lm() function call with
cor(x,y,method="spearman".  It would probably be more informative to
actually plot by the magnitude of the correlation coefficient (all |r| >=
0.20 or something similar) rather than just by those with P <=0.05.

Brian

Brian S. Cade, PhD

U. S. Geological Survey
Fort Collins Science Center
2150 Centre Ave., Bldg. C
Fort Collins, CO  80526-8818

email:  cadeb at usgs.gov <brian_cade at usgs.gov>
tel:  970 226-9326


On Thu, Apr 28, 2016 at 9:35 PM, Hoji, Akihiko <akh22 at pitt.edu> wrote:

> Hi,
>
> A following function was  kindly provided by GGally’s maintainer, Barret
> Schloerke.
>
> function(data, mapping, ...) {
>     p <- ggplot(data = data, mapping = mapping) +
>         geom_point(color = I("blue")) +
>         geom_smooth(method = "lm", color = I("black"), ...) +
>         theme_blank() +
>         theme(panel.border=element_rect(fill=NA, linetype = "solid",
> color="black"))
>
>     lmModel <- eval(substitute(lm(y ~ x, data = data), mapping))
>     fs <- summary(lmModel)$fstatistic
>     pValue <- pf(fs[1], fs[2], fs[3], lower.tail = FALSE)
>
>     if (pValue < 0.05) {
>         p <- p + theme(
>             panel.border = element_rect(
>                 color = "red",
>                 size = 3,
>                 linetype = "solid",
>                 fill = "transparent"
>             )
>         )
>     }
>
>     p
> }
>
> Basically, this function draws red squares  over pairwise corr plots with
> p<0.05.  Now, since I need to use the spearman rank corr, I tried to modify
> the lm function by adding “method=spearman”  but this did not work at al.
> Could anybody suggest the way to add the spearman rank corr function in
> this particular function ?
>
> Thanks.
>
>
>
>
>
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>

	[[alternative HTML version deleted]]



More information about the R-help mailing list