Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- <pre>
- library(outliers)
- grubbsFunction = function(x, p.value) {
- outliers = NULL
- remainedData = NULL
- test = x
- grubbs.result = grubbs.test(test)
- pv = grubbs.result$p.value
- while(pv>p.value) {
- outliers = c(outliers,as.numeric(strsplit(grubbs.result$alternative," ")[[1]][3]))
- test = x[!x %in% outliers]
- grubbs.result = grubbs.test(test)
- pv = grubbs.result$p.value
- }
- remainedData = setdiff(x, outliers)
- df = data.frame(X=x,Outlier=(x %in% outliers))
- return(list(outliers = outliers, remainedData=remainedData, df=df))
- }
- </pre>
- <pre>
- x= c( 0.0000000, 0.0000000 , 0.0469971 , 0.0000000 , 0.0000000 ,14.0940000,
- 0.0000000 , 5.1100000 ,0.0630035, 0.0000000, 34.0150000 ,0.0000000
- , 0.0000000, 9.0930000 , 0.0000000 , 1.1090000)
- </pre>
- <pre>
- grubbs.test(x, type = 10, opposite = FALSE, two.sided = FALSE)
- </pre>
- <pre>
- grubbs.test(x,opposite=FALSE,two.sided=TRUE)
- </pre>
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement