3个或更多参与者的测量技术误差（TEM）

Question

我需要将此公式实现为R函数。令我惊讶的是，没有R套餐具有TEM功能，可供3名或更多参与者使用。有人可以帮帮我吗？

Answer 1

你去吧

# M : data frame of different measurements
tem <- function(M) {
    nrows <- nrow(M)
    ncols <- ncol(M)
    sqrt(sum(apply(M,1,function(x) sum(x^2) - sum(x)^2/ncols))/(nrows*(ncols-1)))
}

这是一个例子

child_data <- data.frame(
    height_a=c(64.50,71.00,58.00,58.00,70.50,69.00,63.00,65.00,62.00,68.00),
    height_b=c(64.00,71.50,59.00,58.00,71.50,67.50,64.00,64.50,62.00,68.00)
)

tem(child_data)
# 0.5477

这是您引用的论文中提供的数字。

然后，例如（三列）

child_data$height_c <- c(63.5,71.3,59.0,58.5,71.4,67.4,64.2,64.50,62.5,67.5)
tem(child_data)
# 0.5