大佬教程收集整理的这篇文章主要介绍了如何重复 Grubbs 检验并显示列的 p 值,大佬教程大佬觉得挺不错的,现在分享给大家,也给大家做个参考。
我是 R 的新手。这可能是一个愚蠢的问题,但我很难过。我试图反复找出一列的异常值。我遵循了这个 How to repeat the Grubbs test and flag the outliers 并得到了预期的结果。但我想显示 p 值以及数据框中的第三列。我尝试了几件事,但我反复得到一个值或只有两个 p 值。 如何显示所有 p 值?
您可以对函数进行以下更改-
library(outliers)
library(ggplot2)
X <- c(152.36,130.38,101.54,96.26,88.03,85.66,83.62,76.53,74.36,73.87,73.36,73.35,68.26,65.25,63.68,63.05,57.53)
grubbs.flag <- function(X) {
outliers <- NULL
test <- x
grubbs.result <- grubbs.test(test)
result <- data.frame(X=X)
i <- 1
result$pv[i] <- grubbs.result$p.value
while(result$pv[i] < 0.05) {
outliers <- c(outliers,as.numeric(strsplit(grubbs.result$alternative," ")[[1]][3]))
test <- x[!x %in% outliers]
grubbs.result <- grubbs.test(test)
i <- i + 1
result$pv[i] <- grubbs.result$p.value
}
return(data.frame(result,Outlier=(x %in% outliers)))
}
grubbs.flag(X)
# X pv Outlier
#1 152.36 0.012614743 TRUE
#2 130.38 0.007648407 TRUE
#3 101.54 0.237267039 falSE
#4 96.26 0.012614743 falSE
#5 88.03 0.012614743 falSE
#6 85.66 0.012614743 falSE
#7 83.62 0.012614743 falSE
#8 76.53 0.012614743 falSE
#9 74.36 0.012614743 falSE
#10 73.87 0.012614743 falSE
#11 73.36 0.012614743 falSE
#12 73.35 0.012614743 falSE
#13 68.26 0.012614743 falSE
#14 65.25 0.012614743 falSE
#15 63.68 0.012614743 falSE
#16 63.05 0.012614743 falSE
#17 57.53 0.012614743 falSE
以上是大佬教程为你收集整理的如何重复 Grubbs 检验并显示列的 p 值全部内容,希望文章能够帮你解决如何重复 Grubbs 检验并显示列的 p 值所遇到的程序开发问题。
如果觉得大佬教程网站内容还不错,欢迎将大佬教程推荐给程序员好友。
本图文内容来源于网友网络收集整理提供,作为学习参考使用,版权属于原作者。
如您有任何意见或建议可联系处理。小编QQ:384754419,请注明来意。