rr1<- matrix(s,400,300)
k=10
a=2
c=a*k#列
#对行进行运算P <- sum(rr1*rr1)
S1 <- apply(rr1,1,function(x){
x <- sum(abs(x))
x <- x*x
return(x)
})
h1<- S1/P
H=rr1[,rank(h) > length(h) -c]#抽取列
e1=rr1-H1
F1=sum(e1*e1)
e=F1/P
请问各位老师,我想给这段代码加上一段循环,使a从2开始,每次加上2,直到50结束,a=a+2,结果能列出每一次的误差e,应该怎么做?谢谢各位。


雷达卡


京公网安备 11010802022788号







