-rw-r--r-- | analysis/year.R | 16 |
1 files changed, 11 insertions, 5 deletions
diff --git a/analysis/year.R b/analysis/year.R index 85d3da1..b7da3c8 100644 --- a/analysis/year.R +++ b/analysis/year.R | |||
@@ -4,16 +4,22 @@ require(hdf5); | |||
4 | 4 | ||
5 | hdf5load("/home/don/exp007/src/influenza.h5", tidy = TRUE); | 5 | hdf5load("/home/don/exp007/src/influenza.h5", tidy = TRUE); |
6 | 6 | ||
7 | d <- influenza.aa.dat; | 7 | A <- influenza.aa.dat; |
8 | |||
9 | B <- influenza.faa; | ||
8 | 10 | ||
9 | # All records for 1918. Based on code from | 11 | # All records for 1918. Based on code from |
10 | # http://wiki.r-project.org/rwiki/doku.php?id=tips:data-frames:select_observations | 12 | # http://wiki.r-project.org/rwiki/doku.php?id=tips:data-frames:select_observations |
11 | p1918 <- d[d$Year == 1918, ] | 13 | C <- A[A$Year == 1918, ] |
12 | 14 | ||
13 | # Countries represented in the 1918 dataset. | 15 | # Countries represented in the 1918 dataset. |
14 | p1918$Country; | 16 | C$Country; |
15 | 17 | ||
16 | # All record with a year value. | 18 | # All record with a year value. |
17 | years <- d[d$Year != 0, ]; | 19 | D <- A[A$Year != 0, ]; |
20 | |||
21 | hist(D$Year); | ||
22 | |||
23 | B[B$GI == 305182, ] | ||
18 | 24 | ||
19 | hist(years$Year); | 25 | B[B$"Protein Type" == "HA", ] |