tornado <- read.csv("sheet/tornado.csv",
header = TRUE)
#1. (a)
sum(tornado$st == "WA" | tornado$st == "MS" )
## [1] 530
#1. (b)
sum(tornado$st == "WA" & tornado$yr > 2012 )
## [1] 7
#1. (c)
a <- tornado$st == "WA"
b <- tornado$yr == 2012:2014
## Warning in tornado$yr == 2012:2014: 长的对象长度不是短的对象长度的整倍数
tornado[ a & b , ]
## yr mo dy date time tz st stf stn f
## 9150 2014 4 27 2014-04-27 18:30:00 3 WA 53 0 EF-0
## 9747 2014 8 13 2014-08-13 19:30:00 3 WA 53 0 EF-0
#1. (d)
Hi <- tornado$st == "HI"
tornado[ Hi ,c(2, 1, 10, 7) ]
## mo yr f st
## 2766 9 2008 EF-0 HI
## 2842 12 2008 EF-0 HI
## 2878 2 2009 EF-0 HI
## 5371 2 2011 EF-0 HI
## 7359 3 2012 EF-0 HI
## 10093 4 2015 EF-0 HI
#1. (e)
orderdate <- order(tornado$date)
orderdate2 <- order(tornado$time, decreasing = T)
orderdate3 <- orderdate & orderdate2
tail (tornado[ orderdate3, ])
## yr mo dy date time tz st stf stn f
## 11136 2015 12 28 2015-12-28 02:53:00 3 LA 22 0 EF-1
## 11137 2015 12 28 2015-12-28 03:20:00 3 LA 22 0 EF-1
## 11138 2015 12 28 2015-12-28 04:46:00 3 AR 5 0 EF-2
## 11139 2015 12 28 2015-12-28 05:43:00 3 MS 28 0 EF-1
## 11140 2015 12 28 2015-12-28 08:30:00 3 FL 12 0 EF-1
## 11141 2015 12 28 2015-12-28 15:58:00 3 NC 37 0 EF-0
# The answer is NC
#1. (g)
sum(tornado$mo== 1)
## [1] 355
sum(tornado$mo== 2)
## [1] 463
sum(tornado$mo== 3)
## [1] 751
sum(tornado$mo== 4)
## [1] 2171
sum(tornado$mo== 5)
## [1] 2500
sum(tornado$mo== 6)
## [1] 1899
sum(tornado$mo== 7)
## [1] 861
sum(tornado$mo== 8)
## [1] 517
sum(tornado$mo== 9)
## [1] 401
sum(tornado$mo== 10)
## [1] 525
sum(tornado$mo== 11)
## [1] 350
sum(tornado$mo== 12)
## [1] 348
airline <- read.csv("sheet/airline.csv",
header = TRUE)
#2.(a)i
i <- airline$DepTime > 2200 & airline$Dest == "BNA"
newdata1 <- na.omit (i)
airline[ newdata1 , c(9,10,17) ]
## UniqueCarrier FlightNum Origin
## 7020 OH 5793 CVG
## 17399 AA 2431 LGA
## 19476 DH 7330 GSP
#2.(a)ii
i <- airline$DepTime > 2200
d <- airline$Origin == "BNA" | airline$Dest == "MEM"
cc <- i & d
newdata2 <- na.omit (cc)
airline[ newdata2 ,c(9,10,17) ]
## UniqueCarrier FlightNum Origin
## 2392 WN 1447 BNA
## 8168 OO 6837 ORD
## 10404 DL 579 PIT
## 19347 DH 6263 ROA
#2.(a)iii
delay <- airline$ArrDelay > 120
newdelay <- na.omit (delay)
airline[ newdelay ,c(9,10,17) ]
## UniqueCarrier FlightNum Origin
## 27 UA 441 ORD
## 65 UA 471 BUF
## 66 UA 471 ORD
## 78 UA 481 ORD
## 86 UA 489 CMH
## 126 UA 518 ORD
## 149 UA 538 ORD
## 154 UA 542 ORD
## 155 UA 543 BOS
## 158 UA 544 SNA
## 165 UA 548 SLC
## 175 UA 558 RNO
## 263 UA 624 SNA
## 266 UA 627 DCA
## 267 UA 627 ORD
## 295 UA 660 ORD
## 296 UA 661 EWR
## 322 UA 684 ORD
## 325 UA 686 MSP
## 356 UA 710 DEN
## 375 UA 726 DEN
## 386 UA 734 IAH
## 401 UA 748 BUR
## 402 UA 749 PDX
## 460 UA 798 SEA
## 461 UA 798 SFO
## 477 UA 825 MHT
## 487 UA 839 JFK
## 505 UA 862 SFO
## 570 UA 935 LAX
## 597 UA 967 SFO
## 637 UA 1011 SEA
## 673 UA 1065 DEN
## 688 UA 1076 SEA
## 694 UA 1081 MHT
## 695 UA 1081 ORD
## 854 UA 1197 ANC
## 865 UA 1206 SEA
## 871 UA 1212 ORD
## 895 UA 1230 ORD
## 904 UA 1237 RSW
## 934 UA 1262 DEN
## 944 UA 1270 DEN
## 955 UA 1279 IAD
## 969 UA 1289 MHT
## 1761 US 1261 CLT
## 2071 US 1130 BDL
## 3352 WN 644 LAS
## 3429 WN 1247 LAX
## 3949 WN 2920 OAK
## 3985 WN 1662 OAK
## 4406 WN 2073 RDU
## 4420 WN 756 RDU
## 4424 WN 2435 RNO
## 4426 WN 289 RNO
## 4763 WN 146 SMF
## 5007 WN 2574 TUL
## 5181 NW 802 SLC
## 5341 NW 1040 HLN
## 5826 NW 1751 DTW
## 5851 NW 1780 DTW
## 5945 NW 1889 DTW
## 7456 OO 4104 SLC
## 7839 OO 6318 PDX
## 7885 OO 6395 BUR
## 8097 OO 6675 DEN
## 8098 OO 6676 FAT
## 8103 OO 6686 FSD
## 8105 OO 6690 DEN
## 8110 OO 6695 CID
## 8111 OO 6696 LNK
## 8115 OO 6701 DEN
## 8117 OO 6703 MEM
## 8261 XE 2962 BTV
## 8268 XE 2657 EWR
## 8562 XE 2199 EWR
## 8571 XE 2965 CRP
## 9148 XE 2754 IAH
## 9178 XE 2101 CLE
## 9314 XE 2660 EWR
## 9358 TZ 109 MDW
## 9428 TZ 271 EWR
## 9441 TZ 279 EWR
## 9444 TZ 284 MDW
## 9467 TZ 387 MDW
## 9564 TZ 925 MDW
## 9641 UA 79 EWR
## 9694 UA 131 ORD
## 9701 UA 138 SFO
## 9794 UA 242 DEN
## 9819 UA 269 SFO
## 9855 UA 306 LAX
## 9905 UA 347 ICT
## 10468 DL 648 IAH
## 10746 DL 987 SAV
## 11149 DL 1529 PHL
## 11273 DL 1651 SLC
## 11344 DL 1726 SAN
## 11685 DL 2121 ATL
## 11797 DL 2452 PBI
## 12114 EV 4258 ATL
## 12914 FL 1910 GPT
## 13262 HP 524 DFW
## 13270 HP 474 DTW
## 13361 HP 548 LAS
## 14241 MQ 3702 ICT
## 14251 MQ 3716 LBB
## 14289 MQ 3759 DFW
## 14297 MQ 3767 DFW
## 14300 MQ 3770 SGF
## 14302 MQ 3772 AMA
## 14303 MQ 3773 DFW
## 14307 MQ 3778 TYS
## 14322 MQ 3798 HOU
## 14327 MQ 3843 DFW
## 14348 MQ 3922 LNK
## 14366 MQ 3940 ORD
## 14370 MQ 3944 BMI
## 14372 MQ 3946 BMI
## 14382 MQ 3956 EVV
## 14384 MQ 3958 ORD
## 14387 MQ 3961 LIT
## 14417 MQ 4005 ORD
## 14436 MQ 4027 ORD
## 14438 MQ 4031 ORD
## 14445 MQ 4038 ORD
## 14448 MQ 4041 ORD
## 14451 MQ 4045 ORD
## 14462 MQ 4057 ORD
## 14502 MQ 4102 PIA
## 14504 MQ 4104 EVV
## 14526 MQ 4127 SYR
## 14530 MQ 4134 GSP
## 14537 MQ 4142 AZO
## 14550 MQ 4159 ORD
## 14555 MQ 4164 PIT
## 14559 MQ 4168 BUF
## 14568 MQ 4177 ORD
## 14571 MQ 4180 MKE
## 14586 MQ 4201 ORD
## 14587 MQ 4202 LSE
## 14603 MQ 4221 ORD
## 14614 MQ 4233 ROC
## 14615 MQ 4234 ORD
## 14618 MQ 4237 ORD
## 14624 MQ 4243 DSM
## 14637 MQ 4257 CLE
## 14638 MQ 4258 DBQ
## 14974 MQ 4840 LGA
## 14987 MQ 4855 CLE
## 15039 MQ 4956 DCA
## 15164 NW 225 DCA
## 15545 AA 44 SFO
## 15667 AA 310 ONT
## 15681 AA 327 DFW
## 15682 AA 328 MIA
## 15696 AA 349 ORD
## 15700 AA 353 LGA
## 15719 AA 371 DFW
## 15720 AA 371 MIA
## 15737 AA 391 DCA
## 15755 AA 413 ORD
## 15778 AA 438 MCI
## 15780 AA 440 DFW
## 15787 AA 453 ORD
## 15790 AA 457 RDU
## 15793 AA 459 DFW
## 15818 AA 484 ORD
## 15891 AA 559 EWR
## 15915 AA 583 PVD
## 15916 AA 584 TUL
## 15925 AA 593 DFW
## 15936 AA 605 LAX
## 15940 AA 608 ORD
## 15946 AA 614 SEA
## 15947 AA 615 DCA
## 15967 AA 645 BDL
## 15970 AA 649 DFW
## 15980 AA 659 DFW
## 15981 AA 661 EWR
## 15988 AA 671 PHL
## 16001 AA 684 MIA
## 16064 AA 762 STL
## 16109 AA 828 PBI
## 16142 AA 872 DEN
## 16209 AA 1024 DFW
## 16240 AA 1057 DFW
## 16246 AA 1062 DFW
## 16271 AA 1086 MFE
## 16302 AA 1117 DTW
## 16366 AA 1179 DFW
## 16384 AA 1198 ORD
## 16437 AA 1244 SAT
## 16466 AA 1269 SJU
## 16516 AA 1320 TUL
## 16564 AA 1365 MSY
## 16581 AA 1386 DFW
## 16582 AA 1387 DFW
## 16636 AA 1442 ORD
## 16645 AA 1452 BNA
## 16662 AA 1468 FLL
## 16682 AA 1487 ORD
## 16692 AA 1496 MSP
## 16717 AA 1521 EWR
## 16748 AA 1548 IAH
## 16761 AA 1564 DFW
## 16762 AA 1565 DCA
## 16781 AA 1581 OKC
## 16843 AA 1648 SNA
## 16859 AA 1668 MIA
## 16880 AA 1690 DFW
## 16891 AA 1700 MIA
## 16892 AA 1701 DTW
## 16917 AA 1724 MIA
## 16958 AA 1766 ORD
## 16964 AA 1771 ORD
## 17014 AA 1820 DFW
## 17015 AA 1820 ELP
## 17100 AA 1908 MIA
## 17104 AA 1913 DFW
## 17140 AA 1946 LAX
## 17163 AA 1966 SEA
## 17185 AA 1987 ORD
## 17187 AA 1988 ORD
## 17188 AA 1989 ATL
## 17191 AA 1992 SJU
## 17193 AA 1994 SFO
## 17678 AS 228 ANC
## 17756 AS 378 SEA
## 18191 B6 184 SAN
## 18210 B6 32 TPA
## 18230 CO 1486 IAH
## 18349 CO 1429 TPA
## 18487 CO 1092 IAH
## 18584 CO 1942 EWR
## 18798 CO 1178 ORD
## 19092 CO 8 EWR
## 19101 CO 1624 MIA
## 19108 CO 1422 MSY
## 19153 DH 6128 BGR
## 19159 DH 6131 CVG
## 19185 DH 6145 CLT
## 19194 DH 6150 CVG
## 19213 DH 6161 ERI
## 19274 DH 6196 CVG
## 19278 DH 6200 SBN
## 19282 DH 6202 CVG
## 19293 DH 6208 TRI
## 19328 DH 6242 BOS
## 19343 DH 6259 RDU
## 19356 DH 6274 LAN
## 19360 DH 6278 TOL
## 19365 DH 6287 PHL
## 19366 DH 6290 CVG
## 19399 DH 7244 BDL
## 19405 DH 7250 IAD
## 19412 DH 7259 IAD
## 19423 DH 7273 BTV
## 19424 DH 7274 IAD
## 19428 DH 7278 ORD
## 19435 DH 7285 IAD
## 19442 DH 7292 BUF
## 19449 DH 7299 IAD
## 19451 DH 7301 EWR
## 19468 DH 7322 CLT
## 19473 DH 7327 IAD
## 19475 DH 7329 IAD
## 19478 DH 7332 ORD
## 19486 DH 7341 ALB
## 19492 DH 7347 IAD
## 19511 DH 7366 CHS
## 19513 DH 7368 BNA
## 19524 DH 7379 BNA
## 19548 DH 7405 BNA
## 19569 DH 7428 ORF
## 19570 DH 7429 CMH
## 19586 DH 7445 GRR
## 19602 DH 7461 TYS
## 19614 DH 7473 CAE
## 19638 DH 7497 ORD
## 19748 DH 7690 CLE
## 19786 DH 7742 CMH
## 19787 DH 7743 ORD
## 19799 DH 7756 ROA
## 19807 DH 7766 TUL
## 19847 DH 7809 JFK
## 19870 DH 7833 BUF
## 19875 DH 7838 ORF
## 19876 DH 7839 IAD
## 19879 DH 7842 ORF
## 19886 DH 7849 BGM
## 19896 DH 7860 BGM
#2.(a)iv
delay2 <- airline$ArrDelay > 120
arrive <- airline$DepDelay <= 0
newdata4 <- delay2 & arrive
newdata5 <- na.omit (newdata4)
airline[ newdata5 ,c(9,10,17) ]
## UniqueCarrier FlightNum Origin
## 9437 TZ 276 SJC
#2.(b)i
o1 <- airline[order(-airline$DepDelay),c(11,16) ]
head (o1)
## TailNum DepDelay
## 5884 N608NW 1080
## 11220 N978DL 324
## 14475 N514MQ 304
## 15132 N352NW 301
## 16220 N5FKAA 297
## 5368 N924RW 296
#2.(b)ii
o2 <- airline[order(airline$DepDelay),c(11,16) ]
head (o2)
## TailNum DepDelay
## 14202 N817MQ -26
## 6293 N933CA -24
## 5109 N986US -22
## 16874 N4XGAA -21
## 5983 N9346 -20
## 6628 N981CA -20
#2.(b)iii
de <- order(airline$Dest)
airlinebydest <- airline[de,]
de2 <- order(-airlinebydest$ArrDelay)
head (airlinebydest[de2, c(11,18,15)])
## TailNum Dest ArrDelay
## 5884 N608NW DTW 1067
## 11220 N978DL CVG 320
## 14475 N514MQ OKC 313
## 15132 N352NW MEM 306
## 5368 N924RW MSP 301
## 16220 N5FKAA ORD 299