Hi, Chung-hong Chan wrote: > Hi, > > You can try with > cor.test(rank(y[1]),rank(y[2])) >
Thanks for this! It didn't solve my problem, but it helped me realize that the formula I was using by hand is invalid for the tie case. I just realized that with R's cor function, the Pearson correlation of the ranks = the Spearman of the original values. I've yet to find the formula for the tied case for Spearman, but at least now I know what the problem is (the formula I was using by hand). Thanks! Ray ______________________________________________ R-help@stat.math.ethz.ch 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.