Fix bug where unclustered nodes were being deleted if other nodes that are in clusters were present. Problem was this line clustervals[clustervals == ""] <- NAcharacter made the non-clustered nodes get dropped when this line got run: clusters <- split(nodeblock, clustervals) Since elsewhere in the code, all NA values are converted to "", keeping non-clustered nodes with the "" value avoids them being dropped by split().