library(rrcov)
data(milk)#rrcov包里的数据
x=scale(milk)
mcd=CovMcd(x)
S=getCorr(mcd);S#S是稳健的相关系数阵。将上面代码运行多次,每次S不一样
键入CovMcd看看其代码,发现有一项seed参数,默认是NULL
然后我尝试着设置参数seed=2
mcd=CovMcd(x,seed=2)
出现下面错误:
Error in covMcd(x = x, raw.only = raw.only, alpha = alpha, nsamp = nsamp, :
invalid 'seed'. Must be compatible with .Random.seed !
这个问题怎么解决,或者如何设置seed使其每次出现的S一样?