手上有组偏态分布资料,对数变换还不行,
试图对5%和95%等百分位数的CI进行确定,
现在试用R,50% 可以用wilcox.test,
5%,这些怎么弄呢,看文献是使用bootstraps,按文献思路是
1.进行bootstraps抽样(有放回,抽样n个)抽B次,
2.对每次抽样计算P5
3.将B个P5作为新序列,估算2.5~97.5%作为范围P5的范围,
觉得这里面的第三步是有问题的。 抽样B个P5还是偏态的,是不是应该用中位数,wilcox.test 来估计其可信限呢。
求大神,是不是有别的方法。求程序。