<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>