setosarersicolorvirginicasetasa 曲 rsicolorM「gimca一、数据可视化1•对于 iris 数据,用每类花(iris$Speciees)的样本数作为高度,制作条形图。> height5<-tablspecies^> barplotCheights)setosaversicolorvirginica2•用每类花的Sepal.Length、Sepal.Width、Petal.Length、Petal.Width 的平均值分别制作条形图,四图同显。> barplotCtapplyCSepal・L^ngth,^口已匚厲已^‘⑴ 已日和),f 口1=4)abarplotCtapplyCSepal.NrcIth’SpE 匚厲已耳,5 已日门) barplot(tappl/(petal・I_E 门gth,5 卩已:彳已3,厲已日门)?匚口1=2)abarplot(tappl/(petal■训rch:h,2pE 匚 b,rriEarQ,©cd=1)osetose■xersicolorVirginia^or-io3•分别制作Sepal.Length、Sepal.Width、Petal.Length、Petal.Width 的直方图(用密度值做代表,设置 prob=T),添加拟合的密度曲线,四图同显。> hist(sepal.Length,prob=T,co"l=53> 11nesCdensltyC^epal.Length^> hl5r(sepal.width,prob=T,col=23> lines(densityC^epal.width^> hist(petal.Length,prob=T,col=3^)> "lines(densityCPetal.Length^ahist(Petal■划彳 cJnh,p「citi=T,匸□1=4)> "lines(densityCPetal■诃■ichzhj)二、中国地图:(Note:首先从网上下载 GIS 数据,解压到 GIS_data 目录http://cos.name/wp-content/uploads/2009/07/chinaprovinceborderdata_tar_gz.zip)setwd('F:/GIS_data')###设置工作目录install.packages('maptools');library(maptools)chinav-readShapePoly('bou2_4p.shp')###获得各省的边界信息HistogramofSepal.LeHiatogrannofSepal.Hist 口 gramofFetal.LengthHistogramofPetal.Widthplot(china)'w:/GiS-data1)instal1.packages'maptool5':11braryCmaptools)1nsnal1.packages'sp':;11braryf'sp')china<-readshapepoly:/bouZ_4p.5hpOp1-oT(china:;>>names(map_data)[1]"AREA""PERIMETER""BOU2_4M_""BOU2_4M_ID""ADCODE93"⑹"ADCODE99""NAME"可以看出 map_data 中有 7 列,对应的字段名如上面显示。>>map_data$AREA#925 个区域单元的面积>>map_data$PERIMETER#925 个区域单元的周长>>map_data$BOU2_4M#有重复数字,特定情况下可作为区域单元 ID>> map_data$ADCOD#93 版 ADCODE 地理编码>> map_data$ADCOD#99 版 ADCODE 地理编码>>map_data$NAME#各区域单元所隶属的省级行政单元的名>>unique(map_data$NA#查看各区域的名称是什么文本...