Please read ?lm! -- where it says:
method:
the method to be used; for fitting, currently only method = "qr" is
supported; method = "model.frame" returns the model frame (the same as
with model = TRUE, see below).
More to the point, your request for a "spearman" method for lm() makes
little or no sense. There *are* rank-based methods for multiple
regression, but that sort of discussion is off topic here. I suggest
you talk with a local statistician as you appear to be out of your
depth statistically; or you might try posting on a statistical site
like stats.stackexchange.com
Cheers,
Bert
Bert Gunter
"The trouble with having an open mind is that people keep coming along
and sticking things into it."
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
On Thu, Apr 28, 2016 at 8:35 PM, Hoji, Akihiko <ak...@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@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.
______________________________________________
R-help@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.
