diff --git a/01-1-2-3-child-heart-survival-times/01-1-child-heart-survival-x.Rmd b/01-1-2-3-child-heart-survival-times/01-1-child-heart-survival-x.Rmd index e7db353..28a0f51 100644 --- a/01-1-2-3-child-heart-survival-times/01-1-child-heart-survival-x.Rmd +++ b/01-1-2-3-child-heart-survival-times/01-1-child-heart-survival-x.Rmd @@ -7,12 +7,23 @@ Data are shown in Table 1.1 (page 23) and are contained in [01-1-child-heart-survival-x.csv](01-1-child-heart-survival-x.csv). The data were originally presented in the [NCHDA 2012-15 report](https://nicor4.nicor.org.uk/chd/an_paeds.nsf/vwContent/Analysis%20Documents?Opendocument), but are best seen on [childrensheartsurgery.info](http://childrensheartsurgery.info/). -```{r figure 1-1} +```{r} library(ggplot2) ThirtyDaySurv <-read.csv("01-1-child-heart-survival-x.csv", header=TRUE) # reads data into ThirtyDaySurv data frame -nhosp=length(ThirtyDaySurv$Hospital) +attach(ThirtyDaySurv) +nhosp=length(Hospital) +``` +First in R base graphics +```{r} +par(mar=c(5,15,4,2)) +barplot(ThirtyDaySurvival,names.arg=Hospital,horiz=T,xlim=c(86,100), xpd=F,las=1, xlab="% surviving 30 days") +``` + + +Now in ggplot2 +```{r} p <- ggplot(ThirtyDaySurv, aes(x=reorder(Hospital,nhosp:1), y= ThirtyDaySurvival, fill=Hospital)) # constructs initial plot object, , starting with top row p <- p + geom_bar(stat = "identity") # assigns bar chart-type p <- p + coord_flip(ylim = c(86,100)) # flips to horizontal bars and limits y-axis @@ -26,4 +37,16 @@ p # draws the plot _Figure 1.1 Bar-chart of 30-day survival rates for thirteen hospitals. The choice of the start of the horizontal axis, here 86%, can have a crucial effect on the impression given by the graphic. If the axis starts at 0%, all the hospitals will look indistinguishable, whereas if we started at 95% the differences would look misleadingly dramatic._ -For other ways of displaying and explaining this data, and more recent results, see [childrensheartsurgery.info](http://childrensheartsurgery.info/). \ No newline at end of file +For other ways of displaying and explaining this data, and more recent results, see [childrensheartsurgery.info](http://childrensheartsurgery.info/), which uses a dot-plot. This may be more appropriate. + +```{r figure 1-1} +p <- ggplot(ThirtyDaySurv, aes(x=reorder(Hospital,nhosp:1), y= ThirtyDaySurvival, fill=Hospital)) # constructs initial plot object, , starting with top row +p <- p + geom_dotplot(binaxis="y",stackdir="center",dotsize=2) # assigns dots chart-type +p <- p + coord_flip(ylim = c(86,100)) # flips to horizontal bars and limits y-axis +p <- p + scale_y_continuous(breaks=seq(86, 100, 2)) # assigns breaks every 2 percent +p <- p + theme(legend.position="none") # removes the legend +p <- p + labs(x="", y="% surviving 30 days") # Adds y-axis label +p # draws the plot +``` + +_Figure 1.1 Bar-chart of 30-day survival rates for thirteen hospitals represented as a dot-plot, in which the non-zero axis is less important as, unlike a bar, it is not directly connected to the data-point._ diff --git a/01-1-2-3-child-heart-survival-times/01-1-child-heart-survival-x.html b/01-1-2-3-child-heart-survival-times/01-1-child-heart-survival-x.html index 55f5758..3d995ee 100644 --- a/01-1-2-3-child-heart-survival-times/01-1-child-heart-survival-x.html +++ b/01-1-2-3-child-heart-survival-times/01-1-child-heart-survival-x.html @@ -379,8 +379,14 @@

Figure 1.1: Survival rates following child heart surgery in thirteen hospita
library(ggplot2)
 
 ThirtyDaySurv <-read.csv("01-1-child-heart-survival-x.csv", header=TRUE) # reads data into ThirtyDaySurv data frame
-nhosp=length(ThirtyDaySurv$Hospital)
-p <- ggplot(ThirtyDaySurv, aes(x=reorder(Hospital,nhosp:1), y= ThirtyDaySurvival, fill=Hospital)) # constructs initial plot object, , starting with top row
+attach(ThirtyDaySurv)
+nhosp=length(Hospital)
+

First in R base graphics

+
par(mar=c(5,15,4,2))
+barplot(ThirtyDaySurvival,names.arg=Hospital,horiz=T,xlim=c(86,100), xpd=F,las=1, xlab="% surviving 30 days")
+

+

Now in ggplot2

+
p <- ggplot(ThirtyDaySurv, aes(x=reorder(Hospital,nhosp:1), y= ThirtyDaySurvival, fill=Hospital)) # constructs initial plot object, , starting with top row
 p <- p + geom_bar(stat = "identity") # assigns bar chart-type
 p <- p + coord_flip(ylim = c(86,100)) # flips to horizontal bars and limits y-axis
 p <- p + scale_y_continuous(breaks=seq(86, 100, 2)) # assigns breaks every 2 percent
@@ -390,7 +396,17 @@ 

Figure 1.1: Survival rates following child heart surgery in thirteen hospita p # draws the plot

Figure 1.1 Bar-chart of 30-day survival rates for thirteen hospitals. The choice of the start of the horizontal axis, here 86%, can have a crucial effect on the impression given by the graphic. If the axis starts at 0%, all the hospitals will look indistinguishable, whereas if we started at 95% the differences would look misleadingly dramatic.

-

For other ways of displaying and explaining this data, and more recent results, see childrensheartsurgery.info.

+

For other ways of displaying and explaining this data, and more recent results, see childrensheartsurgery.info, which uses a dot-plot. This may be more appropriate.

+
p <- ggplot(ThirtyDaySurv, aes(x=reorder(Hospital,nhosp:1), y= ThirtyDaySurvival, fill=Hospital)) # constructs initial plot object, , starting with top row
+p <- p + geom_dotplot(binaxis="y",stackdir="center",dotsize=2) # assigns dots chart-type
+p <- p + coord_flip(ylim = c(86,100)) # flips to horizontal bars and limits y-axis
+p <- p + scale_y_continuous(breaks=seq(86, 100, 2)) # assigns breaks every 2 percent
+p <- p + theme(legend.position="none") # removes the legend
+p <- p +  labs(x="", y="% surviving 30 days") # Adds y-axis label  
+p # draws the plot
+
## `stat_bindot()` using `bins = 30`. Pick better value with `binwidth`.
+

+

Figure 1.1 Bar-chart of 30-day survival rates for thirteen hospitals represented as a dot-plot, in which the non-zero axis is less important as, unlike a bar, it is not directly connected to the data-point.

diff --git a/01-1-2-3-child-heart-survival-times/01-3-child-heart-proportions-x.Rmd b/01-1-2-3-child-heart-survival-times/01-3-child-heart-proportions-x.Rmd index fc5ffd5..5005fd9 100644 --- a/01-1-2-3-child-heart-survival-times/01-3-child-heart-proportions-x.Rmd +++ b/01-1-2-3-child-heart-survival-times/01-3-child-heart-proportions-x.Rmd @@ -16,6 +16,18 @@ library(ggplot2) df <- read.csv("01-1-child-heart-survival-x.csv", header=TRUE) # reads csv into dataframe, df df$Percentage = 100*df$Operations/sum(df$Operations) df$Pos= rank(df$Percentage) +``` +First in R base graphics + +```{r} +par(mar=c(5,15,4,2)) +barplot(df$Percentage,names.arg=df$Hospital,horiz=T, xpd=F,las=1, xlab="Percentage of all operations in 2012-15 \nthat are carried out in each hospital") + +``` + +Now in ggplot2 + +```{r} bp <- ggplot(df, aes(x=reorder(Hospital,-Pos), y=Percentage, fill=Hospital)) #sets initial plot object from the dataframe for Hospitals, reordered by Percentage (descending) as the y-values, colour-filled by Hospital bp <- bp + geom_bar(stat = "identity") + labs(x="Hospital") # makes the plot a bar-chart bp <- bp + coord_flip() # makes it an horizontal bar chart diff --git a/01-1-2-3-child-heart-survival-times/01-3-child-heart-proportions-x.html b/01-1-2-3-child-heart-survival-times/01-3-child-heart-proportions-x.html index 4cb010f..80c1812 100644 --- a/01-1-2-3-child-heart-survival-times/01-3-child-heart-proportions-x.html +++ b/01-1-2-3-child-heart-survival-times/01-3-child-heart-proportions-x.html @@ -380,8 +380,13 @@

Figure 1.3: Percentage of all child heart surgery being carried out in each df <- read.csv("01-1-child-heart-survival-x.csv", header=TRUE) # reads csv into dataframe, df df$Percentage = 100*df$Operations/sum(df$Operations) -df$Pos= rank(df$Percentage) -bp <- ggplot(df, aes(x=reorder(Hospital,-Pos), y=Percentage, fill=Hospital)) #sets initial plot object from the dataframe for Hospitals, reordered by Percentage (descending) as the y-values, colour-filled by Hospital +df$Pos= rank(df$Percentage) +

First in R base graphics

+
par(mar=c(5,15,4,2))
+barplot(df$Percentage,names.arg=df$Hospital,horiz=T, xpd=F,las=1, xlab="Percentage of all operations in 2012-15 \nthat are carried out in each hospital")
+

+

Now in ggplot2

+
bp <- ggplot(df, aes(x=reorder(Hospital,-Pos), y=Percentage, fill=Hospital)) #sets initial plot object from the dataframe for Hospitals, reordered by Percentage (descending) as the y-values, colour-filled by Hospital
 bp <- bp + geom_bar(stat = "identity") + labs(x="Hospital") # makes the plot a bar-chart
 bp <- bp + coord_flip() # makes it an horizontal bar chart
 bp <- bp + scale_y_continuous(breaks=seq(2,16,2)) # breaks every two-count
diff --git a/02-4-reported-partners/02-4-sexual-partners-x.Rmd b/02-4-reported-partners/02-4-sexual-partners-x.Rmd
index 2bf9f95..3af8a30 100644
--- a/02-4-reported-partners/02-4-sexual-partners-x.Rmd
+++ b/02-4-reported-partners/02-4-sexual-partners-x.Rmd
@@ -49,4 +49,4 @@ p # draw the plot
 
 ```
 
-Figure 2.4 Data provided by Natsal-3 based on interviews between 2010 and 2012. The series have been truncated at 50 for reasons of space - the totals go up to 500 for both men and women. Note the clear use of round numbers for ten or more partners, and the tendency for men to report more partners than women.
\ No newline at end of file
+_Figure 2.4 Data provided by Natsal-3 based on interviews between 2010 and 2012. The series have been truncated at 50 for reasons of space - the totals go up to 500 for both men and women. Note the clear use of round numbers for ten or more partners, and the tendency for men to report more partners than women._
\ No newline at end of file
diff --git a/02-4-reported-partners/02-4-sexual-partners-x.html b/02-4-reported-partners/02-4-sexual-partners-x.html
index 97f32b3..61e0690 100644
--- a/02-4-reported-partners/02-4-sexual-partners-x.html
+++ b/02-4-reported-partners/02-4-sexual-partners-x.html
@@ -419,7 +419,7 @@ 

Figure 2.4 (page 55) Bar chart

p # draw the plot
## Warning: Removed 174 rows containing missing values (geom_bar).

-

Figure 2.4 Data provided by Natsal-3 based on interviews between 2010 and 2012. The series have been truncated at 50 for reasons of space - the totals go up to 500 for both men and women. Note the clear use of round numbers for ten or more partners, and the tendency for men to report more partners than women.

+

Figure 2.4 Data provided by Natsal-3 based on interviews between 2010 and 2012. The series have been truncated at 50 for reasons of space - the totals go up to 500 for both men and women. Note the clear use of round numbers for ten or more partners, and the tendency for men to report more partners than women.

diff --git a/02-5-survival-vs-numbers/02-5-child-heart-surgery-x.Rmd b/02-5-survival-vs-numbers/02-5-child-heart-surgery-x.Rmd index 1da7a8f..e6813b7 100755 --- a/02-5-survival-vs-numbers/02-5-child-heart-surgery-x.Rmd +++ b/02-5-survival-vs-numbers/02-5-child-heart-surgery-x.Rmd @@ -11,19 +11,32 @@ Data from 2012-2014 were shown in Table 1.1 (page 23) and are contained in [02-5 ### Figure 2.5 (page 57) Scatterplots -```{r} +Quickly in qplot +```{r} library(ggplot2) # (a) Survival in under-1s, 1991-1995 child.1991 <- read.csv("02-5-child-heart-surgery-1991-x.csv") # read data into dataframe +attach(child.1991) +qplot(Operations,100*Survivors/Operations,xlim = c(0,700),ylim=c(70,100),ylab = "% 30-day survival", label = Hospital, geom=c("point", "text"),hjust=1, vjust=-1,size=2, main = "(a) Survival in under-1s, 1991-1995") + theme(legend.position="none") +# (a) Survival in under-1s, 1991-1995 +#(b) Survival for all children, 2012-2015 +child.2012 <- read.csv("02-5-child-heart-surgery-2012-x.csv") # read data into dataframe all +attach(child.2012) +qplot(Operations,100*Survivors/Operations,xlim = c(0,2000),ylim=c(95,100),ylab = "% 30-day survival", label = Hospital, geom=c("point", "text"),hjust=1, vjust=-1,size=2, main = "(b) Survival for all children, 2012-2015") + theme(legend.position="none") + + +``` + +```{r} p <- ggplot(child.1991, aes(x=Operations, y=100*Survivors/Operations, col=Hospital)) #defines plot axis data fields and colour legend data field p <- p + geom_point(aes(size=1.5)) # defines scatter-type plot p <- p + expand_limits(x = c(0,700),y=c(70,100)) p <- p + scale_size_continuous(name = "Size", guide = FALSE) # turns off otherwise added size legend p <- p + labs(x="Number of operations", y = "% 30-day survival", title="(a) Survival in under-1s, 1991-1995") # Adds title, subtitle, and caption -p <- p + theme(plot.caption=element_text(hjust = 0.5)) # centre justifies the caption p + #(b) Survival for all children, 2012-2015 child.2012 <- read.csv("02-5-child-heart-surgery-2012-x.csv") # read data into dataframe all q <- ggplot(child.2012, aes(x=Operations, y=100*Survivors/Operations,col=Hospital)) #defines plot axis data fields and colour legend data field @@ -35,8 +48,8 @@ q ``` -Figure 2.5 -Scatter-plots of survival rates against number of operations in child heart surgery. For (a) 1991-1995, the Pearson correlation is 0.59 and the rank correlation is 0.85, for (b) 2012-2015, the Pearson correlation is 0.17 and the rank correlation is -0:03. +_Figure 2.5 +Scatter-plots of survival rates against number of operations in child heart surgery. For (a) 1991-1995, the Pearson correlation is 0.59 and the rank correlation is 0.85, for (b) 2012-2015, the Pearson correlation is 0.17 and the rank correlation is -0:03._ ### Correlations in (a) 1991-1995 data diff --git a/02-5-survival-vs-numbers/02-5-child-heart-surgery-x.html b/02-5-survival-vs-numbers/02-5-child-heart-surgery-x.html index 7684e4e..214c53c 100644 --- a/02-5-survival-vs-numbers/02-5-child-heart-surgery-x.html +++ b/02-5-survival-vs-numbers/02-5-child-heart-surgery-x.html @@ -377,16 +377,28 @@

Art of Statistics: Figure 2.5 (page 57) Child heart

Data from 2012-2014 were shown in Table 1.1 (page 23) and are contained in 02-5-child-heart-surgery-2012-x.csv. The data were originally presented in the NCHDA 2012-15 report, but are best seen on childrensheartsurgery.info.

Figure 2.5 (page 57) Scatterplots

+

Quickly in qplot

library(ggplot2)
 # (a) Survival in under-1s, 1991-1995
 child.1991 <- read.csv("02-5-child-heart-surgery-1991-x.csv") # read data into dataframe 
-
-p <- ggplot(child.1991, aes(x=Operations, y=100*Survivors/Operations, col=Hospital)) #defines plot axis data fields and colour legend data field
+attach(child.1991)
+qplot(Operations,100*Survivors/Operations,xlim = c(0,700),ylim=c(70,100),ylab = "% 30-day survival", label = Hospital, geom=c("point", "text"),hjust=1, vjust=-1,size=2, main = "(a) Survival in under-1s, 1991-1995") + theme(legend.position="none")
+

+
# (a) Survival in under-1s, 1991-1995
+#(b) Survival for all children, 2012-2015
+child.2012 <- read.csv("02-5-child-heart-surgery-2012-x.csv") # read data into dataframe all
+attach(child.2012)
+
## The following objects are masked from child.1991:
+## 
+##     Deaths, Hospital, Operations, PercentageDying, Survivors,
+##     ThirtyDaySurvival
+
qplot(Operations,100*Survivors/Operations,xlim = c(0,2000),ylim=c(95,100),ylab = "% 30-day survival", label = Hospital, geom=c("point", "text"),hjust=1, vjust=-1,size=2, main = "(b) Survival for all children, 2012-2015") + theme(legend.position="none")
+

+
p <- ggplot(child.1991, aes(x=Operations, y=100*Survivors/Operations, col=Hospital)) #defines plot axis data fields and colour legend data field
 p <- p + geom_point(aes(size=1.5)) # defines scatter-type plot
 p <- p + expand_limits(x = c(0,700),y=c(70,100))
 p <- p  + scale_size_continuous(name = "Size", guide = FALSE)  # turns off otherwise added size legend
 p <- p +  labs(x="Number of operations", y = "% 30-day survival", title="(a) Survival in under-1s, 1991-1995") # Adds title, subtitle, and caption
-p <- p + theme(plot.caption=element_text(hjust = 0.5)) # centre justifies the caption
 p

#(b) Survival for all children, 2012-2015
@@ -398,12 +410,20 @@ 

Figure 2.5 (page 57) Scatterplots

q <- q + labs(x="Number of operations", y = "% 30-day survival", title="(b) Survival for all children, 2012-2015") # Adds title, subtitle q

-

Figure 2.5 Scatter-plots of survival rates against number of operations in child heart surgery. For (a) 1991-1995, the Pearson correlation is 0.59 and the rank correlation is 0.85, for (b) 2012-2015, the Pearson correlation is 0.17 and the rank correlation is -0:03.

+

Figure 2.5 Scatter-plots of survival rates against number of operations in child heart surgery. For (a) 1991-1995, the Pearson correlation is 0.59 and the rank correlation is 0.85, for (b) 2012-2015, the Pearson correlation is 0.17 and the rank correlation is -0:03.

Correlations in (a) 1991-1995 data

-
attach(child.1991)
-cor.test(Operations, 100*Survivors/Operations,method="pearson")
+
attach(child.1991)
+
## The following objects are masked from child.2012:
+## 
+##     Deaths, Hospital, Operations, PercentageDying, Survivors,
+##     ThirtyDaySurvival
+
## The following objects are masked from child.1991 (pos = 4):
+## 
+##     Deaths, Hospital, Operations, PercentageDying, Survivors,
+##     ThirtyDaySurvival
+
cor.test(Operations, 100*Survivors/Operations,method="pearson")
## 
 ##  Pearson's product-moment correlation
 ## 
@@ -454,7 +474,15 @@ 

Correlations in (a) 1991-1995 data, without Bristol

Correlations in (b) 2012-2014 data

attach(child.2012)
-
## The following objects are masked from child.1991:
+
## The following objects are masked from child.1991 (pos = 3):
+## 
+##     Deaths, Hospital, Operations, PercentageDying, Survivors,
+##     ThirtyDaySurvival
+
## The following objects are masked from child.2012 (pos = 4):
+## 
+##     Deaths, Hospital, Operations, PercentageDying, Survivors,
+##     ThirtyDaySurvival
+
## The following objects are masked from child.1991 (pos = 5):
 ## 
 ##     Deaths, Hospital, Operations, PercentageDying, Survivors,
 ##     ThirtyDaySurvival
diff --git a/02-6-zero-correlations/02-6-zero-correlations-x.Rmd b/02-6-zero-correlations/02-6-zero-correlations-x.Rmd index 2cec963..c9a4310 100644 --- a/02-6-zero-correlations/02-6-zero-correlations-x.Rmd +++ b/02-6-zero-correlations/02-6-zero-correlations-x.Rmd @@ -26,4 +26,4 @@ p ``` -Figure 2.6 Two sets of (fictitious) data-points for which the Pearson correlation coefficients are both 0. This clearly does not mean there is no relationship between the two variables being plotted. +_Figure 2.6 Two sets of (fictitious) data-points for which the Pearson correlation coefficients are both 0. This clearly does not mean there is no relationship between the two variables being plotted._ diff --git a/02-6-zero-correlations/02-6-zero-correlations-x.html b/02-6-zero-correlations/02-6-zero-correlations-x.html index badfae1..626ab4f 100644 --- a/02-6-zero-correlations/02-6-zero-correlations-x.html +++ b/02-6-zero-correlations/02-6-zero-correlations-x.html @@ -393,7 +393,7 @@

Art of Statistics: Figure 2.6 (page 59) Pearson Cor p

-

Figure 2.6 Two sets of (fictitious) data-points for which the Pearson correlation coefficients are both 0. This clearly does not mean there is no relationship between the two variables being plotted.

+

Figure 2.6 Two sets of (fictitious) data-points for which the Pearson correlation coefficients are both 0. This clearly does not mean there is no relationship between the two variables being plotted.

diff --git a/02-7-population-trends/02-7-pop-past-future million in 51.csv b/02-7-population-trends/02-7-pop-past-future million in 51.csv new file mode 100644 index 0000000..7d58858 --- /dev/null +++ b/02-7-population-trends/02-7-pop-past-future million in 51.csv @@ -0,0 +1,67 @@ +Year,China,India,United States of America,Russian Federation,Japan,Indonesia,Germany,Brazil,United Kingdom,Northern Africa,Italy,France,Bangladesh,Nigeria,Pakistan,Ukraine,Mexico,Spain,Viet Nam,Poland,Turkey,Thailand,Egypt,Republic of Korea,Philippines,Ethiopia,Iran (Islamic Republic of),Argentina,Myanmar,Caribbean,Romania,Southern Africa,Canada,South Africa,Oceania,Colombia,Democratic Republic of the Congo,Australia/New Zealand,Dem. People's Republic of Korea,Netherlands,Hungary,Morocco,Algeria,Czechia,Belgium,Nepal,Portugal,Australia,Sri Lanka,Peru,"China, Taiwan Province of China",United Republic of Tanzania,Afghanistan,Greece,Belarus,Bulgaria,Sweden,Austria,Serbia,Kazakhstan,Uzbekistan,Chile,Malaysia,Mozambique,Kenya,Cuba,Iraq,Sudan,Venezuela (Bolivarian Republic of),Uganda,Ghana,Mali,Switzerland,Angola,Cambodia,Yemen,Cameroon,Burkina Faso,Denmark,Madagascar,Finland,Croatia,Tunisia,Ecuador,Georgia,Slovakia,Syrian Arab Republic,Norway,Haiti,Guatemala,Saudi Arabia,Guinea,Bolivia (Plurinational State of),Malawi,Azerbaijan,Ireland,Zimbabwe,Bosnia and Herzegovina,Cote d'Ivoire,Niger,South Sudan,Lithuania,Senegal,Chad,Dominican Republic,Republic of Moldova,Zambia,Burundi,Somalia,Uruguay,Benin,Rwanda,El Salvador,Puerto Rico,Melanesia,"China, Hong Kong SAR",Sierra Leone,New Zealand,Latvia,Kyrgyzstan,Lao People's Democratic Republic,Papua New Guinea,Honduras,Tajikistan,Paraguay,Slovenia,Jamaica,Togo,Armenia,Lebanon,Israel,Central African Republic,Nicaragua,TFYR Macedonia,Albania,Turkmenistan,Eritrea,Libya,Estonia,Singapore +1950,554419.28,376325.2,158804.4,102798.66,82802.08,69543.32,69966.24,53974.73,50616.01,49048.94,46598.6,41879.61,37894.68,37859.74,37542.38,37297.65,28012.56,28069.74,24809.91,24824.02,21408.4,20710.36,20713.33,19211.39,18580.49,18128.03,17119.27,17150.34,17151.51,17075.67,16236.3,15533.18,13733.4,13628.43,12647.78,12340.9,12183.66,10085.34,10549.47,10042.04,9337.72,8985.99,8872.25,8902.62,8637.52,8483.33,8416.97,8177.34,7971.1,7727.73,7623.36,7649.76,7752.12,7668.79,7745,7251,7009.91,6936.44,6732.26,6702.99,6264.06,6188.28,6109.92,6152.24,6076.76,5919.99,5719.19,5733.79,5481.98,5158.19,4980.88,4708.43,4668.09,4548.02,4432.71,4402.32,4307.02,4284.46,4268.27,4083.55,4008.29,3850.3,3605.31,3470.16,3527.01,3436.57,3413.33,3265.27,3221.28,3114.96,3121.34,3093.65,3089.65,2953.87,2927.92,2913.09,2746.85,2661.3,2630.13,2559.7,2582.93,2567.4,2486.64,2502.32,2364.65,2341,2310.44,2308.92,2264.08,2238.51,2255.23,2186.19,2199.9,2218,2165.34,1974,2041.2,1908,1926.57,1740,1682.92,1674.04,1546.72,1531.5,1473.25,1473.09,1402.9,1395.46,1353.51,1334.62,1257.97,1326.65,1294.99,1254.45,1263.17,1210.99,1142.15,1124.52,1101,1022.1 +1951,569611.08,382245.3,160905.04,104306.14,84315.6,70869.61,70299.36,55619.22,50621.42,50237.87,47021.9,42070.92,38705.08,38423.98,37975.28,37815.82,28836.11,28243.64,25364.45,25283.44,21951.21,21263.2,21215.01,19459.43,19246.61,18466.97,17516.58,17507.13,17470.06,17395.84,16512.66,15897.01,14098.68,13953.84,12971.3,12699.95,12428.83,10360.54,10248.5,10167.53,9480.57,9242.8,9039.91,9007.01,8680.17,8658.31,8432.97,8412.6,8108.43,7924.21,7893.06,7847.31,7839.51,7743.25,7717.07,7310.51,7073.23,6931.34,6853.29,6831.32,6476.16,6327.09,6271.49,6248.72,6240.16,6051.29,5901.51,5884.66,5735.17,5308.67,5068.44,4759.04,4726.29,4659.33,4537.54,4473.75,4382.93,4324.49,4308.92,4168.39,4050.49,3886.22,3696.45,3561.15,3555.87,3510.93,3498.53,3300.75,3275.37,3212.33,3198.55,3141.88,3140.34,3008.07,2987.27,2921.11,2829.75,2710.65,2695.12,2638.14,2601.65,2569.61,2545.16,2543.9,2438.41,2381.44,2366.63,2359.12,2307.67,2261.34,2258.48,2250.68,2236.58,2209.7,2203.57,2063.32,2062.44,1947.93,1933.99,1765.56,1723.16,1702.37,1588.56,1584.99,1511.77,1481.63,1436.94,1425.49,1381.87,1358.22,1354.1,1339.84,1331.75,1290.41,1287.5,1231.83,1159.8,1142.97,1113.25,1067.79 +1952,582029.3,388538.62,163324.85,105969.2,85658.68,72309.9,70619.5,57307.48,50686.06,51487.6,47411.94,42365.58,39491.06,39035.15,38483.89,38360.51,29711.63,28442.37,25976.84,25764.93,22508.95,21837.99,21746.81,19829.39,19945.66,18819.7,17933.36,17866.1,17796.4,17734.61,16764.9,16273.35,14482.06,14291.42,13275.57,13064.69,12680.93,10613.77,10049.03,10297.98,9597.9,9528.26,9216.4,9110.14,8727.99,8824.46,8477.26,8620.97,8256.56,8128.52,8163.93,8055.91,7934.98,7814.59,7710.09,7371.97,7127.52,6931.2,6963.4,7041.8,6678.99,6467.64,6450.08,6352.58,6412.36,6180.03,6065.01,6041.33,5989.55,5455.54,5191.22,4811.2,4788.96,4774.56,4656.28,4545.74,4460.32,4367.13,4346.18,4256.77,4096.55,3922.67,3773.36,3655.01,3588.62,3587.66,3591.01,3334.48,3331.58,3311.54,3283.3,3185.91,3193.29,3065.03,3060.01,2924.05,2917.37,2764.71,2770.36,2716.82,2625.13,2578.05,2606.42,2588.91,2517.26,2432.1,2427.99,2403.51,2351.93,2286.27,2264.87,2313.31,2278.31,2212.1,2244.83,2160.1,2084.94,1992.8,1945.96,1795.13,1764.01,1731.23,1631.73,1637.95,1550.76,1491.67,1468.02,1441.91,1418.28,1392.73,1451.49,1353.47,1371.52,1320.98,1316.09,1258.39,1179.54,1163.89,1125.24,1119.93 +1953,592567.86,395160.09,165988.19,107729.37,86869.43,73866.52,70929.26,59030.86,50797.27,52800.14,47765.98,42724.22,40300.36,39685.76,39063.87,38916.26,30634.22,28658.67,26646.17,26261.07,23081.58,22436.53,22309.1,20307.9,20670.54,19184.29,18368.9,18224.15,18135.96,18084.02,17005.72,16662.28,14884.17,14641.02,13572.46,13438.65,12944.39,10857.13,9957.24,10432.3,9692.84,9836.53,9405.44,9205.96,8779.27,8984.31,8537.75,8816.9,8417.12,8340.38,8440.26,8274.75,8038.6,7883,7723.23,7434.61,7175.67,6935.14,7063.16,7315.65,6879.76,6610.73,6640.04,6463.09,6593.11,6304.52,6216.13,6204.1,6244.65,5601.24,5339.18,4864.76,4851.98,4890.47,4783.09,4618.55,4539.44,4413.2,4380.3,4348.72,4144.3,3959.12,3838.51,3752.14,3625.6,3664.92,3690.74,3367.06,3390.05,3413.17,3372.62,3228.38,3248.46,3124.73,3142.81,2921.36,3008.93,2821.31,2852.5,2795.81,2653.1,2591.37,2670.59,2636.26,2600.8,2491.06,2493.51,2445.42,2397.16,2313.21,2274.48,2378.45,2324.56,2221.02,2288.86,2264.11,2108.43,2040.23,1961.66,1828.2,1805.52,1761.05,1676.5,1690.33,1590.5,1502.84,1495.67,1452.09,1461.42,1434.85,1546.24,1367.8,1414.28,1347.46,1348.1,1289.83,1201.12,1187.73,1136.9,1177.58 +1954,601971.22,402077.03,168827.61,109537.79,87980.32,75539.74,71233.13,60784.89,50945.4,54176.15,48084.81,43117.85,41169.07,40370.37,39711.64,39471.67,31599.8,28887.35,27370.7,26764.86,23668.91,23060.43,22901.6,20880.43,21416.1,19560.02,18822.6,18579.04,18492.48,18438.79,17243.82,17064.04,15303.89,15002.73,13870.9,13824.73,13222.52,11099.67,9972.44,10569.69,9768.62,10162.68,9609.51,9290.49,8832.7,9140.06,8604.73,9011.19,8591.03,8559.72,8725.47,8503.22,8150.45,7948.53,7755.06,7497.87,7220.06,6942.62,7153.36,7636.83,7084.38,6756.97,6837.36,6579.63,6782.29,6424.17,6360.44,6373.31,6500.64,5748.21,5504.25,4919.55,4912.89,5004.76,4913.7,4692.42,4620.63,4463.17,4411.66,4444.29,4191.98,3995.14,3894.39,3852.86,3667.06,3741.21,3797.62,3398.91,3450.85,3517.65,3464.55,3271.31,3305.84,3187.17,3233.17,2913.05,3103.9,2878.85,2939.22,2875.27,2685.33,2608.52,2737.83,2685.2,2688.7,2556.59,2562.48,2487.39,2443.63,2342.04,2287.39,2448.83,2374.97,2233.41,2335.45,2374.54,2132.73,2088.48,1980.4,1864.5,1847.78,1792.14,1723.01,1742.35,1631.19,1514.78,1519.79,1461.52,1510.17,1481.96,1635.78,1383.07,1459.85,1370.98,1382.88,1325.43,1224.38,1214.83,1148.21,1239.68 +1955,610834.4,409269.06,171783.84,111355.22,89018.26,77327.79,71537.54,62569.2,51123.71,55615.1,48372.84,43527.76,42121.52,41085.56,40424.3,40019.45,32605.18,29125.35,28147.79,27269.75,24270.59,23711.05,23523.38,21531.33,22179.1,19947.29,19294,18929.35,18867.97,18796.5,17483.94,17479.12,15738.36,15376.83,14176.71,14225.2,13517.51,11347.3,10086.99,10709.63,9828.44,10502.67,9829.72,9361.76,8887.27,9293.49,8670.88,9210.91,8778.45,8786.76,9022.18,8740.99,8270.58,8011.12,7803.54,7561.41,7262.62,6953.41,7234.92,7992,7297.3,6906.84,7039.72,6701.75,6979.93,6539.47,6502.66,6549.3,6758.37,5898.84,5680.41,4975.4,4970.81,5116.07,5045.29,4767.59,4704.23,4517.16,4440.78,4543.54,4238.28,4030.41,3943.53,3957.49,3713.12,3815.4,3911.5,3430.34,3514.07,3625.3,3558.16,3316.14,3365.36,3252.33,3329.39,2899.78,3202.02,2936.31,3029.3,2955.54,2721.54,2628.75,2808.23,2735.21,2780.65,2627.01,2634.46,2531.18,2491.52,2372.57,2303.59,2525.53,2429.33,2247.42,2384.38,2490.1,2157.75,2136.39,2001.57,1904.06,1890.85,1824.71,1771.35,1794.41,1673.01,1527.21,1540.62,1473.84,1563.51,1532.11,1718.96,1399.49,1508.03,1392.45,1419.97,1364.63,1249.19,1245.36,1159.23,1305.09 +1956,619597.55,416728.77,174805.94,113152.38,90004.31,79226.79,71850.98,64387.82,51328.66,57115.14,48638.3,43946.21,43169.93,41830.03,41199.7,40556.53,33647.98,29371.82,28973.87,27769.54,24886.11,24389.52,24172.87,22243.63,22958.36,20347.81,19782.7,19274.53,19262.65,19157.5,17726.63,17908.23,16182.88,15763.82,14492.82,14641.72,13830.44,11602.73,10285.94,10851.92,9875.55,10853.33,10065.82,9420.02,8942.4,9446.01,8731.33,9419.27,8978.77,9021.96,9331.9,8987.98,8399.03,8070.64,7865.98,7625.11,7304.78,6967.63,7308.93,8370.63,7521.59,7060.66,7246.54,6829.09,7186.21,6652.09,6646.64,6732.44,7019.41,6055.52,5863.69,5032.12,5026.67,5224,5176.66,4844.33,4790.71,4574.94,4468.35,4646.58,4282.38,4064.77,3988.48,4066.24,3763.76,3886.75,4032.07,3461.51,3579.73,3736.24,3653.56,3363.66,3427,3320.22,3430.66,2882.85,3303.27,2993.27,3122.63,3037.06,2761.54,2651.63,2881.81,2786.05,2876.43,2700.79,2709.36,2577.74,2540.99,2404.5,2323.06,2607.95,2487.59,2262.42,2435.63,2608.94,2183.52,2183.46,2024.68,1947.22,1934.81,1858.89,1821.54,1847.25,1716.11,1539.79,1558.82,1490.73,1620.59,1584.03,1796.01,1417.22,1558.47,1412.61,1459.09,1406.98,1275.53,1279.41,1170.1,1372.49 +1957,628551.41,424461.11,177851.16,114909.59,90953.95,81230.75,72183.92,66248.95,51559.67,58673.2,48892.97,44375.74,44314.82,42604.51,42036.39,41083.98,34726.72,29628.01,29844.53,28258.52,25514.81,25096.79,24847.82,22999.24,23754.59,20764.45,20288.44,19614.85,19675.04,19524.81,17968.47,18352.23,16630.98,16164.41,14819.08,15075.31,14161.24,11865.65,10547.39,10996.64,9913.23,11212.45,10316.29,9467.55,8997.84,9598.63,8783.55,9635.83,9190.66,9266.08,9655.24,9244.35,8535.81,8126.81,7939.12,7689.05,7347.48,6985.72,7376.69,8764.92,7758.93,7218.59,7458.9,6961.49,7401.42,6764.79,6795.4,6923.12,7285.97,6220.63,6052.16,5089.51,5083,5329.1,5308.14,4922.89,4880.59,4635.96,4495.22,4753.52,4323.92,4098.12,4031.79,4179.32,3818.85,3954.86,4158.98,3492.45,3647.81,3850.43,3751.88,3414.05,3490.68,3390.8,3537.03,2864.14,3407.91,3049.92,3220.16,3120.44,2805.11,2677,2958.54,2837.76,2975.81,2776.55,2787.34,2627.23,2592.13,2437.5,2345.72,2693.82,2549.83,2278.96,2489.14,2728.67,2210.13,2229.83,2049.32,1994.56,1979.73,1894.68,1873.5,1901.78,1760.58,1552.25,1575.36,1512.03,1680.71,1637.13,1868.58,1436.39,1610.82,1432.02,1500.15,1452.17,1303.42,1316.91,1181.02,1440.46 +1958,637853.61,432481.66,180884.46,116615.77,91877.7,83332.05,72548.04,68164.04,51818.62,60285.15,49151.47,44827.61,45546.48,43411.48,42933.59,41606.29,35840.65,29896.93,30754.61,28731.44,26155.91,25833.55,25545.48,23779.37,24570.18,21201.03,20811.06,19951.32,20102.25,19903.82,18202.71,18812.2,17074.87,16579.48,15152.85,15526.39,14508.91,12133.08,10843.98,11144.07,9944.7,11578.56,10578.45,9508.37,9053.64,9752.05,8827.14,9856.95,9412.16,9520.05,9991.99,9510.51,8680.95,8179.27,8019.18,7753.53,7391.25,7008.42,7439.61,9169.45,8009.77,7380.64,7679.27,7098.86,7626,6881.21,6951.16,7121.74,7560.82,6396.46,6245.65,5147.35,5143.67,5432.72,5441.32,5003.53,4974.43,4699.34,4522.33,4864.52,4362.91,4130.49,4075.96,4296.83,3878.13,4019.64,4291.8,3523.05,3718.28,3967.75,3855.08,3466.99,3556.36,3464.05,3649.19,2846.06,3516.4,3106.89,3323.76,3206.36,2852.02,2704.99,3038.33,2890.64,3078.63,2852.99,2868.84,2679.14,2644.99,2471.15,2371.45,2779.55,2616.28,2298.57,2544.92,2846.52,2237.77,2276.13,2075.17,2046.92,2025.67,1931.99,1927.1,1959.14,1806.52,1564.36,1591.51,1535.97,1743.3,1691.36,1939.45,1457.11,1664.61,1451.02,1543.22,1500.02,1332.95,1357.73,1192.25,1507.51 +1959,647555.96,440813.9,183877.87,118266.78,92782.22,85521.98,72955.24,70146.51,52108.97,61946.19,49429.95,45319.11,46847.45,44254.69,43890.98,42130.32,36989.69,30182.78,31698.43,29183.69,26808.62,26600.28,26262.7,24565.35,25408.6,21662,21350.52,20285.53,20540.4,20301.22,18420.45,19289.33,17506.05,17009.99,15489.65,15994.89,14871.77,12400.15,11145.15,11294.69,9973.06,11950.83,10848.97,9547.67,9110.07,9906.69,8863.39,10076.74,9640.84,9784.91,10341.22,9787.03,8834.45,8227.65,8102.2,7818.91,7436.26,7036.55,7499.17,9580.8,8273.56,7546.71,7911.04,7241.23,7860.48,7005.48,7115.51,7328.72,7846.95,6585.11,6445.28,5205.45,5213.4,5536.78,5578.68,5086.54,5072.8,4764.08,4550.68,4979.75,4399.67,4161.96,4123.44,4418.89,3941.23,4081.32,4430.1,3553.13,3791.1,4087.92,3965.82,3521.7,3623.97,3539.97,3768.33,2831.25,3629.36,3165.19,3435.93,3295.58,2902.08,2735.87,3121.09,2945.1,3184.76,2928.99,2954.48,2732.44,2699.6,2505.02,2400.13,2860.69,2687.23,2323.45,2603.08,2959.53,2266.69,2323.41,2102,2105.18,2072.7,1970.69,1982.19,2020.55,1853.94,1575.87,1608.72,1559.67,1807.92,1747.16,2012.15,1479.47,1719.37,1469.84,1588.48,1550.41,1364.28,1401.62,1204.11,1572.14 +1960,657686.14,449480.61,186808.23,119860.29,93673.62,87792.52,73414.43,72207.55,52433.16,63651.89,49742.2,45865.4,48199.75,45137.81,44908.29,42662.15,38174.11,30489.24,32670.63,29611.56,27472.33,27397.18,26996.53,25340.92,26273.03,22151.28,21906.9,20619.08,20986.12,20722.93,18613.94,19784.72,17917.24,17456.86,15825.22,16480.38,15248.25,12662.3,11424.18,11448.82,10001.02,12328.53,11124.89,9590.19,9167.37,10063.01,8893.97,10289.71,9874.48,10061.52,10701.71,10074.51,8996.35,8271.79,8184.63,7885.53,7482.54,7070.77,7556.73,9996,8549.49,7716.63,8157.11,7388.7,8105.44,7141.14,7289.76,7544.49,8146.85,6788.21,6652.29,5263.73,5296.12,5643.18,5722.37,5172.14,5176.27,4829.29,4581.1,5099.37,4434.59,4192.63,4176.27,4545.55,4007.71,4140.13,4573.51,3582.51,3866.16,4210.75,4086.54,3577.41,3693.45,3618.6,3895.4,2822.11,3747.37,3225.67,3558.99,3388.76,2955.15,2769.87,3206.75,3001.59,3294.04,3003.56,3044.85,2786.11,2755.95,2538.65,2431.62,2933.43,2762.9,2355.57,2663.68,3065.19,2297.11,2372.59,2129.59,2170.09,2120.9,2010.68,2038.64,2087.04,1902.88,1586.64,1628.25,1580.51,1874.12,1804.93,2089.89,1503.51,1774.7,1488.67,1636.06,1603.26,1397.49,1448.42,1216.82,1633.09 +1961,668334.9,458494.96,189654.91,121390.41,94560.85,90138.24,73928.05,74351.76,52789.82,65399.33,50095.6,46470.82,49592.8,46062.91,45984.89,43203.64,39394.13,30817.57,33666.77,30012.67,28146.89,28224.2,27744.71,26094.75,27164.62,22671.19,22480.42,20953.08,21438.03,21171.2,18780.2,20299.11,18304.39,17920.67,16157.87,16982.32,15637.73,12917.62,11665.6,11606.4,10030.53,12710.55,11404.86,9638.41,9225.41,10221.76,8919.6,10493.56,10111.65,10350.24,11072.34,10373.4,9166.76,8311.81,8264.14,7953.38,7530.18,7111.14,7613.35,10411.03,8837.35,7890.16,8418.46,7541.33,8361.44,7289.83,7475.35,7769.48,8461.69,7006.63,6866.54,5322.27,5393.41,5753.02,5873.02,5260.5,5285.23,4894.58,4614.03,5223.57,4467.77,4222.5,4235.94,4676.86,4077.07,4196.17,4721.9,3611.15,3943.36,4336.14,4218.88,3633.65,3764.81,3700.02,4030.32,2820.04,3870.76,3288.6,3694.21,3486.3,3011.11,2806.91,3295.29,3060.36,3406.3,3075.92,3140.26,2839.67,2814.1,2571.69,2465.87,2996.1,2843.24,2395.61,2726.8,3162.03,2329.2,2424.07,2157.66,2241.75,2170.34,2051.95,2096.41,2159.12,1953.33,1596.56,1650.81,1597.53,1941.49,1864.61,2174.27,1529.23,1830.4,1507.65,1685.9,1658.36,1432.64,1498.07,1230.5,1689.55 +1962,679732.32,467852.54,192397.74,122842.95,95458.88,92558.01,74489.14,76573.25,53171.33,67188.06,50487.49,47121.38,51030.14,47029.14,47119.36,43749.47,40649.59,31165.08,34684.17,30386.28,28832.81,29081.03,28506.18,26822.59,28081.23,23221.39,23071.43,21287.68,21898.02,21642.64,18924.16,20832.91,18668.58,18401.61,16490.16,17500.17,16041.26,13168.89,11871.71,11766.73,10062.58,13094.82,11690.15,9691.08,9283.54,10384.2,8938.8,10691.69,10352.19,10650.67,11452.48,10683.91,9345.87,8348.33,8340.24,8022.06,7579.57,7156.83,7669.61,10819.46,9138.1,8067.14,8692.82,7699.14,8628.97,7450.4,7674.22,8004.12,8790.59,7240.17,7085.46,5381.37,5503,5866.06,6028.55,5351.8,5399.92,4960.33,4649.35,5352.5,4498.82,4251.42,4303.13,4812.89,4148.73,4249.21,4875.42,3639.21,4022.59,4464.25,4362.86,3690.66,3838.1,3784.44,4171.43,2824.95,3999.42,3353.23,3841.07,3588.16,3069.91,2846.39,3386.86,3121.22,3521.28,3145.55,3240.59,2893.67,2874.19,2603.89,2502.9,3050.6,2927.86,2442.15,2792.35,3250.31,2363.01,2477.2,2185.88,2319.46,2221.12,2094.69,2155.65,2236.56,2005.34,1605.65,1676.25,1612.76,2009.53,1925.28,2264.31,1556.66,1886.56,1527.11,1737.64,1715.41,1469.65,1550.81,1245.08,1741.42 +1963,692280.49,477527.97,195016.21,124193.39,96389.95,95055.67,75080.05,78854.02,53562.81,69020.55,50903.66,47781.4,52532.42,48032.25,48309.32,44285.9,41939.88,31524.39,35722.09,30733.37,29531.34,29967.04,29281.25,27528.21,29016.77,23798.43,23680.43,21621.84,22371.9,22127.04,19059.94,21385.83,19016.66,18899.28,16829.99,18033.55,16461.93,13424.44,12065.47,11928.36,10097.06,13478.23,11985.14,9742.47,9340.42,10552.27,8947.43,10894.27,10597.52,10961.54,11842.14,11005.91,9533.95,8382.61,8414.52,8090.6,7631.46,7205.99,7725.5,11211.86,9454.25,8247.42,8974.08,7862.07,8908.42,7618.35,7888.91,8248.81,9130.35,7487.43,7303.43,5441.61,5618.16,5980.42,6183.75,5446.06,5520.33,5027.82,4686.28,5486.32,4526.75,4279.06,4377.64,4953.73,4222.14,4298.58,5034.65,3667.11,4103.73,4595.51,4516.66,3749.51,3913.4,3872.12,4315.13,2835.02,4132.76,3417.57,3996.94,3693.87,3131.56,2887.07,3481.75,3183.55,3638.63,3212.09,3345.15,2949.93,2936.44,2635.13,2542.86,3102.97,3015.89,2491.28,2860.25,3332.16,2398.41,2530.17,2213.8,2401.46,2273.35,2139.3,2216.71,2318.23,2058.92,1614.02,1703.4,1631.76,2077.58,1984.98,2355.93,1585.76,1943.59,1547.45,1790.53,1773.85,1508.27,1607.17,1260.21,1789.34 +1964,706460.95,487484.54,197491.84,125412.63,97379.85,97638.03,75676.01,81168.65,53945.02,70900.96,51322.65,48402.84,54129.1,49066.06,49551.9,44794.33,43264.27,31885.4,36780.99,31056.32,30244.23,30881.33,30071.1,28220.78,29962.88,24397.02,24308.09,21953.93,22867.74,22610.48,19207.14,21957.35,19359.17,19412.98,17188.13,18581.97,16903.92,13695.74,12282.42,12089.39,10133.08,13857.14,12295.97,9784.36,9394.29,10728.2,8940.09,11115.21,10849.98,11281.02,12241.54,11339.1,9731.36,8416.19,8489.93,8157.72,7686.73,7256.09,7780.69,11577.47,9788.99,8430.84,9253.96,8030.03,9200.16,7787.15,8122.2,8503.99,9476.25,7746.2,7513.29,5503.75,5729.46,6093.32,6331.58,5543.34,5646.32,5098.89,4723.63,5625.16,4550.24,4304.99,4458.61,5099.47,4296.66,4343.46,5200.34,3695.41,4186.64,4730.54,4677.4,3811.66,3990.86,3963.42,4456.69,2847.34,4269.86,3479,4157.97,3802.64,3196.11,2927.38,3580.31,3246.51,3757.96,3275.46,3452.94,3010.86,3001.13,2665.39,2585.97,3161.72,3106.19,2537.64,2930.28,3410.96,2435.2,2580.53,2240.95,2485.43,2327.14,2186.34,2280.05,2402.46,2114.1,1621.88,1730.49,1662.07,2145,2041.21,2443.3,1616.52,2002.12,1569.14,1843.59,1833.06,1548.19,1667.83,1275.44,1834.49 +1965,722562.18,497702.37,199815.54,126483.88,98447,100308.89,76255.93,83498.02,54303.11,72831.67,51728.66,48952.26,55834.04,50127.21,50845.22,45261.94,44623.04,32241.04,37860.01,31357.73,30972.97,31822.8,30875.96,28907.07,30913.93,25013.63,24955.12,22283.39,23391.15,23082.71,19379.57,22547.25,19703.88,19942.3,17571.23,19144.22,17369.88,13989.96,12547.53,12248.34,10169.93,14229.04,12626.95,9811.41,9443.94,10913.72,8914.45,11363.31,11110.83,11607.68,12649.73,11683.53,9938.41,8450.5,8568.45,8222.47,7745.67,7305.16,7835.08,11909,10143.74,8617.08,9526.56,8203.08,9504.7,7951.93,8375.79,8770.1,9824.69,8014.4,7710.55,5568.48,5829.96,6203.3,6467.2,5643.64,5777.83,5174.87,4760.46,5769.22,4568.6,4328.92,4545.34,5250.12,4371.52,4383.45,5373.14,3724.46,4271.13,4869.72,4843.64,3877.81,4070.59,4058.67,4592.61,2859.96,4410.21,3535.64,4321.79,3913.93,3263.64,2966.16,3682.88,3309.57,3878.95,3335.63,3563.41,3077.88,3068.44,2694.54,2632.36,3232.93,3197.86,2577.52,3002.43,3489.21,2473.29,2626.65,2266.94,2569.46,2382.59,2236.21,2346.01,2487.95,2170.86,1629.43,1756.27,1708.63,2211.32,2092.35,2522.53,1648.83,2062.63,1592.43,1896.13,1892.6,1589.18,1733.31,1290.4,1877.81 +1966,740745.64,508161.94,201975.22,127395.88,99595.13,103067.35,76816.34,85837.8,54630.99,74817.55,52115.99,49411.39,57672.99,51217.36,52191.1,45682.31,46011.04,32587.24,38959.33,31637.77,31717.48,32789.1,31697.62,29584.63,31867.56,25641.38,25624.66,22608.75,23944.18,23539.12,19582.33,23154.78,20051.96,20486.44,17982.61,19721.46,17861.88,14310.53,12864.95,12404.19,10208.38,14593.28,12980.27,9821.12,9488.75,11109.88,8866.85,11643.45,11380.68,11941.33,13067.81,12038.9,10152.33,8486.43,8650.32,8284.34,7809.09,7352.24,7888.8,12201.04,10520.88,8806.14,9789.98,8381.46,9822.5,8110.43,8651.16,9047.8,10175.14,8292.78,7890.99,5635.86,5916.91,6309.77,6584.77,5748.59,5915.12,5256.36,4796.49,5918.6,4581.12,4350.81,4638.28,5405.69,4446.82,4418.21,5553.25,3754.36,4357.48,5013.15,5015.2,3948.87,4152.67,4158.12,4721.53,2872.07,4553.43,3586.63,4487.2,4027.76,3334.19,3002.79,3789.21,3372.17,4001.38,3392.09,3676.19,3152.72,3143.84,2722.88,2682.16,3319.08,3290.41,2609.27,3076.56,3568.38,2512.65,2667.08,2291.47,2653.05,2439.2,2289.11,2414.81,2574.48,2229.38,1636.67,1780.26,1774.03,2276.03,2136.64,2591.59,1682.89,2125.24,1617.79,1947.79,1952.14,1631.15,1803.68,1304.99,1919.27 +1967,760771.06,518889.78,203981.93,128164.76,100820.02,105907.4,77352.66,88191.38,54928.53,76856.31,52488.49,49791.94,59620.67,52341.83,53590.93,46060.45,47429.81,32926.31,40074.7,31898.32,32477.96,33778.5,32534.02,30249.46,32826.6,26281.21,26318.12,22932.2,24524.55,23982.58,19810.59,23780.49,20401.08,21045.79,18416.99,20311.37,18378.21,14652.03,13222.69,12556.84,10248.86,14950.8,13354.2,9817.19,9529.04,11316.83,8802.56,11949.36,11657.66,12282.08,13492.5,12406.04,10372.63,8524.56,8734.04,8343.39,7875.91,7397.17,7942.64,12455.71,10917.45,8997.33,10046.17,8565.67,10154.48,8263.55,8947.4,9337.66,10528.05,8580.68,8057.44,5706.2,5991.68,6415,6685.32,5858.64,6058.54,5343.02,4831.75,6073.53,4588.66,4370.86,4737.63,5566.06,4521.71,4448.55,5740.71,3784.87,4445.53,5160.61,5194.85,4023.49,4237.13,4262.01,4843.87,2884.78,4700.04,3632.67,4656.35,4144.4,3407.8,3037.44,3899.24,3434.81,4125.11,3444.85,3791.89,3234.02,3228.5,2750.09,2735.31,3418.32,3383.7,2634.36,3152.83,3648.03,2553.53,2702.68,2314.5,2736.08,2496.92,2344.98,2486.41,2662.23,2289.58,1643.7,1803.06,1855.44,2339.12,2174.85,2652.7,1718.6,2189.88,1644.94,1998.69,2011.76,1674.2,1878.88,1319.25,1958.73 +1968,782008.73,529967.32,205880.42,128836.02,102120.48,108821.56,77841.44,90557.06,55194.46,78934.99,52851.2,50127.14,61579.47,53505.98,55042.4,46409,48894.02,33265.17,41195.84,32146.05,33256.43,34790.95,33377.26,30905.37,33797.04,26946.08,27032.94,23261.28,25128.12,24419.42,20055.92,24426.95,20751.37,21622.59,18861.86,20905.06,18913.2,15001.91,13609.98,12706.79,10289.94,15302.95,13744.39,9808.04,9565.6,11534.26,8739.08,12265.09,11937.61,12629.33,13913.35,12787.49,10604.35,8565.67,8817.94,8399.92,7942.32,7439.58,7998,12682.93,11323.1,9188.82,10297.8,8756.48,10502.25,8413.33,9260.68,9639.84,10882,8872.92,8221.02,5780.84,6056.36,6523.79,6778.72,5971.41,6208.28,5434.04,4866.07,6234.47,4594.02,4389.32,4842.17,5730.91,4593.31,4476.82,5935.86,3815.64,4534.23,5311.62,5387.49,4097.19,4324.06,4370.65,4960.24,2900.18,4851.43,3675.45,4834.28,4263.75,3484.54,3070.69,4013.54,3499.35,4250.03,3495.14,3912.09,3316.23,3313.79,2774.77,2791.59,3527.26,3477.74,2656.3,3231.69,3726.71,2596.57,2736.83,2336.32,2817.34,2556.85,2403.6,2560.73,2750.89,2350.9,1651.09,1825.63,1945.78,2401.14,2210.96,2711.3,1755.34,2256.78,1672.4,2049.15,2071.79,1718.53,1958.91,1333.18,1996.87 +1969,803577.22,541505.08,207733,129472.84,103491.93,111800.09,78254.63,92935.07,55429.64,81035.91,53212.82,50466.39,63417.39,54716.74,56542.43,46746.67,50423.48,33614.19,42309.67,32389.55,34055.36,35826.8,34216.83,31557.58,34787.59,27654.16,27765.24,23605.99,25748.64,24859.64,20305.31,25097.52,21102.23,22219.9,19299.81,21490.95,19458.9,15342.67,14010.34,12855.01,10329.56,15651.92,14144.44,9805.46,9599.65,11761.47,8700.02,12568.3,12214.97,12982.45,14316.71,13186.56,10854.43,8610.47,8899.69,8454.5,8003.19,7479.35,8056.67,12897.41,11723.85,9378.24,10549.23,8954.81,10867.72,8563.19,9585.58,9954.41,11235.49,9162.83,8397.35,5861.41,6114.72,6642.63,6879.18,6083.62,6364.57,5528.17,4899.29,6401.92,4601.03,4406.61,4950.15,5899.85,4658,4506.16,6139.05,3846.18,4622.21,5465.51,5599.63,4164,4413.59,4484.44,5071.93,2921.26,5009.51,3717.47,5027.97,4385.76,3564.47,3103.51,4132.84,3568.38,4376.05,3544.64,4038.92,3391.75,3387.63,2795.05,2850.66,3640.59,3572.71,2680.22,3313.64,3802.1,2642.61,2774.37,2357.4,2895.54,2620.43,2464.55,2637.52,2840.23,2412.57,1659.55,1849.41,2034.91,2462.93,2250.6,2775.16,1792.22,2326.14,1698.14,2099.65,2132.8,1764.34,2043.82,1346.82,2034.53 +1970,824788.46,553578.51,209588.15,130123.43,104925.65,114834.78,78572.98,95326.79,55634.94,83148.81,53578.68,50843.83,65047.77,55981.4,58090.76,47086.76,52029.86,33980.28,43407.29,32635.56,34876.27,36884.91,35046.27,32209.31,35804.73,28415.08,28514.01,23973.06,26381.43,25309.99,20548.91,25794.21,21452.74,22839.45,19718.03,22061.22,20009.94,15661.29,14410.4,13001.94,10366.34,16000.01,14550.03,9817.82,9632.18,11997.93,8701.86,12842.91,12485.76,13341.07,14693.15,13605.53,11126.12,8659.8,8977.64,8507.48,8054.92,7516.24,8119.86,13110,12110.03,9563.87,10803.98,9161.53,11252.49,8715.12,9917.98,10281.7,11587.76,9446.06,8596.98,5949.05,6169.36,6776.38,6994.85,6193.81,6527.64,5624.6,4931.24,6576.31,4612.37,4423.07,5060.4,6072.53,4713.32,4538.84,6350.54,3876.06,4708.64,5621.79,5836.39,4219.77,4505.78,4603.72,5180.03,2950,5175.62,3760.53,5242.4,4510.48,3647.71,3136.58,4257.51,3643.55,4503.11,3594.52,4173.93,3455.61,3444.55,2809.8,2912.34,3754.54,3668.6,2709.84,3399.12,3873.11,2692.26,2818.39,2378.02,2969.86,2688.43,2527.59,2716.66,2930.08,2474.11,1669.58,1875.38,2115.52,2525.07,2297.39,2849.62,1828.71,2398.1,1720.8,2150.6,2195.17,1811.88,2133.53,1360.2,2072.29 +1971,845481.58,566224.81,211461.6,130805.2,106425.27,117922,78786.61,97728.96,55809.75,85269.37,53949.93,51273.47,66424.74,57295.21,59687.14,47433.81,53718.72,34363.92,44485.91,32886.6,35720.57,37964.93,35863.38,32863.7,36851.06,29245.21,29281.27,24366.44,27024.99,25772.8,20783.75,26518.67,21804.27,23482.81,20113.7,22611.99,20562.87,15954.98,14809.52,13147.93,10398.44,16347.2,14960.11,9847.39,9663.32,12243.77,8750.51,13084.16,12747.84,13704.34,15038,14045.82,11417.83,8712.68,9051.22,8558.97,8095.79,7550.34,8187.99,13324.67,12477.06,9745.19,11062.34,9375.14,11657.51,8869.96,10255.9,10621.47,11937.81,9720.4,8827.27,6044.53,6221.64,6927.27,7137.75,6300.55,6697.75,5723.38,4961.76,6757.85,4629.38,4438.81,5172.69,6248.84,4758.15,4575.56,6570.86,3905.27,4793.16,5780.48,6100.99,4263.84,4600.59,4728.7,5284.53,2986.78,5351.2,3805.29,5479.34,4637.83,3734.42,3170.07,4388.46,3726.09,4631.11,3645.52,4317.75,3505.39,3470.32,2818.27,2976.57,3868.34,3765.17,2745.97,3488.29,3937.83,2745.78,2870.82,2398.38,3039.77,2762.27,2592.63,2798.13,3020.39,2535.36,1681.43,1904.02,2185.66,2587.71,2353.56,2936.67,1864.6,2472.66,1739.52,2202.04,2258.96,1861.2,2228.15,1373.31,2110.98 +1972,865686.53,579411.51,213351.77,131514.1,107977.04,121059.51,78901.71,100143.6,55953.32,87408.88,54321.02,51739.84,67597.47,58662.6,61338.26,47783.01,55478.15,34760.52,45549.48,33143.19,36587.23,39061.99,36673.64,33518.39,37925.4,30132.58,30074.3,24782.95,27680.14,26244.32,21011.7,27268.98,22156.21,24148.14,20489.7,23146.8,21120.14,16226.56,15207.77,13291.87,10426.43,16695,15377.09,9891.28,9692.85,12499.43,8839.8,13296.87,13002.28,14072.48,15354.59,14506.62,11721.94,8769.76,9120.96,8608.61,8127.03,7581.12,8260.63,13538.75,12828.63,9922.56,11324.25,9595.76,12083.19,9025.3,10599.85,10974.62,12286.44,9988.38,9083.57,6147.46,6270.12,7094.83,7300.15,6407.3,6875.23,5825.17,4990.55,6946.62,4651.17,4453.99,5287.54,6428.71,4793.55,4615.83,6800.14,3933.73,4876.56,5941.57,6393.89,4298.09,4698.08,4859.61,5385.27,3030.4,5535.87,3851.15,5737.28,4768.08,3824.76,3203.63,4525.11,3815.1,4759.93,3697.15,4469.9,3544.05,3475.02,2821.44,3043.57,3983.7,3861.93,2787.35,3581.07,3998.31,2803.03,2929.69,2418.27,3105.54,2840.84,2659.85,2882.11,3111.26,2596.74,1694.98,1934.83,2247.58,2650.48,2416.74,3033.92,1900.32,2549.77,1754.96,2253.84,2324.01,1912.3,2327.49,1386.07,2150.52 +1973,885145.93,593058.93,215268.69,132250.78,109532.88,124242.3,78936.17,102584.28,56066.83,89592.1,54682.63,52212.31,68691.19,60110.43,63059.48,48127.17,57280.59,35165.56,46604.73,33407.62,37472.3,40164.97,37488.07,34164.54,39026.08,31025.12,30904.27,25213.39,28347.34,26718.43,21233.47,28039.38,22503.24,24829.69,20846.07,23674.5,21689.24,16475.85,15593.35,13431.18,10454.06,17049.17,15804.43,9946.38,9720.76,12764.96,8958.96,13486.44,13252.09,14447.65,15651.32,14985.13,12027.82,8835.14,9187.68,8655.67,8151.6,7607.07,8337.08,13748.32,13173.59,10096.3,11592.7,9827.58,12529.85,9176.05,10951.17,11343.93,12636.97,10256.43,9350.11,6256.19,6311.48,7277.96,7447.29,6523.45,7060.6,5930.48,5017.25,7142.63,4675.46,4469.04,5405.36,6611.92,4822.16,4659.28,7037.85,3960.8,4960.66,6104.53,6714.1,4324.36,4798.51,4996.94,5483.08,3079.3,5727.04,3897.26,6013.86,4902.01,3918.92,3236.73,4664.44,3907.63,4889.44,3747.93,4629.4,3578.49,3506.01,2822.08,3113.68,4102.32,3958.32,2833.13,3677.23,4062.04,2863.74,2989.41,2437.18,3168.29,2919.29,2729.58,2968.99,3203.02,2659.09,1709.98,1966.7,2303.35,2712.78,2480.42,3137.03,1936.84,2629.51,1768.99,2306,2390.21,1965.16,2430.76,1398.47,2189.64 +1974,903613.54,607050.26,217218.39,133008.91,111030.14,127465.23,78917.22,105069.37,56152.18,91852.92,55022.22,52646.2,69884.42,61673.56,64870.83,48455.12,59088.19,35572.48,47661.77,33681.72,38370.24,41259.54,38322.02,34789.94,40149.96,31851.71,31785.5,25644.51,29027.73,27186.06,21451.13,28821.85,22838.19,25519.6,21184.32,24208.02,22280.92,16704.23,15952.08,13562.44,10486.63,17416.96,16247.11,10007.47,9747.02,13040.4,9091.4,13662.05,13501.99,14832.84,15940.9,15477.29,12321.54,8914.18,9252.68,8699.14,8173.97,7626.29,8416.25,13946.82,13525.09,10267.06,11871.23,10076.17,12997.6,9315.37,11312.31,11732.96,12994.03,10533.72,9604.28,6368.35,6341.06,7474.34,7531.42,6661.57,7254.47,6040.04,5041.42,7345.78,4698.78,4484.52,5526.76,6798.21,4847.85,4705.11,7283.18,3985.74,5047.94,6268.71,7059.33,4345.55,4902.17,5141.2,5579.08,3131.18,5920.94,3942.22,6305.29,5040.66,4017.08,3268.64,4802.35,3999.51,5019.47,3795.65,4794.75,3618.59,3627.5,2824.07,3187.41,4226.8,4053.71,2881.75,3776.4,4139.1,2927.47,3042.18,2454.39,3229.69,2990.97,2802.24,3059.25,3296.1,2723.52,1725.96,1998.03,2356.62,2773.75,2535.5,3239.56,1975.52,2711.85,1784.4,2358.46,2457.38,2019.72,2536.89,1410.47,2226.4 +1975,920945.08,621301.72,219205.3,133784.38,112423.06,130724.12,78866.5,107612.1,56211.95,94217.2,55330.69,53010.73,71305.92,63373.57,66787.9,48758.99,60872.4,35975.32,48729.39,33966.57,39277.21,42334.95,39187.7,35386.51,41295.12,32566.82,32730.55,26066.98,29721.97,27641.02,21665.64,29610.74,23156.21,26212.41,21506.86,24756.97,22902.32,16913.84,16274.74,13683.29,10527.27,17803.7,16709.1,10070.14,9771.59,13325.81,9224.06,13831.21,13755.16,15229.95,16232.67,15980.3,12590.29,9009.92,9316.96,8738.36,8197.34,7637.69,8497.17,14130.07,13892.64,10435.53,12162.37,10344.49,13486.63,9438.44,11684.59,12144.14,13360.99,10827.15,9831.41,6482.28,6356.19,7682.48,7522.59,6830.69,7457.36,6154.55,5062.73,7556.03,4718.75,4500.87,5652.48,6987.39,4873.62,4752.54,7535.71,4008.03,5140.36,6433.73,7428.7,4364.51,5009.26,5292.81,5674.14,3184.13,6115.37,3985.1,6608.61,5184.81,4119.44,3298.87,4936.21,4087.95,5149.94,3838.91,4964.83,3671.49,3880.32,2830.17,3265.17,4359.09,4147.53,2931.89,3878.33,4236,2993.88,3082.63,2469.44,3291.08,3051.58,2878.16,3153.26,3390.94,2790.96,1742.62,2027.74,2410.45,2832.76,2575.69,3336.76,2017.37,2796.75,1803.01,2411.23,2525.36,2075.97,2645.14,1422.04,2259.75 +1976,937018.25,635771.73,221239.22,134580.09,113690.11,134010.69,78794.5,110213.08,56247.41,96690.54,55604.29,53295.79,72999.14,65226.23,68813.22,49036.46,62623.76,36373.1,49808.07,34261.87,40189.51,43386.84,40089.03,35948.93,42461.19,33146.89,33737.77,26477.15,30428.03,28081.33,21877.4,30402.64,23456.51,26904.35,21811.92,25323.41,23559.07,17103.19,16554.75,13793.11,10578.46,18210.75,17190.24,10134.72,9794.78,13621.11,9354.49,13994.62,14012.82,15639.9,16528.71,16493.31,12840.3,9124.76,9380.45,8772.54,8223.25,7640.75,8579.63,14296.77,14279.12,10601.84,12468.89,10632.93,13996.7,9544.27,12068.17,12578.41,13739.14,11139.83,10023.47,6596.77,6355.23,7901,7402.87,7034.87,7669.45,6274.04,5081.27,7773.45,4734.44,4518.1,5781.8,7179.4,4900.18,4801.74,7794.66,4027.48,5238.25,6599.21,7818.61,4381.6,5119.83,5454.71,5768.72,3237.59,6308.3,4025.27,6922.98,5334.92,4224.53,3327.26,5064.67,4172.23,5280.72,3876.95,5139.03,3739.66,4289.47,2841.43,3347.17,4499.51,4239.68,2983.58,3982.7,4356.58,3062.96,3108.57,2482.18,3352.62,3098.97,2957.34,3251.16,3487.64,2861.58,1759.73,2055.09,2464.46,2889.58,2598.35,3427.66,2062.41,2884.16,1825.55,2464.34,2594.31,2133.72,2754.7,1433.18,2288.89 +1977,951927.44,650485.03,223324.04,135402.8,114839.12,137322.12,78701.8,112867.87,56262.02,99270.09,55845.07,53516.21,74925.9,67215.81,70946.23,49290.91,64345.88,36763.53,50899.5,34566.89,41108.25,44416.01,41026.48,36481.59,43650.33,33622.39,34810.72,26878.57,31144.32,28509.66,22083.77,31199.33,23741.59,27597.3,22101.73,25905.13,24247.55,17274.47,16796.58,13892.92,10637.07,18636.98,17690.18,10200.92,9816.7,13926.26,9483.63,14152.07,14273.27,16061.32,16826.47,17017.67,13067.54,9255.38,9443,8801.89,8250.7,7636.58,8663.33,14450.53,14681.46,10766.42,12790.55,10936.94,14528.29,9634.68,12460.91,13034.63,14127.79,11470.87,10189.89,6712.4,6341.18,8130.99,7194.28,7271.87,7890.97,6398.94,5097.05,7998.16,4746.93,4536.2,5915.01,7374.23,4927.24,4852.5,8060.65,4044.33,5341.42,6765.52,8231.6,4398.48,5233.68,5627.53,5863.13,3291.17,6501.89,4063.19,7248.83,5490.92,4332.29,3354.06,5189.54,4253.99,5411.87,3910.64,5317.63,3821.19,4827.36,2857.11,3433.44,4647.62,4329.96,3036.75,4089.55,4496.87,3134.8,3122.41,2492.98,3414.32,3135.84,3039.66,3352.84,3586.5,2935.38,1777.32,2080.54,2518.57,2944.38,2606.22,3513.53,2110.46,2973.81,1851.07,2517.87,2664.26,2193.07,2865.64,1443.89,2315.02 +1978,966039.82,665502.28,225449.66,136255.4,115890.28,140665.86,78589.31,115577.67,56263.66,101962.45,56057.03,53696.45,77033.85,69293.55,73194.94,49526.88,66039.49,37137.7,52015.28,34881.38,42039.94,45423.44,42004.66,36997.8,44866.27,34068.32,35972.65,27277.74,31872.23,28930.27,22279.3,32006.9,24014.08,28298.15,22387.77,26502.17,24954.66,17438.25,17015.98,13984.06,10693.83,19081.72,18212.33,10263.11,9836.81,14241.4,9605.59,14309.3,14533.38,16491.08,17125.17,17555.49,13237.73,9392.11,9505.45,8828.01,8277.41,7627.76,8746.98,14597.85,15096.01,10930.78,13123.07,11248.05,15082.99,9711.39,12859.09,13510.42,14525.93,11818.31,10354.5,6831.3,6321.54,8376.15,6955.57,7536.76,8122.53,6530.82,5109.69,8230.22,4758.25,4555.42,6054.91,7571.96,4955.5,4903.2,8336.42,4059.18,5450.55,6933.91,8679.84,4421.13,5350.32,5806.85,5957.93,3343.15,6703.18,4100.35,7585.91,5652.36,4445.83,3379.8,5315.27,4335.65,5543.52,3942.14,5501.45,3913.77,5417.74,2875.97,3523.94,4803.73,4417.52,3089.76,4199.42,4646.18,3209.26,3128.96,2502.51,3477.23,3168.84,3125.03,3458.1,3688.39,3012.83,1795.65,2105.66,2576.47,2997.41,2604.87,3594.49,2162.25,3065.12,1877.69,2571.85,2734.9,2254.45,2979.09,1454.27,2341.43 +1979,979880.62,680915.8,227599.88,137140.56,116878.03,144053.52,78454.85,118342.63,56262.13,104772.87,56246.8,53872.44,79236.78,71391.29,75567.68,49751.26,67709.69,37484.75,53169.67,35204.73,42993.99,46412.31,43027.82,37516.49,46114,34590.23,37252.66,27684.53,32613.89,29349.5,22456.95,32834.45,24278.05,29017.05,22685.18,27113.51,25661.88,17608.54,17235.67,14068.59,10735.99,19543.35,18760.76,10314.15,9854.38,14566.69,9714.16,14473.24,14788.61,16924.75,17422.72,18109.88,13306.7,9521.67,9568.88,8853.19,8300,7617.88,8829.04,14747.56,15516.86,11096.87,13460.2,11554.98,15662.85,9777.29,13257.8,14002.3,14931.74,12178.54,10550.78,6956.58,6306.54,8641.52,6768.72,7821.55,8364.84,6671.66,5118.71,8469.67,4771.35,4576.03,6205.21,7772.65,4985.56,4951.79,8625.69,4072.87,5566.27,7106.15,9179.62,4457.08,5469.12,5986.33,6053.65,3391.58,6921.79,4138.82,7934.28,5818.51,4569.42,3405.18,5448.11,4420.72,5675.93,3974.51,5691.75,4013.31,5953.62,2896.02,3618.53,4968.07,4501.32,3140.56,4313.14,4789.32,3286.18,3135.3,2511.74,3542.52,3207.33,3213.36,3566.67,3794.42,3094.48,1815.1,2132.69,2642.85,3049.11,2602.57,3671.28,2218.58,3157.36,1902.72,2626.29,2805.82,2318.5,3096.73,1464.42,2372.62 +1980,993877.31,696783.52,229763.05,138058.71,117827.36,147490.37,78300.65,121159.76,56265.48,107702.68,56419.28,54070.82,81470.86,73460.72,78068.14,49968.81,69360.87,37796.8,54372.51,35535.11,43975.92,47385.32,44099.14,38050.42,47396.97,35264.9,38668.22,28105.89,33369.71,29772.04,22611.97,33687.79,24537.42,29760.47,23005.07,27737.9,26357.46,17795.73,17472.14,14148.42,10754.77,20019.85,19337.72,10349.27,9869,14902.16,9804.85,14649.11,15035.86,17359.12,17717.37,18683.16,13248.37,9634.5,9633.89,8878.78,8316.34,7609.75,8908.29,14905.89,15939.74,11266.23,13798.13,11848.33,16268.99,9835.18,13653.36,14507.47,15343.92,12549.54,10802.03,7090.13,6303.61,8929.9,6692.11,8120.5,8618.35,6822.84,5123.95,8716.55,4788.24,4598.13,6368.17,7976.45,5017.82,4996.73,8930.77,4086.08,5688.84,7283.46,9740.6,4511.9,5589.58,6163.08,6150.74,3434.95,7164.17,4179.86,8293.68,5988.9,4705.22,3430.76,5592.65,4512.04,5809.27,4009.84,5889.23,4116.82,6359.13,2915.78,3717.17,5140.72,4580.7,3187.67,4431.11,4915.47,3365.44,3146.62,2521.39,3610.98,3258.14,3304.47,3678.29,3905.41,3180.63,1835.84,2163.05,2720.84,3099.75,2605.29,3744.67,2279.82,3249.91,1924.2,2681.24,2876.81,2385.54,3219.47,1474.44,2411.7 +1981,1008000.15,713118.03,231938.96,139002.31,118743.18,150978.84,78126.1,124030.91,56276.32,110752.91,56576.5,54298.05,83721.27,75482.55,80696.95,50179.03,70992.2,38071.26,55627.75,35873.73,44988.36,48337.5,45216.51,38606.82,48715.59,36120.29,40217.63,28543.36,34139.13,30198.59,22739.82,34570.58,24791.14,30532.95,23351.57,28375.99,27039.47,18003.32,17731.23,14223.76,10746.96,20511.6,19943.66,10365.72,9880.59,15249.01,9875.26,14838.49,15273.39,17792.55,18009.44,19277.11,13053.95,9727.72,9700.25,8905.07,8325.49,7604.51,8984.59,15072.36,16363.56,11439.14,14133.84,12133.07,16901.68,9884.21,14046.54,15027.27,15761.8,12930.21,11117.61,7234.3,6315.88,9244.51,6748.19,8434.02,8883.02,6985.16,5125.3,8971.35,4809.81,4621.42,6545.02,8183.19,5051.23,5037.34,9252.85,4098.97,5818.67,7466.49,10366.66,4589.78,5711.6,6327.57,6249.32,3472.89,7431.94,4222.51,8664.06,6164.01,4853.93,3456.14,5750.34,4610.17,5943.59,4048.62,6094.21,4223.2,6604.87,2935.04,3820.13,5315.03,4655.36,3230.47,4553.9,5020.49,3445.28,3164.83,2531.01,3682.67,3323.38,3398.47,3792.94,4020.78,3271.46,1857.85,2197.58,2812.04,3148.09,2615.75,3815.23,2346.8,3342.67,1941.53,2735.33,2947.78,2454.77,3347.78,1483.87,2460.39 +1982,1022253.4,729868.01,234132.63,139964.75,119615.62,154506.27,77944.85,126947.37,56296.24,113916.81,56717.08,54547.69,86007.33,77472.91,83445.86,50379.38,72602.53,38309.71,56931.82,36216.71,46025.36,49267.56,46379.62,39179.55,50068.49,37136.85,41883.33,28993.99,34917.9,30628.34,22842.32,35479.1,25040.95,31330.26,23722.34,29027.16,27717.34,18229.6,18008.56,14295.22,10716.32,21016.82,20575.7,10365.45,9889.82,15607.24,9925.73,15041.07,15502.52,18225.73,18299.16,19891.55,12749.65,9804.17,9767.26,8930.67,8329.43,7602.18,9057.98,15244.05,16790.07,11615.84,14470.63,12409.24,17559.43,9925.62,14438.31,15562.19,16185.89,13323.33,11488.11,7387.66,6341.55,9582.16,6918.1,8764.62,9158.57,7158.26,5123.35,9234.13,4835.1,4645.41,6733.96,8392.94,5085.51,5073.89,9590.23,4111.66,5955.27,7654.82,11048.08,4690.61,5835.18,6484.45,6349.56,3505.36,7721.54,4265.31,9044.47,6344.38,5011.73,3481.13,5920.06,4715.2,6078.82,4089.86,6305.71,4333.39,6716.45,2954.28,3927.71,5489.32,4725.72,3269.25,4680.89,5107.91,3525.4,3188.54,2540.44,3757.18,3401.24,3495.2,3910.66,4140.26,3366.72,1880.7,2235.33,2915.07,3193.69,2632.28,3883.47,2418.84,3435.53,1955.24,2788.31,3019.07,2525.52,3480.45,1492.64,2517.58 +1983,1037123.83,746949.07,236344,140946.86,120438.09,158044.34,77790.08,129882.32,56330.88,117177.06,56838.68,54817.83,88338.24,79462.28,86297.64,50570.03,74196.55,38514.06,58277.39,36550.85,47073.42,50186.2,47594.56,39752.02,51455.03,38285.88,43645.09,29454.74,35697.94,31062.43,22928.93,36403.09,25296.84,32139.71,24113.31,29687.09,28404.88,18473.09,18298.21,14365.39,10669.69,21528.5,21228.29,10354.31,9897.88,15974.42,9958.42,15257.36,15726.8,18660.44,18584.31,20524.67,12389.27,9867.16,9834.42,8953.07,8332.44,7603.1,9127.68,15421.23,17221.21,11797.53,14818.62,12657.71,18239.4,9966.73,14825.79,16107.73,16617.35,13735.27,11895.13,7543.74,6377.1,9931.56,7168.24,9111.1,9445,7340.91,5119.61,9504.28,4862.18,4669.61,6930.39,8606.21,5122.38,5106.99,9938.85,4124.59,6096.69,7847.47,11763.84,4810.5,5959.96,6661.36,6452.08,3531.78,8027.57,4308.11,9432.73,6529.89,5170.56,3506.54,6100.5,4829.09,6214.86,4132.45,6521.54,4448.73,6740.22,2973.46,4039.95,5673.61,4792.9,3304.68,4809.9,5183.91,3608.75,3215.73,2550.93,3834.17,3489.98,3594,4031.35,4265.25,3465.79,1903.55,2273.67,3026.24,3238.59,2651.29,3950.12,2493.14,3527.94,1965.9,2842.62,3091.51,2598.41,3614.69,1501.42,2580.26 +1984,1053210.74,764245.2,238573.86,141950.63,121199.8,161555.58,77704.34,132800.68,56386.23,120510.4,56937.7,55102.74,90732.36,81497.74,89228.95,50750.93,75780.61,38688.12,59653.09,36858.9,48114.11,51108.08,48868.95,40300.4,52873.97,39518.8,45474.71,29920.9,36468.89,31502.16,23013.16,37328.2,25571.67,32943.58,24517.51,30350.09,29121.47,18730.25,18590.14,14437.51,10617.23,22037.61,21893.85,10340.85,9906.37,16347.24,9977.17,15487.21,15951.42,19099.58,18862.03,21173.6,12047.12,9922.43,9901.05,8968.86,8340.44,7607.29,9192.92,15603.72,17659.98,11985.66,15191.63,12853.78,18937.74,10017.06,15205.5,16658.05,17057.79,14174.47,12311.16,7693.67,6416.71,10277.32,7446.02,9472.17,9742.26,7531.24,5116.06,9780.87,4888.2,4693.43,7127.94,8823.75,5164.14,5137.69,10293.05,4138.27,6240.33,8042.9,12484.97,4943.14,6085.5,6895.93,6557.59,3551.65,8342.2,4350.75,9826.06,6720.34,5319.61,3533.48,6289.33,4954.05,6351.57,4174.62,6738.77,4571.29,6747.93,2992.65,4156.82,5881.91,4858.53,3337.86,4937.98,5258.16,3699.47,3243.04,2564.11,3913.06,3586.38,3694.04,4154.89,4397.53,3567.75,1925.31,2308.95,3140.24,3285.6,2667.22,4016.13,2565.8,3619.25,1974.42,2901.59,3166.22,2674.29,3746.72,1511.2,2643.95 +1985,1070863.39,781666.67,240824.12,142971.17,121894.04,165012.2,77717.43,135676.28,56466.13,123896.53,57012.35,55397.48,93199.87,83613.3,92219.49,50920.78,77360.71,38835.88,61049.37,37128.72,49133.88,52041.47,50204.99,40808.57,54323.65,40800.34,47342.7,30388.78,37222.3,31948.03,23103.65,38244.84,25873.95,33730.15,24929.4,31011.69,29883.45,18997.88,18877.24,14513.95,10566.82,22537.38,22565.91,10331.37,9916.58,16723.96,9985.58,15729.64,16179.8,19544.96,19130.09,21837,11783.05,9974.74,9966.15,8975.29,8357.65,7614.87,9253.41,15789.44,18108.3,12181.03,15598.94,12984.41,19651.23,10082.99,15576.4,17210.19,17508.06,14646.62,12716.23,7831.89,6456.22,10609.04,7712.98,9847.9,10050.02,7727.91,5114.24,10063.5,4911.22,4716.12,7321.88,9045.98,5211.44,5166.8,10648.63,4153.11,6384.2,8240.06,13189.12,5084.77,6211.55,7211.11,6666.46,3564.96,8658.86,4392.13,10222.56,6915.93,5450.42,3562.4,6484.74,5091.54,6488.86,4214.92,6955.21,4702.07,6791.72,3011.91,4278.5,6120.11,4923.86,3369.68,5063.19,5337.62,3799.55,3268.24,2580.75,3993.11,3687.9,3794.72,4281.19,4537.79,3671.83,1945.14,2338.64,3252.99,3335.94,2676.58,4082.79,2634.23,3709.09,1981.53,2966.8,3244.02,2753.15,3873.78,1522.4,2705.54 +1986,1090348.06,799181.44,243098.94,144011.51,122519.65,168402.03,77837.5,138499.46,56574.28,127323.03,57061.78,55704.49,95742.43,85818.5,95264.46,51078.79,78934.13,38958.07,62459.56,37354.27,50128.49,52996.47,51607.7,41269.32,55804.07,42120.73,49256.84,30857.24,37957.33,32400.72,23205.29,39144.93,26208.08,34490.42,25348.4,31669.78,30685.82,19276.47,19156.8,14595.76,10519.79,23023.94,23241.27,10327.53,9928.54,17101.14,9983.22,15986.16,16412.71,19996.25,19387.42,22511.24,11601.04,10023.37,10030.07,8972.47,8384.99,7625.39,9306.35,15980.99,18565.48,12384.11,16045.05,13034.39,20378.63,10168.09,15936.38,17757.17,17968.55,15154.52,13104.3,7955.16,6494.09,10921.04,7958.98,10232.73,10368.3,7930.69,5114.66,10352.12,4930.34,4737.95,7509.76,9272.91,5266.83,5194.31,11004.27,4169.32,6527.54,8438.6,13869.01,5229.8,6337.89,7625.31,6778.63,3571.35,8976.21,4435.5,10620.27,7116.74,5565.55,3594.06,6686.16,5243.01,6626.54,4253.17,7170.66,4841.57,6887.37,3031.04,4404.51,6407.67,4988.94,3400.02,5184.59,5424.52,3912.44,3290.3,2602.15,4075.19,3794.04,3895.85,4410.16,4687.28,3777.76,1962.88,2361.72,3364.02,3392.26,2677.28,4148.95,2696.98,3796.92,1987.54,3041.01,3324.46,2837.11,3994.59,1536.03,2763.59 +1987,1111341.73,816792.74,245402.86,145051.66,123083.21,171728.92,78054.09,141273.49,56709.21,130781.87,57089.45,56024.51,98343.81,88101.63,98357.47,51221.15,80503.05,39057.21,63881.3,37538.98,51100.88,53964.41,53066.23,41690.83,57313.31,43493.28,51197.48,31326.47,38673.24,32859.37,23312.07,40032.72,26569.77,35230.25,25775.16,32324.33,31529.82,19565.68,19431.99,14682.65,10475.04,23497.77,23917.9,10327.73,9942.29,17480.92,9971.23,16254.94,16647.95,20451.71,19634.01,23198.53,11502.76,10068.28,10091.63,8960.47,8421.06,7639.35,9352.15,16173.1,19029.88,12594.15,16525.11,13020.86,21119.32,10269.57,16290.15,18302.59,18437.79,15695.41,13481.41,8067.76,6531.89,11218.27,8196.04,10628.59,10696.27,8140.07,5117.28,10647.75,4946.37,4758.15,7692.25,9504.13,5327.23,5219.96,11360.85,4186.97,6670.57,8639.11,14525.66,5381.48,6464.73,8120.09,6893.5,3571.76,9293.28,4478.52,11019.65,7323.97,5666.08,3627.33,6893.9,5408.09,6764.62,4289.1,7385.69,4987.74,7018.11,3049.97,4535.26,6732.13,5053.71,3428.83,5303.2,5516.36,4034.67,3310.74,2626.9,4158.61,3905.16,3997.7,4541.8,4843.95,3885.44,1978.52,2379.28,3474.08,3451.94,2672.17,4216.18,2755.24,3882.94,1992.27,3121.34,3407.32,2924.35,4109.7,1551.21,2819.6 +1988,1132866.34,834489.32,247739.58,146035.58,123594.97,175000.92,78354.02,144001.54,56862.89,134261.53,57103.13,56348.21,100975.32,90450.28,101474.84,51339.58,82083.92,39141.61,65313.71,37691.9,52053.7,54912.33,54547.3,42091.03,58845.21,44932.06,53075.62,31795.52,39362.14,33320.96,23408.42,40925.81,26947.84,35970.54,26206.73,32975.54,32444.16,19860.22,19708.32,14774.04,10435.05,23961.82,24591.49,10331.22,9959.02,17873.67,9956.67,16527.43,16882.19,20909.9,19869.95,23909.95,11540.89,10115.85,10146.63,8936.56,8464.79,7658.84,9397.31,16346.21,19501.23,12809.03,17027.59,13002.55,21871.44,10379.55,16651.81,18866.32,18912.53,16262.53,13854.21,8180.73,6572.27,11513.97,8433.8,11051.5,11031.82,8356.31,5122.44,10952.4,4961.09,4773.68,7871.46,9739.18,5380.91,5244.04,11719.07,4205.97,6813.35,8842.58,15155.22,5554.88,6592.79,8636.94,7010.04,3569.11,9604.3,4508.06,11424.26,7540.25,5741.24,3658.49,7107.98,5584.34,6903.32,4320.93,7600.07,5135.96,7165.3,3069.1,4672.85,7030.18,5119.04,3457.21,5421.46,5609.69,4152.98,3332.79,2649.82,4239.3,4020.3,4100.73,4676.36,5001.11,3994.33,1991.43,2393.53,3581.93,3504.65,2668.59,4291.58,2812.24,3968.45,1995.51,3197.07,3493.89,3006.36,4220.42,1564.19,2876.89 +1989,1153565.89,852270.03,250113.19,146890.54,124068.9,178233.22,78715.61,146691.98,57023.75,137750.97,57113.16,56663.2,103599.23,92844.35,104588.49,51423.53,83697.89,39221.74,66757.4,37827.02,52992.43,55795.11,56006.57,42495.73,60391.87,46458.91,54777.11,32263.56,40014.86,33781.32,23472.58,41849.14,27326.17,36740.88,26639.72,33624.44,33465.44,20153.25,19993.76,14868.66,10402.13,24421.19,25257.67,10336.17,9980.18,18293.51,9948.57,16792.29,17110.71,21368.86,20095.75,24660.58,11777.61,10174.17,10189.62,8897.45,8514.21,7686.57,9450.77,16474.56,19979.13,13025.8,17535.97,13059.61,22633.02,10486.51,17040.19,19475.61,19388.34,16846.09,14233.87,8309.53,6619.16,11827.24,8689.15,11523.27,11372.16,8579.82,5130.36,11268.66,4977.12,4780.63,8050.93,9977.38,5412.15,5266.8,12080.44,4226.12,6956.3,9050.47,15755.94,5770.65,6723.05,9094.67,7126.89,3567.29,9902.54,4506.65,11839.24,7769,5777.5,3682.56,7328.6,5768.09,7042.94,4346.45,7813.81,5280.02,7298.42,3088.98,4820.02,7216.03,5185.94,3486.63,5542.85,5699.22,4249.47,3360.96,2663.92,4311.87,4138.41,4205.65,4814.14,5149.8,4103.91,2000.87,2407.72,3686.37,3536.47,2676.61,4384.36,2872.67,4055.27,1996.87,3253.66,3585.87,3071.77,4328.91,1569.97,2940.26 +1990,1172445.2,870133.48,252529.95,147564.07,124515.56,181436.82,79118.33,149352.15,57183.33,141238.02,57127.12,56960.84,106188.64,95269.99,107678.61,51464.35,85357.87,39306.1,68209.61,37954.55,53921.7,56582.82,57412.22,42923.13,61947.35,48086.52,56226.19,32729.74,40626.25,34237.27,23489.37,42818.44,27692.68,37560.53,27071.46,34271.57,34614.58,20439.6,20293.05,14965.45,10377.65,24879.14,25912.37,10341.17,10006.54,18749.41,9953.33,17041.43,17329.71,21826.66,20311.7,25459.6,12249.11,10248.54,10216.85,8841.37,8567.38,7723.95,9517.67,16540.26,20462.46,13242.13,18038.32,13247.65,23402.51,10582.08,17469.01,20147.59,19861.96,17438.91,14628.26,8465.19,6674.89,12171.44,8973.34,12057.04,11715.22,8811.03,5141.12,11598.63,4996.22,4776.37,8232.8,10218.09,5410.37,5288.45,12446.17,4247.29,7099.73,9263.81,16326.82,6041.09,6856.24,9437.55,7242.77,3569.26,10183.11,4463.42,12267.75,8012.86,5768.48,3696.03,7555.62,5956.86,7183.65,4364.12,8027.25,5415.42,7397.35,3109.99,4978.5,7235.8,5254.98,3517.98,5669.95,5781.46,4312.25,3398.17,2664.43,4372.89,4258.47,4313.06,4955.33,5283.73,4213.74,2006.48,2424.24,3786.94,3538.17,2703.02,4500.48,2939.78,4144.57,1996.23,3281.45,3683.97,3113.31,4436.66,1565.24,3012.95 +1991,1189183.52,888054.88,254974.82,148035.9,124940.02,184615.98,79564.45,151976.58,57339.45,144723.57,57152.52,57237.7,108727.43,97726.32,110730.42,51462.61,87071.51,39403.2,69670.9,38077.21,54840.53,57258.4,58752.39,43378.1,63508.46,49821.08,57375.58,33193.92,41190.16,34687.43,23454.21,43840.47,28044.61,38437.86,27501.31,34916.77,35914.83,20718.04,20609.15,15064.52,10363.56,25336.86,26554.33,10346.75,10039.45,19245.05,9974.2,17271.87,17539.63,22283.13,20518.8,26315.01,12993.66,10341.89,10226.49,8766.89,8625.14,7772.7,9602.94,16540.02,20951.39,13457.24,18529.45,13591.97,24179.6,10663.59,17942.72,20893.63,20332.08,18040.44,15039.51,8652.51,6741.07,12553.45,9286.98,12661.61,12060.73,9050.08,5154.91,11942.82,5019.13,4760,8417.68,10460.99,5370.65,5309.24,12815.22,4269.47,7243.39,9483.27,16867.83,6374.33,6992.52,9641.15,7357.96,3576.38,10443.04,4371.6,12710.01,8272.98,5705.38,3697.45,7789.65,6150.08,7325.62,4373.58,8239.73,5542.05,7455.94,3132.05,5149.5,7051.76,5326.66,3551.59,5803.5,5854.09,4337.24,3446.17,2649.15,4419.63,4380.07,4423.01,5099.95,5400.71,4323.41,2007.88,2443.69,3882.27,3505.25,2752.46,4643.2,3014.62,4236.8,1993.3,3275.43,3789.19,3127.3,4544.29,1548.22,3096.57 +1992,1204003.9,906021.11,257454.27,148318.03,125341.48,187766.09,80046.16,154564.28,57494.54,148195.56,57188.89,57496.97,111221.94,100221.56,113747.14,51417.43,88828.31,39515.51,71130.45,38191.78,55748.88,57837.88,60035.54,43853.45,65075.49,51647.77,58260.74,33655.15,41711.47,35131.66,23373.36,44902.93,28384.35,39360.23,27929.03,35558.68,37346.15,20989.13,20937.4,15165.86,10359.14,25791.49,27181.09,10353.32,10078.1,19773.77,10009.58,17485.95,17740.64,22737.06,20716.68,27219.62,13981.23,10451.42,10219.92,8676.24,8686.74,7830.93,9701.33,16481.11,21443.45,13671.03,19012.72,14071.23,24963.95,10733.36,18458.19,21701.48,20799.08,18652.89,15463.85,8868.26,6815.38,12968.35,9621.5,13325.58,12408.93,9297.11,5171.65,12301.34,5044.93,4732.86,8603.23,10705.67,5298.09,5329.03,13187.09,4292.68,7386.98,9708.54,17378.83,6758.84,7131.71,9729.72,7471.93,3588.59,10682.87,4239.15,13163.02,8549.42,5599.81,3688.1,8029.73,6349.09,7468.55,4375.31,8452.28,5661.14,7488.54,3154.86,5331.8,6701.85,5400.33,3586.88,5942.77,5919.29,4331.33,3503.18,2620.25,4454.19,4502.36,4535.52,5247.84,5502.98,4432.74,2005.51,2465.36,3973.33,3442.81,2821.86,4808.86,3095.81,4331.28,1988.66,3240.59,3899.84,3118.58,4651,1520.81,3189.43 +1993,1217129.13,924057.82,260020.19,148431.37,125717.52,190879.52,80520.09,157132.68,57650.47,151614.72,57226.65,57744.82,113695.14,102761.74,116749.56,51317.64,90600.45,39638.64,72560.43,38295.47,56653.73,58364.89,61275.6,44340.56,66650.25,53532.96,58991.22,34110.92,42209.78,35569.28,23256,45979.22,28711.3,40300.16,28353.14,36195.17,38833.6,21251.94,21265.83,15268.01,10359.59,26237.42,27786.26,10358.97,10118.67,20321.18,10055,17687.48,17928.58,23184.23,20902.88,28149.33,15095.1,10569.42,10200.51,8575.89,8746.78,7892.09,9795.55,16370,21931.69,13882.67,19494.97,14637,25754.11,10794.14,19011.92,22535.94,21263.44,19275.42,15896.43,9105.47,6891.6,13403.73,9968.28,14017.24,12758.88,9552.48,5190.74,12675.46,5071.79,4697.54,8784.89,10951.2,5204.63,5347.21,13564.17,4316.79,7530.71,9938.69,17859.75,7163.24,7273.83,9755.86,7581.76,3605.67,10905.76,4088,13622.73,8842.42,5490.92,3670.65,8272.17,6555.6,7611.47,4369.65,8669.17,5771.4,7519.81,3178.16,5521.76,6299.91,5474,3622.71,6087.29,5983.36,4307.3,3564.46,2582.65,4483.84,4623.28,4651.17,5398.81,5594.11,4541.9,2000.77,2488.78,4064.93,3363.1,2900.85,4987.73,3181.22,4426.58,1984.03,3189.58,4010.79,3099.05,4755.29,1488.06,3287.59 +1994,1228991.93,942204.25,262741.57,148411.86,126063.17,193945.27,80931.04,159705.12,57810.77,154931.79,57251.37,57990.66,116182.27,105355.78,119769.56,51149.54,92349.15,39764.56,73925.08,38383.52,57564.13,58901.67,62495.75,44825.89,68236.23,55431.12,59725.13,34558.12,42712.22,35999.96,23116.18,47032.27,29026.19,41218.9,28772.1,36823.54,40273.7,21506.02,21577.98,15369.12,10358.12,26667.05,28362.25,10361.24,10155.92,20867.13,10103.44,17882.74,18098.35,23619.36,21074.05,29070.62,16172.72,10684.82,10173.36,8474.67,8798.23,7947.21,9861.88,16217.49,22407.16,14091.39,19986.89,15217.04,26548.49,10850.59,19597.24,23347.89,21726.35,19907.63,16330.17,9353.39,6961.26,13841.3,10315.38,14692.69,13109.66,9816.59,5211.33,13066.54,5097.1,4657.8,8956.6,11196.48,5107.67,5362.9,13949.7,4341.62,7674.91,10172.3,18311.09,7544.29,7418.86,9796.98,7683.71,3626.8,11116.95,3948.82,14083.61,9151.76,5431.74,3649.13,8512.17,6772.13,7753.05,4357.28,8896.11,5871.61,7583.95,3201.61,5714.22,6005.1,5544.95,3657.41,6235.97,6055.09,4283.62,3623.28,2543.62,4518.8,4740.38,4770.61,5552.63,5679.83,4651.23,1995.61,2513.05,4163.64,3283.66,2974.64,5166.07,3267.67,4520.73,1981.7,3140.6,4115.1,3085.44,4855,1457.09,3385.35 +1995,1239940,960482.8,265658.85,148288.84,126375.47,196957.85,81240.6,162296.61,57978.32,158115.31,57255.19,58241.79,118706.87,108011.47,122829.15,50905.68,94045.58,39890.5,75198.98,38452.63,58486.38,59491.79,63714.39,45298.61,69835.72,57309.88,60575.64,34994.81,43237.79,36423.24,22965.11,48034.77,29330.07,42088.17,29185.31,37441.98,41595.74,21751.68,21862.3,15467.85,10349.84,27075.23,28904.3,10358.46,10186.31,21396.38,10149.74,18076.74,18247.12,24038.76,21228.51,29960.78,17099.54,10789.38,10142.31,8379.18,8836.42,7990.12,9884.15,16034.96,22863.26,14296.61,20495.6,15759.13,27346.46,10906.04,20208.39,24102.99,22188.67,20550.29,16760.47,9604.45,7018.53,14268.99,10653.56,15320.65,13460.99,10089.88,5232.7,13475.4,5119.01,4616.76,9113.98,11440.58,5020.06,5375.47,14345.49,4367,7819.81,10408.49,18735.84,7871.17,7566.71,9909.09,7775.43,3651.55,11320.35,3843.71,14540.82,9477.33,5459.52,3626.6,8746.61,7000.72,7892.42,4339.08,9137.08,5962.06,7704.89,3224.8,5905.56,5928.08,5611.12,3689.64,6387.97,6140.1,4274.82,3674.94,2508.47,4566.09,4851.92,4894.28,5709.05,5764.71,4760.85,1991.47,2537.44,4274.02,3217.34,3033.39,5333.72,3352.77,4612.23,1983.25,3106.74,4207.84,3090.16,4948.8,1433.02,3478.78 +1996,1249981.44,978893.22,268803.42,148073.92,126654.02,199914.83,81426.34,164913.31,58156.75,161147.69,57229.67,58499.69,121269.65,110732.9,125938.34,50582.46,95687.45,40009.32,76372.72,38503.61,59423.21,60151.47,64933.46,45757.36,71446.11,59155.15,61583.09,35419.68,43793.31,36838.91,22804.17,48978.11,29622.13,42898.52,29593.03,38049.04,42770.54,21989.03,22113.55,15563.26,10332.99,27460.6,29411.42,10350.57,10208.28,21903.38,10192.39,18271.69,18372.12,24441.07,21364.13,30811.85,17822.88,10880.91,10109.07,8291.87,8859.19,8017.86,9855.51,15820.63,23298.91,14497.83,21023.32,16248.23,28147.73,10961.01,20845.89,24786.19,22650.1,21202.12,17185.61,9856.81,7061.17,14682.28,10980.27,15889.45,13812.47,10372.75,5254.86,13902.69,5136.98,4574.89,9256.04,11683.48,4945.57,5384.79,14755.29,4393.21,7965.55,10646.67,19131.58,8132.55,7717.44,10109.79,7855.55,3679.85,11518.26,3780.38,14995.25,9819.96,5591.11,3603.75,8974.08,7241.13,8029.11,4314.7,9394.3,6041.11,7892.39,3248.04,6094.26,6115.17,5671.93,3719.32,6543.4,6241.81,4282.35,3717.35,2478.83,4629.4,4957.18,5022.44,5867.85,5849.54,4870.69,1988.99,2561.99,4398.24,3168.22,3070.96,5488.08,3435.82,4700.78,1989.44,3092.23,4287.34,3116.38,5035.88,1417.74,3566.13 +1997,1259066.97,997405.32,272136.55,147768.37,126902.63,202826.47,81499.87,167545.16,58346.67,164048.57,57184.9,58763.4,123854.64,113522.71,129086.99,50189.31,97281.74,40131.56,77453.34,38538.73,60372.5,60863.51,66151.12,46202.6,73064.76,60976.45,62710.56,35833.97,44371.53,37246.27,22632.81,49868.1,29903.15,43657.02,29997.11,38645.41,43830.15,22219.95,22335.64,15655.48,10309.02,27825.9,29886.84,10338.59,10223.79,22389.8,10232.88,18467.59,18476.51,24827.41,21483.51,31635.25,18381.61,10961.12,10073.06,8211.73,8868.85,8032.88,9783.39,15581.75,23715.49,14694.84,21565.33,16701.35,28954.11,11013.98,21509.29,25410.45,23110.18,21865.93,17608.81,10114.09,7091.77,15088.98,11295.88,16408.95,14165.42,10665.55,5277.71,14347.85,5151.75,4532.4,9384.15,11924.99,4881.25,5391.2,15177.46,4420.26,8111.95,10887.63,19505.58,8337.99,7870.86,10381.86,7925.67,3712.46,11710,3752.43,15445.99,10180.06,5814.01,3579.9,9196.53,7493.25,8163.47,4284.93,9666.58,6112.1,8137.48,3271.01,6281.64,6522.38,5727.76,3746.36,6702.43,6356.47,4305.46,3752.36,2453.3,4705.64,5056.52,5154.91,6028.88,5934.28,4980.34,1987.96,2586.83,4534.55,3133.09,3092.67,5631.54,3517.31,4786.64,1999.6,3093.04,4355.11,3161.35,5117.27,1409.67,3649.14 +1998,1267441.5,1015974.04,275542.6,147381,127127.27,205715.54,81499.23,170170.64,58544.94,166870.57,57153.55,59034.92,126447.97,116385.75,132253.26,49748.58,98821.46,40292.64,78452.9,38558.05,61329.59,61597.28,67378.06,46627.44,74693.7,62794.15,63900.63,36241.59,44959.94,37643.88,22458.08,50709.61,30178.13,44372.11,30401.38,39234.06,44840.53,22449.02,22537.34,15745.65,10280.15,28175.26,30335.73,10323.49,10237.3,22856.31,10272.29,18665.05,18570.7,25199.75,21595.37,32451.71,18864,11030.39,10033.06,8137.56,8870.85,8041.05,9683.61,15349.9,24112.41,14887.76,22113.46,17136.78,29769.8,11064.1,22190.25,26003.54,23569.45,22551.79,18036.49,10380.84,7115.13,15504.32,11597.74,16896.21,14523.57,10968.72,5300.38,14808.79,5164.24,4491.9,9499.4,12163.89,4825.19,5395.25,15602.21,4447.49,8258.48,11133.5,19882.46,8503.3,8026.25,10704.74,7989.9,3750.57,11893.27,3750.49,15884.55,10556.55,6099.92,3555.13,9418.39,7759.26,8296.38,4253.28,9950.22,6186.35,8422.37,3292.14,6470.27,7059.81,5778.71,3769.34,6864.39,6473.84,4353.65,3783.97,2430.55,4786.77,5150.76,5291.18,6192.03,6021.69,5089.31,1987.98,2611.37,4679.02,3108.68,3113.95,5765.21,3597.39,4869.63,2012.06,3103.76,4413.48,3224.22,5195.5,1406.35,3731.44 +1999,1275406.7,1034539.21,278862.28,146919.74,127335.86,208612.56,81481.61,172759.24,58746.73,169688.53,57179.99,59316.11,129029.69,119327.07,135405.58,49290.79,100300.58,40539.72,79391.37,38561.69,62287.33,62306.65,68626.66,47023.88,76335.81,64640.05,65062.66,36648.07,45539.44,38030.11,22288.06,51512.84,30453.88,45058.78,30810.66,39819.28,45898.67,22682.24,22731.99,15835.52,10249.96,28514.8,30765.61,10306.65,10255.16,23305.99,10312.61,18864.11,18669.1,25561.3,21711.57,33291.54,19403.68,11090.27,9986.93,8066.88,8873.1,8051.11,9579.31,15165.57,24489.77,15076.95,22656.29,17584.87,30600.4,11110,22878.16,26607.04,24028.69,23273,18477.61,10663.72,7138.46,15949.77,11883.64,17378.1,14891.89,11282.7,5321.77,15282.52,5175.93,4456.43,9603.74,12398.69,4773.17,5397.71,16016.09,4474,8404.4,11387.2,20294.41,8653.77,8182.71,11044.36,8054.13,3795.7,12064.54,3759.12,16300.23,10947.83,6405.86,3529.18,9645.96,8041.85,8429.11,4224.27,10239.71,6278.94,8720.23,3309.32,6664.1,7593.24,5825.19,3786.58,7028.44,6579.56,4437.8,3818.13,2408.09,4861.13,5241.28,5430.48,6357.22,6114.89,5196.94,1988.32,2634.88,4825.7,3089.02,3156.65,5891.77,3676.51,4949.66,2024.39,3115.58,4466.13,3302.26,5274.16,1403.79,3818.46 +2000,1283198.97,1053050.91,281982.78,146396.51,127533.93,211540.43,81487.76,175287.59,58950.85,172558.87,57293.72,59608.2,131581.24,122352.01,138523.29,48840.07,101719.67,40903.71,80285.56,38550.5,63240.12,62958.02,69905.99,47386.31,77991.57,66537.33,66131.85,37057.45,46095.46,38403.99,22128.13,52286.08,30735.77,45728.32,31229.42,40403.96,47076.39,22924.84,22929.08,15926.19,10221.05,28849.62,31183.66,10289.59,10282.03,23740.91,10355.12,19065.84,18781.94,25914.88,21840.24,34178.04,20093.76,11142.12,9933.61,7997.79,8881.64,8069.28,9487.62,15057.36,24848.92,15262.75,23185.61,18067.69,31450.48,11150.74,23565.41,27250.54,24488.34,24039.27,18938.76,10967.69,7167.25,16440.92,12152.35,17874.73,15274.23,11607.94,5341.19,15766.81,5187.95,4428.07,9699.2,12628.6,4722.06,5399.21,16410.85,4499.37,8549.2,11650.74,20764.31,8808.55,8339.51,11376.17,8122.74,3848.78,12222.25,3766.71,16686.56,11352.97,6700.66,3501.84,9884.05,8342.56,8562.62,4201.09,10531.22,6400.71,9011.48,3321.25,6865.95,8025.7,5867.63,3796.98,7193.99,6663.65,4564.3,3859,2384.16,4920.72,5329.3,5572.22,6524.28,6216.21,5302.7,1988.5,2656.86,4970.37,3069.59,3235.37,6013.74,3754.99,5026.8,2034.82,3121.97,4516.13,3392.8,5355.75,1399.11,3913.96 +2001,1290937.65,1071477.86,284852.39,145814.86,127723.51,214506.5,81535.85,177750.67,59149.34,175508.85,57506.37,59911.25,134107.16,125463.43,141601.44,48402.89,103067.07,41392.1,81139.92,38524.01,64191.47,63543.32,71226.94,47712.79,79665.32,68492.26,67096.41,37471.51,46627.99,38765.02,21984.14,53032.42,31024.73,46385.01,31656.85,40988.91,48394.34,23176.47,23131.81,16018.11,10194.01,29181.83,31592.15,10271.2,10319.02,24161.78,10399.31,19268.53,18913.05,26261.36,21985.28,35117.02,20966.46,11184.4,9872.96,7930.55,8897.79,8097.75,9413.42,15039.97,25188.53,15444.97,23698.91,18588.76,32321.48,11186.54,24251.65,27945.01,24948.48,24854.89,19421.61,11293.26,7202.69,16983.27,12402.47,18390.14,15671.93,11944.59,5358.06,16260.93,5200.48,4408.08,9785.7,12852.76,4672.16,5400.01,16766.9,4523.15,8692.57,11924.95,21303.59,8971.14,8496.38,11695.86,8197.28,3909.29,12366.17,3771.28,17040.15,11771.98,6974.44,3473.62,10134.5,8663.01,8697.13,4185.51,10824.13,6555.83,9290.82,3327.1,7076.73,8329.41,5905.96,3799.91,7360.55,6722.31,4739.15,3907.93,2358.68,4962.59,5414.57,5716.15,6693.06,6327.13,5406.62,1988.35,2677.01,5111.77,3050.66,3359.86,6129.53,3832.2,5100.75,2042.84,3122.41,4564.08,3497.12,5440.57,1391.73,4018.78 +2002,1298646.58,1089807.11,287506.85,145194.88,127902.62,217508.06,81611.87,180151.02,59348.95,178532.17,57801.76,60225.08,136600.67,128666.71,144654.14,47979.59,104355.61,41985.52,81956.5,38484.47,65143.05,64073.16,72590.12,48004.79,81352.06,70497.19,67983.33,37889.37,47140.22,39114.07,21854.68,53751.52,31320.89,47026.17,32094.47,41572.49,49835.76,23438.4,23336.68,16110.36,10167.87,29512.37,31995.05,10252.44,10364.89,24566.34,10443.67,19475.19,19059.3,26601.47,22143.34,36105.81,21979.92,11216.8,9807.13,7866.09,8920.69,8134.41,9352.71,15105.65,25511.88,15623.64,24198.81,19139.66,33214.01,11218,24939.3,28679.57,25408.7,25718.05,19924.52,11638.93,7243.91,17572.65,12634.73,18919.18,16084.89,12293.1,5372.8,16765.12,5213.37,4395.82,9864.33,13072.06,4624.88,5400.15,17087.9,4546.02,8834.73,12208.85,21906.31,9137.35,8653.35,12013.71,8277.02,3976.29,12500.53,3775.81,17366.52,12206,7237.28,3444.75,10396.86,9001.69,8832.29,4176.62,11120.41,6741.57,9564.17,3327.77,7295.39,8536.21,5940.3,3796.03,7528.31,6759.49,4957.22,3963.21,2332.49,4990.69,5497.27,5862.32,6863.16,6447.69,5508.61,1988.26,2695.45,5251.47,3033.9,3522.84,6239.29,3907.61,5171.73,2048.93,3119.03,4610,3614.64,5527.52,1382.73,4131.21 +2003,1306343.91,1108027.85,290027.62,144584.7,128067.79,220545.21,81686.5,182482.15,59580.22,181612.13,58147.04,60550.1,139019,131972.53,147703.4,47580.56,105640.45,42653.41,82747.66,38438.83,66085.8,64554.95,73981.94,48265.03,83031.95,72545.14,68812.71,38309.38,47624.89,39453.54,21728.92,54445.09,31627.81,47648.73,32551.66,42152.15,51390.03,23719.04,23538.54,16200.95,10141.96,29843.94,32403.51,10239.28,10419.03,24950.62,10487.11,19696.97,19215.31,26937.74,22306.34,37149.07,23064.85,11243.57,9740.05,7803.88,8951.42,8175.85,9302.87,15232.33,25830.66,15799.54,24688.7,19716.6,34130.85,11244.89,25627.63,29435.94,25868.52,26624.82,20446.78,12005.13,7291.89,18203.37,12853.12,19462.09,16513.82,12654.62,5386.97,17279.14,5227.1,4388.95,9939.68,13289.6,4579.08,5399.89,17415.27,4570.11,8976.55,12500.48,22556.43,9309.85,8810.42,12336.69,8361.36,4049.55,12633.9,3779.25,17679.36,12656.87,7501.64,3414.35,10670.99,9353.2,8967.76,4171.65,11421.98,6953.11,9836.4,3325.64,7520.56,8680.35,5971.54,3787.27,7698.11,6782.57,5199.55,4022.07,2305.78,5012.27,5579.66,6010.72,7033.82,6576.88,5607.95,1988.97,2712.51,5391.4,3017.81,3701.46,6349.52,3981.67,5240.88,2053.43,3111.01,4655.74,3738.27,5615.95,1372.91,4249.3 +2004,1314007.48,1126135.78,292539.32,144044.61,128213.63,223614.65,81715.66,184738.46,59884.13,184717.78,58494.7,60886.61,141307.49,135393.62,150780.3,47215.97,106995.58,43351.67,83527.68,38396.07,67007.86,65002.23,75381.9,48498.08,84678.49,74624.41,69617.1,38728.7,48073.71,39786.58,21591.16,55114.73,31949.42,48247.4,33040.25,42724.16,53034.22,24028.87,23729.5,16287.18,10115.08,30179.29,32831.1,10239.54,10480.12,25309.45,10528.27,19948.43,19372.54,27273.19,22462.37,38249.98,24118.98,11270.48,9676.9,7743.26,8990.64,8216.81,9258.34,15385.92,26160.33,15973.78,25174.11,20312.71,35074.93,11266.94,26316.61,30186.34,26327.23,27568.44,20986.54,12391.91,7347.25,18865.72,13063.38,20017.07,16959.08,13030.57,5402.76,17803,5242.19,4383.92,10017.6,13509.65,4533.33,5399.42,17806.64,4598.21,9119.18,12796.93,23228.89,9490.23,8967.74,12676.04,8448.88,4128.75,12777.51,3781.29,17997.74,13127.01,7787.66,3381.12,10955.94,9710.04,9103,4166.34,11731.75,7182.45,10116.23,3324.1,7750,8818.44,6000.78,3776.36,7871.22,6802.6,5439.7,4080.44,2278.86,5037.93,5664.61,6161.52,7204.15,6712.84,5703.74,1991.43,2728.78,5534.6,3000.61,3863.27,6468.83,4055.04,5309.7,2057.05,3097.75,4703.4,3858.62,5704.76,1363.57,4370.04 +2005,1321623.49,1144118.67,295129.5,143618.19,128335.77,226712.73,81671.23,186917.36,60286.75,187832.26,58808.48,61233.9,143431.1,138939.48,153909.67,46892.16,108472.23,44042.63,84308.84,38362.67,67903.41,65425.47,76778.15,48708.5,86274.24,76727.08,70421.81,39145.49,48482.61,40115.79,21431.02,55764.05,32287.98,48820.59,33567.64,43285.63,54751.48,24374.26,23904.17,16367.16,10086.47,30521.07,33288.44,10258.18,10546.89,25640.29,10565.72,20238.9,19524.56,27610.41,22602.89,39410.55,25070.8,11301.2,9621.54,7683.76,9038.62,8253.65,9214.99,15541.46,26512.18,16147.06,25659.39,20923.07,36048.29,11284.25,27008.43,30911.91,26784.16,28543.94,21542.01,12798.76,7410.31,19552.54,13270.2,20582.93,17420.8,13421.93,5421.7,18336.72,5258.93,4378.06,10102.48,13735.23,4486.55,5398.96,18294.61,4632.36,9263.4,13096.03,23905.65,9679.75,9125.41,13039.71,8538.61,4212.98,12940.03,3781.53,18336.3,13618.45,8108.88,3344.25,11251.27,10067.01,9237.57,4157.7,12052.16,7423.29,10409.93,3325.61,7982.23,8991.74,6028.96,3765.4,8048.47,6827.76,5658.38,4135.36,2251.98,5075.38,5754.03,6314.71,7373.43,6854.18,5795.49,1996.28,2744.67,5683.27,2981.26,3986.85,6602.97,4127.91,5379.33,2060.27,3079.18,4754.64,3969.01,5792.69,1355.65,4491.04 +2006,1329209.09,1161977.72,297827.36,143322.93,128432.99,229838.2,81540.35,189012.41,60802.8,190943.29,59079.87,61592.88,145368,142614.09,157093.99,46616.39,110092.38,44725.53,85094.62,38341.04,68763.41,65824.16,78159.05,48895.85,87809.42,78850.69,71227.88,39558.89,48846.47,40441.53,21243.92,56390.73,32645.48,49364.58,34138.5,43835.72,56543.01,24759.52,24061.1,16440.22,10055.9,30869.35,33777.92,10298.29,10619.57,25940.62,10599.95,20573.61,19670.15,27949.94,22725.48,40634.95,25893.45,11338.77,9575.04,7624.96,9096.26,8284.74,9173.29,15690.86,26890.26,16319.79,26143.57,21547.46,37052.05,11296.23,27697.91,31607.06,27239.17,29550.66,22113.43,13227.06,7481.41,20262.4,13474.49,21160.53,17899.56,13829.18,5444.3,18880.27,5277.57,4370.7,10196.14,13967.48,4438.9,5398.54,18914.98,4673.07,9409.46,13397.01,24578.3,9881.43,9283.33,13429.26,8630.2,4303.37,13124.27,3779.47,18699.44,14132.06,8468.15,3303.67,11556.76,10421.6,9371.34,4144.73,12383.45,7675.34,10718.32,3331.04,8216.9,9206.58,6056.48,3755.11,8230.17,6860.38,5848.69,4185.92,2225.22,5126.73,5849.36,6470.27,7541.41,7000.56,5882.8,2003.79,2760.28,5837.79,2958.5,4057.35,6755.21,4201.76,5450.21,2063.15,3054.33,4810.11,4066.65,5881.44,1349.3,4611.9 +2007,1336800.51,1179681.24,300595.18,143150.14,128505.25,232989.14,81344.46,191026.64,61414.66,194076.39,59313.51,61960.95,147139.19,146417.02,160332.97,46386.36,111836.35,45393.86,85889.59,38329.59,69597.28,66195.62,79537.25,49062.43,89293.49,81000.41,72031.1,39970.22,49171.59,40763.75,21036.11,57002.23,33019.93,49887.18,34747.5,44374.57,58417.56,25179.82,24203.29,16507.06,10023.89,31225.88,34300.08,10356.59,10697.84,26214.85,10630.12,20946.66,19810.79,28292.72,22833.01,41923.72,26616.79,11380.9,9536.86,7566.95,9163.24,8311.78,9135.21,15841.36,27292.38,16491.69,26625.85,22188.39,38085.91,11303.69,28390.43,32282.53,27691.97,30590.49,22700.21,13675.61,7560.36,20997.69,13676.69,21751.61,18395.39,14252.02,5469.96,19433.52,5298.01,4362.27,10298.09,14205.45,4390.54,5398.33,19632.81,4719.65,9556.89,13700.29,25252.57,10096.73,9441.44,13840.97,8724.3,4398.07,13329.91,3774,19085.94,14668.34,8856.8,3260.09,11873.56,10775.71,9504.35,4128.65,12725.97,7939.57,11038.6,3339.74,8454.79,9447.4,6083.48,3745.26,8415.76,6898.4,6015.42,4233.15,2198.52,5189.72,5949.79,6627.92,7707.97,7152.39,5966.16,2013.54,2775.47,5997.39,2933.06,4086.47,6922.69,4275.8,5522.11,2065.46,3023.91,4870.14,4153.33,5970.36,1344.07,4732.53 +2008,1344415.23,1197146.91,303374.07,143083.19,128550.51,236159.28,81130.94,192979.03,62076.22,197295.99,59502.39,62329.57,148805.81,150347.39,163644.6,46186.43,113661.81,45997.81,86707.8,38325.69,70440.03,66545.76,80953.88,49218.99,90751.86,83184.89,72845.54,40382.39,49479.75,41083.98,20821.08,57623.68,33404.55,50412.13,35379.77,44901.54,60373.61,25620.82,24335.15,16568.73,9991.2,31596.86,34860.72,10423.79,10779.17,26475.86,10652.1,21342.49,19945.83,28641.98,22928.53,43270.14,27294.03,11419.65,9507.33,7510.33,9236.89,8338.45,9099.55,16001.18,27715.83,16661.94,27111.07,22846.76,39148.42,11309.75,29111.42,32955.5,28141.7,31663.9,23298.64,14138.22,7646.11,21759.42,13880.51,22356.39,18907.01,14689.73,5497.8,19996.47,5319.9,4352.43,10407.34,14447.56,4340.51,5398.86,20325.44,4771.41,9705.03,14006.37,25940.77,10323.14,9599.86,14271.23,8821.91,4489.59,13558.47,3763.6,19497.99,15228.53,9263.14,3214.63,12203.96,11133.86,9636.52,4111.44,13082.52,8212.26,11369.28,3350.82,8696.92,9708.17,6110.3,3735.79,8604.13,6940.32,6165.37,4278.33,2171.88,5262.29,6052.19,6787.19,7872.66,7309.73,6047.12,2024.54,2790.12,6161.8,2908.22,4111.05,7097.41,4345.39,5594.51,2067.38,2991.65,4935.76,4232.64,6053.08,1339.71,4851.11 +2009,1352068.09,1214270.13,306076.36,143092.77,128566.66,239340.48,80965.61,194896,62722.61,200688.23,59641.85,62687.52,150454.71,154402.18,167049.58,45994.02,115505.23,46476.07,87565.41,38324.87,71339.19,66881.87,82465.02,49379.21,92220.88,85416.25,73687.57,40799.41,49800.69,41403.98,20617.96,58287.62,33789.83,50970.82,36014.27,45416.18,62409.44,26062.52,24463.02,16626.93,9958.94,31989.9,35465.76,10486.85,10860.29,26741.1,10660.62,21739.04,20075.09,29001.51,23017.06,44664.23,28004.33,11443.83,9486.24,7456.03,9313.48,8370.04,9064.76,16184.16,28155.31,16829.44,27605.38,23524.06,40237.2,11318.6,29894.65,33650.62,28587.32,32771.9,23903.83,14606.6,7737.16,22549.55,14090.21,22974.93,19432.54,15141.1,5526.45,20569.12,5342.65,4341.08,10521.83,14691.28,4287.7,5400.72,20824.89,4827.18,9852.87,14316.21,26661.49,10556.52,9758.75,14714.6,8924.41,4568.06,13810.6,3746.56,19936.37,15813.91,9670.67,3168.77,12550.92,11502.79,9767.76,4096.04,13456.42,8489.03,11707.99,3362.76,8944.71,9977.45,6137.28,3726.33,8793.69,6983.27,6310.26,4323.48,2145.31,5340.78,6152.04,6947.45,8035.02,7472.82,6127.84,2035.42,2804.08,6330.47,2888.58,4183.16,7267.9,4404.23,5666.58,2069.09,2962.64,5007.95,4310.33,6121.05,1335.83,4965.52 +2010,1359755.1,1230980.69,308641.39,143153.87,128551.87,242524.12,80894.79,196796.27,63306.84,204310.19,59729.81,63026.74,152149.1,158578.26,170560.18,45792.5,117318.94,46788.63,88472.51,38323.4,72326.91,67208.81,84107.61,49552.86,93726.62,87702.67,74567.51,41223.89,50155.9,41724.73,20440.35,59016.09,34168.67,51584.66,36635.57,45918.1,64523.26,26490.13,24591.6,16682.92,9927.84,32409.64,36117.64,10536.29,10938.74,27023.14,10652.32,22120.06,20198.35,29373.65,23102.41,46098.59,28803.17,11446.01,9473.07,7404.59,9390.17,8409.95,9029.72,16398.98,28606.29,16993.35,28112.29,24221.41,41350.15,11333.05,30762.7,34385.96,29028.03,33915.13,24512.1,15075.09,7831.97,23369.13,14308.74,23606.78,19970.5,15605.22,5554.84,21151.64,5365.78,4328.15,10639.93,14934.69,4231.66,5404.29,21018.83,4885.88,9999.62,14630.42,27425.68,10794.17,9918.24,15167.1,9032.46,4626.93,14086.32,3722.08,20401.33,16425.58,10067.19,3123.8,12916.23,11887.2,9897.99,4084.48,13850.03,8766.93,12053.22,3374.42,9199.26,10246.84,6164.63,3716.69,8983.31,7025.22,6458.72,4370.06,2118.85,5422.34,6246.27,7108.24,8194.78,7641.63,6209.88,2045.17,2817.21,6502.95,2877.31,4337.14,7425.96,4448.53,5737.72,2070.74,2940.53,5087.21,4390.84,6169.14,1332.1,5074.25 +2011,1367480.26,1247236.03,311051.37,143263.99,128505.4,245707.51,80933.98,198686.69,63811.88,208185.87,59759.93,63343.58,153911.92,162877.08,174184.27,45576.31,119090.02,46909.14,89436.64,38320.95,73409.46,67530.13,85897.56,49744.66,95277.94,90046.76,75491.58,41656.88,50553.03,42047.2,20292.97,59819.9,34538.62,52263.52,37237.88,46406.65,66713.6,26898.16,24722.3,16737,9898.2,32858.82,36819.56,10568.72,11013.08,27327.15,10625.33,22479.8,20315.02,29759.99,23185.01,47570.9,29708.6,11422.81,9468.5,7356.21,9465.89,8459.86,8993.75,16647.38,29068.22,17153.36,28635.13,24939.01,42486.84,11354.65,31727.05,35167.31,29463.29,35093.65,25121.8,15540.99,7930.42,24218.57,14537.89,24252.21,20520.45,16081.9,5582.87,21743.95,5389.24,4313.37,10761.47,15177.36,4171.26,5409.85,20863.99,4947.6,10145.05,14948.92,28238.02,11035.17,10078.34,15627.62,9146.38,4662.58,14386.65,3688.87,20895.31,17064.64,10448.86,3079.88,13300.91,12288.65,10027.1,4077.37,14264.76,9043.51,12404.73,3385.62,9460.8,10516.07,6192.56,3706.98,9172.58,7065.82,6611.69,4418.37,2092.49,5506.6,6333.49,7269.35,8351.6,7815.95,6293.78,2053.44,2829.49,6679.28,2875.58,4588.37,7568.77,4476.15,5807.82,2072.38,2926.66,5174.06,4474.69,6193.5,1328.45,5176.02 +2012,1375198.62,1263065.85,313335.42,143420.6,128426.38,248883.23,81066.23,200560.98,64250.33,212283.78,59733.83,63639.87,155727.05,167297.28,177911.53,45349.33,120828.31,46857.4,90451.88,38317.4,74569.87,67843.98,87813.26,49952.24,96866.64,92444.18,76453.57,42096.74,50986.51,42370.63,20171.26,60689.55,34900.71,52998.21,37824.45,46881.48,68978.68,27289.59,24854.03,16789.1,9869.68,33333.79,37565.85,10586.76,11083.55,27649.93,10581.82,22821.85,20424.56,30158.97,23263.9,49083,30696.96,11378.26,9471.73,7310.38,9540.91,8517.55,8956.98,16921.18,29540.82,17309.75,29170.46,25676.61,43646.63,11382.15,32776.57,35990.19,29893.08,36306.8,25733.05,16006.67,8031.67,25096.15,14776.87,24909.97,21082.38,16571.22,5610.66,22346.57,5412.98,4296.53,10886.67,15419.67,4107.72,5417.07,20420.7,5012.01,10289.21,15271.06,29086.36,11281.47,10239,16097.31,9264.92,4678.12,14710.83,3648.2,21418.6,17731.63,10818.26,3037.25,13703.51,12705.14,10154.95,4073.7,14699.94,9319.71,12763.78,3396.78,9729.16,10788.85,6221.25,3697.58,9361.72,7106.4,6766.1,4467.74,2066.37,5594.11,6415.17,7430.84,8505.65,7995.06,6379.22,2060.39,2840.99,6859.48,2881.92,4916.4,7699.11,4490.42,5877.11,2074.04,2920.04,5267.84,4560.98,6198.26,1324.93,5270.96 +2013,1382793.21,1278562.21,315536.68,143597.23,128312.92,252032.26,81265.14,202408.63,64641.11,216536.07,59668,63919.92,157571.29,171829.3,181712.6,45115.79,122535.97,46697.55,91497.73,38309.45,75787.33,68143.07,89807.43,50169.24,98481.03,94887.72,77435.38,42539.93,51448.2,42691.68,20068.2,61601.24,35255.5,53767.4,38399.4,47342.98,71316.03,27668.26,24985.98,16839.7,9841.7,33824.77,38338.56,10594.48,11151.51,27985.31,10527.67,23150.73,20527.23,30565.72,23339.86,50636.6,31731.69,11321.3,9478.9,7266.14,9615.25,8577.78,8920.22,17207.26,30019.73,17462.98,29706.72,26434.37,44826.85,11412.17,33883.15,36849.92,30317.85,37553.73,26346.25,16477.82,8132.67,25998.34,15022.69,25576.32,21655.72,17072.72,5637.82,22961.15,5436.62,4277.81,11014.56,15661.55,4045.91,5425.13,19809.14,5077.1,10431.78,15596.21,29944.48,11536.62,10400.26,16577.15,9385.29,4681.97,15054.51,3605,21966.31,18426.37,11177.49,2997.33,14120.32,13133.59,10281.3,4071.96,15153.21,9600.19,13132.35,3408.01,10004.45,11065.15,6250.78,3688.72,9551.12,7148.57,6922.08,4517.53,2040.78,5683.81,6494.56,7592.87,8657.79,8177.81,6465.74,2066.12,2851.81,7042.95,2893.51,5276.1,7821.11,4499.65,5945.75,2075.74,2918.98,5366.28,4651,6195.97,1321.56,5360.84 +2014,1390110.39,1293859.29,317718.78,143761.38,128162.87,255131.12,81489.66,204213.13,65015.69,220844.88,59585.67,64190.64,159405.28,176460.5,185546.26,44883.43,124221.6,46521.83,92544.92,38293.06,77030.63,68416.77,91812.57,50385.56,100102.25,97366.77,78411.09,42981.52,51924.18,43005.79,19972.74,62520.02,35604.73,54539.57,38970.39,47791.91,73722.86,28041.37,25116.36,16889.36,9813.34,34318.08,39113.31,10598.53,11219.16,28323.24,10471.17,23474.67,20623.56,30973.35,23413.65,52234.87,32758.02,11264.73,9484.85,7222.15,9689.38,8633.22,8884.71,17487.78,30499.62,17613.8,30228.02,27212.38,46024.25,11439.77,35006.08,37737.91,30738.38,38833.34,26962.56,16962.85,8229.63,26920.47,15270.79,26246.33,22239.9,17585.98,5663.91,23589.8,5459.72,4257.53,11143.91,15903.11,3992.35,5432.84,19203.09,5140.31,10572.47,15923.56,30776.72,11805.51,10562.16,17068.84,9503.77,4686.35,15411.68,3566,22531.35,19148.22,11530.97,2961.85,14546.11,13569.44,10405.84,4069.94,15620.97,9891.79,13513.13,3419.55,10286.71,11345.36,6281.19,3680.72,9741.44,7194.56,7079.16,4566.7,2016.13,5774.57,6576.4,7755.79,8809.22,8362.75,6552.58,2070.85,2862.09,7228.92,2906.22,5603.28,7941.33,4515.39,6014,2077.5,2920.78,5466.24,4746.05,6204.11,1318.36,5448.34 +2015,1397028.55,1309053.98,319929.16,143888,127974.96,258162.11,81707.79,205962.11,65397.08,225135.66,59504.21,64457.2,161200.89,181181.74,189380.51,44657.7,125890.95,46397.66,93571.57,38265.23,78271.47,68657.6,93778.17,50593.66,101716.36,99873.03,79360.49,43417.77,52403.67,43309.61,19876.62,63419.64,35949.71,55291.23,39542.98,48228.7,76196.62,28414.09,25243.92,16938.5,9783.93,34803.32,39871.53,10603.76,11287.94,28656.28,10418.47,23799.56,20714.04,31376.67,23485.76,53879.96,33736.49,11217.8,9485.77,7177.4,9763.57,8678.66,8851.28,17749.65,30976.02,17762.68,30723.16,28010.69,47236.26,11461.43,36115.65,38647.8,31155.13,40144.87,27582.82,17467.91,8319.77,27859.31,15517.64,26916.21,22834.52,18110.62,5688.7,24234.09,5481.97,4236.02,11273.66,16144.37,3951.52,5439.32,18734.99,5199.84,10711.06,16252.43,31557.14,12091.53,10724.71,17573.61,9617.48,4700.11,15777.45,3535.96,23108.47,19896.97,11882.14,2931.93,14976.99,14009.41,10528.39,4065.98,16100.59,10199.27,13908.13,3431.55,10575.95,11629.55,6312.48,3673.73,9933.15,7245.7,7237.03,4614.53,1992.66,5865.4,6663.97,7919.83,8960.83,8548.65,6639.12,2074.79,2871.93,7416.8,2916.95,5851.48,8064.55,4546.1,6082.04,2079.31,2923.35,5565.28,4846.98,6234.96,1315.32,5535.26 \ No newline at end of file diff --git a/02-7-population-trends/02-7-pop-past-future-2.7-necessities.csv b/02-7-population-trends/02-7-pop-past-future-2.7-necessities.csv new file mode 100644 index 0000000..c6ad2e8 --- /dev/null +++ b/02-7-population-trends/02-7-pop-past-future-2.7-necessities.csv @@ -0,0 +1,67 @@ +Year,Europe,World,Africa,Asia,Oceania,NorthernAmerica,LatinAmericaandCarribean,,,,,,, +1950,549,2536,229,1404,13,173,169,,,,,,, +1951,554,2584,233,1435,13,175,173,,,,,,, +1952,560,2631,238,1463,13,178,178,,,,,,, +1953,566,2677,243,1491,14,181,183,,,,,,, +1954,571,2724,248,1518,14,184,188,,,,,,, +1955,577,2772,254,1546,14,188,193,,,,,,, +1956,583,2821,260,1575,14,191,198,,,,,,, +1957,589,2872,266,1604,15,195,204,,,,,,, +1958,594,2924,272,1635,15,198,209,,,,,,, +1959,600,2978,278,1667,15,201,215,,,,,,, +1960,606,3033,285,1700,16,205,221,,,,,,, +1961,612,3090,292,1735,16,208,227,,,,,,, +1962,618,3149,299,1771,16,211,234,,,,,,, +1963,624,3210,307,1808,17,214,240,,,,,,, +1964,630,3274,314,1848,17,217,247,,,,,,, +1965,635,3340,322,1891,18,220,253,,,,,,, +1966,640,3408,331,1937,18,222,260,,,,,,, +1967,645,3479,339,1985,18,224,267,,,,,,, +1968,649,3552,348,2035,19,227,274,,,,,,, +1969,653,3626,357,2086,19,229,281,,,,,,, +1970,657,3701,366,2138,20,231,288,,,,,,, +1971,662,3776,376,2189,20,233,295,,,,,,, +1972,666,3852,386,2241,20,236,303,,,,,,, +1973,670,3928,396,2293,21,238,310,,,,,,, +1974,674,4003,407,2344,21,240,318,,,,,,, +1975,678,4079,418,2394,22,242,325,,,,,,, +1976,681,4154,429,2444,22,245,333,,,,,,, +1977,685,4229,441,2493,22,247,341,,,,,,, +1978,688,4304,454,2542,22,250,348,,,,,,, +1979,691,4381,467,2592,23,252,356,,,,,,, +1980,694,4458,480,2642,23,254,364,,,,,,, +1981,697,4538,494,2694,23,257,372,,,,,,, +1982,700,4619,508,2747,24,259,380,,,,,,, +1983,703,4702,522,2802,24,262,389,,,,,,, +1984,705,4786,537,2858,25,264,397,,,,,,, +1985,708,4874,553,2916,25,267,405,,,,,,, +1986,711,4964,569,2976,25,269,413,,,,,,, +1987,714,5056,585,3038,26,272,421,,,,,,, +1988,717,5149,601,3100,26,275,430,,,,,,, +1989,719,5241,618,3162,27,278,438,,,,,,, +1990,722,5331,635,3221,27,280,446,,,,,,, +1991,724,5419,652,3279,28,283,454,,,,,,, +1992,725,5504,669,3334,28,286,462,,,,,,, +1993,727,5588,687,3387,28,289,471,,,,,,, +1994,728,5670,705,3439,29,292,479,,,,,,, +1995,728,5751,723,3489,29,295,487,,,,,,, +1996,728,5832,741,3539,30,299,495,,,,,,, +1997,728,5911,760,3588,30,302,503,,,,,,, +1998,728,5989,779,3636,30,306,511,,,,,,, +1999,727,6067,798,3683,31,309,518,,,,,,, +2000,727,6145,818,3730,31,313,526,,,,,,, +2001,727,6223,838,3777,32,316,533,,,,,,, +2002,728,6302,859,3824,32,319,540,,,,,,, +2003,728,6381,880,3871,33,322,547,,,,,,, +2004,729,6461,902,3918,33,325,555,,,,,,, +2005,730,6542,925,3964,34,328,562,,,,,,, +2006,732,6624,948,4011,34,331,569,,,,,,, +2007,733,6706,972,4057,35,334,576,,,,,,, +2008,734,6790,997,4103,35,337,583,,,,,,, +2009,736,6874,1023,4149,36,340,590,,,,,,, +2010,737,6958,1049,4194,37,343,598,,,,,,, +2011,738,7043,1077,4240,37,346,605,,,,,,, +2012,739,7128,1105,4286,38,348,612,,,,,,, +2013,740,7213,1134,4331,38,351,619,,,,,,, +2014,740,7298,1164,4376,39,353,626,,,,,,, +2015,741,7383,1194,4420,40,356,632,,,,,,, diff --git a/02-7-population-trends/02-7-population-history-x.Rmd b/02-7-population-trends/02-7-population-history-x.Rmd new file mode 100644 index 0000000..8b05377 --- /dev/null +++ b/02-7-population-trends/02-7-population-history-x.Rmd @@ -0,0 +1,59 @@ +--- + title: "Art of Statistics: Figure 2.7 (pages 62) World Population Trends" + #output: md_document + output: html_document +--- + +Data is from +[United Nations Population Division: World Population Prospects 2017](https://population.un.org/wpp/Download/Standard/Population/), specifically [Total Population - Both Sexes](https://population.un.org/wpp/DVD/Files/1_Indicators%20(Standard)/EXCEL_FILES/1_Population/WPP2017_POP_F01_1_TOTAL_POPULATION_BOTH_SEXES.xlsx) + +```{r} +library(ggplot2) + +popTrends<-read.csv("02-7-pop-past-future-2.7-necessities.csv",header=TRUE) # reads data into popTrends data frame + +p <- ggplot(popTrends, aes(x=Year, y)) # constructs initial plot object, p +p <- p + geom_line(aes(y = Europe, col = "Europe"), size=1.5) # adds a y-series +p <- p + geom_line(aes(y = World, col = "World"), size=1.5) # adds a y-series +p <- p + geom_line(aes(y = Africa, col = "Africa"), size=1.5) # adds a y-series +p <- p + geom_line(aes(y = Oceania, col = "Oceania"), size=1.5) # adds a y-series +p <- p + geom_line(aes(y = Asia, col = "Asia"), size=1.5) # adds a y-series +p <- p + geom_line(aes(y = LatinAmericaandCarribean, col = "S America"), size=1.5) # adds a y-series +p <- p + geom_line(aes(y = NorthernAmerica, col = "N America"), size=1.5) # adds a y-series +p <- p + labs(subtitle="(a) On a standard scale", y="Millions", x="Year") # Adds title, subtitle + +# Some of the colours rather pale to see in palette Accent, Set1 seems better +p <- p + scale_colour_brewer(palette = "Set1") +#p <- p + scale_colour_brewer(palette = "Accent") + +p <- p + scale_size_continuous(name = "Size", guide = FALSE) # turns off unwanted size legend +p <- p + theme(legend.position="bottom", legend.box = "horizontal") # positions the colour legend +p +``` + +_Figure 2.7(a) . Total population for the world and continents between 1950 and 2015, both sexes combined._ + +```{r} +# now for the logarithmic version +millIn51 <- read.csv("02-7-pop-past-future million in 51.csv",header=TRUE) # reads data into millIn51 data frame + +logbreaks <- c(1, 5, 10, 50, 100, 500, 1000, 5000, 10000) +q <- p + labs(subtitle="(b) On a logarithmic scale", y="Millions", x="Year") # reword title, subtitle, and caption + +# try trick from https://www.biostars.org/p/234142/ using aes_string? + +plotIt <- function(q, index){ #recursive function to draw lines for countries that had at least a million population in 1951 + if(index==1){return(q)}else{ #end recursion when index reaches 1 + q <- q + geom_line(aes(y=millIn51[,index]/1000), col = "grey76", size=0.5) # or add the y-series + plotIt(q, index-1) # and resursively call next index + } +} +q <- plotIt(q, 141) # call recursive function for the 140 countries in the CSV (recursion ends when index reaches 1, the x=Year field) +q <- q + geom_line(aes(y = Oceania, col = "Oceania"), size=1.5) # for visibility put Oceania back on top of the sea of grey +q <- q + scale_y_log10(breaks=logbreaks) # change to logarithmic y-scale +q <- q + theme(legend.position="bottom", legend.box = "horizontal") +q + +``` + +_Figure 2.7(b) . Total population for the world and continents between 1950 and 2015, both sexes combined, on a logarithmic scale, with the trend-lines for individual countries that had a population of at least one million in 1951._ \ No newline at end of file diff --git a/02-7-population-trends/02-7-population-history-x.html b/02-7-population-trends/02-7-population-history-x.html new file mode 100644 index 0000000..a7f0288 --- /dev/null +++ b/02-7-population-trends/02-7-population-history-x.html @@ -0,0 +1,451 @@ + + + + + + + + + + + + + +Art of Statistics: Figure 2.7 (pages 62) World Population Trends + + + + + + + + + + + + + + + + + + + + + + + + + + + + +
+ + + + + + + + + + + + + + + + +

Data is from United Nations Population Division: World Population Prospects 2017, specifically Total Population - Both Sexes

+
library(ggplot2)
+
+popTrends<-read.csv("02-7-pop-past-future-2.7-necessities.csv",header=TRUE) # reads data into popTrends data frame
+
+p <- ggplot(popTrends, aes(x=Year, y)) # constructs initial plot object, p
+p <- p + geom_line(aes(y = Europe, col = "Europe"), size=1.5) # adds a y-series
+p <- p + geom_line(aes(y = World, col = "World"), size=1.5) # adds a y-series
+p <- p + geom_line(aes(y = Africa, col = "Africa"), size=1.5) # adds a y-series
+p <- p + geom_line(aes(y = Oceania, col = "Oceania"), size=1.5) # adds a y-series
+p <- p + geom_line(aes(y = Asia, col = "Asia"), size=1.5) # adds a y-series
+p <- p + geom_line(aes(y = LatinAmericaandCarribean, col = "S America"), size=1.5) # adds a y-series
+p <- p + geom_line(aes(y = NorthernAmerica, col = "N America"), size=1.5) # adds a y-series
+p <- p +  labs(subtitle="(a) On a standard scale", y="Millions", x="Year") # Adds title, subtitle
+
+# Some of the colours rather pale to see in palette Accent, Set1 seems better
+p <- p + scale_colour_brewer(palette = "Set1")
+#p <- p + scale_colour_brewer(palette = "Accent")
+
+p <- p  + scale_size_continuous(name = "Size", guide = FALSE)  # turns off unwanted size legend
+p <- p + theme(legend.position="bottom", legend.box = "horizontal") # positions the colour legend
+p
+

+

Figure 2.7(a) . Total population for the world and continents between 1950 and 2015, both sexes combined.

+
# now for the logarithmic version
+millIn51 <- read.csv("02-7-pop-past-future million in 51.csv",header=TRUE) # reads data into millIn51 data frame
+
+logbreaks <- c(1, 5, 10, 50, 100, 500, 1000, 5000, 10000)
+q <- p +  labs(subtitle="(b) On a logarithmic scale", y="Millions", x="Year") # reword title, subtitle, and caption
+
+# try trick from https://www.biostars.org/p/234142/  using aes_string?
+
+plotIt <- function(q, index){ #recursive function to draw lines for countries that had at least a million population in 1951
+    if(index==1){return(q)}else{ #end recursion when index reaches 1
+      q <- q + geom_line(aes(y=millIn51[,index]/1000), col = "grey76", size=0.5) # or add the y-series
+      plotIt(q, index-1) # and resursively call next index
+    }  
+}
+q <- plotIt(q, 141) # call recursive function for the 140 countries in the CSV (recursion ends when index reaches 1, the x=Year field)
+q <- q + geom_line(aes(y = Oceania, col = "Oceania"), size=1.5) # for visibility put Oceania back on top of the sea of grey
+q <- q + scale_y_log10(breaks=logbreaks) # change to logarithmic y-scale
+q <- q + theme(legend.position="bottom", legend.box = "horizontal")  
+q
+

+

Figure 2.7(b) . Total population for the world and continents between 1950 and 2015, both sexes combined, on a logarithmic scale, with the trend-lines for individual countries that had a population of at least one million in 1951.

+ + + + +
+ + + + + + + + diff --git a/08-1-chevalier-de-mere/08-1-mere-simulation-x.Rmd b/08-1-chevalier-de-mere/08-1-mere-simulation-x.Rmd index 6a4d350..3cb0b52 100644 --- a/08-1-chevalier-de-mere/08-1-mere-simulation-x.Rmd +++ b/08-1-chevalier-de-mere/08-1-mere-simulation-x.Rmd @@ -1,7 +1,8 @@ --- - title: "Art of Statistics: Figure 8-1 (page 206) Simulation of two dice games played by the Chevaler de Mere" + title: "Art of Statistics: Figure 8-1 (page 206) Simulation of two dice games played by the Chevalier de Mere" #output: md_document output: html_document + #output: pdf_document --- The data is simulated within the code, and so will be different every time - this code will not produce the plots in the book, as I forgot to save the random-number seed. @@ -15,7 +16,34 @@ sixes2<-rbinom(nsims,24,1/36) # simulate the number of six-pairs in 24 throws of game1=cumsum(sixes1>0)/(1:nsims) # proportion of 'successes' game2=cumsum(sixes2>0)/(1:nsims) df=cbind.data.frame(n,game1,game2) +``` +First in R base graphics +```{r} +par(mfrow=c(2,1),mar=c(2,4,1,2)) +nsim1=100 +plot(1:nsim1, game1[1:nsim1], type="n",xlab="number of games",ylab="Proportion of wins", + ylim=c(0,1.0),axes=F, main="A simulation of the Chevalier playing his two games") +axis(1) #,at=(seq(0,100,by=5))) +axis(2, las=1 ) +lines(1:nsim1,game1[1:nsim1],lty=1, lwd=5,col="blue") +lines(1:nsim1,game2[1:nsim1],lty=1, lwd=5,col="red") +lines(c(1,nsim1),c(0.5,0.5),lty=3,lwd=1) +legend(20,1,bty="n",col=c("blue","red"),legend=c("Game 1: one die","Game 2: two dice"),lwd=5, cex=0.8) + +plot(nsim1:nsims, type="n",xlab="number of games",ylab="Proportion of wins", + ylim=c(0.4,0.6),axes=F, main="") +axis(1) #,at=(seq(0,100,by=5))) +axis(2, las=1 ) +lines(nsim1:nsims,game1[nsim1:nsims],lty=1, lwd=5,col="blue") +lines(nsim1:nsims,game2[nsim1:nsims],lty=1, lwd=5,col="red") +lines(c(1,nsims),c(0.5,0.5),lty=3,lwd=1) +``` + +_Figure 8.1. A computer simulation of 10,000 repeats of two games. In Game 1, you win if you throw a six in at most 4 throws of a fair die; in Game 2, you win if you throw a double-six in at most 24 throws of two fair dice. Over the first 100 games of each (upper chart) it may sometimes look like there is a higher chance of winning Game 2, but after thousands of plays (lower chart) it becomes apparent that Game 1 provides a marginally better bet._ + +Now in ggplot +```{r} library(ggplot2) P <- list() # prepare list for the plots for(i in 1:2) # 2 runs with parameter i @@ -28,7 +56,7 @@ for(i in 1:2) # 2 runs with parameter i if (i==1) {p <- p + geom_text(aes(x = 80, y = 1, label = "Game 1: one die", color = "Game 1", hjust=0)) p <- p + geom_text(aes(x = 80, y = 0.9, label = "Game 2: two dice", color = "Game 2", hjust=0))} p <- p +labs(x="Number of games", y="Proportion of wins") # Adds title, subtitle, and caption - p <- p + scale_colour_brewer(palette = "Accent") # sets the colour palette + p <- p + scale_colour_brewer(palette = "Set1") # sets the colour palette p <- p + theme(legend.position="none")#, legend.box = "horizontal") # removes the legend p <- p + theme(plot.caption=element_text(hjust = 0.5)) # centre justifies any caption P <- c(P, list(p)) # add plot to the plot-list @@ -42,4 +70,3 @@ for(i in 1:2) # 2 runs with parameter i ``` -Figure 8.1. A computer simulation of 10,000 repeats of two games. In Game 1, you win if you throw a six in at most 4 throws of a fair die; in Game 2, you win if you throw a double- six in at most 24 throws of two fair dice. Over the first 100 games of each (upper chart) it looks like there is a higher chance of winning Game 2, but after thousands of plays (lower chart) it becomes apparent that Game 1 provides a marginally better bet. \ No newline at end of file diff --git a/08-1-chevalier-de-mere/08-1-mere-simulation-x.html b/08-1-chevalier-de-mere/08-1-mere-simulation-x.html index 5e03e12..08c193d 100644 --- a/08-1-chevalier-de-mere/08-1-mere-simulation-x.html +++ b/08-1-chevalier-de-mere/08-1-mere-simulation-x.html @@ -11,7 +11,7 @@ -Art of Statistics: Figure 8-1 (page 206) Simulation of two dice games played by the Chevaler de Mere +Art of Statistics: Figure 8-1 (page 206) Simulation of two dice games played by the Chevalier de Mere