从零开始的 R 中值函数

R median function from scratch

我是 R 初学者,我尝试从头开始制作中值函数。

这是我的代码:

mymedian <- function(x) {
  len <- length(x)
  sorted <- sort(x)
  if (len %% 2 == 0) {
    med1 <- sorted[len / 2]
    med2 <- sorted[(len + 1) %/% 2]
    result <- sorted[med1 + med2 / 2]
    return(result)
  } else {
    result <- sorted[(len + 1)/2]
    return(result)
  }
}

我得到“NA”输出。我找不到问题所在。

主要问题是您试图用非整数(例如 168.5)索引您的 sorted 向量。将您的函数与此进行比较:

mymedian <- function(x){
  len <- length(x)
  sorted <-sort(x)
  if(len%%2==0){
    i <- len/2
    med1<-sorted[i]
    med2 <- sorted[i+1]
    result <- sum(med1,med2)/2
    return(result)
  }else{
    result<-sorted[(len+1)/2]
    return(result)
  }
}