Use<-c("Used","Available")
par(mfrow = c(2,2),mar = c(3,4,3,2))
boxplot(habitat$dfir~Use[habitat$use],col = "darkgreen",
main = "Fir basal area measurement\n in used and available habitat",
ylab = "Basal Area Measurement")
boxplot(habitat$junip~Use[habitat$use],col = "blue",
main = "Juniper basal area measurement\n in used and available habitat",
ylab = "Basal Area Measurement")
boxplot(habitat$aspen~Use[habitat$use],col = "orange",
main = "Aspen basal area measurement\n in used and available habitat",
ylab = "Basal Area Measurement")
boxplot(habitat$other~Use[habitat$use],col = "brown",
main = "Other basal area measurement\n in used and available habitat",
ylab = "Basal Area Measurement")
Run the code above in your browser using DataLab