NSW

R/format.data/NSW.R

@@ -13,12 +13,24 @@ data.nswbrc$Date.of.measure <- as.vector(sapply(data.nswbrc$Date.of.measure, fun
     unlist(strsplit(x, "/"))[3]}))  ## Extract the years only
 
 
+## write a function that create a 10*10m quadrat per plot
+fun.quadrat <- function(x,y,max.x,max.y){
+vec.x <- seq(0,max.x,by=10)
+vec.y <- seq(0,max.y,by=10)
+x.cut <- cut(x,breaks=vec.x)
+y.cut <- cut(y,breaks=vec.y)
 
-data.nswbrt <- read.csv("../../data/raw/NSW/NSW_data_BRtreatments.csv", header = TRUE, 
+return(unclass(factor(paste(x.cut,y.cut))))
+
+}
+
+
+
+data.nswbrt <- read.csv("data/raw/NSW/NSW_data_BRtreatments.csv", header = TRUE, 
     stringsAsFactors = FALSE)
 data.nswbrt$Date.of.measure <- as.vector(sapply(data.nswbrt$Date.of.measure, function(x) {
     unlist(strsplit(x, "/"))[3]}))  ## Extract the years only
-data.nswbs1 <- read.csv("../../data/raw/NSW/NSW_data_BS1.csv", header = TRUE, stringsAsFactors = FALSE)
+data.nswbs1 <- read.csv("data/raw/NSW/NSW_data_BS1.csv", header = TRUE, stringsAsFactors = FALSE)
 
 data.nswbs1$Date.of.measure <- as.character(format(as.Date(data.nswbs1$Date.of.measure, 
     format = "%d-%b-%y"), format = "%d/%b/%Y"))