On Jul 12, 2011, at 7:27 AM, Mitra, Sumona wrote:
Dear all,
I am new to programming in R.
You see to think there is a "++" operation in R. That is not so.
I deal with microarray data,which is a data frame object type. I
need to carry out a few statistical procedures on this, one of them
being the pearson corelation. I need to do this between each row
which is a gene. So the desired result is a square matrix with the
pearson corelation value between each row. So the first column would
be (1,1)=0,(1,2),(1,3) and so on.
I do not understand what that means. You should offer either a minimal
dataset or at the very least the results of str(a).
I uploaded the data frame as "a":-
If by that you mean you made a failed effort at attaching the data in
a file, then you need to read the Posting Guide for what the server
will accept as a file type.
a <- read.csv("a.csv", header= TRUE, row.names=1)
and then I started the script:-
pearson <- function(a){
r <- matrix[x,y]
for(x in as.vector(a[,1], mode="double")){
I do not see a need for as.vector here or at any point later. a[,1]
would already be a vector and if it is not numeric to begin with, then
you are going to get junk.
x++{
for(y in as.vector(a[2,], mode="double")){
y <- x+1
x++
{
r <- (cor.test(as.vector(as.matrix(a)[x,], mode="double"),
as.vector(as.matrix(a)[y,], mode="double"))$p.value)
}
}
}
r[x,y]==r[y,x]
}
return(r)
}
However whenever I run it,I get the error:-
pearson(a)
Error in matrix[x, y] : object of type 'closure' is not subsettable
Please help!
Best Regards
Sumona Mitra
______________________________________________
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.
David Winsemius, MD
West Hartford, CT
______________________________________________
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.