2

在 R 中，使用库中的'geyser'数据集MASS，我们正在构建一个非参数公差区间 cont。90% 的间歇泉喷发 w/ 99% CI

请在下面找到我正在使用的方程式：

no_param_pi <- function(x, conf.level = 0.90) {
  n <- length(x)
  x <- sort(x)
  j <- max(floor((n + 1) * (1 - conf.level) / 2), 1)
  conf.level <- (n + 1 - 2 * j)/(n + 1)
  interval <- c(x[j], x[n + 1 - j])
  attr(interval, "conf.level") <- conf.level
  interval
}

这是我正在使用的代码，其中我遗漏了一些东西......：

library(MASS)
library(tolerance)
head(geyser)

eruptions = subset(geyser, select = waiting:duration, increasing)

nptol.int(eruptions, alpha = 1 - 0.99, P = 0.9)

寻找两者

4

1 回答 1

0

你能试试这几行代码吗：

eruptions <- geyser
nptol.int(eruptions$duration, alpha = 1 - 0.99, P = 0.9, side = 1)
nptol.int(eruptions$duration, alpha = 1 - 0.99, P = 0.9, side = 2)
nptol.int(eruptions$waiting, alpha = 1 - 0.99, P = 0.9, side = 1)
nptol.int(eruptions$waiting, alpha = 1 - 0.99, P = 0.9, side = 2)

这是你要找的吗？你有很多变化，所以我建议你查看帮助下的文档。?nptol.int

于 2022-02-10T00:12:50.203 回答