Thank Dennis!!!.....it works....
regards
Gaurav Kumar
www.gauravkumar.org
PhD Student, Chemistry and Biomolecular Sciences, Macquarie , Sydney, Australia.
MS (Computational Biology), NCBS-TIFR, Bangalore, India.
--- On Wed, 20/1/10, Dennis Murphy <djmu...@gmail.com> wrote:
From: Dennis Murphy <djmu...@gmail.com>
Subject: Re: [R] R help for read.table()
To: gau...@gauravkumar.org
Cc: r-help@r-project.org
Date: Wednesday, 20 January, 2010, 4:43 PM
R does vectorized arithmetic. Your loop was superfluous. Observe:
> x <- read.table(textConnection("
+ S:C 2.011085038928
+ S:A 21.496800549900762
+ S:J 0.183181039138149
+ P:E 9.641984304606304"))
> x
V1 V2
1 S:C 2.0110850
2 S:A 21.4968005
3 S:J 0.1831810
4 P:E 9.6419843
> transform(x, pval = 2 * (1 - pnorm(abs(V2))))
V1 V2 pval
1 S:C 2.0110850 0.04431648
2 S:A 21.4968005 0.00000000
3 S:J 0.1831810 0.85465597
4 P:E 9.6419843 0.00000000
HTH,
Dennis
On Wed, Jan 20, 2010 at 3:04 AM, Gaurav Kumar <gau...@gauravkumar.org> wrote:
I've a table containing two columns
seperated by space, as shown below.
S:C 2.011085038928
S:A 21.496800549900762
S:J 0.183181039138149
P:E 9.641984304606304
I'm reading this table inside a loop
but unable to access the first column as file[i,1].
file <-
read.table("data.txt",header=FALSE)
for ( i in 1:nrow(file))
{
cat(file[i,1],":")
pvalue <-
2*(1-pnorm(abs(file[i,2])))
cat(pvalue,"\n")
}
can anyone help me how to read the
first column.
regards
Gaurav Kumar
www.gauravkumar.org
PhD Student, Chemistry and Biomolecular Sciences, Macquarie , Sydney, Australia.
MS (Computational Biology), NCBS-TIFR, Bangalore, India.
[[alternative HTML version deleted]]
______________________________________________
R-help@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.
[[alternative HTML version deleted]]
______________________________________________
R-help@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.
