R/savedata.R

In RNED/agriTutorial: Tutorial Analysis of Some Agricultural Experiments

rice=data.frame(yield=c(
 3.320, 6.101, 5.355, 3.766, 5.096, 7.442, 4.660, 6.573,  7.018, 3.188, 5.595, 6.706, 3.625, 6.357, 8.592, 5.232, 7.016, 8.480,
 5.468, 5.442, 8.452, 5.759, 6.398, 8.662, 6.215, 6.953,  9.112, 4.246, 6.209, 8.042, 5.255, 6.992, 9.080, 6.829, 7.565, 9.660,
 3.132, 6.860, 9.314, 5.389, 6.857, 9.224, 5.217, 7.254, 10.360, 3.864, 5.122, 5.536, 4.311, 4.873, 6.462, 5.915, 5.495, 8.020,
 4.752, 6.780, 6.546, 4.809, 5.925, 7.646, 5.170, 7.442,  9.942, 5.788, 5.988, 6.698, 6.130, 6.533, 8.526, 7.106, 6.914, 9.140,
 4.842, 6.768, 7.414, 5.742, 7.856, 9.016, 5.869, 7.626,  8.966, 4.375, 6.894, 8.508, 4.315, 6.974, 9.680, 5.389, 7.812, 9.896,
 4.507, 4.815, 5.244, 4.875, 4.166, 5.584, 5.400, 4.225,  7.642, 4.756, 5.390, 7.092, 5.295, 5.163, 7.212, 6.046, 4.478, 8.714,
 4.422, 6.509, 8.650, 5.308, 6.569, 8.514, 6.318, 7.991,  9.320, 4.863, 5.779, 6.902, 5.345, 6.164, 7.778, 6.011, 7.362, 9.128,
 4.678, 6.573, 8.032, 5.896, 7.422, 9.294, 7.309, 8.950,  9.712))
rice$nrate=rep(rep(c(0,50,80,110,140),each=9),3)
rice$nitrogen=as.factor(rice$nrate)
rice$Replicate = gl(3, 45,length=135, labels=c("R1","R2","R3"))
rice$management = gl(3, 3,length=135, labels=c("Minimum","Optimum","Intensive"))
rice$variety = gl(3, 1,length=135, labels=c("V1","V2","V3"))
rice$Main = gl(5, 9, length=135)
rice$Sub =  gl(3, 3, length=135)
save(rice, file = "data/rice.RData")


beet=data.frame(yield=c(9.9, 7.8, 10.7, 20.3, 22.6, 23.9, 27.5, 30.3, 29.2, 31.4, 27.2, 33.4, 28.1, 25.7, 31.9))
beet$nrate=rep(c(0,35,70,105,140),each=3)
beet$Replicate = gl(3,1,length=15)
save(beet, file = "data/beet.RData")

 greenrice=data.frame(uptake=c(
1.033, 0.506, 0.692, 0.254, 0.154, 0.111, 0.099, 0.089, 0.250, 0.595, 0.503, 1.089, 0.428, 0.248, 0.484, 0.507, 1.297, 0.493, 0.836, 0.321, 0.142,
0.102, 0.103, 0.142, 0.757, 0.754, 0.613, 0.373, 0.234, 0.498, 0.397, 0.453, 0.493, 0.457, 0.267, 0.477, 0.373, 0.534, 1.007, 0.739, 0.548, 0.588,
0.349, 1.034, 0.133, 0.093, 0.098, 0.138, 0.152, 0.129, 0.084, 0.141, 0.587, 0.305, 0.372, 0.619, 0.677, 0.537, 0.974, 0.327, 0.367, 0.625, 0.831,
0.713))
 greenrice$W=c(
 10,10,10,10,40,40,40,40,0,0,0,0,20,20,20,20,0,0,0,0,40,40,40,40,10,10,10,10,20,20,20,20,20,20,20,20,0,0,0,0,10,10,
 10, 10, 40, 40, 40, 40, 40, 40, 40, 40, 20, 20, 20, 20,  0,  0,  0,  0, 10, 10, 10, 10)
greenrice$N=c(
 270, 90, 180, 0, 90, 180, 0, 270, 0, 180, 90, 270, 90, 0, 180, 270, 270, 90, 180, 0, 270, 180, 0, 90, 270, 180, 90, 0, 0, 270, 90,
 180, 90, 180, 0, 270, 0, 90, 270, 180, 180, 90, 0, 270, 90, 0, 180, 270, 180,  90, 0, 270, 90, 0, 180, 270, 270, 90, 180, 0, 0, 90,
 270, 180)
greenrice$Replicate = gl(4,16)
greenrice$Main = gl(4,4,length=64)
save(greenrice, file = "data/greenrice.RData")

sorghum=data.frame(y=c(
       5.00,4.84,4.02,3.75,3.13,
       4.42,4.30,3.67,3.23,2.83,
       4.42,4.10,3.46,3.09,2.82,
       4.01,3.89,3.21,2.89,2.56,
       3.36,3.10,2.67,2.47,2.16,
       5.82,5.60,5.05,4.72,4.46,
       5.73,5.59,5.00,4.65,4.42,
       5.31,5.19,4.86,4.44,4.22,
       4.92,4.66,4.56,4.16,3.99,
       3.96,3.86,3.50,3.13,2.95,
       5.65,5.97,5.27,5.07,4.52,
       5.39,5.49,5.08,4.87,4.32,
       5.15,5.28,4.93,4.67,4.15,
       4.50,4.89,4.74,4.49,4.10,
       3.75,3.74,3.55,3.28,3.00,
       5.86,5.60,5.37,5.00,4.37,
       5.82,5.55,5.29,4.95,4.07,
       5.26,5.06,4.76,4.48,3.94,
       4.87,4.75,4.55,4.33,3.83,
       3.96,3.76,3.56,3.18,2.96))
sorghum$variety =  gl(4,25,length=100)
sorghum$Replicate= gl(5,5,length=100)
sorghum$factweek = gl(5,1,length=100)
sorghum$factplot = gl(20,5,length=100)
sorghum$varweek = rep(1:5,20)
sorghum$varblock=rep(rep(1:5,each=5),4)
sorghum
save(sorghum, file = "data/sorghum.RData")

 turnip=data.frame(log_yield=c(
0.00,1.39,2.40,2.56,2.48,0.81,1.66,2.33,2.56,2.53,0.69,1.50,2.01,2.56,
2.30,0.92,1.56,1.95,1.87,2.01,0.69,1.79,2.60,2.89,2.76,1.25,2.22,2.63,
2.69,2.98,0.69,1.95,2.37,3.04,2.71,1.01,1.95,2.17,2.40,2.48,1.18,2.22,
2.67,2.71,3.14,1.01,2.40,2.98,3.21,2.89,1.75,2.63,2.56,2.83,2.80,0.92,
2.22,2.46,2.53,2.56))
turnip$Replicate = gl(3,20,length=60)
turnip$Density = gl(5,1,length=60)
turnip$Spacing= gl(4,5,length=60)
turnip$rowspacing = rep(rep(c(4,8,16,32),each=5),3)
turnip$density = rep(c(0.5,2.0,8.0,20.0,32.0),12)
turnip
save(turnip, file = "data/turnip.RData")