| -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); hdf5load("/home/don/exp007/src/influenza.h5", tidy = TRUE); -d <- influenza.aa.dat; +A <- influenza.aa.dat; + +B <- influenza.faa; # All records for 1918. Based on code from # http://wiki.r-project.org/rwiki/doku.php?id=tips:data-frames:select_observations -p1918 <- d[d$Year == 1918, ] +C <- A[A$Year == 1918, ] # Countries represented in the 1918 dataset. -p1918$Country; +C$Country; # All record with a year value. -years <- d[d$Year != 0, ]; +D <- A[A$Year != 0, ]; + +hist(D$Year); + +B[B$GI == 305182, ] -hist(years$Year); +B[B$"Protein Type" == "HA", ] |
