我需要使用箱线图可视化圆值的分布。 我的值分布在一天中的不同时段(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()
函数使用自定义中位数和四分位数函数?