我有 6 个基于分位数函数创建的变量(我在“0.25、.5 和 .75”的默认数量之上创建了不同的分位数。我想创建一个类似于以下输出的汇总表,但使用我的变量。任何帮助将不胜感激。
我使用 cbind 将所有变量合并到一个表中,但格式已关闭。
ten_quantile = quantile(tt$HPTACC, probs = .1)
twenty_quantile = quantile(tt$HPTACC, probs = .2)
thirty_quantile = quantile(tt$HPTACC, probs = .3)
med = quantile(tt$HPTACC, probs = .5)
sixty_quantile = quantile(tt$HPTACC, probs = .6)
upper = quantile(tt$HPTACC, probs = .75)
rr <- cbind(ten_quantile, twenty_quantile, thirty_quantile, med, sixty_quantile, upper)
rr
ten_quantile twenty_quantile thirty_quantile med sixty_quantile upper
10% 7172.1 17137.72 22099.3 29542 32302 35864.5
但是我想有这样的东西,但有额外的分位数:
Sfttime <- quantile(tt$HPTACC)
## Create quantile soft time into df ##
SummaryTable <- data.frame(
Sfttime)
Sfttime
0% 0.00
10% #
20% #
25% 20072.25
35% #
50% 29542.00
60% #
75% 35864.50
100% 889368.00
我认为您可以使用quantile()函数。您可以配置probs参数中的值。我将使用mtcars作为示例,但它应该适用于您的数据
QuantData <- quantile(mtcars$drat)
QuantData <- as.data.frame(QuantData)
Result <- data.frame(Quantile = row.names(QuantData),
Value = QuantData[,1])
QuantData <- quantile(mtcars$drat,
probs = c(0,0.1,0.20,0.25,0.35,0.5,0.6,0.75,1))
QuantData <- as.data.frame(QuantData)
Result <- data.frame(Quantile = row.names(QuantData),
Value = QuantData[,1])
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句