请教一个R语言报错问题

目前在用R语言做关于copula函数可靠性分析中，
遇到报错
Error in svd(x, nu = 0L, nv = 0L) : 'x'里有无穷值或遗漏值
详细报错是：
1: In dfunc[i] <- (f.dy - f.y)/0.001 :
number of items to replace is not a multiple of replacement length
原程序如下
#4#以下函数计算标准独立空间中的梯度
GRAD<-function(func,y){
dy<-c(0,0,0,0,0,0)
f.dy<-c(0,0,0,0,0,0)
f.y<-c(0,0,0,0,0,0)
dfunc<-c(0,0,0,0,0,0)
for(i in 1:6){
dy<-y
dy[i]<-y[i]+1e-3
f.dy<-func(dy)
f.y<-func(y)
dfunc[i]<-(f.dy-f.y)/1e-3
}
return(dfunc)
}
恳请指导一下。万分感谢！