程序代码如下:
t<-0;f<-0;g<-0;k<-0;p<-0;h<-0;l<-0;m<-0;n<-0;o<-0;u<-0;y<-0;r<-0
for (i in 1:186)
{
for (j in 1:25)
{
if(is.na(dusio1m[i,j])==TRUE) {t<-(t+1)}
else if(as.numeric(dusio1m[i,j])<0.92000) {if(is.na(rus1m1[i,j])==TRUE) {o<-(o+1)}
else if(as.numeric(rus1m1[i,j])>0) {f<-(f+1)}
else {h<-(h+1)}}