Displaying 1 result from an estimated 1 matches for "file1_gene_cor".
2005 Jun 10
0
Top N correlations from 'cor' for very large datasets being run many times
...ader=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_corrs=corrs[corrs_ordered[1:n]] #Get top n correlations - these would be added to some data structure and used for subsequent analysis
}
[[alternativ...