corresponding pdf

oxplot(Di.change ~ Fertilizer, data = BS, col = c("brown", "red"))
ggplot(data = BS, aes(x = Fertilizer, y = Di.change)) + geom_boxplot()
7.5
2 4 6 8
Di.change
5.0
2.5
F
NF
F
Fertilizer
NF
(d) Use the tapply command to find the numeric summaries of the diameter changes for the two levels
of fertilization.
with(data = BS, tapply(Di.change, list(Fertilizer), summary))
$F
Min. 1st Qu. Median Mean 3rd Qu. Max.
2.91 4.32 4.76 5.27 6.52 8.92
$NF
Min. 1st Qu. Median Mean 3rd Qu. Max.
1.02 1.92 2.71 2.72 3.16 5.71
with(data = BS, tapply(Di.change, list(Fertilizer), sd))
F NF
1.383 1.101
(e) Create a scatter plot of the height changes against the diameter changes and describe the relationship.
plot(Ht.change ~ Di.change, data = BS, cex = 0.5, pch = 19, col = "blue")
p