Loading required package: vcd library(DAAG) data(rainforest) names(rainforest) rainforest summary(rainforest) str(rainforest) ?rainforest Making packages.html . done + 1 dim(rainforest) attach(rainforest) hist(dbh) cor(rainforest[,1:6]) cor(rainforest[,1:6],use="pair") library(vmv) tablemissing(rainforest, sortby="both") dev.new() hist(dbh) dev.new() hist(root) dev.new() hist(rootsk) help.start() apply(rainforest[,1:6],2,hist) dev.new() par(mfrow=c(3,2)) apply(rainforest[,1:6],2,hist) median(rainforest[,1:6]) mean(rainforest[,1:6]) mean(rainforest[,1:6],na.rm=TRUE) median(rainforest[,1:6],na.rm=TRUE) barplot(table(species)) dev.new() barplot(table(species)) dev.new() plot(wood) hist(wood) dev.new() stripchart(wood) stripchart(wood,pch=19) dev.new();tablemissing(rainforest, sortby="both") dev.new();tablemissing(rainforest, sortby="variable") dev.new() splom(rainforest[,1:6]) library(YaleToolkit) gpairs(rainforest) dev.new() dev.new() plot(dbh,wood,col=species,pch=19) dev.new();plot(dbh,wood,col=species,pch=19,log="xy") ggplot(rainforest,aes(x=dbh,y=wood))+geom_point()+facet_grid(.~species) ggplot(rainforest,aes(x=log(dbh),y=log(wood)))+geom_point()+facet_grid(.~species) dev.new();plot(species,wood)