spat_data <- read.table("BaldDistAndArea_2011.txt", header=T)
unique(spat_data$Bald)
bald_names <- c(1,  29,  32,  42,  50,  57,  59,  62,  67,  87,  88,
                91,  93, 103)
bald <- c("1a", "29a", "32", "42", "50", "57", "58",
          "62a","67", "87", "88","92", "93", "26a")
mat_dis <- matrix(rep(0,length(bald)*length(bald)),length(bald),length(bald))

for (i in 1:length(bald)){
  for ( j in 1: length(bald)){
    if(i!=j){
    mat_dis[i,j] <- spat_data$NEAR_DIST.m.[spat_data$Bald==bald[i]  &
                                           spat_data$Bald_near==bald[j]]   
    }
  }
}

colnames(mat_dis) <- bald_names
write.table(round(mat_dis,0),file="Bald_mat_dis.txt", sep = " ",
            col.names = TRUE)