以下是我的程序,
> mydata<-read.csv("E:/data/cangshan/rtest2.csv")
> clean<-mydata[-which(is.na(mydata)),]
> inside<-clean[clean$range==1,]
> outside<-clean[clean$range==0,]
> diff.slope<-vector()
> for (i in 1:nrow(inside)){
+ diff.slope<-cbind(diff.slope,abs(inside$slope-outside$slope))
+ }
> diff.slope<-t(diff.slope)
> diff.dem<-vector()
> for (i in 1:nrow(inside)){
+ diff.dem<-rbind(diff.dem,abs(inside$dem-outside$dem))
+ }
> diff.settlement<-vector()
> for (i in 1:nrow(inside)){
+ diff.settlement<-rbind(diff.settlement,abs(inside$settlement-outside$settlement))
+ }
> diff.road<-vector()
> for (i in 1:nrow(inside)){
+ diff.road<-rbind(diff.road,abs(inside$road-outside$road))
+ }
> whole.matrix<-diff.slope+diff.dem+diff.settlement+diff.road
> diff.sort<-sort(whole.matrix)
> pair.index<-vector()
> for ( j in 1:length(diff.sort)){
+ pair.index<-rbind(pair.index,which(whole.matrix==diff.sort[j],arr.ind = T))
+ }
>
> a<-pair.index[!duplicated(pair.index[,1]),]
运行到这里,便提示Error in pair.index[, 1] : subscript out of bounds。