#Afternoon - so, we decided/stuart pointed out that We only need 680 positions/site, because the seeds will be in between the plug points. So I'm doing the whole thing again but trying to get 680 points, or thereabouts. cc <- read.csv("Plug Datasheet_Verification_04July2010_a.csv") str(cc) Cleaning up cc: cc$tray wantextra<-c("19", "20","21", "22","23", "24", "25", "26", "27", "28") cc cc<-cc[!is.na(cc$tray),] #removed the extra na rows cc<-cc[!cc$tray %in% wantextra,] cc #I have to re-randomize the trays, and make sure that there are 3 of each kind in each plot: sample(1:9,3) sample(10:18,3) #bm = Bob Mahoney: #rune = Runestone: #hlh = Hegg Lake Hill sample(1:9) #[1] 2 9 7 8 3 1 5 4 6 sample(10:18) #[1] 17 16 15 18 13 12 10 14 11 #so: wantbm<-c("2", "9", "7", "17","16", "15") wantrune<-c("8","3", "1", "18", "13", "12") wanthlh<-c("5", "4", "6", "10", "14", "11") cc$site<-"bm" str(cc) cc[cc$tray %in% wantrune,"site"] cc[cc$tray %in% wantrune,"site"]<-"rune" cc[cc$tray %in% wanthlh,"site"] cc[cc$tray %in% wanthlh,"site"]<-"hlh" head(cc) tail(cc) #so cc now has a site column, and the trays have been randomized into each site. table(cc$site, cc$status) # This has the number of plugs/site #I'm going to stick w/ 8 pos/meter, plus 1 pos at the end of each row. I'll stake the 0 and the last positions to mark my rows. Maybe w/ metal disks. 8*8.625*10 #this would work, 690 positions + 10, so total of 700 positions. 8 pos/meter * 8.625 meters/row * 10 rows + the 10 "0" rows #plots will have 5 sets of double rows (10cm apart), each double will be 1m apart. (8*8)+(6) #number of pos. per row (70) 70*10 700*3 # = 2100 = the total number of positions. #stake either end of the row, so we won't be actually planting positions 0 and 8.625 1.1*5 #Plot will be 5.5 meters long, and 8.625 meters wide #10 rows per site, 30 total rows rep(1:30, each=70) df<-data.frame(pos=rep(seq(0,8.625, by= .125), 30), row=rep(1:30, each=70)) str(df) dim(df) 2100-(60) 2040/3 #= 680 = perfect, the max number of envelope positions is 671, so this will be great! df$site<-"bm" wantrune2<-11:20 wanthlh2<-21:30 str(df) df[df$row %in% wantrune2,"site"] df[df$row %in% wantrune2,"site"]<-"rune" df[df$row %in% wanthlh2,"site"] df[df$row %in% wanthlh2,"site"]<-"hlh" str(df) head(df) tail(df) table(df$row, df$site) write.csv(df, "KG_positions_03July.csv", row.names= FALSE)