proc multtest用PFDR计算q value，输出的p值和原始p值一样是怎么回事？

在SAS里用两种方法做FDR校正

1. data aa;
   
2.    input Set_ID $ Raw_P @@;
   
3.    datalines;
   
4. test55        0.046040529
   
5. test56        0.567874575
   
6. test57        0.559127236
   
7. ……
   
8. run;
   
9. proc multtest inpvalues(Raw_P)=aa
   
10. fdr PFDR;
    
11. id Set_ID;
    
12. run;

复制代码

但PFDR出来的结果跟原始的p值一样，不明白是什么问题（但有些数据集又可以输出校正后的p值），请问有没有大神知道。
另外一个问题：SAS里计算q value的code（也就是PFDR）和R的qvalue  得出的结果不一样，参考的都是Storey提出的“正假发现率”  ，那么两种软件计算的结果为什么不同呢……

应该是要将P值进行从小到大排序再行计算把