有个40X5的矩阵,每列是40个用户测同一样本的结果,开始时用grubbs.test敲5遍命令检测离群值,之后发现了apply和sapply函数可以批量操作,但是有个问题:怎么在apply(data, margin, fun)里的fun参数里加参数。
现在用下面的函数跑的结果是无法在grubbs.test里写参数,只能用默认参数检测1个离群值:
apply (data, 2, grubbs.test)
像下面的写法就报错:
apply (data, 2, grubbs.test (type = 11, opposite = FALSE, two.sided = TRUE) )
而我的数据需要在高值和低值两头都检测离群值:
grubbs.test(data[,1], type = 11, opposite = FALSE, two.sided = TRUE)
看grubbs.test函数说明,传入的数据只能是个向量。。。。