[R] Repetitive correlation test

[Amarjit Singh Sethi]

Dear Burns
Replacing comma through tilde has not helped; I got the message:

Error in cor.test.formula(frm, alternative = "two.sided", method = "pearson") : 
'formula' missing or invalid.

As you know, the syntax for test for corelation coefficient b/w variables x1 & x2  is as:

cor=cor.test(x1,x2,alternative = "two.sided", method = "pearson")

That is, the variables essentially need be separated by a comma. So please help.

--- On Sat, 25/10/08, Patrick Burns <pburns at pburns.seanet.com> wrote:

> From: Patrick Burns <pburns at pburns.seanet.com>
> Subject: Re: [R] Repetitive correlation test
> To: set_alt at yahoo.co.in
> Date: Saturday, 25 October, 2008, 2:39 PM
> I think you mean to paste a tilde rather than a
> comma in the command that creates 'frm'.
> 
> Patrick Burns
> patrick at burns-stat.com
> +44 (0)20 8525 0696
> http://www.burns-stat.com
> (home of S Poetry and "A Guide for the Unwilling S
> User")
> 
> Amarjit Singh Sethi wrote:
> > Dear all,
> > Through the following code, I wanted to perform
> correlation test repetedly (through loop) on different
> combinations of variables of a data set.
> >
> > Code:
> >
> >
> x=read.table("sample.txt",header=T,sep="\t")
> > out="corout.txt"
> > sink(out)
> > nm = names(x)
> > print(nm)
> > nvr=3
> > # nvr=Total no. of variables in the input data file
> > for (i in 1:(nvr-1))
> > {
> > for(j in (i+1):nvr) 
> > {
> > frm=
> as.formula(paste(nm[i+1],",",nm[j+1],sep=""))
>            ($$)
> >
> crl=cor.test(frm,alternative="two.sided",method="pearson")
> 
> > smr=summary(crl)
> > print(smr)
> > }
> > }
> > sink()
> >
> > Sample Data:
> >
> > S No	V1	V2	V3
> > 1	15	10	4
> > 2	6	4	7
> > 3	10	5	2
> > 4	8	6	6
> >
> > But the code does not work; there seems to be some
> problem in the statement marked ($$) above. The message I
> get is:
> >
> > Error in parse(text = x) : unexpected ',' in
> "V1,"
> >
> > Kindly help to get rid of the bug
> >
> > Regards
> >
> > ajss
> >
> >
> >       Unlimited freedom, unlimited storage. Get it
> now, on
> http://help.yahoo.com/l/in/yahoo/mail/yahoomail/tools/tools-08.html/
> >
> > ______________________________________________
> > 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.
> >
> >
> >


      Download prohibited? No problem. CHAT from any browser, without download. Go to http://in.webmessenger.yahoo.com/