关于rest工具箱里的reho.m

Submitted by Lighthawk on Mon, 11/21/2011 - 21:14

老师、学长好~ 请教个问题
在使用 REST_V1.5_101101 这个工具包后，自己根据定义写了一个计算reho的代码。即对每个体素的所有时间点直接升序sort，然后计算KCC值。没有加入任何优化。但计算出来的结果，每个体素的reho值普遍比工具箱里算的值相对偏小。模糊可以看到，自己算的图像黑白的区域与工具箱的结果是很相似的。
回头看工具箱里的代码时，发现reho.m计算rank时取了平均，注释写的是"Compute mean of tied ranks"。
我的问题是：
1、tied ranks是什么意思啊？
2、为什么不直接用[result, Index] = sort(matrix) 里的Index，而自己重新写了一个rank呢？
3、reho.m里的处理方法比原定义的方法有什么优点吗？
谢谢老师、学长回复哈~~

Submitted by YAN Chao-Gan on Thu, 11/24/2011 - 05:45

没错，% Compute mean of tied ranks就是指“排序中有等值的情况出现，此时取序的均值”
当时为了提高ReHo计算效率，把for循环变成矩阵运算，为了优化这一段费老大劲了。