r - R中的四分位变异系数
问题描述
我正在尝试计算每 3 次连续测量的 R 的四分位变异系数。
cqv_versatile(
Q,
na.rm = TRUE,
digits = 3)
Q
49
44
34
33
37
48
20
48
37
42
44
35
40
有人知道如何计算cqv_versatile
这个连续间隔吗?
解决方案
这是使用基数 R 的一种方法,将hst
变量更改为所需的子长度(为了清楚起见,第一个系数是根据索引 1、2、3 的值估计的;第二个系数是根据索引为 4、5 的值估计的,6)
Q=as.numeric(read.table(text="49 44 34 33 37 48 20 48 37 42 44 35 40"))
hst=3
sapply(seq(1,length(Q),hst),function(x){
tmp=x:(pmin(x+hst-1,length(Q)))
q1=quantile(Q[tmp],0.25)
q3=quantile(Q[tmp],0.75)
as.numeric((q3-q1)/(q3+q1))
})
[1] 0.08771930 0.09677419 0.19718310 0.05521472 0.00000000
正如评论中提到的,你应该在估计系数时使用更多的数据。
如果您想根据滚动窗口进行估计,那么索引 1,2,3 -> 2,3,4 -> ... 只需替换seq(1,length(Q),hst)
为seq(1,length(Q),1)
.
推荐阅读
- pandas - 字典到 plot.bar 有两个值,一个在 y 轴上,另一个在 x 轴上,带有 dict.key
- powershell - PowerShell - 将父文件夹名称添加到文件名
- python - Pandas:循环遍历许多列并在执行 value_counts 时将不同的映射应用于列
- django - Django Pagination 不工作并显示每个页面中的所有项目
- python - Centos 7 - Dlib 被编译为使用 SSE41 指令,但这些在您的机器上不可用
- javascript - tensorflow.js 向模型添加许多图像样本会填充视频卡内存并崩溃
- awk - 如何使用 awk 在文件增长时连续读取文件?
- vba - 阅读文档中的每个单词并将其替换为单词的第一个字母
- javascript - nodemon 错误:错误:找不到模块“vue”
- django - django - 将文件上传到现有表单上的模型