数据分析展现与语言

stock.data<- AdjustedClose")

3

4WhauVolcano")

5type="n",xlab="MetresWest",

6z=volcano,color.palette=terrain.colors,plot.title=title(main="TheTopographyofxlab="MetersNorth",ylab="Metersplot.axes={axis(1,seq(100,800,by=100))axis(2,seq(100,600,by=100))},key.axes=axis(4,seq(90,190,by=10)))

7y=10*1:ncol(volcano),z=volcano,color.palette=terrain.colors,plot.title=title(main="TheTopographyofMaungaWhau",xlab="MetersNorth",ylab="Metersplot.axes={axis(1,seq(100,800,by=axis(2,seq(100,600,by=key.axes=axis(4,seq(90,190,by=10)))

8z<-2*x<-10*y<-10*(1:ncol(z))zlim<-range(z)zlen<-zlim[2]-zlim[1]+1colorlut<-terrain.colors(zlen)col<-colorlut[z-zlim[1]+1]rgl.surface(x,y,z,color=col,

9

for(iin1:length(rownames(USArrests))){for(jin1:length(x$names)){}}sd<-data.frame(col=colors,}widths=c(8,1),heights=c(8,1))

#ColorScalepar(mar=image(x=1,y=0:length(breaks),z=t(matrix(breaks))*1.001,xlab="",ylab="",xaxt="n")map("state",boundary=map("state",col="white",add=TRUE)title("MurderRatesbyUSStatein1973\n(arrestsper100,000residents)",line=2)

map("county","new

map('italy',fill=TRUE,col=

load(url("A_adm1.RData"))administrativeregions")

zoom=10,destfile="London.png",maptype=" PlotOnStaticMap(london,lat=air$lat,lon=air$lon,

destfile pe="sa PlotOnStaticMap(london,lat=air$lat,lon=

73.99867),zoom=14,destfile="Manhattan.png",maptype="hybrid");

GetMap.OSM(lonR=c(-74.67102,-74.63943),latR=c(40.33804,40.3556),scale=7500,destfile=

'sKeyholeMarkupLanguage(KML)formatcitiesage="rgdal")[1],"cities")

sfdata)plot(sfdata,col="orange",

shapefilessfdata<-read.shapefile(sf)write.shapefile(sfdata,"newsf")

main="Relationshipbetweencardistanceandxlab="Speed( perhour)",ylab="Distancetravelled(

main="Relationshipbetweencardistanceandspeed",xlab="Speed( perhour)",ylab="Distancetravelled(

main="Relationshipbetweencardistanceandxlab="Speed( perhour)",ylab="Distancetravelled(

#plotcommandhere#plotcommandhere

for(iin1:3)

for(iin1:3)

hipbetween",PM[10],"andconcentrations(",mu*g^-concentrations(",mu*g^-

f(x)==frac(1,sqrt(2*pi)*sigma)~ mtext(expression(paste("where",mu,"isthemeanofthedistributionand",sigma,"thestandarddeviation.")),

plot(1:200,type="n",main="FontsunderWindows",axes=FALSE,xlab="",text(0,180,"Arial\n(family=\"sans\",font=1)",text(0,140,"ArialBold\n(family=\"sans\",text(0,100,"ArialItalic\n(family=\"sans\",font=3)",text(0,60,"ArialBoldItalic\n(family=\"sans\",text(70,180,"Times\n(family=\"serif\",font=1)",text(70,140,"TimesBold\n(family=\"serif\",

text(70,100,"TimesItalic\n(family=\"serif\",font=3)",text(70,60,"TimesBoldItalic\n(family=\"serif\",font=4)",text(130,180,"CourierNew\n(family=\"mono\",text(130,140,"CourierNewBold\n(family=\"mono\",font=2)",text(130,100,"CourierNewItalic\n(family=\"mono\",text(130,60,"CourierNewBoldItalic\n(family=\"mono\",

==lines(inter,col="black",lwd=0.8) pal<-colorRampPalette(c("#f2f2f2","black"))colors<-pal(100)fsub<-flights[flights$airline=="AA",]t<-max(fsub$cnt)for(jin1:length(fsub$airline))air1<-airports[airports$iata==fsub[j,]$airport1,]air2<-airports[airports$iata==inter<-gcIntermediate(c(air1[1,]$long,air1[1,]$lat),c(air2[1,]$long,air2[1,]$lat),n=100,addStartEnd=TRUE)colindex<-round((fsub[j,]$cnt/ t)*length(colors))}

pal<-colorRampPalette(c("#f2f2f2","red"))colors<-pal(100)fsub<-flights[flights$airline=="AA",]fsub<-fsub[order(fsub$cnt),]t<-for(jin1:length(fsub$airline))colindex<-round((fsub[j,]$cnt/ t)*length(colors)) #Uniquecarriers<-#colors<-for(iin1:length(carriers))pdf(paste("carrier",carriers[i],".pdf",sep=""),width=11,

fsub<-flights[flights$airline==carriers[i],]fsub<-fsub[order(fsub$cnt),]t<-for(jin1:length(fsub$airline))air1<-airports[airports$iata==fsub[j,]$airport1,]air2<-airports[airports$iata==inter<-gcIntermediate(c(air1[1,]$long,air1[1,]$lat),c(air2[1,]$long,air2[1,]$lat),n=100,colindex<-r