如何计算 h 点（在 R 中）答案

【问题标题】：how to calculate h-point (in R)如何计算 h 点（在 R 中）
【发布时间】：2021-03-25 23:01:51
【问题描述】：

我正在尝试编写一个函数来计算 h 点。该函数是在秩频率数据帧上定义的。考虑以下data.frame：

DATA <-data.frame(frequency=c(64,58,54,32,29,29,25,17,17,15,12,12,10), rank=c(seq(1, 13)))

h点的公式是：

if {有一个 r = f(r), h-point = r } else { h-point = f(i)j-f(j)i / j-i+f(i)-f(j)其中 f(i) 和 f(j) 是第 i 个和第 j 个秩的对应频率，i 和 j 是 if(j) 的相邻秩。

这是我到目前为止所做的：

h_point <- function(data){
  x <- seq(nrow(data))
  f_x <- data[["frequency"]][x]
  h <- which(x == f_x)
  if(length(h)>1) h
  else{
    i <- which(x < f_x)
    j <- which(x > f_x)
    s <- which(outer(i,j,"-") == -1, TRUE)
    i <- i[s[,1]]
    j <- j[s[,2]]
    cat("i: ",i, "j: ", j,"\n")
    f_x[i]*j - f_x[j]*i / (i-j + f_x[i]-f_x[j])
  }
}

在 DATA 中，h 点是 12——因为 x = f_x。然而，

h_point(DATA)
i:   j:   
numeric(0)

我在这里做错了什么？

【问题讨论】：

这不是您之前的问题how to calculate h-point 的重复吗？您是否逐步通过该函数来检查中间结果？

标签： r function if-statement text-mining

【解决方案1】：

我看过你之前的帖子how to calculate h-point，但必须说我不太遵循你计算h点的方法。

根据我找到的h点的定义

_{参考：https://www.researchgate.net/figure/The-definition-of-the-h-point-cf-Popescu-Altmann-2006-25_fig1_281010850}

我认为更简单的方法是使用approxfun创建函数频率（排名），然后使用uniroot找到h点：

get_h_point <- function(DATA) {
    fn_interp <- approxfun(DATA$rank, DATA$frequency)
    fn_root <- function(x) fn_interp(x) - x
    uniroot(fn_root, range(DATA$rank))$root
}

get_h_point(DATA)
#[1] 12

【讨论】：