基本绘图系统.R

hist(airquality$Wind)

hist(airquality$Temp,xlab = 'Temp')

boxplot(airquality$Wind)

boxplot(airquality$Wind,xlab='wind',ylab='speed(mph)')

boxplot(Wind~Month,airquality,xlab='Month',ylab='wind')

plot(airquality$Wind,airquality$Temp)

with(airquality,plot(Wind,Temp,main = 'wind and temp in NYC'))

with(airquality,plot(Wind,Temp,main = 'wind and temp in NYC',type='n'))
with(subset(airquality,Month==9),
     points(Wind,Temp,col='red'))


with(subset(airquality,Month==5),
     points(Wind,Temp,col='blue'))

with(subset(airquality,Month %in% c(6,7,8)),
     points(Wind,Temp,col='black'))


fit=lm(Temp~Wind,airquality)

abline(fit,lwd=2)

legend('topright',pch=1,
       col=c('red','blue','black'),
       legend = c('sep','may','other'))