Advertisement
Guest User

Untitled

a guest
Apr 26th, 2015
196
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.47 KB | None | 0 0
  1. attach(ufa)
  2.  
  3. ufa_subset<-ufa[c("X","Y","FACILITYID","VICINITY","TBOX_L","TBOX_W","TBOX_STAT","SCI_NM","CMMN_NM","DATE_PLANT","DBH")]
  4.  
  5. ufa_subset$FACILITYID <- gsub("-","",ufa_subset$FACILITYID)
  6. ufa_subset$TBOX_STAT <- ufa_subset$TBOX_STAT=="Plant"
  7. ufa_subset$DATE_PLANT <- sub("T.*","",ufa_subset$DATE_PLANT)
  8.  
  9. ufa_subset$Other <- as.character(lapply(strsplit(as.character(SCI_NM), split=" x "), "[", 1))
  10. ufa_subset$Cultivar <- as.character(lapply(strsplit(as.character(SCI_NM), split=" x "), "[", 2))
  11. ufa_subset$Genus <- sub("^(\\w+)\\s?(.*)$","\\2",ufa_subset$Other)
  12. ufa_subset$Species <- sub("^(\\w+)\\s?(.*)$","\\1",ufa_subset$Other)
  13.  
  14. ufa_subset$Genus <- sub("^\\(See Notes\\)$","NA",ufa_subset$Genus)
  15. ufa_subset$Species <- sub("^Other$","NA",ufa_subset$Species)
  16.  
  17. ufa_subset$Other <- NULL
  18. ufa_subset$SCI_NM <- NULL
  19.  
  20. names(ufa_subset)[names(ufa_subset)=="X"] <- "Point X"
  21. names(ufa_subset)[names(ufa_subset)=="Y"] <- "Point Y"
  22. names(ufa_subset)[names(ufa_subset)=="FACILITYID"] <- "External ID Number"
  23. names(ufa_subset)[names(ufa_subset)=="VICINITY"] <- "Street Address"
  24. names(ufa_subset)[names(ufa_subset)=="TBOX_L"] <- "Plot Length"
  25. names(ufa_subset)[names(ufa_subset)=="TBOX_W"] <- "Plot Width"
  26. names(ufa_subset)[names(ufa_subset)=="TBOX_STAT"] <- "Tree Present"
  27. names(ufa_subset)[names(ufa_subset)=="CMMN_NM"] <- "Common Name"
  28. names(ufa_subset)[names(ufa_subset)=="DATE_PLANT"] <- "Date Planted"
  29. names(ufa_subset)[names(ufa_subset)=="DBH"] <- "Diameter"
  30.  
  31. write.csv(ufa_subset,'~/ufa_subset.csv')
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement