作者Edster (Edster)
看板R_Language
标题[问题] quantile, 如何将预设值全部改掉呢？
时间Mon Nov 16 22:24:18 2015
我de了数天的bug,
终於找到程式结果不如预期的原因(时间就在debug中度过).
是quantile非常吃样本个数, 虽然原先就知道, 但以为这是个小事
哪知当样本数小於20跟大於20时.
同一组样本算出来的95%却落差非常多.
举例来说, 我的资料如下.
x17 = c(-95.290065, -94.676170, -94.042493, -93.766582, -93.281851,
-92.621728, -92.275779, -91.985500, -86.475105, -80.372520,
-80.040033, -78.674849, -68.573205, -60.486220, -57.897714,
-4.204705, 119.349514)
x31 = c(-95.747361, -94.438780, -94.198810, -93.922347, -93.881117,
-93.788013, -93.658992, -93.609898, -93.599487, -93.599094,
-93.533554, -93.517688, -93.478389, -93.475430, -93.364654,
-93.048528, -92.872143, -92.533137, -90.225752, -89.077365,
-84.827016, -84.166512, -69.243811, -68.747692, -65.920400,
-62.867256, -56.184730, -48.477508, 9.279569, 10.153318,
264.201379)
x17 = rep(x17, 17); x31 = rep(x31, 31)
quantile(x17, 0.95, type=6); quantile(x17, 0.95, type=7)
quantile(x31, 0.95, type=6); quantile(x31, 0.95, type=7)
quantile的default是type=7 (R的预设值)
现在我想让这个type在这整只程式里全部改成type=6
(matlab & spss 都用这个, 也是我比较希望的算法.)
请问有什麽方法可以让我电脑里的关於quantile的type全改成type=6
一个个改也是可以, 只是这只程式用到太多次quantile, 有的还写在aggregate里面.
怕会漏掉.
感谢.
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 111.241.16.98
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1447683861.A.A9F.html
※ 编辑: Edster (111.241.16.98), 11/16/2015 22:27:38
1F：→ Edster: 贴错资料, 已修改. 11/16 22:28
2F：→ Edster: 最後我用取代 na.rm=TRUE 成na.rm=TRUE, type=6 11/16 22:34
3F：→ Edster: 花了比编辑这篇文章还少的时间解决这件事. 11/16 22:35
4F：→ Edster: 不过我还是想知道如何变成R中的预设值. 11/16 22:35
5F：→ Wush978: 在.Rprofile里面把quantile这个函数改掉？ 11/16 22:43
6F：→ Wush978: 不过如果你呼叫quantile的函数是来自套件的话，就改不到 11/16 22:44