Re: [R] Correlation matrix one side with significance

From: Henrique Dallazuanna <wwwhsd_at_gmail.com>
Date: Thu, 06 Mar 2008 16:17:26 -0300

Please provide a example of what you want

On 06/03/2008, Martin Kaffanke <technik_at_roomandspace.com> wrote:
>
> Am Mittwoch, den 05.03.2008, 14:38 -0300 schrieb Henrique Dallazuanna:
>
> > Try this:
> >
> > On 05/03/2008, Martin Kaffanke <technik_at_roomandspace.com> wrote:
> > > Hi there!
> > >
> > > In my case,
> > >
> > > cor(d[1:20])
> > >
> > > makes me a good correlation matrix.
> > >
> > > Now I'd like to have it one sided, means only the left bottom side to be
> > > printed (the others are the same) and I'd like to have * where the
> > > p-value is lower than 0.05 and ** lower than 0.01.
> > >
> > > How can I do this?
> >
> > d <- matrix(rexp(16, 2), 4)
> > corr <- cor(d)
> > sign <- symnum(cor(d), cutpoints=c(0.05, 0.01), corr = T,
> > symbols=c("***", "**", "*"), abbr=T, diag=F)
> >
> > noquote(mapply(function(x, y)paste(x, format(y, dig=3), sep=''),
> > as.data.frame(unclass(sign)), as.data.frame(corr)))
>
>
> Seems that we mark the value itself, but not the p-value.
>
> So lets say, in a way I have to get the lower left half of a
>
> cor(el[1:20])
>
> Then I need to calc all the values with a cor.test() to see for the
> p-value. And the p-value should be lower than .05 or .01 -> this should
> make the * to the value.
>
> Thanks,
>
> Martin
>
>

-- 
Henrique Dallazuanna
Curitiba-Paraná-Brasil
25° 25' 40" S 49° 16' 22" O

______________________________________________
R-help_at_r-project.org mailing list
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.
Received on Thu 06 Mar 2008 - 19:24:01 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 06 Mar 2008 - 20:30:19 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive