Displaying 1 result from an estimated 1 matches for "cor_ij".
Did you mean:
car_id
2005 Jun 10
0
Top N correlations from 'cor' for very large datasets being run many times
...all possible correlations for each row and then store only those that I want to keep. - too slow
file1 = read.table("test.txt", header=F, quote="", sep="\t", comment.char="", as.is=1)
for (i in 1:15000){
corrs = vector(length=15000)
for (j in 1:15000){
cor_ij = cor(as.numeric(file1[i,]), as.numeric(file1[j,]), method = "pearson", use = "pairwise.complete.obs")
corrs[j]=file1_gene_cor
}
corrs[i]=NA; #Set correlation to self as NA
corrs_ordered=order(corrs,decreasing=TRUE) #Order correlations from largest to smallest
top_cor...