Hi Lindsey,
Below are the plot treatment assignments. I included some notes about how the decisions
were made. I also included the R script I wrote to do this.
PlotNo Treatment
[1,] "1" "0.5 yr 2"
[2,] "2" "0.50 yr 1"
[3,] "3" "0.25 yr 1"
[4,] "4" "Control"
[5,] "5" "0.75 yr 1"
[6,] "6" "0.25 yr 1"
[7,] "7" "Control"
[8,] "8" "0.5 yr 2"
[9,] "9" "0.75 yr 1"
[10,] "10" "0.5 yr 1"
Notes:
Plots were divided into two blocks based on position (Plots 1-5 and Plots 6-10).
each block will get one of the two replicate treatments assigned to it.
Treatments are as follows:
Control
0.25" in yr 1
0.50" in yr 1
0.75" in yr 1
0.50" in yr 1 and 2
Plots 4 and 7 were selected as controls because they contained
few sugar maple trees. Since the controls are untreated, we can
use some nearby sugar maple trees from outside the plot.
Plot 1 was assigned a 0.5" treatment in year two because it is closer to
Hubbard Brook than plots 2-5, which will make it easier to spray. Ian said
he would prefer to have the one he has to spray twice closer to the stream.
Plot 6 was selected as a 0.25" treatment because it is odd and the 0.25"
treatment was considered less important than the other treatments.
At plot 6 the lysimeter chemistry is different, the soils are rocky, and
the vegetation is not ideal (i.e., no red maple, not many trees).
R script used to assign treatments to plots:
PlotNo <- c(1,2,3,4,5,6,7,8,9,10)
Treatment <- c("0.5 yr 2","","","Control","","0.25 yr 1","Control","","","")
Assign <- cbind(PlotNo,Treatment)
Block_1_to_5 <- sample(c("0.25 yr 1","0.50 yr 1","0.75 yr 1"), 3, replace=F)
Block_6_to_10 <- sample(c("0.5 yr 1","0.5 yr 2","0.75 yr 1"), 3, replace=F)
Assign[2, "Treatment"] <- Block_1_to_5[1]
Assign[3, "Treatment"] <- Block_1_to_5[2]
Assign[5, "Treatment"] <- Block_1_to_5[3]
Assign[8, "Treatment"] <- Block_6_to_10[1]
Assign[9, "Treatment"] <- Block_6_to_10[2]
Assign[10, "Treatment"] <- Block_6_to_10[3]
Assign
|
