将自定义中位数和四分位数函数用于ggplot2 boxplot()

时间:2019-02-05 13:51:34

标签: r ggplot2 boxplot

我需要使用箱线图可视化圆值的分布。 我的值分布在一天中的不同时段(0小时= 24小时):

value <- c(18.8,1.2,4.4,3.8,3.4,0.2,2.8,2.3,4.8,7.8,21.5,1.5,22.2,23.4,1.6,19.9,0.5,20.5,21.1,22.2,18.7,17.1,16.7,16.4,16.3,18.1,15.7,18.9,16.7,15.2)
condition <- c(rep("A", 10), rep("B",10), rep("C",10))
df <- cbind.data.frame(condition, value)

我想使用circular R包来计算圆中值以及第一和第三圆四分位数:

library(circular)

get_circular_median <- function(vector){
  circular <- circular(vector, units=c("hours"), rotation=c("clock"), modulo=c("2pi"))
  circular_median <- (median.circular(circular))[[1]]
  return(circular_median)
}

get_third_circular_quartile <- function(vector){
  circular <- circular(vector, units=c("hours"), rotation=c("clock"), modulo=c("2pi"))
  third_quartile <- quantile.circular(circular, probs=c(0.25))[[1]]
  return(third_quartile)
}

get_first_circular_quartile <- function(vector){
  circular <- circular(vector, units=c("hours"), rotation=c("clock"), modulo=c("2pi"))
  first_quartile <- quantile.circular(circular, probs=c(0.75))[[1]]
  return(first_quartile)
}

是否可以对ggplot2的boxplot()函数使用自定义中位数和四分位数函数?

0 个答案:

没有答案