ECON 413
Data Structures
Erol Taymaz
Department of Economics
Middle East Technical University
Topics
- Data structures in R
- Data frames and indexing
- data.table package
R script
- Open a new file (R Script)
Save the file
file name: 004_studentID.R- Load the libraries
library(ggplot2)
library(data.table)
Refresh our memory
- Key operators
- <- is the assignment operator (creates an object)
- = parameter assignment for functions
- == test equality (returns a logical variable)
## [1] 5
# The following also works
b = 8
b
## [1] 8
# Duplicate the object
d <- b
d
## [1] 8
# If you use the same name, it will be replaced WITHOUT AND WARNING
a <- 10
a
## [1] 10
## [1] "a" "b" "d"
- List all objects with their description
## a : num 10
## b : num 8
## d : num 8
## [1] "b" "d"
## character(0)
- # means comment. Anything after the # symbol is ignored by R
Data modes
- Numeric
- Character
- Logical
- Factor
- Date
a <- 5
b <- "a"
c <- TRUE
d <- T
e <- "1"
# e + 1
class(e)
## [1] "character"
## [1] FALSE
f <- as.numeric(e)
class(f)
## [1] "numeric"
## [1] 2
Data types
- Vectors: A single column of data (even a column of length 1)
- Matrices: Two-dimensiona data (all rows and columns of all the same type)
- Arrays: Multi-dimensional data (all rows and columns of all the same type)
- Dataframes: Similar to matrics, but the columns can have different data types
- Dataframes are the most common data structure in R
- Datatables: Similar to dataframes but uses a different syntax
- Lists: are arbitrary combinations of disparate object types in R
Creating vectors
- : the sequence operator
- c() function
- rep() function
- seq() function
- sample() function
- rnorm() function
# : creates an index
aa <- 1:3
aa
## [1] 1 2 3
## [1] 100 101 102 103
## [1] -3 -2 -1 0 1 2 3
# c() combines values into a vector or a list
bb <- c(-1, 0, 1)
bb
## [1] -1 0 1
## [1] -1 0 1
## [1] -1 0 1 -1 0 1
bb <- c(1, 2, "a", "b", 5)
bb
## [1] "1" "2" "a" "b" "5"
## [1] FALSE
## [1] TRUE
# Why?
bb <- c(1, 2, "a", "b", TRUE)
bb
## [1] "1" "2" "a" "b" "TRUE"
## [1] 1 0 1
## [1] 1 1 2 2 3 3
## [1] 1 2 3 1 2 3
## [1] 1 2 3 1 2
rep(1:3, each = 2, times = 3)
## [1] 1 1 2 2 3 3 1 1 2 2 3 3 1 1 2 2 3 3
rep(1:3, each = 2, len = 7)
## [1] 1 1 2 2 3 3 1
## [1] -1 0 1 2 3 4
## [1] -1 1 3
seq(-1, 4, length.out = 3)
## [1] -1.0 1.5 4.0
# seq(-1, 4, by = 2, length.out = 3)
set.seed(123) aa <- rnorm(10) aa
set.seed(123) bb <- rnorm(10, mean = 0, sd = 1) bb
mean(rnorm(1000000, mean = 1, sd = 1)) ```
aa <- c(10:1)
aa
# 3rd element of aa
aa[3]
# 3rd and 6th elements
aa[c(3, 6)]
# Get even numbers
aa[seq(1, 10, by = 2)]
rm(list=ls())
Vector operations
# Simple operations
aa <- c(1:10)
bb <- sample(aa, 10)
cc <- 5
dd <- aa * bb
ee <- aa > bb
ff <- aa > cc
max(aa)
## [1] 10
## [1] 10
## [1] 8 2 9 5 6 6 7 10 9 10
# aa * ee
# Changing elements
aa[5] <- 50
aa
## [1] 1 2 3 4 50 6 7 8 9 10
aa[1:3] <- 10 * aa[1:3]
aa
## [1] 10 20 30 4 50 6 7 8 9 10
## [1] 10 20 30 4 50 6 7 8 9 10 NA NA NA NA 15
# Lag variables
aa <- c(1:10)
lagaa <- c(NA, aa[1:length(aa)-1])
# Growth rates
graa <- 100*((aa/lagaa) - 1)
graa
## [1] NA 100.00000 50.00000 33.33333 25.00000 20.00000 16.66667
## [8] 14.28571 12.50000 11.11111
## [1] NA 100.0 50.0 33.3 25.0 20.0 16.7 14.3 12.5 11.1
# Why vectors?
# Time
aa <- sample(100, 1000000, replace = TRUE)
bb <- sample(100, 1000000, replace = TRUE)
sumAB <- 0
system.time(for (i in c(1:length(aa))) {
sumAB <- sumAB + (aa[i] < bb[i])
})
## user system elapsed
## 0.226 0.008 0.234
## [1] 495116
system.time(sumBA <- sum(aa < bb))
## user system elapsed
## 0.004 0.000 0.005
## [1] 495116
Factors
- Factors are categorical variables
aa <- c("small", "large", "medium")
aa
## [1] "small" "large" "medium"
bb <- sample(3, 10, replace = TRUE)
bb
## [1] 1 2 3 3 1 2 1 1 3 3
## [1] "small" "large" "medium" "medium" "small" "large" "small"
## [8] "small" "medium" "medium"
## [1] "character"
## [1] small large medium medium small large small small medium medium
## Levels: large medium small
## [1] 3 1 2 2 3 1 3 3 2 2
## [1] NA NA NA NA NA NA NA NA NA NA
aaof <- ordered(aaf, levels=c("small", "medium", "large"))
aaof
## [1] small large medium medium small large small small medium medium
## Levels: small < medium < large
## [1] 3 1 2 2 3 1 3 3 2 2
## [1] 1 3 2 2 1 3 1 1 2 2
# Converting characters to numeric values
aa <- factor(c("3", "11", "2", "23", "313", "2"))
aa
## [1] 3 11 2 23 313 2
## Levels: 11 2 23 3 313
aan <- as.numeric(aa)
aan
## [1] 4 1 2 3 5 2
aan <- as.numeric(as.character(aa))
aan
## [1] 3 11 2 23 313 2
Data frames
- nxm table
- columns: variables (same data type)
- rows: observations
- Each column is a vector
- An nxm data frame is an m-element list
- Each column is an element of the list
Creating data frames
a = c(1:3)
b = c(T, F, T)
c = c("1", "One", "Two")
uu <- data.frame(a, b, c)
a <- a + 10
sum(a)
## [1] 36
## [1] 6
aa <- data.frame(a = c(1:3),
b = c(T, F, T),
c = c("1", "One", "Two"))
bb <- data.frame(c(1:3),
c(T, F, T),
c("1", "One", "Two"))
class(aa)
## [1] "data.frame"
## a b c
## 1 1 TRUE 1
## 2 2 FALSE One
## 3 3 TRUE Two
## c.1.3. c.T..F..T. c..1....One....Two..
## 1 1 TRUE 1
## 2 2 FALSE One
## 3 3 TRUE Two
names(bb)<-c("a", "b", "c")
bb
## a b c
## 1 1 TRUE 1
## 2 2 FALSE One
## 3 3 TRUE Two
aa <- data.frame(a = c(1:3), b = 1.5)
aa
## a b
## 1 1 1.5
## 2 2 1.5
## 3 3 1.5
aa <- data.frame(a = c(1:3), b = c(1:6))
aa
## a b
## 1 1 1
## 2 2 2
## 3 3 3
## 4 1 4
## 5 2 5
## 6 3 6
# aa <- data.frame(a = c(1:3), b = c(1:5))
# Empty data frame
aa <- data.frame(a = as.numeric(), b = as.character())
aa
## [1] a b
## <0 rows> (or 0-length row.names)
## [1] 0 2
Indexing data frames
- Rows, columns, and cells in a dataframe is accessed through indexing
df <- data.frame(a = c(1,2,3,4,5), b=c(4, 3, 6, 1, 1),
c=c("a", "a", "b", "c", "c"))
dfc <- data.frame(a = c(1,2,3,4,5), b=c(4, 3, 6, 1, 1),
c=c("a", "a", "b", "c", "c"),
stringsAsFactors = FALSE)
str(df)
## 'data.frame': 5 obs. of 3 variables:
## $ a: num 1 2 3 4 5
## $ b: num 4 3 6 1 1
## $ c: Factor w/ 3 levels "a","b","c": 1 1 2 3 3
## 'data.frame': 5 obs. of 3 variables:
## $ a: num 1 2 3 4 5
## $ b: num 4 3 6 1 1
## $ c: chr "a" "a" "b" "c" ...
# Variables (columns)
df[, 3]
## [1] a a b c c
## Levels: a b c
## c
## 1 a
## 2 a
## 3 b
## 4 c
## 5 c
## [1] a a b c c
## Levels: a b c
## c
## 1 a
## 2 a
## 3 b
## 4 c
## 5 c
## [1] a a b c c
## Levels: a b c
## [1] a a b c c
## Levels: a b c
## [1] a
## Levels: a b c
## [1] a
## Levels: a b c
## a b c
## 2 2 3 a
## b c
## 3 6 b
## 4 1 c
## 5 1 c
## [1] b c
## Levels: a b c
aa <- matrix(c(1, 2, 2, 1), nrow = 2, byrow = T)
aa
## [,1] [,2]
## [1,] 1 2
## [2,] 2 1
## [1] "4" "2"
Adding new variables into data frame
df$as <- df$c == "a"
df$ab <- df$a * df$b
df$sel1[df$as] <- 5
df$sel1[!df$as] <- 0
df$sel2 <- ifelse(df$c == "a", 5, 0)
df <- within(df, sel3 <- ifelse(c == "a", 5, 0))
Combining two data frames
dd <- data.frame(d = c(10, 20, 30, 40, 50), e=c("a", "a", "b", "d", "e"))
f <- c(1:5)
# Column append
dd
## d e
## 1 10 a
## 2 20 a
## 3 30 b
## 4 40 d
## 5 50 e
## d e f
## 1 10 a 1
## 2 20 a 2
## 3 30 b 3
## 4 40 d 4
## 5 50 e 5
# Row append
ee <- data.frame(e = "z", d = 3)
dr <- rbind(dd, ee)
dr
## d e
## 1 10 a
## 2 20 a
## 3 30 b
## 4 40 d
## 5 50 e
## 6 3 z
Using variables in a data frame
aa <- data.frame(a = c(1:1000))
# Generate a random walk
set.seed(123)
aa$b <- rnorm(1000)
aa$c <- cumsum(aa$b)
# Plot the random walk
ggplot(aa, aes(x = a, y = c)) + geom_path()
# Generate another random walk
aa$d <- rnorm(1000)
aa$e <- cumsum(aa$d)
# Plot the 2D random walk
ggplot(aa, aes(x = c, y = e)) + geom_path()
# Any correlation
m1 <- lm(c ~ e, data = aa)
summary(m1)
##
## Call:
## lm(formula = c ~ e, data = aa)
##
## Residuals:
## Min 1Q Median 3Q Max
## -16.7193 -4.3469 -0.2624 5.2217 14.1622
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 2.3850 0.3652 6.53 1.05e-10 ***
## e 0.3168 0.0173 18.31 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 6.309 on 998 degrees of freedom
## Multiple R-squared: 0.2515, Adjusted R-squared: 0.2508
## F-statistic: 335.4 on 1 and 998 DF, p-value: < 2.2e-16
# Calculate growth rate of e
aa$gre <- c(NA, 100 * (aa$e[2:1000] / aa$e[1:999]) - 1)
ggplot(aa, aes(x = c(1:1000), y = gre)) + geom_path()
ggplot(aa, aes(x = c(1:1000), y = e)) + geom_path()
Data operations
- Create a data set (data table)
- Filter a data table down to the rows and columns of interest
- Clean its values with arithmetic and string operations
- Summarize groups of rows (conditional on something)
- Compute new columns (variables) based on existing columns
- Join a data table with others, for further analysis
- Plot one variable vs. another (in many different ways)
- Model one variable as a function of others (for example, linear regression)
Data table operations
- DT[ i, j, by ]
- Take DT
- subset rows using i,
- then calculate j
- grouped by by
Data table operations
- Create a firm level data set
# Empty the global environment
rm(list=ls())
# Parameters
nf <- 1000
nt <- 20
ns <- 10
nobs <- nf*nt
kelas <- 0.1
lelas <- 0.2
melas <- 0.75
id <- rep(c(1:nf), each = nt)
year <- rep(c(1:nt), times = nf)
sect <- rep(c(1:ns), each = nobs / ns)
kg <- rnorm(nobs, 0.01, 0.02)
lg <- (kg/10) + rnorm(nobs, 0.01, 0.01)
mg <- (lg/10) + rnorm(nobs, 0.005, 0.01)
tc <- rnorm(nobs, 0.01, 0.01)
fs <- max(0.1, rep(rnorm(nf, 2, 0.5), each = nt))
re <- max(0.1, rep(rnorm(nf, 0, 0.5), each = nt))
ee <- rnorm(nobs, 0, 0.1)
# Data table
fdata <- data.table(id, year, sect, kg, lg, mg, tc, fs, re, ee)
fdata[, lk := fs*(1 + cumsum(kg)), by = id]
fdata[, ll := fs*(1 + cumsum(lg)), by = id]
fdata[, lm := fs*(1 + cumsum(mg)), by = id]
fdata[, lA := fs*(1 + cumsum(tc)), by = id]
fdata[, fe := max(0.1, rep(rnorm(nf, 0, 0.5)) + mean(lk)/5 + mean(ll)/5), by = id]
fdata[, lq := re + (lA + kelas*lk + lelas*ll + melas*lm) + ee]
rm(nf, nt, ns, nobs, kelas, lelas, melas, id, year, sect, kg, lg, mg,
tc, fs, re, ee)
Analyze the data set
# Basic statistics
class(fdata)
## [1] "data.table" "data.frame"
## [1] "id" "year" "sect" "kg" "lg" "mg" "tc" "fs" "re" "ee"
## [11] "lk" "ll" "lm" "lA" "fe" "lq"
## id year sect kg lg mg tc
## 1: 1 1 1 -0.0002320744 -0.022823040 -0.0059726355 0.004636073
## 2: 1 2 1 0.0147387576 0.012889209 0.0007184099 0.021840268
## 3: 1 3 1 -0.0008317834 -0.003767392 0.0144581771 0.009797892
## 4: 1 4 1 0.0343845529 0.006160969 0.0025308890 0.007914230
## 5: 1 5 1 0.0134827175 0.010668339 0.0073971701 0.015954826
## 6: 1 6 1 -0.0023053663 0.014908510 0.0077967036 0.004953235
## fs re ee lk ll lm lA
## 1: 3.863068 1.549448 0.13122076 3.862172 3.774901 3.839995 3.880978
## 2: 3.863068 1.549448 0.04699733 3.919108 3.824693 3.842771 3.965348
## 3: 3.863068 1.549448 0.15445384 3.915895 3.810139 3.898624 4.003198
## 4: 3.863068 1.549448 0.08272013 4.048725 3.833940 3.908401 4.033771
## 5: 3.863068 1.549448 -0.04270201 4.100810 3.875152 3.936976 4.095406
## 6: 3.863068 1.549448 -0.17411472 4.091904 3.932745 3.967096 4.114540
## fe lq
## 1: 3.068317 9.582840
## 2: 3.068317 9.600721
## 3: 3.068317 9.784685
## 4: 3.068317 9.768900
## 5: 3.068317 9.739995
## 6: 3.068317 9.660935
## Classes 'data.table' and 'data.frame': 20000 obs. of 16 variables:
## $ id : int 1 1 1 1 1 1 1 1 1 1 ...
## $ year: int 1 2 3 4 5 6 7 8 9 10 ...
## $ sect: int 1 1 1 1 1 1 1 1 1 1 ...
## $ kg : num -0.000232 0.014739 -0.000832 0.034385 0.013483 ...
## $ lg : num -0.02282 0.01289 -0.00377 0.00616 0.01067 ...
## $ mg : num -0.005973 0.000718 0.014458 0.002531 0.007397 ...
## $ tc : num 0.00464 0.02184 0.0098 0.00791 0.01595 ...
## $ fs : num 3.86 3.86 3.86 3.86 3.86 ...
## $ re : num 1.55 1.55 1.55 1.55 1.55 ...
## $ ee : num 0.1312 0.047 0.1545 0.0827 -0.0427 ...
## $ lk : num 3.86 3.92 3.92 4.05 4.1 ...
## $ ll : num 3.77 3.82 3.81 3.83 3.88 ...
## $ lm : num 3.84 3.84 3.9 3.91 3.94 ...
## $ lA : num 3.88 3.97 4 4.03 4.1 ...
## $ fe : num 3.07 3.07 3.07 3.07 3.07 ...
## $ lq : num 9.58 9.6 9.78 9.77 9.74 ...
## - attr(*, ".internal.selfref")=<externalptr>
## id year sect kg
## Min. : 1.0 Min. : 1.00 Min. : 1.0 Min. :-0.066906
## 1st Qu.: 250.8 1st Qu.: 5.75 1st Qu.: 3.0 1st Qu.:-0.003536
## Median : 500.5 Median :10.50 Median : 5.5 Median : 0.009693
## Mean : 500.5 Mean :10.50 Mean : 5.5 Mean : 0.009808
## 3rd Qu.: 750.2 3rd Qu.:15.25 3rd Qu.: 8.0 3rd Qu.: 0.023310
## Max. :1000.0 Max. :20.00 Max. :10.0 Max. : 0.088428
## lg mg tc
## Min. :-0.030734 Min. :-0.0312275 Min. :-0.029682
## 1st Qu.: 0.004022 1st Qu.:-0.0006562 1st Qu.: 0.003516
## Median : 0.010873 Median : 0.0061980 Median : 0.010156
## Mean : 0.010931 Mean : 0.0061835 Mean : 0.010122
## 3rd Qu.: 0.017872 3rd Qu.: 0.0129005 3rd Qu.: 0.016868
## Max. : 0.051246 Max. : 0.0437439 Max. : 0.048689
## fs re ee lk
## Min. :3.863 Min. :1.549 Min. :-0.4132090 Min. :3.324
## 1st Qu.:3.863 1st Qu.:1.549 1st Qu.:-0.0671504 1st Qu.:4.004
## Median :3.863 Median :1.549 Median :-0.0004495 Median :4.214
## Mean :3.863 Mean :1.549 Mean :-0.0005325 Mean :4.261
## 3rd Qu.:3.863 3rd Qu.:1.549 3rd Qu.: 0.0662047 3rd Qu.:4.479
## Max. :3.863 Max. :1.549 Max. : 0.4206090 Max. :5.566
## ll lm lA fe
## Min. :3.741 Min. :3.645 Min. :3.726 Min. :2.862
## 1st Qu.:4.076 1st Qu.:3.962 1st Qu.:4.061 1st Qu.:3.202
## Median :4.286 Median :4.093 Median :4.259 Median :3.310
## Mean :4.305 Mean :4.113 Mean :4.274 Mean :3.330
## 3rd Qu.:4.510 3rd Qu.:4.238 3rd Qu.:4.464 3rd Qu.:3.435
## Max. :5.208 Max. :4.881 Max. :5.201 Max. :4.195
## lq
## Min. : 9.124
## 1st Qu.: 9.831
## Median :10.177
## Mean :10.196
## 3rd Qu.:10.539
## Max. :11.698
# Check outliers
ggplot(fdata, aes(x = lq)) + geom_density() + labs(title = "Output")
ggplot(fdata, aes(x = lk)) + geom_density() + labs(title = "Capital")
ggplot(fdata, aes(x = ll)) + geom_density() + labs(title = "Labor")
ggplot(fdata, aes(x = lm)) + geom_density() + labs(title = "Inputs")
# Check correlations
ggplot(fdata, aes(x = lk, y = lq)) + geom_point() +
labs(x = "Capital", y = "Output")
ggplot(fdata, aes(x = ll, y = lq)) + geom_point() +
labs(x = "Labor", y = "Output")
ggplot(fdata, aes(x = lm, y = lq)) + geom_point() +
labs(x = "Inputs", y = "Output")
Deleting variables
# One variable
names(fdata)
## [1] "id" "year" "sect" "kg" "lg" "mg" "tc" "fs" "re" "ee"
## [11] "lk" "ll" "lm" "lA" "fe" "lq"
fdata[, kg := NULL]
names(fdata)
## [1] "id" "year" "sect" "lg" "mg" "tc" "fs" "re" "ee" "lk"
## [11] "ll" "lm" "lA" "fe" "lq"
# Many variables
fdata[, c("lg", "mg", "tc", "fs", "fe", "ee") := NULL]
names(fdata)
## [1] "id" "year" "sect" "re" "lk" "ll" "lm" "lA" "lq"
Subsetting rows
## id year sect re lk ll lm lA lq
## 1: 1 3 1 1.549448 3.915895 3.810139 3.898624 4.003198 9.784685
## 2: 1 4 1 1.549448 4.048725 3.833940 3.908401 4.033771 9.768900
## 3: 1 5 1 1.549448 4.100810 3.875152 3.936976 4.095406 9.739995
## id year sect re lk ll lm lA lq
## 1: 1 3 1 1.549448 3.915895 3.810139 3.898624 4.003198 9.784685
## 2: 1 4 1 1.549448 4.048725 3.833940 3.908401 4.033771 9.768900
## 3: 1 5 1 1.549448 4.100810 3.875152 3.936976 4.095406 9.739995
# Rows other than 3-5
fdata[-(3:5)]
## id year sect re lk ll lm lA
## 1: 1 1 1 1.549448 3.862172 3.774901 3.839995 3.880978
## 2: 1 2 1 1.549448 3.919108 3.824693 3.842771 3.965348
## 3: 1 6 1 1.549448 4.091904 3.932745 3.967096 4.114540
## 4: 1 7 1 1.549448 3.990932 3.992303 4.025372 4.091948
## 5: 1 8 1 1.549448 3.979831 4.035224 4.067952 4.101471
## ---
## 19993: 1000 16 10 1.549448 4.226719 4.150365 4.233986 4.389690
## 19994: 1000 17 10 1.549448 4.290319 4.176399 4.264723 4.341451
## 19995: 1000 18 10 1.549448 4.230260 4.141397 4.290881 4.371437
## 19996: 1000 19 10 1.549448 4.217664 4.239751 4.397182 4.361404
## 19997: 1000 20 10 1.549448 4.156285 4.269480 4.404235 4.381244
## lq
## 1: 9.582840
## 2: 9.600721
## 3: 9.660935
## 4: 9.674646
## 5: 9.817215
## ---
## 19993: 10.304541
## 19994: 10.431146
## 19995: 10.230929
## 19996: 10.529978
## 19997: 10.484083
## id year sect re lk ll lm lA
## 1: 1 1 1 1.549448 3.862172 3.774901 3.839995 3.880978
## 2: 1 2 1 1.549448 3.919108 3.824693 3.842771 3.965348
## 3: 1 6 1 1.549448 4.091904 3.932745 3.967096 4.114540
## 4: 1 7 1 1.549448 3.990932 3.992303 4.025372 4.091948
## 5: 1 8 1 1.549448 3.979831 4.035224 4.067952 4.101471
## ---
## 19993: 1000 16 10 1.549448 4.226719 4.150365 4.233986 4.389690
## 19994: 1000 17 10 1.549448 4.290319 4.176399 4.264723 4.341451
## 19995: 1000 18 10 1.549448 4.230260 4.141397 4.290881 4.371437
## 19996: 1000 19 10 1.549448 4.217664 4.239751 4.397182 4.361404
## 19997: 1000 20 10 1.549448 4.156285 4.269480 4.404235 4.381244
## lq
## 1: 9.582840
## 2: 9.600721
## 3: 9.660935
## 4: 9.674646
## 5: 9.817215
## ---
## 19993: 10.304541
## 19994: 10.431146
## 19995: 10.230929
## 19996: 10.529978
## 19997: 10.484083
# Rows conditional on variable values
fdata[lq >= 10]
## id year sect re lk ll lm lA
## 1: 1 9 1 1.549448 4.176540 4.082325 4.095582 4.157613
## 2: 1 10 1 1.549448 4.171845 4.105331 4.115125 4.252260
## 3: 1 11 1 1.549448 4.145885 4.137243 4.094832 4.277519
## 4: 1 12 1 1.549448 4.234963 4.160947 4.185049 4.271940
## 5: 1 13 1 1.549448 4.307689 4.215517 4.193729 4.366108
## ---
## 12574: 1000 16 10 1.549448 4.226719 4.150365 4.233986 4.389690
## 12575: 1000 17 10 1.549448 4.290319 4.176399 4.264723 4.341451
## 12576: 1000 18 10 1.549448 4.230260 4.141397 4.290881 4.371437
## 12577: 1000 19 10 1.549448 4.217664 4.239751 4.397182 4.361404
## 12578: 1000 20 10 1.549448 4.156285 4.269480 4.404235 4.381244
## lq
## 1: 10.10790
## 2: 10.21832
## 3: 10.06598
## 4: 10.36843
## 5: 10.18135
## ---
## 12574: 10.30454
## 12575: 10.43115
## 12576: 10.23093
## 12577: 10.52998
## 12578: 10.48408
fdata[lq >= 10 & ll <= 4]
## id year sect re lk ll lm lA lq
## 1: 22 8 1 1.549448 3.964360 3.943044 4.166672 4.147934 10.04763
## 2: 22 9 1 1.549448 3.907169 3.959647 4.138709 4.318437 10.18808
## 3: 22 10 1 1.549448 3.994934 3.984031 4.178632 4.338461 10.14727
## 4: 44 3 1 1.549448 3.929529 3.921664 3.957233 4.074655 10.01993
## 5: 44 10 1 1.549448 4.292750 3.986652 3.939943 4.250642 10.07869
## 6: 48 6 1 1.549448 4.337953 3.926515 3.911822 4.158060 10.05402
## 7: 60 5 1 1.549448 4.286924 3.980216 3.957640 4.127286 10.03818
## 8: 66 5 1 1.549448 4.103973 3.964526 3.989679 4.168325 10.03675
## 9: 81 5 1 1.549448 4.160035 3.983833 4.047617 4.100158 10.03234
## 10: 84 6 1 1.549448 4.090370 3.937021 3.904897 4.279514 10.09745
## 11: 86 7 1 1.549448 4.108229 3.958409 4.072148 4.156292 10.01233
## 12: 86 8 1 1.549448 4.127602 3.992879 4.134535 4.203489 10.11454
## 13: 86 9 1 1.549448 4.150056 3.998646 4.111365 4.266858 10.11756
## 14: 104 5 2 1.549448 4.087900 3.984045 3.970657 4.205672 10.02807
## 15: 105 3 2 1.549448 4.079946 3.901249 3.917645 4.184081 10.04988
## 16: 105 5 2 1.549448 4.105789 3.983337 3.911483 4.336507 10.17413
## 17: 117 6 2 1.549448 4.184468 3.914455 4.027427 4.221266 10.16526
## 18: 171 9 2 1.549448 4.393681 3.990453 3.874322 4.244626 10.02542
## 19: 179 5 2 1.549448 4.274190 3.988259 3.916190 4.081211 10.06478
## 20: 187 5 2 1.549448 3.928935 3.971678 4.212971 4.086787 10.04109
## 21: 196 5 2 1.549448 4.167278 3.912802 3.941541 4.053841 10.06467
## 22: 196 7 2 1.549448 4.196675 3.982085 4.031395 4.192934 10.10295
## 23: 203 6 3 1.549448 4.127320 3.943892 3.999382 4.149666 10.06831
## 24: 207 5 3 1.549448 3.533799 3.971118 3.919580 4.245862 10.04473
## 25: 208 6 3 1.549448 4.203470 3.993528 4.119346 4.285253 10.17616
## 26: 232 5 3 1.549448 4.102305 3.971926 4.080606 4.114355 10.05200
## 27: 238 4 3 1.549448 4.127385 3.971595 4.059243 4.235015 10.02852
## 28: 241 9 3 1.549448 4.038185 3.992956 4.097272 4.169077 10.05565
## 29: 274 6 3 1.549448 4.285715 3.989212 4.086154 4.045542 10.12868
## 30: 288 7 3 1.549448 4.205684 3.992251 4.109665 4.182581 10.16742
## 31: 288 8 3 1.549448 4.123643 3.996831 4.165512 4.197117 10.14770
## 32: 304 9 4 1.549448 3.888055 3.931938 4.042809 4.153552 10.12966
## 33: 322 4 4 1.549448 4.311193 3.973509 4.058342 4.145852 10.07326
## 34: 332 6 4 1.549448 3.986982 3.944099 4.168934 4.166780 10.04141
## 35: 343 4 4 1.549448 4.194759 3.971668 3.907479 4.087669 10.02019
## 36: 353 8 4 1.549448 4.386497 3.797304 4.029919 4.120017 10.05302
## 37: 353 10 4 1.549448 4.498942 3.919748 4.087537 4.232995 10.20960
## 38: 353 11 4 1.549448 4.450754 3.903496 4.109351 4.283214 10.24865
## 39: 353 12 4 1.549448 4.436245 3.857895 4.161753 4.376189 10.27744
## 40: 353 13 4 1.549448 4.695506 3.919274 4.116723 4.435363 10.25172
## 41: 353 14 4 1.549448 4.837503 3.994337 4.102823 4.531332 10.31149
## 42: 360 11 4 1.549448 3.891386 3.944686 4.071107 4.266197 10.09632
## 43: 360 13 4 1.549448 3.979733 3.981024 4.095422 4.361224 10.21422
## 44: 372 9 4 1.549448 4.187286 3.977012 4.156904 4.317882 10.17376
## 45: 376 2 4 1.549448 3.902858 3.950702 4.064942 3.958760 10.02872
## 46: 376 8 4 1.549448 4.132459 3.999163 4.303971 4.320557 10.39004
## 47: 488 5 5 1.549448 3.872855 3.988658 4.055062 4.125596 10.04974
## 48: 503 7 6 1.549448 4.172103 3.988974 4.088575 4.253803 10.23107
## 49: 514 7 6 1.549448 3.800607 3.976381 4.107767 4.314769 10.16383
## 50: 514 10 6 1.549448 3.398345 3.998618 4.087823 4.400269 10.05965
## 51: 515 6 6 1.549448 3.991894 3.971719 4.033537 4.100363 10.11305
## 52: 548 6 6 1.549448 3.493123 3.989201 3.966229 4.205552 10.04471
## 53: 549 4 6 1.549448 4.056349 3.940114 4.063638 4.199362 10.08180
## 54: 555 4 6 1.549448 3.764762 3.925230 4.072521 4.224980 10.07619
## 55: 555 5 6 1.549448 3.823802 3.875611 4.060127 4.254434 10.08968
## 56: 555 6 6 1.549448 3.942906 3.893384 4.029976 4.267406 10.00381
## 57: 555 7 6 1.549448 4.013744 3.908506 4.099076 4.264481 10.03565
## 58: 555 8 6 1.549448 4.079353 3.951379 4.115003 4.300655 10.25782
## 59: 555 9 6 1.549448 4.064608 3.978550 4.159558 4.325170 10.13528
## 60: 576 11 6 1.549448 4.133885 3.961912 4.027718 4.325095 10.27078
## 61: 576 12 6 1.549448 4.185800 3.985690 3.975917 4.351113 10.27484
## 62: 580 4 6 1.549448 4.052231 3.971976 4.037273 4.137004 10.01279
## 63: 590 6 6 1.549448 3.885280 3.783660 4.218220 4.151171 10.14796
## 64: 590 7 6 1.549448 3.966494 3.810581 4.219157 4.217682 10.03372
## 65: 590 8 6 1.549448 4.025476 3.872732 4.333581 4.291607 10.07089
## 66: 590 9 6 1.549448 4.068039 3.925001 4.319601 4.303562 10.33646
## 67: 590 10 6 1.549448 4.110593 3.981840 4.417890 4.334186 10.33488
## 68: 597 5 6 1.549448 4.081260 3.949865 4.095587 4.211193 10.20527
## 69: 597 6 6 1.549448 4.139364 3.965978 4.156614 4.227594 10.00824
## 70: 597 7 6 1.549448 4.170292 3.963587 4.214163 4.337503 10.16262
## 71: 598 5 6 1.549448 4.186374 3.985959 4.094650 4.140344 10.02744
## 72: 607 7 7 1.549448 4.003811 3.996909 4.035484 4.152300 10.04161
## 73: 611 8 7 1.549448 4.038371 3.979344 4.205507 4.208096 10.22450
## 74: 653 5 7 1.549448 4.040855 3.999637 4.076984 4.156147 10.00904
## 75: 670 7 7 1.549448 4.034499 3.958361 4.100971 4.160245 10.04954
## 76: 677 6 7 1.549448 4.429054 3.999156 3.990160 4.194015 10.11680
## 77: 694 9 7 1.549448 4.111514 3.914892 4.017000 4.211639 10.01480
## 78: 708 6 8 1.549448 3.814376 3.972530 4.078417 4.035947 10.00603
## 79: 721 4 8 1.549448 4.301520 3.997013 4.029182 4.065889 10.01151
## 80: 722 6 8 1.549448 3.944978 3.884962 3.916787 4.179817 10.01132
## 81: 736 7 8 1.549448 4.112922 3.953057 4.053268 4.173261 10.10585
## 82: 758 8 8 1.549448 4.252038 3.992350 4.081153 4.066778 10.04613
## 83: 841 5 9 1.549448 4.250160 3.990124 4.127973 4.061800 10.03974
## 84: 841 6 9 1.549448 4.367890 3.991815 4.136932 4.083888 10.06158
## 85: 861 7 9 1.549448 3.973755 3.962595 4.097687 4.236036 10.05534
## 86: 870 7 9 1.549448 4.088852 3.998658 3.804062 4.322108 10.00264
## 87: 870 8 9 1.549448 4.135867 3.978269 3.799667 4.279405 10.02226
## 88: 876 5 9 1.549448 4.163170 3.896405 3.957041 4.243617 10.05893
## 89: 876 6 9 1.549448 4.135281 3.911980 3.987460 4.334178 10.13110
## 90: 876 7 9 1.549448 4.056646 3.974623 4.047332 4.365960 10.20689
## 91: 889 6 9 1.549448 4.250195 3.974654 4.075458 4.248495 10.05386
## 92: 928 6 10 1.549448 3.905643 3.877204 3.919821 4.216733 10.00391
## 93: 957 3 10 1.549448 4.126499 3.969568 3.959378 4.068321 10.01070
## 94: 965 9 10 1.549448 4.016721 3.980972 3.870732 4.259614 10.01143
## 95: 967 5 10 1.549448 4.116947 3.939725 4.081508 4.250174 10.20110
## 96: 967 7 10 1.549448 4.091198 3.895581 4.222321 4.302242 10.27767
## 97: 967 8 10 1.549448 4.159804 3.955221 4.269998 4.360279 10.30976
## 98: 1000 6 10 1.549448 3.790789 3.947072 4.004013 4.111599 10.06579
## 99: 1000 7 10 1.549448 3.822544 3.969823 4.039003 4.132749 10.03942
## id year sect re lk ll lm lA lq
fdata[lq >= 10 && ll <= 4]
## Empty data.table (0 rows) of 9 cols: id,year,sect,re,lk,ll...
fdata[lq >= 10 | ll <= 4]
## id year sect re lk ll lm lA
## 1: 1 1 1 1.549448 3.862172 3.774901 3.839995 3.880978
## 2: 1 2 1 1.549448 3.919108 3.824693 3.842771 3.965348
## 3: 1 3 1 1.549448 3.915895 3.810139 3.898624 4.003198
## 4: 1 4 1 1.549448 4.048725 3.833940 3.908401 4.033771
## 5: 1 5 1 1.549448 4.100810 3.875152 3.936976 4.095406
## ---
## 15645: 1000 16 10 1.549448 4.226719 4.150365 4.233986 4.389690
## 15646: 1000 17 10 1.549448 4.290319 4.176399 4.264723 4.341451
## 15647: 1000 18 10 1.549448 4.230260 4.141397 4.290881 4.371437
## 15648: 1000 19 10 1.549448 4.217664 4.239751 4.397182 4.361404
## 15649: 1000 20 10 1.549448 4.156285 4.269480 4.404235 4.381244
## lq
## 1: 9.582840
## 2: 9.600721
## 3: 9.784685
## 4: 9.768900
## 5: 9.739995
## ---
## 15645: 10.304541
## 15646: 10.431146
## 15647: 10.230929
## 15648: 10.529978
## 15649: 10.484083
# All rows where id == 1
fdata[id == 1]
## id year sect re lk ll lm lA lq
## 1: 1 1 1 1.549448 3.862172 3.774901 3.839995 3.880978 9.582840
## 2: 1 2 1 1.549448 3.919108 3.824693 3.842771 3.965348 9.600721
## 3: 1 3 1 1.549448 3.915895 3.810139 3.898624 4.003198 9.784685
## 4: 1 4 1 1.549448 4.048725 3.833940 3.908401 4.033771 9.768900
## 5: 1 5 1 1.549448 4.100810 3.875152 3.936976 4.095406 9.739995
## 6: 1 6 1 1.549448 4.091904 3.932745 3.967096 4.114540 9.660935
## 7: 1 7 1 1.549448 3.990932 3.992303 4.025372 4.091948 9.674646
## 8: 1 8 1 1.549448 3.979831 4.035224 4.067952 4.101471 9.817215
## 9: 1 9 1 1.549448 4.176540 4.082325 4.095582 4.157613 10.107900
## 10: 1 10 1 1.549448 4.171845 4.105331 4.115125 4.252260 10.218320
## 11: 1 11 1 1.549448 4.145885 4.137243 4.094832 4.277519 10.065978
## 12: 1 12 1 1.549448 4.234963 4.160947 4.185049 4.271940 10.368427
## 13: 1 13 1 1.549448 4.307689 4.215517 4.193729 4.366108 10.181351
## 14: 1 14 1 1.549448 4.404392 4.255034 4.242622 4.382836 10.384899
## 15: 1 15 1 1.549448 4.421546 4.344362 4.316797 4.432053 10.594560
## 16: 1 16 1 1.549448 4.546914 4.359401 4.330609 4.486893 10.689910
## 17: 1 17 1 1.549448 4.494947 4.326558 4.309247 4.536981 10.516520
## 18: 1 18 1 1.549448 4.529802 4.452499 4.365506 4.592035 10.783428
## 19: 1 19 1 1.549448 4.514030 4.441219 4.367828 4.617422 10.670061
## 20: 1 20 1 1.549448 4.605257 4.494579 4.411256 4.681515 10.925218
setkey(fdata, id, year)
fdata[list(1)]
## id year sect re lk ll lm lA lq
## 1: 1 1 1 1.549448 3.862172 3.774901 3.839995 3.880978 9.582840
## 2: 1 2 1 1.549448 3.919108 3.824693 3.842771 3.965348 9.600721
## 3: 1 3 1 1.549448 3.915895 3.810139 3.898624 4.003198 9.784685
## 4: 1 4 1 1.549448 4.048725 3.833940 3.908401 4.033771 9.768900
## 5: 1 5 1 1.549448 4.100810 3.875152 3.936976 4.095406 9.739995
## 6: 1 6 1 1.549448 4.091904 3.932745 3.967096 4.114540 9.660935
## 7: 1 7 1 1.549448 3.990932 3.992303 4.025372 4.091948 9.674646
## 8: 1 8 1 1.549448 3.979831 4.035224 4.067952 4.101471 9.817215
## 9: 1 9 1 1.549448 4.176540 4.082325 4.095582 4.157613 10.107900
## 10: 1 10 1 1.549448 4.171845 4.105331 4.115125 4.252260 10.218320
## 11: 1 11 1 1.549448 4.145885 4.137243 4.094832 4.277519 10.065978
## 12: 1 12 1 1.549448 4.234963 4.160947 4.185049 4.271940 10.368427
## 13: 1 13 1 1.549448 4.307689 4.215517 4.193729 4.366108 10.181351
## 14: 1 14 1 1.549448 4.404392 4.255034 4.242622 4.382836 10.384899
## 15: 1 15 1 1.549448 4.421546 4.344362 4.316797 4.432053 10.594560
## 16: 1 16 1 1.549448 4.546914 4.359401 4.330609 4.486893 10.689910
## 17: 1 17 1 1.549448 4.494947 4.326558 4.309247 4.536981 10.516520
## 18: 1 18 1 1.549448 4.529802 4.452499 4.365506 4.592035 10.783428
## 19: 1 19 1 1.549448 4.514030 4.441219 4.367828 4.617422 10.670061
## 20: 1 20 1 1.549448 4.605257 4.494579 4.411256 4.681515 10.925218
fdata[list(c(1:2), c(10,20)), list(id, year, sect, lq)]
## id year sect lq
## 1: 1 10 1 10.21832
## 2: 2 20 1 10.98291
fdata[(id %in% c(1,2)) & (year %in% c(10, 20))]
## id year sect re lk ll lm lA lq
## 1: 1 10 1 1.549448 4.171845 4.105331 4.115125 4.252260 10.21832
## 2: 1 20 1 1.549448 4.605257 4.494579 4.411256 4.681515 10.92522
## 3: 2 10 1 1.549448 4.451089 4.262958 4.041126 4.245163 10.05635
## 4: 2 20 1 1.549448 4.918969 4.646606 4.350668 4.648884 10.98291
fdata[list(c(1:2))][year %in% c(10, 20)]
## id year sect re lk ll lm lA lq
## 1: 1 10 1 1.549448 4.171845 4.105331 4.115125 4.252260 10.21832
## 2: 1 20 1 1.549448 4.605257 4.494579 4.411256 4.681515 10.92522
## 3: 2 10 1 1.549448 4.451089 4.262958 4.041126 4.245163 10.05635
## 4: 2 20 1 1.549448 4.918969 4.646606 4.350668 4.648884 10.98291
sM <- data.table(id=c(1,1,2,2), year=c(10,20,10,20))
sM
## id year
## 1: 1 10
## 2: 1 20
## 3: 2 10
## 4: 2 20
## id year sect re lk ll lm lA lq
## 1: 1 10 1 1.549448 4.171845 4.105331 4.115125 4.252260 10.21832
## 2: 1 20 1 1.549448 4.605257 4.494579 4.411256 4.681515 10.92522
## 3: 2 10 1 1.549448 4.451089 4.262958 4.041126 4.245163 10.05635
## 4: 2 20 1 1.549448 4.918969 4.646606 4.350668 4.648884 10.98291
Selecting and computing columns
## [1] 1 1 1 1 1 1 1 1 1 1 1 1 1
## [14] 1 1 1 1 1 1 1 2 2 2 2 2 2
## [27] 2 2 2 2 2 2 2 2 2 2 2 2 2
## [40] 2 3 3 3 3 3 3 3 3 3 3 3 3
## [53] 3 3 3 3 3 3 3 3 4 4 4 4 4
## [66] 4 4 4 4 4 4 4 4 4 4 4 4 4
## [79] 4 4 5 5 5 5 5 5 5 5 5 5 5
## [92] 5 5 5 5 5 5 5 5 5 6 6 6 6
## [105] 6 6 6 6 6 6 6 6 6 6 6 6 6
## [118] 6 6 6 7 7 7 7 7 7 7 7 7 7
## [131] 7 7 7 7 7 7 7 7 7 7 8 8 8
## [144] 8 8 8 8 8 8 8 8 8 8 8 8 8
## [157] 8 8 8 8 9 9 9 9 9 9 9 9 9
## [170] 9 9 9 9 9 9 9 9 9 9 9 10 10
## [183] 10 10 10 10 10 10 10 10 10 10 10 10 10
## [196] 10 10 10 10 10 11 11 11 11 11 11 11 11
## [209] 11 11 11 11 11 11 11 11 11 11 11 11 12
## [222] 12 12 12 12 12 12 12 12 12 12 12 12 12
## [235] 12 12 12 12 12 12 13 13 13 13 13 13 13
## [248] 13 13 13 13 13 13 13 13 13 13 13 13 13
## [261] 14 14 14 14 14 14 14 14 14 14 14 14 14
## [274] 14 14 14 14 14 14 14 15 15 15 15 15 15
## [287] 15 15 15 15 15 15 15 15 15 15 15 15 15
## [300] 15 16 16 16 16 16 16 16 16 16 16 16 16
## [313] 16 16 16 16 16 16 16 16 17 17 17 17 17
## [326] 17 17 17 17 17 17 17 17 17 17 17 17 17
## [339] 17 17 18 18 18 18 18 18 18 18 18 18 18
## [352] 18 18 18 18 18 18 18 18 18 19 19 19 19
## [365] 19 19 19 19 19 19 19 19 19 19 19 19 19
## [378] 19 19 19 20 20 20 20 20 20 20 20 20 20
## [391] 20 20 20 20 20 20 20 20 20 20 21 21 21
## [404] 21 21 21 21 21 21 21 21 21 21 21 21 21
## [417] 21 21 21 21 22 22 22 22 22 22 22 22 22
## [430] 22 22 22 22 22 22 22 22 22 22 22 23 23
## [443] 23 23 23 23 23 23 23 23 23 23 23 23 23
## [456] 23 23 23 23 23 24 24 24 24 24 24 24 24
## [469] 24 24 24 24 24 24 24 24 24 24 24 24 25
## [482] 25 25 25 25 25 25 25 25 25 25 25 25 25
## [495] 25 25 25 25 25 25 26 26 26 26 26 26 26
## [508] 26 26 26 26 26 26 26 26 26 26 26 26 26
## [521] 27 27 27 27 27 27 27 27 27 27 27 27 27
## [534] 27 27 27 27 27 27 27 28 28 28 28 28 28
## [547] 28 28 28 28 28 28 28 28 28 28 28 28 28
## [560] 28 29 29 29 29 29 29 29 29 29 29 29 29
## [573] 29 29 29 29 29 29 29 29 30 30 30 30 30
## [586] 30 30 30 30 30 30 30 30 30 30 30 30 30
## [599] 30 30 31 31 31 31 31 31 31 31 31 31 31
## [612] 31 31 31 31 31 31 31 31 31 32 32 32 32
## [625] 32 32 32 32 32 32 32 32 32 32 32 32 32
## [638] 32 32 32 33 33 33 33 33 33 33 33 33 33
## [651] 33 33 33 33 33 33 33 33 33 33 34 34 34
## [664] 34 34 34 34 34 34 34 34 34 34 34 34 34
## [677] 34 34 34 34 35 35 35 35 35 35 35 35 35
## [690] 35 35 35 35 35 35 35 35 35 35 35 36 36
## [703] 36 36 36 36 36 36 36 36 36 36 36 36 36
## [716] 36 36 36 36 36 37 37 37 37 37 37 37 37
## [729] 37 37 37 37 37 37 37 37 37 37 37 37 38
## [742] 38 38 38 38 38 38 38 38 38 38 38 38 38
## [755] 38 38 38 38 38 38 39 39 39 39 39 39 39
## [768] 39 39 39 39 39 39 39 39 39 39 39 39 39
## [781] 40 40 40 40 40 40 40 40 40 40 40 40 40
## [794] 40 40 40 40 40 40 40 41 41 41 41 41 41
## [807] 41 41 41 41 41 41 41 41 41 41 41 41 41
## [820] 41 42 42 42 42 42 42 42 42 42 42 42 42
## [833] 42 42 42 42 42 42 42 42 43 43 43 43 43
## [846] 43 43 43 43 43 43 43 43 43 43 43 43 43
## [859] 43 43 44 44 44 44 44 44 44 44 44 44 44
## [872] 44 44 44 44 44 44 44 44 44 45 45 45 45
## [885] 45 45 45 45 45 45 45 45 45 45 45 45 45
## [898] 45 45 45 46 46 46 46 46 46 46 46 46 46
## [911] 46 46 46 46 46 46 46 46 46 46 47 47 47
## [924] 47 47 47 47 47 47 47 47 47 47 47 47 47
## [937] 47 47 47 47 48 48 48 48 48 48 48 48 48
## [950] 48 48 48 48 48 48 48 48 48 48 48 49 49
## [963] 49 49 49 49 49 49 49 49 49 49 49 49 49
## [976] 49 49 49 49 49 50 50 50 50 50 50 50 50
## [989] 50 50 50 50 50 50 50 50 50 50 50 50 51
## [1002] 51 51 51 51 51 51 51 51 51 51 51 51 51
## [1015] 51 51 51 51 51 51 52 52 52 52 52 52 52
## [1028] 52 52 52 52 52 52 52 52 52 52 52 52 52
## [1041] 53 53 53 53 53 53 53 53 53 53 53 53 53
## [1054] 53 53 53 53 53 53 53 54 54 54 54 54 54
## [1067] 54 54 54 54 54 54 54 54 54 54 54 54 54
## [1080] 54 55 55 55 55 55 55 55 55 55 55 55 55
## [1093] 55 55 55 55 55 55 55 55 56 56 56 56 56
## [1106] 56 56 56 56 56 56 56 56 56 56 56 56 56
## [1119] 56 56 57 57 57 57 57 57 57 57 57 57 57
## [1132] 57 57 57 57 57 57 57 57 57 58 58 58 58
## [1145] 58 58 58 58 58 58 58 58 58 58 58 58 58
## [1158] 58 58 58 59 59 59 59 59 59 59 59 59 59
## [1171] 59 59 59 59 59 59 59 59 59 59 60 60 60
## [1184] 60 60 60 60 60 60 60 60 60 60 60 60 60
## [1197] 60 60 60 60 61 61 61 61 61 61 61 61 61
## [1210] 61 61 61 61 61 61 61 61 61 61 61 62 62
## [1223] 62 62 62 62 62 62 62 62 62 62 62 62 62
## [1236] 62 62 62 62 62 63 63 63 63 63 63 63 63
## [1249] 63 63 63 63 63 63 63 63 63 63 63 63 64
## [1262] 64 64 64 64 64 64 64 64 64 64 64 64 64
## [1275] 64 64 64 64 64 64 65 65 65 65 65 65 65
## [1288] 65 65 65 65 65 65 65 65 65 65 65 65 65
## [1301] 66 66 66 66 66 66 66 66 66 66 66 66 66
## [1314] 66 66 66 66 66 66 66 67 67 67 67 67 67
## [1327] 67 67 67 67 67 67 67 67 67 67 67 67 67
## [1340] 67 68 68 68 68 68 68 68 68 68 68 68 68
## [1353] 68 68 68 68 68 68 68 68 69 69 69 69 69
## [1366] 69 69 69 69 69 69 69 69 69 69 69 69 69
## [1379] 69 69 70 70 70 70 70 70 70 70 70 70 70
## [1392] 70 70 70 70 70 70 70 70 70 71 71 71 71
## [1405] 71 71 71 71 71 71 71 71 71 71 71 71 71
## [1418] 71 71 71 72 72 72 72 72 72 72 72 72 72
## [1431] 72 72 72 72 72 72 72 72 72 72 73 73 73
## [1444] 73 73 73 73 73 73 73 73 73 73 73 73 73
## [1457] 73 73 73 73 74 74 74 74 74 74 74 74 74
## [1470] 74 74 74 74 74 74 74 74 74 74 74 75 75
## [1483] 75 75 75 75 75 75 75 75 75 75 75 75 75
## [1496] 75 75 75 75 75 76 76 76 76 76 76 76 76
## [1509] 76 76 76 76 76 76 76 76 76 76 76 76 77
## [1522] 77 77 77 77 77 77 77 77 77 77 77 77 77
## [1535] 77 77 77 77 77 77 78 78 78 78 78 78 78
## [1548] 78 78 78 78 78 78 78 78 78 78 78 78 78
## [1561] 79 79 79 79 79 79 79 79 79 79 79 79 79
## [1574] 79 79 79 79 79 79 79 80 80 80 80 80 80
## [1587] 80 80 80 80 80 80 80 80 80 80 80 80 80
## [1600] 80 81 81 81 81 81 81 81 81 81 81 81 81
## [1613] 81 81 81 81 81 81 81 81 82 82 82 82 82
## [1626] 82 82 82 82 82 82 82 82 82 82 82 82 82
## [1639] 82 82 83 83 83 83 83 83 83 83 83 83 83
## [1652] 83 83 83 83 83 83 83 83 83 84 84 84 84
## [1665] 84 84 84 84 84 84 84 84 84 84 84 84 84
## [1678] 84 84 84 85 85 85 85 85 85 85 85 85 85
## [1691] 85 85 85 85 85 85 85 85 85 85 86 86 86
## [1704] 86 86 86 86 86 86 86 86 86 86 86 86 86
## [1717] 86 86 86 86 87 87 87 87 87 87 87 87 87
## [1730] 87 87 87 87 87 87 87 87 87 87 87 88 88
## [1743] 88 88 88 88 88 88 88 88 88 88 88 88 88
## [1756] 88 88 88 88 88 89 89 89 89 89 89 89 89
## [1769] 89 89 89 89 89 89 89 89 89 89 89 89 90
## [1782] 90 90 90 90 90 90 90 90 90 90 90 90 90
## [1795] 90 90 90 90 90 90 91 91 91 91 91 91 91
## [1808] 91 91 91 91 91 91 91 91 91 91 91 91 91
## [1821] 92 92 92 92 92 92 92 92 92 92 92 92 92
## [1834] 92 92 92 92 92 92 92 93 93 93 93 93 93
## [1847] 93 93 93 93 93 93 93 93 93 93 93 93 93
## [1860] 93 94 94 94 94 94 94 94 94 94 94 94 94
## [1873] 94 94 94 94 94 94 94 94 95 95 95 95 95
## [1886] 95 95 95 95 95 95 95 95 95 95 95 95 95
## [1899] 95 95 96 96 96 96 96 96 96 96 96 96 96
## [1912] 96 96 96 96 96 96 96 96 96 97 97 97 97
## [1925] 97 97 97 97 97 97 97 97 97 97 97 97 97
## [1938] 97 97 97 98 98 98 98 98 98 98 98 98 98
## [1951] 98 98 98 98 98 98 98 98 98 98 99 99 99
## [1964] 99 99 99 99 99 99 99 99 99 99 99 99 99
## [1977] 99 99 99 99 100 100 100 100 100 100 100 100 100
## [1990] 100 100 100 100 100 100 100 100 100 100 100 101 101
## [2003] 101 101 101 101 101 101 101 101 101 101 101 101 101
## [2016] 101 101 101 101 101 102 102 102 102 102 102 102 102
## [2029] 102 102 102 102 102 102 102 102 102 102 102 102 103
## [2042] 103 103 103 103 103 103 103 103 103 103 103 103 103
## [2055] 103 103 103 103 103 103 104 104 104 104 104 104 104
## [2068] 104 104 104 104 104 104 104 104 104 104 104 104 104
## [2081] 105 105 105 105 105 105 105 105 105 105 105 105 105
## [2094] 105 105 105 105 105 105 105 106 106 106 106 106 106
## [2107] 106 106 106 106 106 106 106 106 106 106 106 106 106
## [2120] 106 107 107 107 107 107 107 107 107 107 107 107 107
## [2133] 107 107 107 107 107 107 107 107 108 108 108 108 108
## [2146] 108 108 108 108 108 108 108 108 108 108 108 108 108
## [2159] 108 108 109 109 109 109 109 109 109 109 109 109 109
## [2172] 109 109 109 109 109 109 109 109 109 110 110 110 110
## [2185] 110 110 110 110 110 110 110 110 110 110 110 110 110
## [2198] 110 110 110 111 111 111 111 111 111 111 111 111 111
## [2211] 111 111 111 111 111 111 111 111 111 111 112 112 112
## [2224] 112 112 112 112 112 112 112 112 112 112 112 112 112
## [2237] 112 112 112 112 113 113 113 113 113 113 113 113 113
## [2250] 113 113 113 113 113 113 113 113 113 113 113 114 114
## [2263] 114 114 114 114 114 114 114 114 114 114 114 114 114
## [2276] 114 114 114 114 114 115 115 115 115 115 115 115 115
## [2289] 115 115 115 115 115 115 115 115 115 115 115 115 116
## [2302] 116 116 116 116 116 116 116 116 116 116 116 116 116
## [2315] 116 116 116 116 116 116 117 117 117 117 117 117 117
## [2328] 117 117 117 117 117 117 117 117 117 117 117 117 117
## [2341] 118 118 118 118 118 118 118 118 118 118 118 118 118
## [2354] 118 118 118 118 118 118 118 119 119 119 119 119 119
## [2367] 119 119 119 119 119 119 119 119 119 119 119 119 119
## [2380] 119 120 120 120 120 120 120 120 120 120 120 120 120
## [2393] 120 120 120 120 120 120 120 120 121 121 121 121 121
## [2406] 121 121 121 121 121 121 121 121 121 121 121 121 121
## [2419] 121 121 122 122 122 122 122 122 122 122 122 122 122
## [2432] 122 122 122 122 122 122 122 122 122 123 123 123 123
## [2445] 123 123 123 123 123 123 123 123 123 123 123 123 123
## [2458] 123 123 123 124 124 124 124 124 124 124 124 124 124
## [2471] 124 124 124 124 124 124 124 124 124 124 125 125 125
## [2484] 125 125 125 125 125 125 125 125 125 125 125 125 125
## [2497] 125 125 125 125 126 126 126 126 126 126 126 126 126
## [2510] 126 126 126 126 126 126 126 126 126 126 126 127 127
## [2523] 127 127 127 127 127 127 127 127 127 127 127 127 127
## [2536] 127 127 127 127 127 128 128 128 128 128 128 128 128
## [2549] 128 128 128 128 128 128 128 128 128 128 128 128 129
## [2562] 129 129 129 129 129 129 129 129 129 129 129 129 129
## [2575] 129 129 129 129 129 129 130 130 130 130 130 130 130
## [2588] 130 130 130 130 130 130 130 130 130 130 130 130 130
## [2601] 131 131 131 131 131 131 131 131 131 131 131 131 131
## [2614] 131 131 131 131 131 131 131 132 132 132 132 132 132
## [2627] 132 132 132 132 132 132 132 132 132 132 132 132 132
## [2640] 132 133 133 133 133 133 133 133 133 133 133 133 133
## [2653] 133 133 133 133 133 133 133 133 134 134 134 134 134
## [2666] 134 134 134 134 134 134 134 134 134 134 134 134 134
## [2679] 134 134 135 135 135 135 135 135 135 135 135 135 135
## [2692] 135 135 135 135 135 135 135 135 135 136 136 136 136
## [2705] 136 136 136 136 136 136 136 136 136 136 136 136 136
## [2718] 136 136 136 137 137 137 137 137 137 137 137 137 137
## [2731] 137 137 137 137 137 137 137 137 137 137 138 138 138
## [2744] 138 138 138 138 138 138 138 138 138 138 138 138 138
## [2757] 138 138 138 138 139 139 139 139 139 139 139 139 139
## [2770] 139 139 139 139 139 139 139 139 139 139 139 140 140
## [2783] 140 140 140 140 140 140 140 140 140 140 140 140 140
## [2796] 140 140 140 140 140 141 141 141 141 141 141 141 141
## [2809] 141 141 141 141 141 141 141 141 141 141 141 141 142
## [2822] 142 142 142 142 142 142 142 142 142 142 142 142 142
## [2835] 142 142 142 142 142 142 143 143 143 143 143 143 143
## [2848] 143 143 143 143 143 143 143 143 143 143 143 143 143
## [2861] 144 144 144 144 144 144 144 144 144 144 144 144 144
## [2874] 144 144 144 144 144 144 144 145 145 145 145 145 145
## [2887] 145 145 145 145 145 145 145 145 145 145 145 145 145
## [2900] 145 146 146 146 146 146 146 146 146 146 146 146 146
## [2913] 146 146 146 146 146 146 146 146 147 147 147 147 147
## [2926] 147 147 147 147 147 147 147 147 147 147 147 147 147
## [2939] 147 147 148 148 148 148 148 148 148 148 148 148 148
## [2952] 148 148 148 148 148 148 148 148 148 149 149 149 149
## [2965] 149 149 149 149 149 149 149 149 149 149 149 149 149
## [2978] 149 149 149 150 150 150 150 150 150 150 150 150 150
## [2991] 150 150 150 150 150 150 150 150 150 150 151 151 151
## [3004] 151 151 151 151 151 151 151 151 151 151 151 151 151
## [3017] 151 151 151 151 152 152 152 152 152 152 152 152 152
## [3030] 152 152 152 152 152 152 152 152 152 152 152 153 153
## [3043] 153 153 153 153 153 153 153 153 153 153 153 153 153
## [3056] 153 153 153 153 153 154 154 154 154 154 154 154 154
## [3069] 154 154 154 154 154 154 154 154 154 154 154 154 155
## [3082] 155 155 155 155 155 155 155 155 155 155 155 155 155
## [3095] 155 155 155 155 155 155 156 156 156 156 156 156 156
## [3108] 156 156 156 156 156 156 156 156 156 156 156 156 156
## [3121] 157 157 157 157 157 157 157 157 157 157 157 157 157
## [3134] 157 157 157 157 157 157 157 158 158 158 158 158 158
## [3147] 158 158 158 158 158 158 158 158 158 158 158 158 158
## [3160] 158 159 159 159 159 159 159 159 159 159 159 159 159
## [3173] 159 159 159 159 159 159 159 159 160 160 160 160 160
## [3186] 160 160 160 160 160 160 160 160 160 160 160 160 160
## [3199] 160 160 161 161 161 161 161 161 161 161 161 161 161
## [3212] 161 161 161 161 161 161 161 161 161 162 162 162 162
## [3225] 162 162 162 162 162 162 162 162 162 162 162 162 162
## [3238] 162 162 162 163 163 163 163 163 163 163 163 163 163
## [3251] 163 163 163 163 163 163 163 163 163 163 164 164 164
## [3264] 164 164 164 164 164 164 164 164 164 164 164 164 164
## [3277] 164 164 164 164 165 165 165 165 165 165 165 165 165
## [3290] 165 165 165 165 165 165 165 165 165 165 165 166 166
## [3303] 166 166 166 166 166 166 166 166 166 166 166 166 166
## [3316] 166 166 166 166 166 167 167 167 167 167 167 167 167
## [3329] 167 167 167 167 167 167 167 167 167 167 167 167 168
## [3342] 168 168 168 168 168 168 168 168 168 168 168 168 168
## [3355] 168 168 168 168 168 168 169 169 169 169 169 169 169
## [3368] 169 169 169 169 169 169 169 169 169 169 169 169 169
## [3381] 170 170 170 170 170 170 170 170 170 170 170 170 170
## [3394] 170 170 170 170 170 170 170 171 171 171 171 171 171
## [3407] 171 171 171 171 171 171 171 171 171 171 171 171 171
## [3420] 171 172 172 172 172 172 172 172 172 172 172 172 172
## [3433] 172 172 172 172 172 172 172 172 173 173 173 173 173
## [3446] 173 173 173 173 173 173 173 173 173 173 173 173 173
## [3459] 173 173 174 174 174 174 174 174 174 174 174 174 174
## [3472] 174 174 174 174 174 174 174 174 174 175 175 175 175
## [3485] 175 175 175 175 175 175 175 175 175 175 175 175 175
## [3498] 175 175 175 176 176 176 176 176 176 176 176 176 176
## [3511] 176 176 176 176 176 176 176 176 176 176 177 177 177
## [3524] 177 177 177 177 177 177 177 177 177 177 177 177 177
## [3537] 177 177 177 177 178 178 178 178 178 178 178 178 178
## [3550] 178 178 178 178 178 178 178 178 178 178 178 179 179
## [3563] 179 179 179 179 179 179 179 179 179 179 179 179 179
## [3576] 179 179 179 179 179 180 180 180 180 180 180 180 180
## [3589] 180 180 180 180 180 180 180 180 180 180 180 180 181
## [3602] 181 181 181 181 181 181 181 181 181 181 181 181 181
## [3615] 181 181 181 181 181 181 182 182 182 182 182 182 182
## [3628] 182 182 182 182 182 182 182 182 182 182 182 182 182
## [3641] 183 183 183 183 183 183 183 183 183 183 183 183 183
## [3654] 183 183 183 183 183 183 183 184 184 184 184 184 184
## [3667] 184 184 184 184 184 184 184 184 184 184 184 184 184
## [3680] 184 185 185 185 185 185 185 185 185 185 185 185 185
## [3693] 185 185 185 185 185 185 185 185 186 186 186 186 186
## [3706] 186 186 186 186 186 186 186 186 186 186 186 186 186
## [3719] 186 186 187 187 187 187 187 187 187 187 187 187 187
## [3732] 187 187 187 187 187 187 187 187 187 188 188 188 188
## [3745] 188 188 188 188 188 188 188 188 188 188 188 188 188
## [3758] 188 188 188 189 189 189 189 189 189 189 189 189 189
## [3771] 189 189 189 189 189 189 189 189 189 189 190 190 190
## [3784] 190 190 190 190 190 190 190 190 190 190 190 190 190
## [3797] 190 190 190 190 191 191 191 191 191 191 191 191 191
## [3810] 191 191 191 191 191 191 191 191 191 191 191 192 192
## [3823] 192 192 192 192 192 192 192 192 192 192 192 192 192
## [3836] 192 192 192 192 192 193 193 193 193 193 193 193 193
## [3849] 193 193 193 193 193 193 193 193 193 193 193 193 194
## [3862] 194 194 194 194 194 194 194 194 194 194 194 194 194
## [3875] 194 194 194 194 194 194 195 195 195 195 195 195 195
## [3888] 195 195 195 195 195 195 195 195 195 195 195 195 195
## [3901] 196 196 196 196 196 196 196 196 196 196 196 196 196
## [3914] 196 196 196 196 196 196 196 197 197 197 197 197 197
## [3927] 197 197 197 197 197 197 197 197 197 197 197 197 197
## [3940] 197 198 198 198 198 198 198 198 198 198 198 198 198
## [3953] 198 198 198 198 198 198 198 198 199 199 199 199 199
## [3966] 199 199 199 199 199 199 199 199 199 199 199 199 199
## [3979] 199 199 200 200 200 200 200 200 200 200 200 200 200
## [3992] 200 200 200 200 200 200 200 200 200 201 201 201 201
## [4005] 201 201 201 201 201 201 201 201 201 201 201 201 201
## [4018] 201 201 201 202 202 202 202 202 202 202 202 202 202
## [4031] 202 202 202 202 202 202 202 202 202 202 203 203 203
## [4044] 203 203 203 203 203 203 203 203 203 203 203 203 203
## [4057] 203 203 203 203 204 204 204 204 204 204 204 204 204
## [4070] 204 204 204 204 204 204 204 204 204 204 204 205 205
## [4083] 205 205 205 205 205 205 205 205 205 205 205 205 205
## [4096] 205 205 205 205 205 206 206 206 206 206 206 206 206
## [4109] 206 206 206 206 206 206 206 206 206 206 206 206 207
## [4122] 207 207 207 207 207 207 207 207 207 207 207 207 207
## [4135] 207 207 207 207 207 207 208 208 208 208 208 208 208
## [4148] 208 208 208 208 208 208 208 208 208 208 208 208 208
## [4161] 209 209 209 209 209 209 209 209 209 209 209 209 209
## [4174] 209 209 209 209 209 209 209 210 210 210 210 210 210
## [4187] 210 210 210 210 210 210 210 210 210 210 210 210 210
## [4200] 210 211 211 211 211 211 211 211 211 211 211 211 211
## [4213] 211 211 211 211 211 211 211 211 212 212 212 212 212
## [4226] 212 212 212 212 212 212 212 212 212 212 212 212 212
## [4239] 212 212 213 213 213 213 213 213 213 213 213 213 213
## [4252] 213 213 213 213 213 213 213 213 213 214 214 214 214
## [4265] 214 214 214 214 214 214 214 214 214 214 214 214 214
## [4278] 214 214 214 215 215 215 215 215 215 215 215 215 215
## [4291] 215 215 215 215 215 215 215 215 215 215 216 216 216
## [4304] 216 216 216 216 216 216 216 216 216 216 216 216 216
## [4317] 216 216 216 216 217 217 217 217 217 217 217 217 217
## [4330] 217 217 217 217 217 217 217 217 217 217 217 218 218
## [4343] 218 218 218 218 218 218 218 218 218 218 218 218 218
## [4356] 218 218 218 218 218 219 219 219 219 219 219 219 219
## [4369] 219 219 219 219 219 219 219 219 219 219 219 219 220
## [4382] 220 220 220 220 220 220 220 220 220 220 220 220 220
## [4395] 220 220 220 220 220 220 221 221 221 221 221 221 221
## [4408] 221 221 221 221 221 221 221 221 221 221 221 221 221
## [4421] 222 222 222 222 222 222 222 222 222 222 222 222 222
## [4434] 222 222 222 222 222 222 222 223 223 223 223 223 223
## [4447] 223 223 223 223 223 223 223 223 223 223 223 223 223
## [4460] 223 224 224 224 224 224 224 224 224 224 224 224 224
## [4473] 224 224 224 224 224 224 224 224 225 225 225 225 225
## [4486] 225 225 225 225 225 225 225 225 225 225 225 225 225
## [4499] 225 225 226 226 226 226 226 226 226 226 226 226 226
## [4512] 226 226 226 226 226 226 226 226 226 227 227 227 227
## [4525] 227 227 227 227 227 227 227 227 227 227 227 227 227
## [4538] 227 227 227 228 228 228 228 228 228 228 228 228 228
## [4551] 228 228 228 228 228 228 228 228 228 228 229 229 229
## [4564] 229 229 229 229 229 229 229 229 229 229 229 229 229
## [4577] 229 229 229 229 230 230 230 230 230 230 230 230 230
## [4590] 230 230 230 230 230 230 230 230 230 230 230 231 231
## [4603] 231 231 231 231 231 231 231 231 231 231 231 231 231
## [4616] 231 231 231 231 231 232 232 232 232 232 232 232 232
## [4629] 232 232 232 232 232 232 232 232 232 232 232 232 233
## [4642] 233 233 233 233 233 233 233 233 233 233 233 233 233
## [4655] 233 233 233 233 233 233 234 234 234 234 234 234 234
## [4668] 234 234 234 234 234 234 234 234 234 234 234 234 234
## [4681] 235 235 235 235 235 235 235 235 235 235 235 235 235
## [4694] 235 235 235 235 235 235 235 236 236 236 236 236 236
## [4707] 236 236 236 236 236 236 236 236 236 236 236 236 236
## [4720] 236 237 237 237 237 237 237 237 237 237 237 237 237
## [4733] 237 237 237 237 237 237 237 237 238 238 238 238 238
## [4746] 238 238 238 238 238 238 238 238 238 238 238 238 238
## [4759] 238 238 239 239 239 239 239 239 239 239 239 239 239
## [4772] 239 239 239 239 239 239 239 239 239 240 240 240 240
## [4785] 240 240 240 240 240 240 240 240 240 240 240 240 240
## [4798] 240 240 240 241 241 241 241 241 241 241 241 241 241
## [4811] 241 241 241 241 241 241 241 241 241 241 242 242 242
## [4824] 242 242 242 242 242 242 242 242 242 242 242 242 242
## [4837] 242 242 242 242 243 243 243 243 243 243 243 243 243
## [4850] 243 243 243 243 243 243 243 243 243 243 243 244 244
## [4863] 244 244 244 244 244 244 244 244 244 244 244 244 244
## [4876] 244 244 244 244 244 245 245 245 245 245 245 245 245
## [4889] 245 245 245 245 245 245 245 245 245 245 245 245 246
## [4902] 246 246 246 246 246 246 246 246 246 246 246 246 246
## [4915] 246 246 246 246 246 246 247 247 247 247 247 247 247
## [4928] 247 247 247 247 247 247 247 247 247 247 247 247 247
## [4941] 248 248 248 248 248 248 248 248 248 248 248 248 248
## [4954] 248 248 248 248 248 248 248 249 249 249 249 249 249
## [4967] 249 249 249 249 249 249 249 249 249 249 249 249 249
## [4980] 249 250 250 250 250 250 250 250 250 250 250 250 250
## [4993] 250 250 250 250 250 250 250 250 251 251 251 251 251
## [5006] 251 251 251 251 251 251 251 251 251 251 251 251 251
## [5019] 251 251 252 252 252 252 252 252 252 252 252 252 252
## [5032] 252 252 252 252 252 252 252 252 252 253 253 253 253
## [5045] 253 253 253 253 253 253 253 253 253 253 253 253 253
## [5058] 253 253 253 254 254 254 254 254 254 254 254 254 254
## [5071] 254 254 254 254 254 254 254 254 254 254 255 255 255
## [5084] 255 255 255 255 255 255 255 255 255 255 255 255 255
## [5097] 255 255 255 255 256 256 256 256 256 256 256 256 256
## [5110] 256 256 256 256 256 256 256 256 256 256 256 257 257
## [5123] 257 257 257 257 257 257 257 257 257 257 257 257 257
## [5136] 257 257 257 257 257 258 258 258 258 258 258 258 258
## [5149] 258 258 258 258 258 258 258 258 258 258 258 258 259
## [5162] 259 259 259 259 259 259 259 259 259 259 259 259 259
## [5175] 259 259 259 259 259 259 260 260 260 260 260 260 260
## [5188] 260 260 260 260 260 260 260 260 260 260 260 260 260
## [5201] 261 261 261 261 261 261 261 261 261 261 261 261 261
## [5214] 261 261 261 261 261 261 261 262 262 262 262 262 262
## [5227] 262 262 262 262 262 262 262 262 262 262 262 262 262
## [5240] 262 263 263 263 263 263 263 263 263 263 263 263 263
## [5253] 263 263 263 263 263 263 263 263 264 264 264 264 264
## [5266] 264 264 264 264 264 264 264 264 264 264 264 264 264
## [5279] 264 264 265 265 265 265 265 265 265 265 265 265 265
## [5292] 265 265 265 265 265 265 265 265 265 266 266 266 266
## [5305] 266 266 266 266 266 266 266 266 266 266 266 266 266
## [5318] 266 266 266 267 267 267 267 267 267 267 267 267 267
## [5331] 267 267 267 267 267 267 267 267 267 267 268 268 268
## [5344] 268 268 268 268 268 268 268 268 268 268 268 268 268
## [5357] 268 268 268 268 269 269 269 269 269 269 269 269 269
## [5370] 269 269 269 269 269 269 269 269 269 269 269 270 270
## [5383] 270 270 270 270 270 270 270 270 270 270 270 270 270
## [5396] 270 270 270 270 270 271 271 271 271 271 271 271 271
## [5409] 271 271 271 271 271 271 271 271 271 271 271 271 272
## [5422] 272 272 272 272 272 272 272 272 272 272 272 272 272
## [5435] 272 272 272 272 272 272 273 273 273 273 273 273 273
## [5448] 273 273 273 273 273 273 273 273 273 273 273 273 273
## [5461] 274 274 274 274 274 274 274 274 274 274 274 274 274
## [5474] 274 274 274 274 274 274 274 275 275 275 275 275 275
## [5487] 275 275 275 275 275 275 275 275 275 275 275 275 275
## [5500] 275 276 276 276 276 276 276 276 276 276 276 276 276
## [5513] 276 276 276 276 276 276 276 276 277 277 277 277 277
## [5526] 277 277 277 277 277 277 277 277 277 277 277 277 277
## [5539] 277 277 278 278 278 278 278 278 278 278 278 278 278
## [5552] 278 278 278 278 278 278 278 278 278 279 279 279 279
## [5565] 279 279 279 279 279 279 279 279 279 279 279 279 279
## [5578] 279 279 279 280 280 280 280 280 280 280 280 280 280
## [5591] 280 280 280 280 280 280 280 280 280 280 281 281 281
## [5604] 281 281 281 281 281 281 281 281 281 281 281 281 281
## [5617] 281 281 281 281 282 282 282 282 282 282 282 282 282
## [5630] 282 282 282 282 282 282 282 282 282 282 282 283 283
## [5643] 283 283 283 283 283 283 283 283 283 283 283 283 283
## [5656] 283 283 283 283 283 284 284 284 284 284 284 284 284
## [5669] 284 284 284 284 284 284 284 284 284 284 284 284 285
## [5682] 285 285 285 285 285 285 285 285 285 285 285 285 285
## [5695] 285 285 285 285 285 285 286 286 286 286 286 286 286
## [5708] 286 286 286 286 286 286 286 286 286 286 286 286 286
## [5721] 287 287 287 287 287 287 287 287 287 287 287 287 287
## [5734] 287 287 287 287 287 287 287 288 288 288 288 288 288
## [5747] 288 288 288 288 288 288 288 288 288 288 288 288 288
## [5760] 288 289 289 289 289 289 289 289 289 289 289 289 289
## [5773] 289 289 289 289 289 289 289 289 290 290 290 290 290
## [5786] 290 290 290 290 290 290 290 290 290 290 290 290 290
## [5799] 290 290 291 291 291 291 291 291 291 291 291 291 291
## [5812] 291 291 291 291 291 291 291 291 291 292 292 292 292
## [5825] 292 292 292 292 292 292 292 292 292 292 292 292 292
## [5838] 292 292 292 293 293 293 293 293 293 293 293 293 293
## [5851] 293 293 293 293 293 293 293 293 293 293 294 294 294
## [5864] 294 294 294 294 294 294 294 294 294 294 294 294 294
## [5877] 294 294 294 294 295 295 295 295 295 295 295 295 295
## [5890] 295 295 295 295 295 295 295 295 295 295 295 296 296
## [5903] 296 296 296 296 296 296 296 296 296 296 296 296 296
## [5916] 296 296 296 296 296 297 297 297 297 297 297 297 297
## [5929] 297 297 297 297 297 297 297 297 297 297 297 297 298
## [5942] 298 298 298 298 298 298 298 298 298 298 298 298 298
## [5955] 298 298 298 298 298 298 299 299 299 299 299 299 299
## [5968] 299 299 299 299 299 299 299 299 299 299 299 299 299
## [5981] 300 300 300 300 300 300 300 300 300 300 300 300 300
## [5994] 300 300 300 300 300 300 300 301 301 301 301 301 301
## [6007] 301 301 301 301 301 301 301 301 301 301 301 301 301
## [6020] 301 302 302 302 302 302 302 302 302 302 302 302 302
## [6033] 302 302 302 302 302 302 302 302 303 303 303 303 303
## [6046] 303 303 303 303 303 303 303 303 303 303 303 303 303
## [6059] 303 303 304 304 304 304 304 304 304 304 304 304 304
## [6072] 304 304 304 304 304 304 304 304 304 305 305 305 305
## [6085] 305 305 305 305 305 305 305 305 305 305 305 305 305
## [6098] 305 305 305 306 306 306 306 306 306 306 306 306 306
## [6111] 306 306 306 306 306 306 306 306 306 306 307 307 307
## [6124] 307 307 307 307 307 307 307 307 307 307 307 307 307
## [6137] 307 307 307 307 308 308 308 308 308 308 308 308 308
## [6150] 308 308 308 308 308 308 308 308 308 308 308 309 309
## [6163] 309 309 309 309 309 309 309 309 309 309 309 309 309
## [6176] 309 309 309 309 309 310 310 310 310 310 310 310 310
## [6189] 310 310 310 310 310 310 310 310 310 310 310 310 311
## [6202] 311 311 311 311 311 311 311 311 311 311 311 311 311
## [6215] 311 311 311 311 311 311 312 312 312 312 312 312 312
## [6228] 312 312 312 312 312 312 312 312 312 312 312 312 312
## [6241] 313 313 313 313 313 313 313 313 313 313 313 313 313
## [6254] 313 313 313 313 313 313 313 314 314 314 314 314 314
## [6267] 314 314 314 314 314 314 314 314 314 314 314 314 314
## [6280] 314 315 315 315 315 315 315 315 315 315 315 315 315
## [6293] 315 315 315 315 315 315 315 315 316 316 316 316 316
## [6306] 316 316 316 316 316 316 316 316 316 316 316 316 316
## [6319] 316 316 317 317 317 317 317 317 317 317 317 317 317
## [6332] 317 317 317 317 317 317 317 317 317 318 318 318 318
## [6345] 318 318 318 318 318 318 318 318 318 318 318 318 318
## [6358] 318 318 318 319 319 319 319 319 319 319 319 319 319
## [6371] 319 319 319 319 319 319 319 319 319 319 320 320 320
## [6384] 320 320 320 320 320 320 320 320 320 320 320 320 320
## [6397] 320 320 320 320 321 321 321 321 321 321 321 321 321
## [6410] 321 321 321 321 321 321 321 321 321 321 321 322 322
## [6423] 322 322 322 322 322 322 322 322 322 322 322 322 322
## [6436] 322 322 322 322 322 323 323 323 323 323 323 323 323
## [6449] 323 323 323 323 323 323 323 323 323 323 323 323 324
## [6462] 324 324 324 324 324 324 324 324 324 324 324 324 324
## [6475] 324 324 324 324 324 324 325 325 325 325 325 325 325
## [6488] 325 325 325 325 325 325 325 325 325 325 325 325 325
## [6501] 326 326 326 326 326 326 326 326 326 326 326 326 326
## [6514] 326 326 326 326 326 326 326 327 327 327 327 327 327
## [6527] 327 327 327 327 327 327 327 327 327 327 327 327 327
## [6540] 327 328 328 328 328 328 328 328 328 328 328 328 328
## [6553] 328 328 328 328 328 328 328 328 329 329 329 329 329
## [6566] 329 329 329 329 329 329 329 329 329 329 329 329 329
## [6579] 329 329 330 330 330 330 330 330 330 330 330 330 330
## [6592] 330 330 330 330 330 330 330 330 330 331 331 331 331
## [6605] 331 331 331 331 331 331 331 331 331 331 331 331 331
## [6618] 331 331 331 332 332 332 332 332 332 332 332 332 332
## [6631] 332 332 332 332 332 332 332 332 332 332 333 333 333
## [6644] 333 333 333 333 333 333 333 333 333 333 333 333 333
## [6657] 333 333 333 333 334 334 334 334 334 334 334 334 334
## [6670] 334 334 334 334 334 334 334 334 334 334 334 335 335
## [6683] 335 335 335 335 335 335 335 335 335 335 335 335 335
## [6696] 335 335 335 335 335 336 336 336 336 336 336 336 336
## [6709] 336 336 336 336 336 336 336 336 336 336 336 336 337
## [6722] 337 337 337 337 337 337 337 337 337 337 337 337 337
## [6735] 337 337 337 337 337 337 338 338 338 338 338 338 338
## [6748] 338 338 338 338 338 338 338 338 338 338 338 338 338
## [6761] 339 339 339 339 339 339 339 339 339 339 339 339 339
## [6774] 339 339 339 339 339 339 339 340 340 340 340 340 340
## [6787] 340 340 340 340 340 340 340 340 340 340 340 340 340
## [6800] 340 341 341 341 341 341 341 341 341 341 341 341 341
## [6813] 341 341 341 341 341 341 341 341 342 342 342 342 342
## [6826] 342 342 342 342 342 342 342 342 342 342 342 342 342
## [6839] 342 342 343 343 343 343 343 343 343 343 343 343 343
## [6852] 343 343 343 343 343 343 343 343 343 344 344 344 344
## [6865] 344 344 344 344 344 344 344 344 344 344 344 344 344
## [6878] 344 344 344 345 345 345 345 345 345 345 345 345 345
## [6891] 345 345 345 345 345 345 345 345 345 345 346 346 346
## [6904] 346 346 346 346 346 346 346 346 346 346 346 346 346
## [6917] 346 346 346 346 347 347 347 347 347 347 347 347 347
## [6930] 347 347 347 347 347 347 347 347 347 347 347 348 348
## [6943] 348 348 348 348 348 348 348 348 348 348 348 348 348
## [6956] 348 348 348 348 348 349 349 349 349 349 349 349 349
## [6969] 349 349 349 349 349 349 349 349 349 349 349 349 350
## [6982] 350 350 350 350 350 350 350 350 350 350 350 350 350
## [6995] 350 350 350 350 350 350 351 351 351 351 351 351 351
## [7008] 351 351 351 351 351 351 351 351 351 351 351 351 351
## [7021] 352 352 352 352 352 352 352 352 352 352 352 352 352
## [7034] 352 352 352 352 352 352 352 353 353 353 353 353 353
## [7047] 353 353 353 353 353 353 353 353 353 353 353 353 353
## [7060] 353 354 354 354 354 354 354 354 354 354 354 354 354
## [7073] 354 354 354 354 354 354 354 354 355 355 355 355 355
## [7086] 355 355 355 355 355 355 355 355 355 355 355 355 355
## [7099] 355 355 356 356 356 356 356 356 356 356 356 356 356
## [7112] 356 356 356 356 356 356 356 356 356 357 357 357 357
## [7125] 357 357 357 357 357 357 357 357 357 357 357 357 357
## [7138] 357 357 357 358 358 358 358 358 358 358 358 358 358
## [7151] 358 358 358 358 358 358 358 358 358 358 359 359 359
## [7164] 359 359 359 359 359 359 359 359 359 359 359 359 359
## [7177] 359 359 359 359 360 360 360 360 360 360 360 360 360
## [7190] 360 360 360 360 360 360 360 360 360 360 360 361 361
## [7203] 361 361 361 361 361 361 361 361 361 361 361 361 361
## [7216] 361 361 361 361 361 362 362 362 362 362 362 362 362
## [7229] 362 362 362 362 362 362 362 362 362 362 362 362 363
## [7242] 363 363 363 363 363 363 363 363 363 363 363 363 363
## [7255] 363 363 363 363 363 363 364 364 364 364 364 364 364
## [7268] 364 364 364 364 364 364 364 364 364 364 364 364 364
## [7281] 365 365 365 365 365 365 365 365 365 365 365 365 365
## [7294] 365 365 365 365 365 365 365 366 366 366 366 366 366
## [7307] 366 366 366 366 366 366 366 366 366 366 366 366 366
## [7320] 366 367 367 367 367 367 367 367 367 367 367 367 367
## [7333] 367 367 367 367 367 367 367 367 368 368 368 368 368
## [7346] 368 368 368 368 368 368 368 368 368 368 368 368 368
## [7359] 368 368 369 369 369 369 369 369 369 369 369 369 369
## [7372] 369 369 369 369 369 369 369 369 369 370 370 370 370
## [7385] 370 370 370 370 370 370 370 370 370 370 370 370 370
## [7398] 370 370 370 371 371 371 371 371 371 371 371 371 371
## [7411] 371 371 371 371 371 371 371 371 371 371 372 372 372
## [7424] 372 372 372 372 372 372 372 372 372 372 372 372 372
## [7437] 372 372 372 372 373 373 373 373 373 373 373 373 373
## [7450] 373 373 373 373 373 373 373 373 373 373 373 374 374
## [7463] 374 374 374 374 374 374 374 374 374 374 374 374 374
## [7476] 374 374 374 374 374 375 375 375 375 375 375 375 375
## [7489] 375 375 375 375 375 375 375 375 375 375 375 375 376
## [7502] 376 376 376 376 376 376 376 376 376 376 376 376 376
## [7515] 376 376 376 376 376 376 377 377 377 377 377 377 377
## [7528] 377 377 377 377 377 377 377 377 377 377 377 377 377
## [7541] 378 378 378 378 378 378 378 378 378 378 378 378 378
## [7554] 378 378 378 378 378 378 378 379 379 379 379 379 379
## [7567] 379 379 379 379 379 379 379 379 379 379 379 379 379
## [7580] 379 380 380 380 380 380 380 380 380 380 380 380 380
## [7593] 380 380 380 380 380 380 380 380 381 381 381 381 381
## [7606] 381 381 381 381 381 381 381 381 381 381 381 381 381
## [7619] 381 381 382 382 382 382 382 382 382 382 382 382 382
## [7632] 382 382 382 382 382 382 382 382 382 383 383 383 383
## [7645] 383 383 383 383 383 383 383 383 383 383 383 383 383
## [7658] 383 383 383 384 384 384 384 384 384 384 384 384 384
## [7671] 384 384 384 384 384 384 384 384 384 384 385 385 385
## [7684] 385 385 385 385 385 385 385 385 385 385 385 385 385
## [7697] 385 385 385 385 386 386 386 386 386 386 386 386 386
## [7710] 386 386 386 386 386 386 386 386 386 386 386 387 387
## [7723] 387 387 387 387 387 387 387 387 387 387 387 387 387
## [7736] 387 387 387 387 387 388 388 388 388 388 388 388 388
## [7749] 388 388 388 388 388 388 388 388 388 388 388 388 389
## [7762] 389 389 389 389 389 389 389 389 389 389 389 389 389
## [7775] 389 389 389 389 389 389 390 390 390 390 390 390 390
## [7788] 390 390 390 390 390 390 390 390 390 390 390 390 390
## [7801] 391 391 391 391 391 391 391 391 391 391 391 391 391
## [7814] 391 391 391 391 391 391 391 392 392 392 392 392 392
## [7827] 392 392 392 392 392 392 392 392 392 392 392 392 392
## [7840] 392 393 393 393 393 393 393 393 393 393 393 393 393
## [7853] 393 393 393 393 393 393 393 393 394 394 394 394 394
## [7866] 394 394 394 394 394 394 394 394 394 394 394 394 394
## [7879] 394 394 395 395 395 395 395 395 395 395 395 395 395
## [7892] 395 395 395 395 395 395 395 395 395 396 396 396 396
## [7905] 396 396 396 396 396 396 396 396 396 396 396 396 396
## [7918] 396 396 396 397 397 397 397 397 397 397 397 397 397
## [7931] 397 397 397 397 397 397 397 397 397 397 398 398 398
## [7944] 398 398 398 398 398 398 398 398 398 398 398 398 398
## [7957] 398 398 398 398 399 399 399 399 399 399 399 399 399
## [7970] 399 399 399 399 399 399 399 399 399 399 399 400 400
## [7983] 400 400 400 400 400 400 400 400 400 400 400 400 400
## [7996] 400 400 400 400 400 401 401 401 401 401 401 401 401
## [8009] 401 401 401 401 401 401 401 401 401 401 401 401 402
## [8022] 402 402 402 402 402 402 402 402 402 402 402 402 402
## [8035] 402 402 402 402 402 402 403 403 403 403 403 403 403
## [8048] 403 403 403 403 403 403 403 403 403 403 403 403 403
## [8061] 404 404 404 404 404 404 404 404 404 404 404 404 404
## [8074] 404 404 404 404 404 404 404 405 405 405 405 405 405
## [8087] 405 405 405 405 405 405 405 405 405 405 405 405 405
## [8100] 405 406 406 406 406 406 406 406 406 406 406 406 406
## [8113] 406 406 406 406 406 406 406 406 407 407 407 407 407
## [8126] 407 407 407 407 407 407 407 407 407 407 407 407 407
## [8139] 407 407 408 408 408 408 408 408 408 408 408 408 408
## [8152] 408 408 408 408 408 408 408 408 408 409 409 409 409
## [8165] 409 409 409 409 409 409 409 409 409 409 409 409 409
## [8178] 409 409 409 410 410 410 410 410 410 410 410 410 410
## [8191] 410 410 410 410 410 410 410 410 410 410 411 411 411
## [8204] 411 411 411 411 411 411 411 411 411 411 411 411 411
## [8217] 411 411 411 411 412 412 412 412 412 412 412 412 412
## [8230] 412 412 412 412 412 412 412 412 412 412 412 413 413
## [8243] 413 413 413 413 413 413 413 413 413 413 413 413 413
## [8256] 413 413 413 413 413 414 414 414 414 414 414 414 414
## [8269] 414 414 414 414 414 414 414 414 414 414 414 414 415
## [8282] 415 415 415 415 415 415 415 415 415 415 415 415 415
## [8295] 415 415 415 415 415 415 416 416 416 416 416 416 416
## [8308] 416 416 416 416 416 416 416 416 416 416 416 416 416
## [8321] 417 417 417 417 417 417 417 417 417 417 417 417 417
## [8334] 417 417 417 417 417 417 417 418 418 418 418 418 418
## [8347] 418 418 418 418 418 418 418 418 418 418 418 418 418
## [8360] 418 419 419 419 419 419 419 419 419 419 419 419 419
## [8373] 419 419 419 419 419 419 419 419 420 420 420 420 420
## [8386] 420 420 420 420 420 420 420 420 420 420 420 420 420
## [8399] 420 420 421 421 421 421 421 421 421 421 421 421 421
## [8412] 421 421 421 421 421 421 421 421 421 422 422 422 422
## [8425] 422 422 422 422 422 422 422 422 422 422 422 422 422
## [8438] 422 422 422 423 423 423 423 423 423 423 423 423 423
## [8451] 423 423 423 423 423 423 423 423 423 423 424 424 424
## [8464] 424 424 424 424 424 424 424 424 424 424 424 424 424
## [8477] 424 424 424 424 425 425 425 425 425 425 425 425 425
## [8490] 425 425 425 425 425 425 425 425 425 425 425 426 426
## [8503] 426 426 426 426 426 426 426 426 426 426 426 426 426
## [8516] 426 426 426 426 426 427 427 427 427 427 427 427 427
## [8529] 427 427 427 427 427 427 427 427 427 427 427 427 428
## [8542] 428 428 428 428 428 428 428 428 428 428 428 428 428
## [8555] 428 428 428 428 428 428 429 429 429 429 429 429 429
## [8568] 429 429 429 429 429 429 429 429 429 429 429 429 429
## [8581] 430 430 430 430 430 430 430 430 430 430 430 430 430
## [8594] 430 430 430 430 430 430 430 431 431 431 431 431 431
## [8607] 431 431 431 431 431 431 431 431 431 431 431 431 431
## [8620] 431 432 432 432 432 432 432 432 432 432 432 432 432
## [8633] 432 432 432 432 432 432 432 432 433 433 433 433 433
## [8646] 433 433 433 433 433 433 433 433 433 433 433 433 433
## [8659] 433 433 434 434 434 434 434 434 434 434 434 434 434
## [8672] 434 434 434 434 434 434 434 434 434 435 435 435 435
## [8685] 435 435 435 435 435 435 435 435 435 435 435 435 435
## [8698] 435 435 435 436 436 436 436 436 436 436 436 436 436
## [8711] 436 436 436 436 436 436 436 436 436 436 437 437 437
## [8724] 437 437 437 437 437 437 437 437 437 437 437 437 437
## [8737] 437 437 437 437 438 438 438 438 438 438 438 438 438
## [8750] 438 438 438 438 438 438 438 438 438 438 438 439 439
## [8763] 439 439 439 439 439 439 439 439 439 439 439 439 439
## [8776] 439 439 439 439 439 440 440 440 440 440 440 440 440
## [8789] 440 440 440 440 440 440 440 440 440 440 440 440 441
## [8802] 441 441 441 441 441 441 441 441 441 441 441 441 441
## [8815] 441 441 441 441 441 441 442 442 442 442 442 442 442
## [8828] 442 442 442 442 442 442 442 442 442 442 442 442 442
## [8841] 443 443 443 443 443 443 443 443 443 443 443 443 443
## [8854] 443 443 443 443 443 443 443 444 444 444 444 444 444
## [8867] 444 444 444 444 444 444 444 444 444 444 444 444 444
## [8880] 444 445 445 445 445 445 445 445 445 445 445 445 445
## [8893] 445 445 445 445 445 445 445 445 446 446 446 446 446
## [8906] 446 446 446 446 446 446 446 446 446 446 446 446 446
## [8919] 446 446 447 447 447 447 447 447 447 447 447 447 447
## [8932] 447 447 447 447 447 447 447 447 447 448 448 448 448
## [8945] 448 448 448 448 448 448 448 448 448 448 448 448 448
## [8958] 448 448 448 449 449 449 449 449 449 449 449 449 449
## [8971] 449 449 449 449 449 449 449 449 449 449 450 450 450
## [8984] 450 450 450 450 450 450 450 450 450 450 450 450 450
## [8997] 450 450 450 450 451 451 451 451 451 451 451 451 451
## [9010] 451 451 451 451 451 451 451 451 451 451 451 452 452
## [9023] 452 452 452 452 452 452 452 452 452 452 452 452 452
## [9036] 452 452 452 452 452 453 453 453 453 453 453 453 453
## [9049] 453 453 453 453 453 453 453 453 453 453 453 453 454
## [9062] 454 454 454 454 454 454 454 454 454 454 454 454 454
## [9075] 454 454 454 454 454 454 455 455 455 455 455 455 455
## [9088] 455 455 455 455 455 455 455 455 455 455 455 455 455
## [9101] 456 456 456 456 456 456 456 456 456 456 456 456 456
## [9114] 456 456 456 456 456 456 456 457 457 457 457 457 457
## [9127] 457 457 457 457 457 457 457 457 457 457 457 457 457
## [9140] 457 458 458 458 458 458 458 458 458 458 458 458 458
## [9153] 458 458 458 458 458 458 458 458 459 459 459 459 459
## [9166] 459 459 459 459 459 459 459 459 459 459 459 459 459
## [9179] 459 459 460 460 460 460 460 460 460 460 460 460 460
## [9192] 460 460 460 460 460 460 460 460 460 461 461 461 461
## [9205] 461 461 461 461 461 461 461 461 461 461 461 461 461
## [9218] 461 461 461 462 462 462 462 462 462 462 462 462 462
## [9231] 462 462 462 462 462 462 462 462 462 462 463 463 463
## [9244] 463 463 463 463 463 463 463 463 463 463 463 463 463
## [9257] 463 463 463 463 464 464 464 464 464 464 464 464 464
## [9270] 464 464 464 464 464 464 464 464 464 464 464 465 465
## [9283] 465 465 465 465 465 465 465 465 465 465 465 465 465
## [9296] 465 465 465 465 465 466 466 466 466 466 466 466 466
## [9309] 466 466 466 466 466 466 466 466 466 466 466 466 467
## [9322] 467 467 467 467 467 467 467 467 467 467 467 467 467
## [9335] 467 467 467 467 467 467 468 468 468 468 468 468 468
## [9348] 468 468 468 468 468 468 468 468 468 468 468 468 468
## [9361] 469 469 469 469 469 469 469 469 469 469 469 469 469
## [9374] 469 469 469 469 469 469 469 470 470 470 470 470 470
## [9387] 470 470 470 470 470 470 470 470 470 470 470 470 470
## [9400] 470 471 471 471 471 471 471 471 471 471 471 471 471
## [9413] 471 471 471 471 471 471 471 471 472 472 472 472 472
## [9426] 472 472 472 472 472 472 472 472 472 472 472 472 472
## [9439] 472 472 473 473 473 473 473 473 473 473 473 473 473
## [9452] 473 473 473 473 473 473 473 473 473 474 474 474 474
## [9465] 474 474 474 474 474 474 474 474 474 474 474 474 474
## [9478] 474 474 474 475 475 475 475 475 475 475 475 475 475
## [9491] 475 475 475 475 475 475 475 475 475 475 476 476 476
## [9504] 476 476 476 476 476 476 476 476 476 476 476 476 476
## [9517] 476 476 476 476 477 477 477 477 477 477 477 477 477
## [9530] 477 477 477 477 477 477 477 477 477 477 477 478 478
## [9543] 478 478 478 478 478 478 478 478 478 478 478 478 478
## [9556] 478 478 478 478 478 479 479 479 479 479 479 479 479
## [9569] 479 479 479 479 479 479 479 479 479 479 479 479 480
## [9582] 480 480 480 480 480 480 480 480 480 480 480 480 480
## [9595] 480 480 480 480 480 480 481 481 481 481 481 481 481
## [9608] 481 481 481 481 481 481 481 481 481 481 481 481 481
## [9621] 482 482 482 482 482 482 482 482 482 482 482 482 482
## [9634] 482 482 482 482 482 482 482 483 483 483 483 483 483
## [9647] 483 483 483 483 483 483 483 483 483 483 483 483 483
## [9660] 483 484 484 484 484 484 484 484 484 484 484 484 484
## [9673] 484 484 484 484 484 484 484 484 485 485 485 485 485
## [9686] 485 485 485 485 485 485 485 485 485 485 485 485 485
## [9699] 485 485 486 486 486 486 486 486 486 486 486 486 486
## [9712] 486 486 486 486 486 486 486 486 486 487 487 487 487
## [9725] 487 487 487 487 487 487 487 487 487 487 487 487 487
## [9738] 487 487 487 488 488 488 488 488 488 488 488 488 488
## [9751] 488 488 488 488 488 488 488 488 488 488 489 489 489
## [9764] 489 489 489 489 489 489 489 489 489 489 489 489 489
## [9777] 489 489 489 489 490 490 490 490 490 490 490 490 490
## [9790] 490 490 490 490 490 490 490 490 490 490 490 491 491
## [9803] 491 491 491 491 491 491 491 491 491 491 491 491 491
## [9816] 491 491 491 491 491 492 492 492 492 492 492 492 492
## [9829] 492 492 492 492 492 492 492 492 492 492 492 492 493
## [9842] 493 493 493 493 493 493 493 493 493 493 493 493 493
## [9855] 493 493 493 493 493 493 494 494 494 494 494 494 494
## [9868] 494 494 494 494 494 494 494 494 494 494 494 494 494
## [9881] 495 495 495 495 495 495 495 495 495 495 495 495 495
## [9894] 495 495 495 495 495 495 495 496 496 496 496 496 496
## [9907] 496 496 496 496 496 496 496 496 496 496 496 496 496
## [9920] 496 497 497 497 497 497 497 497 497 497 497 497 497
## [9933] 497 497 497 497 497 497 497 497 498 498 498 498 498
## [9946] 498 498 498 498 498 498 498 498 498 498 498 498 498
## [9959] 498 498 499 499 499 499 499 499 499 499 499 499 499
## [9972] 499 499 499 499 499 499 499 499 499 500 500 500 500
## [9985] 500 500 500 500 500 500 500 500 500 500 500 500 500
## [9998] 500 500 500 501 501 501 501 501 501 501 501 501 501
## [10011] 501 501 501 501 501 501 501 501 501 501 502 502 502
## [10024] 502 502 502 502 502 502 502 502 502 502 502 502 502
## [10037] 502 502 502 502 503 503 503 503 503 503 503 503 503
## [10050] 503 503 503 503 503 503 503 503 503 503 503 504 504
## [10063] 504 504 504 504 504 504 504 504 504 504 504 504 504
## [10076] 504 504 504 504 504 505 505 505 505 505 505 505 505
## [10089] 505 505 505 505 505 505 505 505 505 505 505 505 506
## [10102] 506 506 506 506 506 506 506 506 506 506 506 506 506
## [10115] 506 506 506 506 506 506 507 507 507 507 507 507 507
## [10128] 507 507 507 507 507 507 507 507 507 507 507 507 507
## [10141] 508 508 508 508 508 508 508 508 508 508 508 508 508
## [10154] 508 508 508 508 508 508 508 509 509 509 509 509 509
## [10167] 509 509 509 509 509 509 509 509 509 509 509 509 509
## [10180] 509 510 510 510 510 510 510 510 510 510 510 510 510
## [10193] 510 510 510 510 510 510 510 510 511 511 511 511 511
## [10206] 511 511 511 511 511 511 511 511 511 511 511 511 511
## [10219] 511 511 512 512 512 512 512 512 512 512 512 512 512
## [10232] 512 512 512 512 512 512 512 512 512 513 513 513 513
## [10245] 513 513 513 513 513 513 513 513 513 513 513 513 513
## [10258] 513 513 513 514 514 514 514 514 514 514 514 514 514
## [10271] 514 514 514 514 514 514 514 514 514 514 515 515 515
## [10284] 515 515 515 515 515 515 515 515 515 515 515 515 515
## [10297] 515 515 515 515 516 516 516 516 516 516 516 516 516
## [10310] 516 516 516 516 516 516 516 516 516 516 516 517 517
## [10323] 517 517 517 517 517 517 517 517 517 517 517 517 517
## [10336] 517 517 517 517 517 518 518 518 518 518 518 518 518
## [10349] 518 518 518 518 518 518 518 518 518 518 518 518 519
## [10362] 519 519 519 519 519 519 519 519 519 519 519 519 519
## [10375] 519 519 519 519 519 519 520 520 520 520 520 520 520
## [10388] 520 520 520 520 520 520 520 520 520 520 520 520 520
## [10401] 521 521 521 521 521 521 521 521 521 521 521 521 521
## [10414] 521 521 521 521 521 521 521 522 522 522 522 522 522
## [10427] 522 522 522 522 522 522 522 522 522 522 522 522 522
## [10440] 522 523 523 523 523 523 523 523 523 523 523 523 523
## [10453] 523 523 523 523 523 523 523 523 524 524 524 524 524
## [10466] 524 524 524 524 524 524 524 524 524 524 524 524 524
## [10479] 524 524 525 525 525 525 525 525 525 525 525 525 525
## [10492] 525 525 525 525 525 525 525 525 525 526 526 526 526
## [10505] 526 526 526 526 526 526 526 526 526 526 526 526 526
## [10518] 526 526 526 527 527 527 527 527 527 527 527 527 527
## [10531] 527 527 527 527 527 527 527 527 527 527 528 528 528
## [10544] 528 528 528 528 528 528 528 528 528 528 528 528 528
## [10557] 528 528 528 528 529 529 529 529 529 529 529 529 529
## [10570] 529 529 529 529 529 529 529 529 529 529 529 530 530
## [10583] 530 530 530 530 530 530 530 530 530 530 530 530 530
## [10596] 530 530 530 530 530 531 531 531 531 531 531 531 531
## [10609] 531 531 531 531 531 531 531 531 531 531 531 531 532
## [10622] 532 532 532 532 532 532 532 532 532 532 532 532 532
## [10635] 532 532 532 532 532 532 533 533 533 533 533 533 533
## [10648] 533 533 533 533 533 533 533 533 533 533 533 533 533
## [10661] 534 534 534 534 534 534 534 534 534 534 534 534 534
## [10674] 534 534 534 534 534 534 534 535 535 535 535 535 535
## [10687] 535 535 535 535 535 535 535 535 535 535 535 535 535
## [10700] 535 536 536 536 536 536 536 536 536 536 536 536 536
## [10713] 536 536 536 536 536 536 536 536 537 537 537 537 537
## [10726] 537 537 537 537 537 537 537 537 537 537 537 537 537
## [10739] 537 537 538 538 538 538 538 538 538 538 538 538 538
## [10752] 538 538 538 538 538 538 538 538 538 539 539 539 539
## [10765] 539 539 539 539 539 539 539 539 539 539 539 539 539
## [10778] 539 539 539 540 540 540 540 540 540 540 540 540 540
## [10791] 540 540 540 540 540 540 540 540 540 540 541 541 541
## [10804] 541 541 541 541 541 541 541 541 541 541 541 541 541
## [10817] 541 541 541 541 542 542 542 542 542 542 542 542 542
## [10830] 542 542 542 542 542 542 542 542 542 542 542 543 543
## [10843] 543 543 543 543 543 543 543 543 543 543 543 543 543
## [10856] 543 543 543 543 543 544 544 544 544 544 544 544 544
## [10869] 544 544 544 544 544 544 544 544 544 544 544 544 545
## [10882] 545 545 545 545 545 545 545 545 545 545 545 545 545
## [10895] 545 545 545 545 545 545 546 546 546 546 546 546 546
## [10908] 546 546 546 546 546 546 546 546 546 546 546 546 546
## [10921] 547 547 547 547 547 547 547 547 547 547 547 547 547
## [10934] 547 547 547 547 547 547 547 548 548 548 548 548 548
## [10947] 548 548 548 548 548 548 548 548 548 548 548 548 548
## [10960] 548 549 549 549 549 549 549 549 549 549 549 549 549
## [10973] 549 549 549 549 549 549 549 549 550 550 550 550 550
## [10986] 550 550 550 550 550 550 550 550 550 550 550 550 550
## [10999] 550 550 551 551 551 551 551 551 551 551 551 551 551
## [11012] 551 551 551 551 551 551 551 551 551 552 552 552 552
## [11025] 552 552 552 552 552 552 552 552 552 552 552 552 552
## [11038] 552 552 552 553 553 553 553 553 553 553 553 553 553
## [11051] 553 553 553 553 553 553 553 553 553 553 554 554 554
## [11064] 554 554 554 554 554 554 554 554 554 554 554 554 554
## [11077] 554 554 554 554 555 555 555 555 555 555 555 555 555
## [11090] 555 555 555 555 555 555 555 555 555 555 555 556 556
## [11103] 556 556 556 556 556 556 556 556 556 556 556 556 556
## [11116] 556 556 556 556 556 557 557 557 557 557 557 557 557
## [11129] 557 557 557 557 557 557 557 557 557 557 557 557 558
## [11142] 558 558 558 558 558 558 558 558 558 558 558 558 558
## [11155] 558 558 558 558 558 558 559 559 559 559 559 559 559
## [11168] 559 559 559 559 559 559 559 559 559 559 559 559 559
## [11181] 560 560 560 560 560 560 560 560 560 560 560 560 560
## [11194] 560 560 560 560 560 560 560 561 561 561 561 561 561
## [11207] 561 561 561 561 561 561 561 561 561 561 561 561 561
## [11220] 561 562 562 562 562 562 562 562 562 562 562 562 562
## [11233] 562 562 562 562 562 562 562 562 563 563 563 563 563
## [11246] 563 563 563 563 563 563 563 563 563 563 563 563 563
## [11259] 563 563 564 564 564 564 564 564 564 564 564 564 564
## [11272] 564 564 564 564 564 564 564 564 564 565 565 565 565
## [11285] 565 565 565 565 565 565 565 565 565 565 565 565 565
## [11298] 565 565 565 566 566 566 566 566 566 566 566 566 566
## [11311] 566 566 566 566 566 566 566 566 566 566 567 567 567
## [11324] 567 567 567 567 567 567 567 567 567 567 567 567 567
## [11337] 567 567 567 567 568 568 568 568 568 568 568 568 568
## [11350] 568 568 568 568 568 568 568 568 568 568 568 569 569
## [11363] 569 569 569 569 569 569 569 569 569 569 569 569 569
## [11376] 569 569 569 569 569 570 570 570 570 570 570 570 570
## [11389] 570 570 570 570 570 570 570 570 570 570 570 570 571
## [11402] 571 571 571 571 571 571 571 571 571 571 571 571 571
## [11415] 571 571 571 571 571 571 572 572 572 572 572 572 572
## [11428] 572 572 572 572 572 572 572 572 572 572 572 572 572
## [11441] 573 573 573 573 573 573 573 573 573 573 573 573 573
## [11454] 573 573 573 573 573 573 573 574 574 574 574 574 574
## [11467] 574 574 574 574 574 574 574 574 574 574 574 574 574
## [11480] 574 575 575 575 575 575 575 575 575 575 575 575 575
## [11493] 575 575 575 575 575 575 575 575 576 576 576 576 576
## [11506] 576 576 576 576 576 576 576 576 576 576 576 576 576
## [11519] 576 576 577 577 577 577 577 577 577 577 577 577 577
## [11532] 577 577 577 577 577 577 577 577 577 578 578 578 578
## [11545] 578 578 578 578 578 578 578 578 578 578 578 578 578
## [11558] 578 578 578 579 579 579 579 579 579 579 579 579 579
## [11571] 579 579 579 579 579 579 579 579 579 579 580 580 580
## [11584] 580 580 580 580 580 580 580 580 580 580 580 580 580
## [11597] 580 580 580 580 581 581 581 581 581 581 581 581 581
## [11610] 581 581 581 581 581 581 581 581 581 581 581 582 582
## [11623] 582 582 582 582 582 582 582 582 582 582 582 582 582
## [11636] 582 582 582 582 582 583 583 583 583 583 583 583 583
## [11649] 583 583 583 583 583 583 583 583 583 583 583 583 584
## [11662] 584 584 584 584 584 584 584 584 584 584 584 584 584
## [11675] 584 584 584 584 584 584 585 585 585 585 585 585 585
## [11688] 585 585 585 585 585 585 585 585 585 585 585 585 585
## [11701] 586 586 586 586 586 586 586 586 586 586 586 586 586
## [11714] 586 586 586 586 586 586 586 587 587 587 587 587 587
## [11727] 587 587 587 587 587 587 587 587 587 587 587 587 587
## [11740] 587 588 588 588 588 588 588 588 588 588 588 588 588
## [11753] 588 588 588 588 588 588 588 588 589 589 589 589 589
## [11766] 589 589 589 589 589 589 589 589 589 589 589 589 589
## [11779] 589 589 590 590 590 590 590 590 590 590 590 590 590
## [11792] 590 590 590 590 590 590 590 590 590 591 591 591 591
## [11805] 591 591 591 591 591 591 591 591 591 591 591 591 591
## [11818] 591 591 591 592 592 592 592 592 592 592 592 592 592
## [11831] 592 592 592 592 592 592 592 592 592 592 593 593 593
## [11844] 593 593 593 593 593 593 593 593 593 593 593 593 593
## [11857] 593 593 593 593 594 594 594 594 594 594 594 594 594
## [11870] 594 594 594 594 594 594 594 594 594 594 594 595 595
## [11883] 595 595 595 595 595 595 595 595 595 595 595 595 595
## [11896] 595 595 595 595 595 596 596 596 596 596 596 596 596
## [11909] 596 596 596 596 596 596 596 596 596 596 596 596 597
## [11922] 597 597 597 597 597 597 597 597 597 597 597 597 597
## [11935] 597 597 597 597 597 597 598 598 598 598 598 598 598
## [11948] 598 598 598 598 598 598 598 598 598 598 598 598 598
## [11961] 599 599 599 599 599 599 599 599 599 599 599 599 599
## [11974] 599 599 599 599 599 599 599 600 600 600 600 600 600
## [11987] 600 600 600 600 600 600 600 600 600 600 600 600 600
## [12000] 600 601 601 601 601 601 601 601 601 601 601 601 601
## [12013] 601 601 601 601 601 601 601 601 602 602 602 602 602
## [12026] 602 602 602 602 602 602 602 602 602 602 602 602 602
## [12039] 602 602 603 603 603 603 603 603 603 603 603 603 603
## [12052] 603 603 603 603 603 603 603 603 603 604 604 604 604
## [12065] 604 604 604 604 604 604 604 604 604 604 604 604 604
## [12078] 604 604 604 605 605 605 605 605 605 605 605 605 605
## [12091] 605 605 605 605 605 605 605 605 605 605 606 606 606
## [12104] 606 606 606 606 606 606 606 606 606 606 606 606 606
## [12117] 606 606 606 606 607 607 607 607 607 607 607 607 607
## [12130] 607 607 607 607 607 607 607 607 607 607 607 608 608
## [12143] 608 608 608 608 608 608 608 608 608 608 608 608 608
## [12156] 608 608 608 608 608 609 609 609 609 609 609 609 609
## [12169] 609 609 609 609 609 609 609 609 609 609 609 609 610
## [12182] 610 610 610 610 610 610 610 610 610 610 610 610 610
## [12195] 610 610 610 610 610 610 611 611 611 611 611 611 611
## [12208] 611 611 611 611 611 611 611 611 611 611 611 611 611
## [12221] 612 612 612 612 612 612 612 612 612 612 612 612 612
## [12234] 612 612 612 612 612 612 612 613 613 613 613 613 613
## [12247] 613 613 613 613 613 613 613 613 613 613 613 613 613
## [12260] 613 614 614 614 614 614 614 614 614 614 614 614 614
## [12273] 614 614 614 614 614 614 614 614 615 615 615 615 615
## [12286] 615 615 615 615 615 615 615 615 615 615 615 615 615
## [12299] 615 615 616 616 616 616 616 616 616 616 616 616 616
## [12312] 616 616 616 616 616 616 616 616 616 617 617 617 617
## [12325] 617 617 617 617 617 617 617 617 617 617 617 617 617
## [12338] 617 617 617 618 618 618 618 618 618 618 618 618 618
## [12351] 618 618 618 618 618 618 618 618 618 618 619 619 619
## [12364] 619 619 619 619 619 619 619 619 619 619 619 619 619
## [12377] 619 619 619 619 620 620 620 620 620 620 620 620 620
## [12390] 620 620 620 620 620 620 620 620 620 620 620 621 621
## [12403] 621 621 621 621 621 621 621 621 621 621 621 621 621
## [12416] 621 621 621 621 621 622 622 622 622 622 622 622 622
## [12429] 622 622 622 622 622 622 622 622 622 622 622 622 623
## [12442] 623 623 623 623 623 623 623 623 623 623 623 623 623
## [12455] 623 623 623 623 623 623 624 624 624 624 624 624 624
## [12468] 624 624 624 624 624 624 624 624 624 624 624 624 624
## [12481] 625 625 625 625 625 625 625 625 625 625 625 625 625
## [12494] 625 625 625 625 625 625 625 626 626 626 626 626 626
## [12507] 626 626 626 626 626 626 626 626 626 626 626 626 626
## [12520] 626 627 627 627 627 627 627 627 627 627 627 627 627
## [12533] 627 627 627 627 627 627 627 627 628 628 628 628 628
## [12546] 628 628 628 628 628 628 628 628 628 628 628 628 628
## [12559] 628 628 629 629 629 629 629 629 629 629 629 629 629
## [12572] 629 629 629 629 629 629 629 629 629 630 630 630 630
## [12585] 630 630 630 630 630 630 630 630 630 630 630 630 630
## [12598] 630 630 630 631 631 631 631 631 631 631 631 631 631
## [12611] 631 631 631 631 631 631 631 631 631 631 632 632 632
## [12624] 632 632 632 632 632 632 632 632 632 632 632 632 632
## [12637] 632 632 632 632 633 633 633 633 633 633 633 633 633
## [12650] 633 633 633 633 633 633 633 633 633 633 633 634 634
## [12663] 634 634 634 634 634 634 634 634 634 634 634 634 634
## [12676] 634 634 634 634 634 635 635 635 635 635 635 635 635
## [12689] 635 635 635 635 635 635 635 635 635 635 635 635 636
## [12702] 636 636 636 636 636 636 636 636 636 636 636 636 636
## [12715] 636 636 636 636 636 636 637 637 637 637 637 637 637
## [12728] 637 637 637 637 637 637 637 637 637 637 637 637 637
## [12741] 638 638 638 638 638 638 638 638 638 638 638 638 638
## [12754] 638 638 638 638 638 638 638 639 639 639 639 639 639
## [12767] 639 639 639 639 639 639 639 639 639 639 639 639 639
## [12780] 639 640 640 640 640 640 640 640 640 640 640 640 640
## [12793] 640 640 640 640 640 640 640 640 641 641 641 641 641
## [12806] 641 641 641 641 641 641 641 641 641 641 641 641 641
## [12819] 641 641 642 642 642 642 642 642 642 642 642 642 642
## [12832] 642 642 642 642 642 642 642 642 642 643 643 643 643
## [12845] 643 643 643 643 643 643 643 643 643 643 643 643 643
## [12858] 643 643 643 644 644 644 644 644 644 644 644 644 644
## [12871] 644 644 644 644 644 644 644 644 644 644 645 645 645
## [12884] 645 645 645 645 645 645 645 645 645 645 645 645 645
## [12897] 645 645 645 645 646 646 646 646 646 646 646 646 646
## [12910] 646 646 646 646 646 646 646 646 646 646 646 647 647
## [12923] 647 647 647 647 647 647 647 647 647 647 647 647 647
## [12936] 647 647 647 647 647 648 648 648 648 648 648 648 648
## [12949] 648 648 648 648 648 648 648 648 648 648 648 648 649
## [12962] 649 649 649 649 649 649 649 649 649 649 649 649 649
## [12975] 649 649 649 649 649 649 650 650 650 650 650 650 650
## [12988] 650 650 650 650 650 650 650 650 650 650 650 650 650
## [13001] 651 651 651 651 651 651 651 651 651 651 651 651 651
## [13014] 651 651 651 651 651 651 651 652 652 652 652 652 652
## [13027] 652 652 652 652 652 652 652 652 652 652 652 652 652
## [13040] 652 653 653 653 653 653 653 653 653 653 653 653 653
## [13053] 653 653 653 653 653 653 653 653 654 654 654 654 654
## [13066] 654 654 654 654 654 654 654 654 654 654 654 654 654
## [13079] 654 654 655 655 655 655 655 655 655 655 655 655 655
## [13092] 655 655 655 655 655 655 655 655 655 656 656 656 656
## [13105] 656 656 656 656 656 656 656 656 656 656 656 656 656
## [13118] 656 656 656 657 657 657 657 657 657 657 657 657 657
## [13131] 657 657 657 657 657 657 657 657 657 657 658 658 658
## [13144] 658 658 658 658 658 658 658 658 658 658 658 658 658
## [13157] 658 658 658 658 659 659 659 659 659 659 659 659 659
## [13170] 659 659 659 659 659 659 659 659 659 659 659 660 660
## [13183] 660 660 660 660 660 660 660 660 660 660 660 660 660
## [13196] 660 660 660 660 660 661 661 661 661 661 661 661 661
## [13209] 661 661 661 661 661 661 661 661 661 661 661 661 662
## [13222] 662 662 662 662 662 662 662 662 662 662 662 662 662
## [13235] 662 662 662 662 662 662 663 663 663 663 663 663 663
## [13248] 663 663 663 663 663 663 663 663 663 663 663 663 663
## [13261] 664 664 664 664 664 664 664 664 664 664 664 664 664
## [13274] 664 664 664 664 664 664 664 665 665 665 665 665 665
## [13287] 665 665 665 665 665 665 665 665 665 665 665 665 665
## [13300] 665 666 666 666 666 666 666 666 666 666 666 666 666
## [13313] 666 666 666 666 666 666 666 666 667 667 667 667 667
## [13326] 667 667 667 667 667 667 667 667 667 667 667 667 667
## [13339] 667 667 668 668 668 668 668 668 668 668 668 668 668
## [13352] 668 668 668 668 668 668 668 668 668 669 669 669 669
## [13365] 669 669 669 669 669 669 669 669 669 669 669 669 669
## [13378] 669 669 669 670 670 670 670 670 670 670 670 670 670
## [13391] 670 670 670 670 670 670 670 670 670 670 671 671 671
## [13404] 671 671 671 671 671 671 671 671 671 671 671 671 671
## [13417] 671 671 671 671 672 672 672 672 672 672 672 672 672
## [13430] 672 672 672 672 672 672 672 672 672 672 672 673 673
## [13443] 673 673 673 673 673 673 673 673 673 673 673 673 673
## [13456] 673 673 673 673 673 674 674 674 674 674 674 674 674
## [13469] 674 674 674 674 674 674 674 674 674 674 674 674 675
## [13482] 675 675 675 675 675 675 675 675 675 675 675 675 675
## [13495] 675 675 675 675 675 675 676 676 676 676 676 676 676
## [13508] 676 676 676 676 676 676 676 676 676 676 676 676 676
## [13521] 677 677 677 677 677 677 677 677 677 677 677 677 677
## [13534] 677 677 677 677 677 677 677 678 678 678 678 678 678
## [13547] 678 678 678 678 678 678 678 678 678 678 678 678 678
## [13560] 678 679 679 679 679 679 679 679 679 679 679 679 679
## [13573] 679 679 679 679 679 679 679 679 680 680 680 680 680
## [13586] 680 680 680 680 680 680 680 680 680 680 680 680 680
## [13599] 680 680 681 681 681 681 681 681 681 681 681 681 681
## [13612] 681 681 681 681 681 681 681 681 681 682 682 682 682
## [13625] 682 682 682 682 682 682 682 682 682 682 682 682 682
## [13638] 682 682 682 683 683 683 683 683 683 683 683 683 683
## [13651] 683 683 683 683 683 683 683 683 683 683 684 684 684
## [13664] 684 684 684 684 684 684 684 684 684 684 684 684 684
## [13677] 684 684 684 684 685 685 685 685 685 685 685 685 685
## [13690] 685 685 685 685 685 685 685 685 685 685 685 686 686
## [13703] 686 686 686 686 686 686 686 686 686 686 686 686 686
## [13716] 686 686 686 686 686 687 687 687 687 687 687 687 687
## [13729] 687 687 687 687 687 687 687 687 687 687 687 687 688
## [13742] 688 688 688 688 688 688 688 688 688 688 688 688 688
## [13755] 688 688 688 688 688 688 689 689 689 689 689 689 689
## [13768] 689 689 689 689 689 689 689 689 689 689 689 689 689
## [13781] 690 690 690 690 690 690 690 690 690 690 690 690 690
## [13794] 690 690 690 690 690 690 690 691 691 691 691 691 691
## [13807] 691 691 691 691 691 691 691 691 691 691 691 691 691
## [13820] 691 692 692 692 692 692 692 692 692 692 692 692 692
## [13833] 692 692 692 692 692 692 692 692 693 693 693 693 693
## [13846] 693 693 693 693 693 693 693 693 693 693 693 693 693
## [13859] 693 693 694 694 694 694 694 694 694 694 694 694 694
## [13872] 694 694 694 694 694 694 694 694 694 695 695 695 695
## [13885] 695 695 695 695 695 695 695 695 695 695 695 695 695
## [13898] 695 695 695 696 696 696 696 696 696 696 696 696 696
## [13911] 696 696 696 696 696 696 696 696 696 696 697 697 697
## [13924] 697 697 697 697 697 697 697 697 697 697 697 697 697
## [13937] 697 697 697 697 698 698 698 698 698 698 698 698 698
## [13950] 698 698 698 698 698 698 698 698 698 698 698 699 699
## [13963] 699 699 699 699 699 699 699 699 699 699 699 699 699
## [13976] 699 699 699 699 699 700 700 700 700 700 700 700 700
## [13989] 700 700 700 700 700 700 700 700 700 700 700 700 701
## [14002] 701 701 701 701 701 701 701 701 701 701 701 701 701
## [14015] 701 701 701 701 701 701 702 702 702 702 702 702 702
## [14028] 702 702 702 702 702 702 702 702 702 702 702 702 702
## [14041] 703 703 703 703 703 703 703 703 703 703 703 703 703
## [14054] 703 703 703 703 703 703 703 704 704 704 704 704 704
## [14067] 704 704 704 704 704 704 704 704 704 704 704 704 704
## [14080] 704 705 705 705 705 705 705 705 705 705 705 705 705
## [14093] 705 705 705 705 705 705 705 705 706 706 706 706 706
## [14106] 706 706 706 706 706 706 706 706 706 706 706 706 706
## [14119] 706 706 707 707 707 707 707 707 707 707 707 707 707
## [14132] 707 707 707 707 707 707 707 707 707 708 708 708 708
## [14145] 708 708 708 708 708 708 708 708 708 708 708 708 708
## [14158] 708 708 708 709 709 709 709 709 709 709 709 709 709
## [14171] 709 709 709 709 709 709 709 709 709 709 710 710 710
## [14184] 710 710 710 710 710 710 710 710 710 710 710 710 710
## [14197] 710 710 710 710 711 711 711 711 711 711 711 711 711
## [14210] 711 711 711 711 711 711 711 711 711 711 711 712 712
## [14223] 712 712 712 712 712 712 712 712 712 712 712 712 712
## [14236] 712 712 712 712 712 713 713 713 713 713 713 713 713
## [14249] 713 713 713 713 713 713 713 713 713 713 713 713 714
## [14262] 714 714 714 714 714 714 714 714 714 714 714 714 714
## [14275] 714 714 714 714 714 714 715 715 715 715 715 715 715
## [14288] 715 715 715 715 715 715 715 715 715 715 715 715 715
## [14301] 716 716 716 716 716 716 716 716 716 716 716 716 716
## [14314] 716 716 716 716 716 716 716 717 717 717 717 717 717
## [14327] 717 717 717 717 717 717 717 717 717 717 717 717 717
## [14340] 717 718 718 718 718 718 718 718 718 718 718 718 718
## [14353] 718 718 718 718 718 718 718 718 719 719 719 719 719
## [14366] 719 719 719 719 719 719 719 719 719 719 719 719 719
## [14379] 719 719 720 720 720 720 720 720 720 720 720 720 720
## [14392] 720 720 720 720 720 720 720 720 720 721 721 721 721
## [14405] 721 721 721 721 721 721 721 721 721 721 721 721 721
## [14418] 721 721 721 722 722 722 722 722 722 722 722 722 722
## [14431] 722 722 722 722 722 722 722 722 722 722 723 723 723
## [14444] 723 723 723 723 723 723 723 723 723 723 723 723 723
## [14457] 723 723 723 723 724 724 724 724 724 724 724 724 724
## [14470] 724 724 724 724 724 724 724 724 724 724 724 725 725
## [14483] 725 725 725 725 725 725 725 725 725 725 725 725 725
## [14496] 725 725 725 725 725 726 726 726 726 726 726 726 726
## [14509] 726 726 726 726 726 726 726 726 726 726 726 726 727
## [14522] 727 727 727 727 727 727 727 727 727 727 727 727 727
## [14535] 727 727 727 727 727 727 728 728 728 728 728 728 728
## [14548] 728 728 728 728 728 728 728 728 728 728 728 728 728
## [14561] 729 729 729 729 729 729 729 729 729 729 729 729 729
## [14574] 729 729 729 729 729 729 729 730 730 730 730 730 730
## [14587] 730 730 730 730 730 730 730 730 730 730 730 730 730
## [14600] 730 731 731 731 731 731 731 731 731 731 731 731 731
## [14613] 731 731 731 731 731 731 731 731 732 732 732 732 732
## [14626] 732 732 732 732 732 732 732 732 732 732 732 732 732
## [14639] 732 732 733 733 733 733 733 733 733 733 733 733 733
## [14652] 733 733 733 733 733 733 733 733 733 734 734 734 734
## [14665] 734 734 734 734 734 734 734 734 734 734 734 734 734
## [14678] 734 734 734 735 735 735 735 735 735 735 735 735 735
## [14691] 735 735 735 735 735 735 735 735 735 735 736 736 736
## [14704] 736 736 736 736 736 736 736 736 736 736 736 736 736
## [14717] 736 736 736 736 737 737 737 737 737 737 737 737 737
## [14730] 737 737 737 737 737 737 737 737 737 737 737 738 738
## [14743] 738 738 738 738 738 738 738 738 738 738 738 738 738
## [14756] 738 738 738 738 738 739 739 739 739 739 739 739 739
## [14769] 739 739 739 739 739 739 739 739 739 739 739 739 740
## [14782] 740 740 740 740 740 740 740 740 740 740 740 740 740
## [14795] 740 740 740 740 740 740 741 741 741 741 741 741 741
## [14808] 741 741 741 741 741 741 741 741 741 741 741 741 741
## [14821] 742 742 742 742 742 742 742 742 742 742 742 742 742
## [14834] 742 742 742 742 742 742 742 743 743 743 743 743 743
## [14847] 743 743 743 743 743 743 743 743 743 743 743 743 743
## [14860] 743 744 744 744 744 744 744 744 744 744 744 744 744
## [14873] 744 744 744 744 744 744 744 744 745 745 745 745 745
## [14886] 745 745 745 745 745 745 745 745 745 745 745 745 745
## [14899] 745 745 746 746 746 746 746 746 746 746 746 746 746
## [14912] 746 746 746 746 746 746 746 746 746 747 747 747 747
## [14925] 747 747 747 747 747 747 747 747 747 747 747 747 747
## [14938] 747 747 747 748 748 748 748 748 748 748 748 748 748
## [14951] 748 748 748 748 748 748 748 748 748 748 749 749 749
## [14964] 749 749 749 749 749 749 749 749 749 749 749 749 749
## [14977] 749 749 749 749 750 750 750 750 750 750 750 750 750
## [14990] 750 750 750 750 750 750 750 750 750 750 750 751 751
## [15003] 751 751 751 751 751 751 751 751 751 751 751 751 751
## [15016] 751 751 751 751 751 752 752 752 752 752 752 752 752
## [15029] 752 752 752 752 752 752 752 752 752 752 752 752 753
## [15042] 753 753 753 753 753 753 753 753 753 753 753 753 753
## [15055] 753 753 753 753 753 753 754 754 754 754 754 754 754
## [15068] 754 754 754 754 754 754 754 754 754 754 754 754 754
## [15081] 755 755 755 755 755 755 755 755 755 755 755 755 755
## [15094] 755 755 755 755 755 755 755 756 756 756 756 756 756
## [15107] 756 756 756 756 756 756 756 756 756 756 756 756 756
## [15120] 756 757 757 757 757 757 757 757 757 757 757 757 757
## [15133] 757 757 757 757 757 757 757 757 758 758 758 758 758
## [15146] 758 758 758 758 758 758 758 758 758 758 758 758 758
## [15159] 758 758 759 759 759 759 759 759 759 759 759 759 759
## [15172] 759 759 759 759 759 759 759 759 759 760 760 760 760
## [15185] 760 760 760 760 760 760 760 760 760 760 760 760 760
## [15198] 760 760 760 761 761 761 761 761 761 761 761 761 761
## [15211] 761 761 761 761 761 761 761 761 761 761 762 762 762
## [15224] 762 762 762 762 762 762 762 762 762 762 762 762 762
## [15237] 762 762 762 762 763 763 763 763 763 763 763 763 763
## [15250] 763 763 763 763 763 763 763 763 763 763 763 764 764
## [15263] 764 764 764 764 764 764 764 764 764 764 764 764 764
## [15276] 764 764 764 764 764 765 765 765 765 765 765 765 765
## [15289] 765 765 765 765 765 765 765 765 765 765 765 765 766
## [15302] 766 766 766 766 766 766 766 766 766 766 766 766 766
## [15315] 766 766 766 766 766 766 767 767 767 767 767 767 767
## [15328] 767 767 767 767 767 767 767 767 767 767 767 767 767
## [15341] 768 768 768 768 768 768 768 768 768 768 768 768 768
## [15354] 768 768 768 768 768 768 768 769 769 769 769 769 769
## [15367] 769 769 769 769 769 769 769 769 769 769 769 769 769
## [15380] 769 770 770 770 770 770 770 770 770 770 770 770 770
## [15393] 770 770 770 770 770 770 770 770 771 771 771 771 771
## [15406] 771 771 771 771 771 771 771 771 771 771 771 771 771
## [15419] 771 771 772 772 772 772 772 772 772 772 772 772 772
## [15432] 772 772 772 772 772 772 772 772 772 773 773 773 773
## [15445] 773 773 773 773 773 773 773 773 773 773 773 773 773
## [15458] 773 773 773 774 774 774 774 774 774 774 774 774 774
## [15471] 774 774 774 774 774 774 774 774 774 774 775 775 775
## [15484] 775 775 775 775 775 775 775 775 775 775 775 775 775
## [15497] 775 775 775 775 776 776 776 776 776 776 776 776 776
## [15510] 776 776 776 776 776 776 776 776 776 776 776 777 777
## [15523] 777 777 777 777 777 777 777 777 777 777 777 777 777
## [15536] 777 777 777 777 777 778 778 778 778 778 778 778 778
## [15549] 778 778 778 778 778 778 778 778 778 778 778 778 779
## [15562] 779 779 779 779 779 779 779 779 779 779 779 779 779
## [15575] 779 779 779 779 779 779 780 780 780 780 780 780 780
## [15588] 780 780 780 780 780 780 780 780 780 780 780 780 780
## [15601] 781 781 781 781 781 781 781 781 781 781 781 781 781
## [15614] 781 781 781 781 781 781 781 782 782 782 782 782 782
## [15627] 782 782 782 782 782 782 782 782 782 782 782 782 782
## [15640] 782 783 783 783 783 783 783 783 783 783 783 783 783
## [15653] 783 783 783 783 783 783 783 783 784 784 784 784 784
## [15666] 784 784 784 784 784 784 784 784 784 784 784 784 784
## [15679] 784 784 785 785 785 785 785 785 785 785 785 785 785
## [15692] 785 785 785 785 785 785 785 785 785 786 786 786 786
## [15705] 786 786 786 786 786 786 786 786 786 786 786 786 786
## [15718] 786 786 786 787 787 787 787 787 787 787 787 787 787
## [15731] 787 787 787 787 787 787 787 787 787 787 788 788 788
## [15744] 788 788 788 788 788 788 788 788 788 788 788 788 788
## [15757] 788 788 788 788 789 789 789 789 789 789 789 789 789
## [15770] 789 789 789 789 789 789 789 789 789 789 789 790 790
## [15783] 790 790 790 790 790 790 790 790 790 790 790 790 790
## [15796] 790 790 790 790 790 791 791 791 791 791 791 791 791
## [15809] 791 791 791 791 791 791 791 791 791 791 791 791 792
## [15822] 792 792 792 792 792 792 792 792 792 792 792 792 792
## [15835] 792 792 792 792 792 792 793 793 793 793 793 793 793
## [15848] 793 793 793 793 793 793 793 793 793 793 793 793 793
## [15861] 794 794 794 794 794 794 794 794 794 794 794 794 794
## [15874] 794 794 794 794 794 794 794 795 795 795 795 795 795
## [15887] 795 795 795 795 795 795 795 795 795 795 795 795 795
## [15900] 795 796 796 796 796 796 796 796 796 796 796 796 796
## [15913] 796 796 796 796 796 796 796 796 797 797 797 797 797
## [15926] 797 797 797 797 797 797 797 797 797 797 797 797 797
## [15939] 797 797 798 798 798 798 798 798 798 798 798 798 798
## [15952] 798 798 798 798 798 798 798 798 798 799 799 799 799
## [15965] 799 799 799 799 799 799 799 799 799 799 799 799 799
## [15978] 799 799 799 800 800 800 800 800 800 800 800 800 800
## [15991] 800 800 800 800 800 800 800 800 800 800 801 801 801
## [16004] 801 801 801 801 801 801 801 801 801 801 801 801 801
## [16017] 801 801 801 801 802 802 802 802 802 802 802 802 802
## [16030] 802 802 802 802 802 802 802 802 802 802 802 803 803
## [16043] 803 803 803 803 803 803 803 803 803 803 803 803 803
## [16056] 803 803 803 803 803 804 804 804 804 804 804 804 804
## [16069] 804 804 804 804 804 804 804 804 804 804 804 804 805
## [16082] 805 805 805 805 805 805 805 805 805 805 805 805 805
## [16095] 805 805 805 805 805 805 806 806 806 806 806 806 806
## [16108] 806 806 806 806 806 806 806 806 806 806 806 806 806
## [16121] 807 807 807 807 807 807 807 807 807 807 807 807 807
## [16134] 807 807 807 807 807 807 807 808 808 808 808 808 808
## [16147] 808 808 808 808 808 808 808 808 808 808 808 808 808
## [16160] 808 809 809 809 809 809 809 809 809 809 809 809 809
## [16173] 809 809 809 809 809 809 809 809 810 810 810 810 810
## [16186] 810 810 810 810 810 810 810 810 810 810 810 810 810
## [16199] 810 810 811 811 811 811 811 811 811 811 811 811 811
## [16212] 811 811 811 811 811 811 811 811 811 812 812 812 812
## [16225] 812 812 812 812 812 812 812 812 812 812 812 812 812
## [16238] 812 812 812 813 813 813 813 813 813 813 813 813 813
## [16251] 813 813 813 813 813 813 813 813 813 813 814 814 814
## [16264] 814 814 814 814 814 814 814 814 814 814 814 814 814
## [16277] 814 814 814 814 815 815 815 815 815 815 815 815 815
## [16290] 815 815 815 815 815 815 815 815 815 815 815 816 816
## [16303] 816 816 816 816 816 816 816 816 816 816 816 816 816
## [16316] 816 816 816 816 816 817 817 817 817 817 817 817 817
## [16329] 817 817 817 817 817 817 817 817 817 817 817 817 818
## [16342] 818 818 818 818 818 818 818 818 818 818 818 818 818
## [16355] 818 818 818 818 818 818 819 819 819 819 819 819 819
## [16368] 819 819 819 819 819 819 819 819 819 819 819 819 819
## [16381] 820 820 820 820 820 820 820 820 820 820 820 820 820
## [16394] 820 820 820 820 820 820 820 821 821 821 821 821 821
## [16407] 821 821 821 821 821 821 821 821 821 821 821 821 821
## [16420] 821 822 822 822 822 822 822 822 822 822 822 822 822
## [16433] 822 822 822 822 822 822 822 822 823 823 823 823 823
## [16446] 823 823 823 823 823 823 823 823 823 823 823 823 823
## [16459] 823 823 824 824 824 824 824 824 824 824 824 824 824
## [16472] 824 824 824 824 824 824 824 824 824 825 825 825 825
## [16485] 825 825 825 825 825 825 825 825 825 825 825 825 825
## [16498] 825 825 825 826 826 826 826 826 826 826 826 826 826
## [16511] 826 826 826 826 826 826 826 826 826 826 827 827 827
## [16524] 827 827 827 827 827 827 827 827 827 827 827 827 827
## [16537] 827 827 827 827 828 828 828 828 828 828 828 828 828
## [16550] 828 828 828 828 828 828 828 828 828 828 828 829 829
## [16563] 829 829 829 829 829 829 829 829 829 829 829 829 829
## [16576] 829 829 829 829 829 830 830 830 830 830 830 830 830
## [16589] 830 830 830 830 830 830 830 830 830 830 830 830 831
## [16602] 831 831 831 831 831 831 831 831 831 831 831 831 831
## [16615] 831 831 831 831 831 831 832 832 832 832 832 832 832
## [16628] 832 832 832 832 832 832 832 832 832 832 832 832 832
## [16641] 833 833 833 833 833 833 833 833 833 833 833 833 833
## [16654] 833 833 833 833 833 833 833 834 834 834 834 834 834
## [16667] 834 834 834 834 834 834 834 834 834 834 834 834 834
## [16680] 834 835 835 835 835 835 835 835 835 835 835 835 835
## [16693] 835 835 835 835 835 835 835 835 836 836 836 836 836
## [16706] 836 836 836 836 836 836 836 836 836 836 836 836 836
## [16719] 836 836 837 837 837 837 837 837 837 837 837 837 837
## [16732] 837 837 837 837 837 837 837 837 837 838 838 838 838
## [16745] 838 838 838 838 838 838 838 838 838 838 838 838 838
## [16758] 838 838 838 839 839 839 839 839 839 839 839 839 839
## [16771] 839 839 839 839 839 839 839 839 839 839 840 840 840
## [16784] 840 840 840 840 840 840 840 840 840 840 840 840 840
## [16797] 840 840 840 840 841 841 841 841 841 841 841 841 841
## [16810] 841 841 841 841 841 841 841 841 841 841 841 842 842
## [16823] 842 842 842 842 842 842 842 842 842 842 842 842 842
## [16836] 842 842 842 842 842 843 843 843 843 843 843 843 843
## [16849] 843 843 843 843 843 843 843 843 843 843 843 843 844
## [16862] 844 844 844 844 844 844 844 844 844 844 844 844 844
## [16875] 844 844 844 844 844 844 845 845 845 845 845 845 845
## [16888] 845 845 845 845 845 845 845 845 845 845 845 845 845
## [16901] 846 846 846 846 846 846 846 846 846 846 846 846 846
## [16914] 846 846 846 846 846 846 846 847 847 847 847 847 847
## [16927] 847 847 847 847 847 847 847 847 847 847 847 847 847
## [16940] 847 848 848 848 848 848 848 848 848 848 848 848 848
## [16953] 848 848 848 848 848 848 848 848 849 849 849 849 849
## [16966] 849 849 849 849 849 849 849 849 849 849 849 849 849
## [16979] 849 849 850 850 850 850 850 850 850 850 850 850 850
## [16992] 850 850 850 850 850 850 850 850 850 851 851 851 851
## [17005] 851 851 851 851 851 851 851 851 851 851 851 851 851
## [17018] 851 851 851 852 852 852 852 852 852 852 852 852 852
## [17031] 852 852 852 852 852 852 852 852 852 852 853 853 853
## [17044] 853 853 853 853 853 853 853 853 853 853 853 853 853
## [17057] 853 853 853 853 854 854 854 854 854 854 854 854 854
## [17070] 854 854 854 854 854 854 854 854 854 854 854 855 855
## [17083] 855 855 855 855 855 855 855 855 855 855 855 855 855
## [17096] 855 855 855 855 855 856 856 856 856 856 856 856 856
## [17109] 856 856 856 856 856 856 856 856 856 856 856 856 857
## [17122] 857 857 857 857 857 857 857 857 857 857 857 857 857
## [17135] 857 857 857 857 857 857 858 858 858 858 858 858 858
## [17148] 858 858 858 858 858 858 858 858 858 858 858 858 858
## [17161] 859 859 859 859 859 859 859 859 859 859 859 859 859
## [17174] 859 859 859 859 859 859 859 860 860 860 860 860 860
## [17187] 860 860 860 860 860 860 860 860 860 860 860 860 860
## [17200] 860 861 861 861 861 861 861 861 861 861 861 861 861
## [17213] 861 861 861 861 861 861 861 861 862 862 862 862 862
## [17226] 862 862 862 862 862 862 862 862 862 862 862 862 862
## [17239] 862 862 863 863 863 863 863 863 863 863 863 863 863
## [17252] 863 863 863 863 863 863 863 863 863 864 864 864 864
## [17265] 864 864 864 864 864 864 864 864 864 864 864 864 864
## [17278] 864 864 864 865 865 865 865 865 865 865 865 865 865
## [17291] 865 865 865 865 865 865 865 865 865 865 866 866 866
## [17304] 866 866 866 866 866 866 866 866 866 866 866 866 866
## [17317] 866 866 866 866 867 867 867 867 867 867 867 867 867
## [17330] 867 867 867 867 867 867 867 867 867 867 867 868 868
## [17343] 868 868 868 868 868 868 868 868 868 868 868 868 868
## [17356] 868 868 868 868 868 869 869 869 869 869 869 869 869
## [17369] 869 869 869 869 869 869 869 869 869 869 869 869 870
## [17382] 870 870 870 870 870 870 870 870 870 870 870 870 870
## [17395] 870 870 870 870 870 870 871 871 871 871 871 871 871
## [17408] 871 871 871 871 871 871 871 871 871 871 871 871 871
## [17421] 872 872 872 872 872 872 872 872 872 872 872 872 872
## [17434] 872 872 872 872 872 872 872 873 873 873 873 873 873
## [17447] 873 873 873 873 873 873 873 873 873 873 873 873 873
## [17460] 873 874 874 874 874 874 874 874 874 874 874 874 874
## [17473] 874 874 874 874 874 874 874 874 875 875 875 875 875
## [17486] 875 875 875 875 875 875 875 875 875 875 875 875 875
## [17499] 875 875 876 876 876 876 876 876 876 876 876 876 876
## [17512] 876 876 876 876 876 876 876 876 876 877 877 877 877
## [17525] 877 877 877 877 877 877 877 877 877 877 877 877 877
## [17538] 877 877 877 878 878 878 878 878 878 878 878 878 878
## [17551] 878 878 878 878 878 878 878 878 878 878 879 879 879
## [17564] 879 879 879 879 879 879 879 879 879 879 879 879 879
## [17577] 879 879 879 879 880 880 880 880 880 880 880 880 880
## [17590] 880 880 880 880 880 880 880 880 880 880 880 881 881
## [17603] 881 881 881 881 881 881 881 881 881 881 881 881 881
## [17616] 881 881 881 881 881 882 882 882 882 882 882 882 882
## [17629] 882 882 882 882 882 882 882 882 882 882 882 882 883
## [17642] 883 883 883 883 883 883 883 883 883 883 883 883 883
## [17655] 883 883 883 883 883 883 884 884 884 884 884 884 884
## [17668] 884 884 884 884 884 884 884 884 884 884 884 884 884
## [17681] 885 885 885 885 885 885 885 885 885 885 885 885 885
## [17694] 885 885 885 885 885 885 885 886 886 886 886 886 886
## [17707] 886 886 886 886 886 886 886 886 886 886 886 886 886
## [17720] 886 887 887 887 887 887 887 887 887 887 887 887 887
## [17733] 887 887 887 887 887 887 887 887 888 888 888 888 888
## [17746] 888 888 888 888 888 888 888 888 888 888 888 888 888
## [17759] 888 888 889 889 889 889 889 889 889 889 889 889 889
## [17772] 889 889 889 889 889 889 889 889 889 890 890 890 890
## [17785] 890 890 890 890 890 890 890 890 890 890 890 890 890
## [17798] 890 890 890 891 891 891 891 891 891 891 891 891 891
## [17811] 891 891 891 891 891 891 891 891 891 891 892 892 892
## [17824] 892 892 892 892 892 892 892 892 892 892 892 892 892
## [17837] 892 892 892 892 893 893 893 893 893 893 893 893 893
## [17850] 893 893 893 893 893 893 893 893 893 893 893 894 894
## [17863] 894 894 894 894 894 894 894 894 894 894 894 894 894
## [17876] 894 894 894 894 894 895 895 895 895 895 895 895 895
## [17889] 895 895 895 895 895 895 895 895 895 895 895 895 896
## [17902] 896 896 896 896 896 896 896 896 896 896 896 896 896
## [17915] 896 896 896 896 896 896 897 897 897 897 897 897 897
## [17928] 897 897 897 897 897 897 897 897 897 897 897 897 897
## [17941] 898 898 898 898 898 898 898 898 898 898 898 898 898
## [17954] 898 898 898 898 898 898 898 899 899 899 899 899 899
## [17967] 899 899 899 899 899 899 899 899 899 899 899 899 899
## [17980] 899 900 900 900 900 900 900 900 900 900 900 900 900
## [17993] 900 900 900 900 900 900 900 900 901 901 901 901 901
## [18006] 901 901 901 901 901 901 901 901 901 901 901 901 901
## [18019] 901 901 902 902 902 902 902 902 902 902 902 902 902
## [18032] 902 902 902 902 902 902 902 902 902 903 903 903 903
## [18045] 903 903 903 903 903 903 903 903 903 903 903 903 903
## [18058] 903 903 903 904 904 904 904 904 904 904 904 904 904
## [18071] 904 904 904 904 904 904 904 904 904 904 905 905 905
## [18084] 905 905 905 905 905 905 905 905 905 905 905 905 905
## [18097] 905 905 905 905 906 906 906 906 906 906 906 906 906
## [18110] 906 906 906 906 906 906 906 906 906 906 906 907 907
## [18123] 907 907 907 907 907 907 907 907 907 907 907 907 907
## [18136] 907 907 907 907 907 908 908 908 908 908 908 908 908
## [18149] 908 908 908 908 908 908 908 908 908 908 908 908 909
## [18162] 909 909 909 909 909 909 909 909 909 909 909 909 909
## [18175] 909 909 909 909 909 909 910 910 910 910 910 910 910
## [18188] 910 910 910 910 910 910 910 910 910 910 910 910 910
## [18201] 911 911 911 911 911 911 911 911 911 911 911 911 911
## [18214] 911 911 911 911 911 911 911 912 912 912 912 912 912
## [18227] 912 912 912 912 912 912 912 912 912 912 912 912 912
## [18240] 912 913 913 913 913 913 913 913 913 913 913 913 913
## [18253] 913 913 913 913 913 913 913 913 914 914 914 914 914
## [18266] 914 914 914 914 914 914 914 914 914 914 914 914 914
## [18279] 914 914 915 915 915 915 915 915 915 915 915 915 915
## [18292] 915 915 915 915 915 915 915 915 915 916 916 916 916
## [18305] 916 916 916 916 916 916 916 916 916 916 916 916 916
## [18318] 916 916 916 917 917 917 917 917 917 917 917 917 917
## [18331] 917 917 917 917 917 917 917 917 917 917 918 918 918
## [18344] 918 918 918 918 918 918 918 918 918 918 918 918 918
## [18357] 918 918 918 918 919 919 919 919 919 919 919 919 919
## [18370] 919 919 919 919 919 919 919 919 919 919 919 920 920
## [18383] 920 920 920 920 920 920 920 920 920 920 920 920 920
## [18396] 920 920 920 920 920 921 921 921 921 921 921 921 921
## [18409] 921 921 921 921 921 921 921 921 921 921 921 921 922
## [18422] 922 922 922 922 922 922 922 922 922 922 922 922 922
## [18435] 922 922 922 922 922 922 923 923 923 923 923 923 923
## [18448] 923 923 923 923 923 923 923 923 923 923 923 923 923
## [18461] 924 924 924 924 924 924 924 924 924 924 924 924 924
## [18474] 924 924 924 924 924 924 924 925 925 925 925 925 925
## [18487] 925 925 925 925 925 925 925 925 925 925 925 925 925
## [18500] 925 926 926 926 926 926 926 926 926 926 926 926 926
## [18513] 926 926 926 926 926 926 926 926 927 927 927 927 927
## [18526] 927 927 927 927 927 927 927 927 927 927 927 927 927
## [18539] 927 927 928 928 928 928 928 928 928 928 928 928 928
## [18552] 928 928 928 928 928 928 928 928 928 929 929 929 929
## [18565] 929 929 929 929 929 929 929 929 929 929 929 929 929
## [18578] 929 929 929 930 930 930 930 930 930 930 930 930 930
## [18591] 930 930 930 930 930 930 930 930 930 930 931 931 931
## [18604] 931 931 931 931 931 931 931 931 931 931 931 931 931
## [18617] 931 931 931 931 932 932 932 932 932 932 932 932 932
## [18630] 932 932 932 932 932 932 932 932 932 932 932 933 933
## [18643] 933 933 933 933 933 933 933 933 933 933 933 933 933
## [18656] 933 933 933 933 933 934 934 934 934 934 934 934 934
## [18669] 934 934 934 934 934 934 934 934 934 934 934 934 935
## [18682] 935 935 935 935 935 935 935 935 935 935 935 935 935
## [18695] 935 935 935 935 935 935 936 936 936 936 936 936 936
## [18708] 936 936 936 936 936 936 936 936 936 936 936 936 936
## [18721] 937 937 937 937 937 937 937 937 937 937 937 937 937
## [18734] 937 937 937 937 937 937 937 938 938 938 938 938 938
## [18747] 938 938 938 938 938 938 938 938 938 938 938 938 938
## [18760] 938 939 939 939 939 939 939 939 939 939 939 939 939
## [18773] 939 939 939 939 939 939 939 939 940 940 940 940 940
## [18786] 940 940 940 940 940 940 940 940 940 940 940 940 940
## [18799] 940 940 941 941 941 941 941 941 941 941 941 941 941
## [18812] 941 941 941 941 941 941 941 941 941 942 942 942 942
## [18825] 942 942 942 942 942 942 942 942 942 942 942 942 942
## [18838] 942 942 942 943 943 943 943 943 943 943 943 943 943
## [18851] 943 943 943 943 943 943 943 943 943 943 944 944 944
## [18864] 944 944 944 944 944 944 944 944 944 944 944 944 944
## [18877] 944 944 944 944 945 945 945 945 945 945 945 945 945
## [18890] 945 945 945 945 945 945 945 945 945 945 945 946 946
## [18903] 946 946 946 946 946 946 946 946 946 946 946 946 946
## [18916] 946 946 946 946 946 947 947 947 947 947 947 947 947
## [18929] 947 947 947 947 947 947 947 947 947 947 947 947 948
## [18942] 948 948 948 948 948 948 948 948 948 948 948 948 948
## [18955] 948 948 948 948 948 948 949 949 949 949 949 949 949
## [18968] 949 949 949 949 949 949 949 949 949 949 949 949 949
## [18981] 950 950 950 950 950 950 950 950 950 950 950 950 950
## [18994] 950 950 950 950 950 950 950 951 951 951 951 951 951
## [19007] 951 951 951 951 951 951 951 951 951 951 951 951 951
## [19020] 951 952 952 952 952 952 952 952 952 952 952 952 952
## [19033] 952 952 952 952 952 952 952 952 953 953 953 953 953
## [19046] 953 953 953 953 953 953 953 953 953 953 953 953 953
## [19059] 953 953 954 954 954 954 954 954 954 954 954 954 954
## [19072] 954 954 954 954 954 954 954 954 954 955 955 955 955
## [19085] 955 955 955 955 955 955 955 955 955 955 955 955 955
## [19098] 955 955 955 956 956 956 956 956 956 956 956 956 956
## [19111] 956 956 956 956 956 956 956 956 956 956 957 957 957
## [19124] 957 957 957 957 957 957 957 957 957 957 957 957 957
## [19137] 957 957 957 957 958 958 958 958 958 958 958 958 958
## [19150] 958 958 958 958 958 958 958 958 958 958 958 959 959
## [19163] 959 959 959 959 959 959 959 959 959 959 959 959 959
## [19176] 959 959 959 959 959 960 960 960 960 960 960 960 960
## [19189] 960 960 960 960 960 960 960 960 960 960 960 960 961
## [19202] 961 961 961 961 961 961 961 961 961 961 961 961 961
## [19215] 961 961 961 961 961 961 962 962 962 962 962 962 962
## [19228] 962 962 962 962 962 962 962 962 962 962 962 962 962
## [19241] 963 963 963 963 963 963 963 963 963 963 963 963 963
## [19254] 963 963 963 963 963 963 963 964 964 964 964 964 964
## [19267] 964 964 964 964 964 964 964 964 964 964 964 964 964
## [19280] 964 965 965 965 965 965 965 965 965 965 965 965 965
## [19293] 965 965 965 965 965 965 965 965 966 966 966 966 966
## [19306] 966 966 966 966 966 966 966 966 966 966 966 966 966
## [19319] 966 966 967 967 967 967 967 967 967 967 967 967 967
## [19332] 967 967 967 967 967 967 967 967 967 968 968 968 968
## [19345] 968 968 968 968 968 968 968 968 968 968 968 968 968
## [19358] 968 968 968 969 969 969 969 969 969 969 969 969 969
## [19371] 969 969 969 969 969 969 969 969 969 969 970 970 970
## [19384] 970 970 970 970 970 970 970 970 970 970 970 970 970
## [19397] 970 970 970 970 971 971 971 971 971 971 971 971 971
## [19410] 971 971 971 971 971 971 971 971 971 971 971 972 972
## [19423] 972 972 972 972 972 972 972 972 972 972 972 972 972
## [19436] 972 972 972 972 972 973 973 973 973 973 973 973 973
## [19449] 973 973 973 973 973 973 973 973 973 973 973 973 974
## [19462] 974 974 974 974 974 974 974 974 974 974 974 974 974
## [19475] 974 974 974 974 974 974 975 975 975 975 975 975 975
## [19488] 975 975 975 975 975 975 975 975 975 975 975 975 975
## [19501] 976 976 976 976 976 976 976 976 976 976 976 976 976
## [19514] 976 976 976 976 976 976 976 977 977 977 977 977 977
## [19527] 977 977 977 977 977 977 977 977 977 977 977 977 977
## [19540] 977 978 978 978 978 978 978 978 978 978 978 978 978
## [19553] 978 978 978 978 978 978 978 978 979 979 979 979 979
## [19566] 979 979 979 979 979 979 979 979 979 979 979 979 979
## [19579] 979 979 980 980 980 980 980 980 980 980 980 980 980
## [19592] 980 980 980 980 980 980 980 980 980 981 981 981 981
## [19605] 981 981 981 981 981 981 981 981 981 981 981 981 981
## [19618] 981 981 981 982 982 982 982 982 982 982 982 982 982
## [19631] 982 982 982 982 982 982 982 982 982 982 983 983 983
## [19644] 983 983 983 983 983 983 983 983 983 983 983 983 983
## [19657] 983 983 983 983 984 984 984 984 984 984 984 984 984
## [19670] 984 984 984 984 984 984 984 984 984 984 984 985 985
## [19683] 985 985 985 985 985 985 985 985 985 985 985 985 985
## [19696] 985 985 985 985 985 986 986 986 986 986 986 986 986
## [19709] 986 986 986 986 986 986 986 986 986 986 986 986 987
## [19722] 987 987 987 987 987 987 987 987 987 987 987 987 987
## [19735] 987 987 987 987 987 987 988 988 988 988 988 988 988
## [19748] 988 988 988 988 988 988 988 988 988 988 988 988 988
## [19761] 989 989 989 989 989 989 989 989 989 989 989 989 989
## [19774] 989 989 989 989 989 989 989 990 990 990 990 990 990
## [19787] 990 990 990 990 990 990 990 990 990 990 990 990 990
## [19800] 990 991 991 991 991 991 991 991 991 991 991 991 991
## [19813] 991 991 991 991 991 991 991 991 992 992 992 992 992
## [19826] 992 992 992 992 992 992 992 992 992 992 992 992 992
## [19839] 992 992 993 993 993 993 993 993 993 993 993 993 993
## [19852] 993 993 993 993 993 993 993 993 993 994 994 994 994
## [19865] 994 994 994 994 994 994 994 994 994 994 994 994 994
## [19878] 994 994 994 995 995 995 995 995 995 995 995 995 995
## [19891] 995 995 995 995 995 995 995 995 995 995 996 996 996
## [19904] 996 996 996 996 996 996 996 996 996 996 996 996 996
## [19917] 996 996 996 996 997 997 997 997 997 997 997 997 997
## [19930] 997 997 997 997 997 997 997 997 997 997 997 998 998
## [19943] 998 998 998 998 998 998 998 998 998 998 998 998 998
## [19956] 998 998 998 998 998 999 999 999 999 999 999 999 999
## [19969] 999 999 999 999 999 999 999 999 999 999 999 999 1000
## [19982] 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
## [19995] 1000 1000 1000 1000 1000 1000
## id
## 1: 1
## 2: 1
## 3: 1
## 4: 1
## 5: 1
## ---
## 19996: 1000
## 19997: 1000
## 19998: 1000
## 19999: 1000
## 20000: 1000
## id
## 1: 1
## 2: 1
## 3: 1
## 4: 1
## 5: 1
## ---
## 19996: 1000
## 19997: 1000
## 19998: 1000
## 19999: 1000
## 20000: 1000
# Columns id, year and lq
fdata[, list(id, year, lq)]
## id year lq
## 1: 1 1 9.582840
## 2: 1 2 9.600721
## 3: 1 3 9.784685
## 4: 1 4 9.768900
## 5: 1 5 9.739995
## ---
## 19996: 1000 16 10.304541
## 19997: 1000 17 10.431146
## 19998: 1000 18 10.230929
## 19999: 1000 19 10.529978
## 20000: 1000 20 10.484083
## id year lq
## 1: 1 1 9.582840
## 2: 1 2 9.600721
## 3: 1 3 9.784685
## 4: 1 4 9.768900
## 5: 1 5 9.739995
## ---
## 19996: 1000 16 10.304541
## 19997: 1000 17 10.431146
## 19998: 1000 18 10.230929
## 19999: 1000 19 10.529978
## 20000: 1000 20 10.484083
# Computing columns
fdata[, mean(lq)]
## [1] 10.19558
## V1
## 1: 10.19558
fdata[, list(mean(lq), sd(lq))]
## V1 V2
## 1: 10.19558 0.4421356
fdata[, list(mean(lq), mean(ll))]
## V1 V2
## 1: 10.19558 4.304935
# Assigning names to computed columns
fdata[, list(mean_lq = mean(lq), sd_lq = sd(lq))]
## mean_lq sd_lq
## 1: 10.19558 0.4421356
Selecting rows and columns
# Rows 3-5 of column id
fdata[3:5, id]
## [1] 1 1 1
# Rows 3-5 of columns id, year and lq
fdata[3:5, .(id, year, lq)]
## id year lq
## 1: 1 3 9.784685
## 2: 1 4 9.768900
## 3: 1 5 9.739995
# id, year, lq and lk for missing lk's
fdata[is.na(lk), .(id, year, lq, lk)]
## Empty data.table (0 rows) of 4 cols: id,year,lq,lk
# First 10 rows of id, year, lq and lk for missing lk's
fdata[is.na(lk), .(id, year, lq, lk)][1:10]
## id year lq lk
## 1: NA NA NA NA
## 2: NA NA NA NA
## 3: NA NA NA NA
## 4: NA NA NA NA
## 5: NA NA NA NA
## 6: NA NA NA NA
## 7: NA NA NA NA
## 8: NA NA NA NA
## 9: NA NA NA NA
## 10: NA NA NA NA
Rows and columns by groups
# Doing j by group
sectq <- fdata[, list(mean.lq = mean(lq), sd.lq = sd(lq)), by = sect]
sectq
## sect mean.lq sd.lq
## 1: 1 10.19695 0.4469915
## 2: 2 10.18742 0.4377487
## 3: 3 10.20416 0.4448391
## 4: 4 10.19486 0.4443234
## 5: 5 10.21022 0.4410429
## 6: 6 10.19225 0.4353711
## 7: 7 10.18851 0.4326625
## 8: 8 10.17703 0.4432725
## 9: 9 10.20973 0.4429429
## 10: 10 10.19466 0.4517431
# Doing j by several groups
sectq <- fdata[, list(mean.lq = mean(lq), sd.lq = sd(lq)),
by = list(sect, year)]
sectq
## sect year mean.lq sd.lq
## 1: 1 1 9.531483 0.1094064
## 2: 1 2 9.610779 0.1267184
## 3: 1 3 9.659092 0.1335480
## 4: 1 4 9.736916 0.1424732
## 5: 1 5 9.816365 0.1543937
## ---
## 196: 10 16 10.593747 0.2533653
## 197: 10 17 10.641686 0.2507292
## 198: 10 18 10.712511 0.2542738
## 199: 10 19 10.783951 0.2658066
## 200: 10 20 10.854523 0.2803954
# Call functions in by
fdata[, list(mean.lq = mean(lq), sd.lq = sd(lq)), by = list(year>5)]
## year mean.lq sd.lq
## 1: FALSE 9.674969 0.1624902
## 2: TRUE 10.369114 0.3624671
fdata[, list(mean.lq = mean(lq), sd.lq = sd(lq)), by = list(sect, year>5)]
## sect year mean.lq sd.lq
## 1: 1 FALSE 9.670927 0.1661829
## 2: 1 TRUE 10.372284 0.3663270
## 3: 2 FALSE 9.663978 0.1597828
## 4: 2 TRUE 10.361894 0.3538207
## 5: 3 FALSE 9.675822 0.1600963
## 6: 3 TRUE 10.380271 0.3622282
## 7: 4 FALSE 9.676876 0.1681837
## 8: 4 TRUE 10.367524 0.3667846
## 9: 5 FALSE 9.684569 0.1638380
## 10: 5 TRUE 10.385431 0.3571861
## 11: 6 FALSE 9.684314 0.1621445
## 12: 6 TRUE 10.361562 0.3595526
## 13: 7 FALSE 9.680890 0.1576551
## 14: 7 TRUE 10.357717 0.3560383
## 15: 8 FALSE 9.654526 0.1622511
## 16: 8 TRUE 10.351193 0.3631124
## 17: 9 FALSE 9.688940 0.1645010
## 18: 9 TRUE 10.383332 0.3633289
## 19: 10 FALSE 9.668854 0.1583017
## 20: 10 TRUE 10.369933 0.3752853
# Grouping only on a subset
fdata[1:5, list(mean.lq = mean(lq), sd.lq = sd(lq)), by= year]
## year mean.lq sd.lq
## 1: 1 9.582840 NA
## 2: 2 9.600721 NA
## 3: 3 9.784685 NA
## 4: 4 9.768900 NA
## 5: 5 9.739995 NA
# Using .N to get the total number of observations
fdata[!is.na(lk), .N, by = sect]
## sect N
## 1: 1 2000
## 2: 2 2000
## 3: 3 2000
## 4: 4 2000
## 5: 5 2000
## 6: 6 2000
## 7: 7 2000
## 8: 8 2000
## 9: 9 2000
## 10: 10 2000
fdata[!is.na(lk), .N, keyby = sect]
## sect N
## 1: 1 2000
## 2: 2 2000
## 3: 3 2000
## 4: 4 2000
## 5: 5 2000
## 6: 6 2000
## 7: 7 2000
## 8: 8 2000
## 9: 9 2000
## 10: 10 2000
fdata[!is.na(lk), .N, by = sect][order(-sect)]
## sect N
## 1: 10 2000
## 2: 9 2000
## 3: 8 2000
## 4: 7 2000
## 5: 6 2000
## 6: 5 2000
## 7: 4 2000
## 8: 3 2000
## 9: 2 2000
## 10: 1 2000
Fast row subset
setkey(fdata, sect, year)
fdata[1, mean(lq)]
## [1] 9.58284
fdata[1, mean(lq), by = year]
## year V1
## 1: 1 9.58284
## [1] 9.659092
fdata[.(1, 3), mean(lq), by = year]
## year V1
## 1: 3 9.659092
# ad hoc subset
fdata[.(1), mean(lq), on = "sect"]
## [1] 10.19695
fdata[.(1, 3), mean(lq), on = c("sect", "year")]
## [1] 9.659092
Creating/updating variables
## id year sect re lk ll lm lA lq
## 1: 1 1 1 1.549448 3.862172 3.774901 3.839995 3.880978 9.582840
## 2: 2 1 1 1.549448 3.911743 3.884606 3.802311 3.913071 9.352522
## 3: 3 1 1 1.549448 3.905252 3.934949 3.917053 3.931591 9.486415
## 4: 4 1 1 1.549448 3.987684 3.848051 3.882126 3.944045 9.446215
## 5: 5 1 1 1.549448 3.805230 3.862148 3.877565 3.899033 9.631169
## 6: 6 1 1 1.549448 3.886235 3.901880 3.880069 3.889287 9.484189
fdata[, Q1 := exp(lq)]
fdata[, large := ifelse(lq > 10, 1, 0)]
fdata[, large2 := lq > 10]
fdata[, extralarge := lq > 10 & ll > 4]
head(fdata)
## id year sect re lk ll lm lA lq
## 1: 1 1 1 1.549448 3.862172 3.774901 3.839995 3.880978 9.582840
## 2: 2 1 1 1.549448 3.911743 3.884606 3.802311 3.913071 9.352522
## 3: 3 1 1 1.549448 3.905252 3.934949 3.917053 3.931591 9.486415
## 4: 4 1 1 1.549448 3.987684 3.848051 3.882126 3.944045 9.446215
## 5: 5 1 1 1.549448 3.805230 3.862148 3.877565 3.899033 9.631169
## 6: 6 1 1 1.549448 3.886235 3.901880 3.880069 3.889287 9.484189
## Q1 large large2 extralarge
## 1: 14513.58 0 FALSE FALSE
## 2: 11527.86 0 FALSE FALSE
## 3: 13179.46 0 FALSE FALSE
## 4: 12660.16 0 FALSE FALSE
## 5: 15232.24 0 FALSE FALSE
## 6: 13150.16 0 FALSE FALSE
# Be careful: DT <- DT[...] is redundant. The variables are added into DT
fdata[, c("K", "L", "M") := list(exp(lk), exp(ll), exp(lm))]
head(fdata)
## id year sect re lk ll lm lA lq
## 1: 1 1 1 1.549448 3.862172 3.774901 3.839995 3.880978 9.582840
## 2: 2 1 1 1.549448 3.911743 3.884606 3.802311 3.913071 9.352522
## 3: 3 1 1 1.549448 3.905252 3.934949 3.917053 3.931591 9.486415
## 4: 4 1 1 1.549448 3.987684 3.848051 3.882126 3.944045 9.446215
## 5: 5 1 1 1.549448 3.805230 3.862148 3.877565 3.899033 9.631169
## 6: 6 1 1 1.549448 3.886235 3.901880 3.880069 3.889287 9.484189
## Q1 large large2 extralarge K L M
## 1: 14513.58 0 FALSE FALSE 47.56854 43.59320 46.52526
## 2: 11527.86 0 FALSE FALSE 49.98602 48.64775 44.80460
## 3: 13179.46 0 FALSE FALSE 49.66261 51.15955 50.25213
## 4: 12660.16 0 FALSE FALSE 53.92985 46.90155 48.52730
## 5: 15232.24 0 FALSE FALSE 44.93560 47.56742 48.30644
## 6: 13150.16 0 FALSE FALSE 48.72709 49.49542 48.42757
Changing variable names
## [1] "id" "year" "sect" "re" "lk"
## [6] "ll" "lm" "lA" "lq" "Q1"
## [11] "large" "large2" "extralarge" "K" "L"
## [16] "M"
setnames(fdata, "Q1", "Q")
names(fdata)
## [1] "id" "year" "sect" "re" "lk"
## [6] "ll" "lm" "lA" "lq" "Q"
## [11] "large" "large2" "extralarge" "K" "L"
## [16] "M"
# Sorting
fdata[order(-year)]
## id year sect re lk ll lm lA
## 1: 1 20 1 1.549448 4.605257 4.494579 4.411256 4.681515
## 2: 2 20 1 1.549448 4.918969 4.646606 4.350668 4.648884
## 3: 3 20 1 1.549448 4.649388 4.585628 4.236517 4.569557
## 4: 4 20 1 1.549448 4.168683 4.558516 4.348127 4.623504
## 5: 5 20 1 1.549448 3.921117 4.825542 4.371843 4.818680
## ---
## 19996: 996 1 10 1.549448 3.873763 3.884820 3.923239 3.883875
## 19997: 997 1 10 1.549448 3.967849 3.877679 3.923888 3.900817
## 19998: 998 1 10 1.549448 3.862811 3.917578 3.900807 3.879504
## 19999: 999 1 10 1.549448 3.996398 3.846305 3.945201 3.956596
## 20000: 1000 1 10 1.549448 3.869087 3.891409 3.873159 3.899274
## lq Q large large2 extralarge K L
## 1: 10.925218 55559.96 1 TRUE TRUE 100.00864 89.53051
## 2: 10.982912 58859.69 1 TRUE TRUE 136.86145 104.23067
## 3: 10.855711 51829.28 1 TRUE TRUE 104.52101 98.06476
## 4: 10.749249 46595.02 1 TRUE TRUE 64.63028 95.44170
## 5: 11.184850 72030.86 1 TRUE TRUE 50.45676 124.65396
## ---
## 19996: 9.472405 12996.10 0 FALSE FALSE 48.12314 48.65817
## 19997: 9.718873 16628.49 0 FALSE FALSE 52.87067 48.31194
## 19998: 9.682799 16039.32 0 FALSE FALSE 47.59895 50.27854
## 19999: 9.631278 15233.89 0 FALSE FALSE 54.40182 46.81975
## 20000: 9.691597 16181.07 0 FALSE FALSE 47.89866 48.97987
## M
## 1: 82.37285
## 2: 77.53023
## 3: 69.16650
## 4: 77.33347
## 5: 79.18941
## ---
## 19996: 50.56396
## 19997: 50.59678
## 19998: 49.44236
## 19999: 51.68673
## 20000: 48.09405
Chaining
fdata[, list(sum(Q), round(sum(L, na.rm = TRUE), 0)), by = year]
## year V1 V2
## 1: 1 13963060 49630
## 2: 2 14969435 51850
## 3: 3 15989231 54175
## 4: 4 17245769 56498
## 5: 5 18474463 59034
## 6: 6 19931278 61661
## 7: 7 21287052 64330
## 8: 8 22782065 67136
## 9: 9 24493255 69946
## 10: 10 26229945 72826
## 11: 11 28153923 76123
## 12: 12 30280745 79615
## 13: 13 32537075 82977
## 14: 14 34913483 86832
## 15: 15 37470938 90596
## 16: 16 40262567 94478
## 17: 17 43038900 98502
## 18: 18 46234144 103017
## 19: 19 49869149 107413
## 20: 20 53412670 112377
fdata[, list(sum(Q), round(sum(L, na.rm = TRUE), 0)), by = year][order(-year)]
## year V1 V2
## 1: 20 53412670 112377
## 2: 19 49869149 107413
## 3: 18 46234144 103017
## 4: 17 43038900 98502
## 5: 16 40262567 94478
## 6: 15 37470938 90596
## 7: 14 34913483 86832
## 8: 13 32537075 82977
## 9: 12 30280745 79615
## 10: 11 28153923 76123
## 11: 10 26229945 72826
## 12: 9 24493255 69946
## 13: 8 22782065 67136
## 14: 7 21287052 64330
## 15: 6 19931278 61661
## 16: 5 18474463 59034
## 17: 4 17245769 56498
## 18: 3 15989231 54175
## 19: 2 14969435 51850
## 20: 1 13963060 49630
Melting and casting
- Melting: Transforming the dataset into long format
- Casting: Transforming the dataset into wide format
# Sectoral prices
nt <- 10
set.seed(123)
sPrice <- data.table(
year = c(1:nt),
p1 = 1 + cumsum(rnorm(nt, 0.02, 0.04)),
p2 = 1 + cumsum(rnorm(nt, 0.03, 0.05)),
p3 = 1 + cumsum(rnorm(nt, 0.04, 0.02)),
p4 = 1 + cumsum(rnorm(nt, 0.05, 0.03)),
p5 = 1 + cumsum(rnorm(nt, 0.02, 0.04)),
p6 = 1 + cumsum(rnorm(nt, 0.03, 0.06)),
p7 = 1 + cumsum(rnorm(nt, 0.04, 0.02)),
p8 = 1 + cumsum(rnorm(nt, 0.05, 0.07)),
p9 = 1 + cumsum(rnorm(nt, 0.05, 0.04)),
p10 = 1 + cumsum(rnorm(nt, 0.05, 0.06)))
head(sPrice)
## year p1 p2 p3 p4 p5 p6 p7
## 1: 1 0.997581 1.091204 1.018644 1.062794 0.9922117 1.045199 1.047593
## 2: 2 1.008374 1.139195 1.054284 1.103942 1.0038950 1.073486 1.077546
## 3: 3 1.090722 1.189233 1.073764 1.180796 0.9732792 1.100914 1.110882
## 4: 4 1.113543 1.224767 1.099186 1.257140 1.0800374 1.213030 1.130511
## 5: 5 1.138714 1.226975 1.126685 1.331787 1.1483559 1.229484 1.149075
## 6: 6 1.227317 1.346321 1.132951 1.402446 1.1234316 1.350472 1.195145
## p8 p9 p10
## 1: 1.015628 1.050231 1.109610
## 2: 0.903986 1.115642 1.192514
## 3: 1.024388 1.150815 1.256838
## 4: 1.024744 1.226590 1.269164
## 5: 1.026583 1.267771 1.400803
## 6: 1.148373 1.331042 1.414787
sPriceM <- melt(sPrice, id.vars=1, variable.name = "sectName", value.name = "price")
head(sPriceM)
## year sectName price
## 1: 1 p1 0.997581
## 2: 2 p1 1.008374
## 3: 3 p1 1.090722
## 4: 4 p1 1.113543
## 5: 5 p1 1.138714
## 6: 6 p1 1.227317
sPriceM$sect <- as.numeric(substring(sPriceM$sectName, 2))
head(sPriceM)
## year sectName price sect
## 1: 1 p1 0.997581 1
## 2: 2 p1 1.008374 1
## 3: 3 p1 1.090722 1
## 4: 4 p1 1.113543 1
## 5: 5 p1 1.138714 1
## 6: 6 p1 1.227317 1
sPriceM[, sectName := NULL]
ggplot(sPriceM, aes(x=year, y=price, group=sect, color=as.factor(sect))) + geom_line() +
ggtitle("Price indices")
setkey(sPriceM, sect, year)
setkey(fdata, sect, year)
fdataM <- fdata[sPriceM]
head(fdataM)
## id year sect re lk ll lm lA lq
## 1: 1 1 1 1.549448 3.862172 3.774901 3.839995 3.880978 9.582840
## 2: 2 1 1 1.549448 3.911743 3.884606 3.802311 3.913071 9.352522
## 3: 3 1 1 1.549448 3.905252 3.934949 3.917053 3.931591 9.486415
## 4: 4 1 1 1.549448 3.987684 3.848051 3.882126 3.944045 9.446215
## 5: 5 1 1 1.549448 3.805230 3.862148 3.877565 3.899033 9.631169
## 6: 6 1 1 1.549448 3.886235 3.901880 3.880069 3.889287 9.484189
## Q large large2 extralarge K L M price
## 1: 14513.58 0 FALSE FALSE 47.56854 43.59320 46.52526 0.997581
## 2: 11527.86 0 FALSE FALSE 49.98602 48.64775 44.80460 0.997581
## 3: 13179.46 0 FALSE FALSE 49.66261 51.15955 50.25213 0.997581
## 4: 12660.16 0 FALSE FALSE 53.92985 46.90155 48.52730 0.997581
## 5: 15232.24 0 FALSE FALSE 44.93560 47.56742 48.30644 0.997581
## 6: 13150.16 0 FALSE FALSE 48.72709 49.49542 48.42757 0.997581
Merging datasets
- Create sector level variables
sdata <- fdata[, list(sQ = sum(Q, na.rm=T), sL = sum(L, na.rm=T), sM = sum(M, na.rm=T)),
by = list(sect, year)]
summary(sdata)
## sect year sQ sL
## Min. : 1.0 Min. : 1.00 Min. :1382622 Min. : 4942
## 1st Qu.: 3.0 1st Qu.: 5.75 1st Qu.:1933496 1st Qu.: 6055
## Median : 5.5 Median :10.50 Median :2709047 Median : 7468
## Mean : 5.5 Mean :10.50 Mean :2957696 Mean : 7695
## 3rd Qu.: 8.0 3rd Qu.:15.25 3rd Qu.:3873070 3rd Qu.: 9239
## Max. :10.0 Max. :20.00 Max. :5424697 Max. :11435
## sM
## Min. :4858
## 1st Qu.:5439
## Median :6181
## Mean :6227
## 3rd Qu.:6948
## Max. :8107
## [1] 200 5
# Set keys for datasets
setkey(fdata, sect, year)
setkey(sdata, sect, year)
# Merge it with the firm level data
hdata <- fdata[sdata]
# That's all...
dim(hdata)
## [1] 20000 19
## [1] 20000 16
identical(hdata$ll, fdata$ll)
## [1] TRUE
# Shorter way
fdata[, c("sQ", "sL", "sM") := list(sum(Q, na.rm=T), sum(L, na.rm=T), sum(M, na.rm=T)),
by = list(sect, year)]
identical(hdata, fdata)
## [1] FALSE
# Find deviations from sectoral geometric averages
fdata[, lqdev := lq - mean(lq, na.rm=T), by = sect]
fdata[, lkdev := lk - mean(lk, na.rm=T), by = sect]
fdata[, lldev := ll - mean(ll, na.rm=T), by = sect]
fdata[, lmdev := lm - mean(lm, na.rm=T), by = sect]
# Find the growth rate of output
setkey(fdata, id, year)
fdata[, grq := lq - shift(lq, 1), by = id]
names(fdata)
## [1] "id" "year" "sect" "re" "lk"
## [6] "ll" "lm" "lA" "lq" "Q"
## [11] "large" "large2" "extralarge" "K" "L"
## [16] "M" "sQ" "sL" "sM" "lqdev"
## [21] "lkdev" "lldev" "lmdev" "grq"
fdata[, mean(grq, na.rm = TRUE)]
## [1] 0.06946501
Merging datasets
DT <- data.table(x=rep(c("b","a","c"),each=3), y=c(1,3,6), v=1:9)
X <- data.table(x=c("c","b"), v=8:7, foo=c(4,2))
DT[X, on="x"] # right join
## x y v i.v foo
## 1: c 1 7 8 4
## 2: c 3 8 8 4
## 3: c 6 9 8 4
## 4: b 1 1 7 2
## 5: b 3 2 7 2
## 6: b 6 3 7 2
X[DT, on="x"] # left join
## x v foo y i.v
## 1: b 7 2 1 1
## 2: b 7 2 3 2
## 3: b 7 2 6 3
## 4: a NA NA 1 4
## 5: a NA NA 3 5
## 6: a NA NA 6 6
## 7: c 8 4 1 7
## 8: c 8 4 3 8
## 9: c 8 4 6 9
DT[X, on="x", nomatch=0] # inner join
## x y v i.v foo
## 1: c 1 7 8 4
## 2: c 3 8 8 4
## 3: c 6 9 8 4
## 4: b 1 1 7 2
## 5: b 3 2 7 2
## 6: b 6 3 7 2
DT[!X, on="x"] # not join
## x y v
## 1: a 1 4
## 2: a 3 5
## 3: a 6 6
DT[X, on=c(y="v")] # join DT$y to X$v
## x y v i.x foo
## 1: <NA> 8 NA c 4
## 2: <NA> 7 NA b 2
## x y v i.x foo
## 1: <NA> 8 NA c 4
## 2: <NA> 7 NA b 2
DT[X, on=.(y<=foo)] # non-equi join
## x y v i.x i.v
## 1: b 4 1 c 8
## 2: b 4 2 c 8
## 3: a 4 4 c 8
## 4: a 4 5 c 8
## 5: c 4 7 c 8
## 6: c 4 8 c 8
## 7: b 2 1 b 7
## 8: a 2 4 b 7
## 9: c 2 7 b 7
DT[X, on="y<=foo"] # same
## x y v i.x i.v
## 1: b 4 1 c 8
## 2: b 4 2 c 8
## 3: a 4 4 c 8
## 4: a 4 5 c 8
## 5: c 4 7 c 8
## 6: c 4 8 c 8
## 7: b 2 1 b 7
## 8: a 2 4 b 7
## 9: c 2 7 b 7
DT[X, on=c("y<=foo")] # same
## x y v i.x i.v
## 1: b 4 1 c 8
## 2: b 4 2 c 8
## 3: a 4 4 c 8
## 4: a 4 5 c 8
## 5: c 4 7 c 8
## 6: c 4 8 c 8
## 7: b 2 1 b 7
## 8: a 2 4 b 7
## 9: c 2 7 b 7
DT[X, on=.(y>=foo)] # non-equi join
## x y v i.x i.v
## 1: b 4 3 c 8
## 2: a 4 6 c 8
## 3: c 4 9 c 8
## 4: b 2 2 b 7
## 5: b 2 3 b 7
## 6: a 2 5 b 7
## 7: a 2 6 b 7
## 8: c 2 8 b 7
## 9: c 2 9 b 7
DT[X, on=.(x, y<=foo)] # non-equi join
## x y v i.v
## 1: c 4 7 8
## 2: c 4 8 8
## 3: b 2 1 7
DT[X, .(x,y,x.y,v), on=.(x, y>=foo)] # select x's join columns as well
## x y x.y v
## 1: c 4 6 9
## 2: b 2 3 2
## 3: b 2 6 3
DT[X, on="x", mult="first"] # first row of each group
## x y v i.v foo
## 1: c 1 7 8 4
## 2: b 1 1 7 2
DT[X, on="x", mult="last"] # last row of each group
## x y v i.v foo
## 1: c 6 9 8 4
## 2: b 6 3 7 2
Exercises
- Get the Penn World Table 9 dataset
- Find the url of the data base
- Use the haven library, and read.stata function
- Make it data table
- Check the structure of the dataset
- Calculate GDP per capita for all countries
- Use the rgdpo and pop variables
- Check the structure of the dataset
- Drop Bermuda (BMU) from the sample
- The countrycode variable keeps country codes
- Find the average GDP per capita for small and large countries for 2010
- Small means population less than 1 million
- Find the average GDP per capita for poor and rich countries for all years after 2010
- Poor means GDP per capita less than 10,000 USD
- Create a GDP per capita data set for Turkey, US, Korea and China
- Find GDP per capita relative to the US level for Turkey, Korea and China
- Plot relative GDP per capita for these countries
Exercises
# Load libraries
library(haven)
library(data.table)
library(ggplot2)
# Get the data
pwt <- read_stata("http://www.rug.nl/ggdc/docs/pwt90.dta")
pwt <- data.table(pwt)
# Structure
# str(pwt)
# GDP per capita
pwt[, gdpc := rgdpo / pop]
str(pwt)
## Classes 'data.table' and 'data.frame': 11830 obs. of 48 variables:
## $ countrycode : chr "ABW" "ABW" "ABW" "ABW" ...
## ..- attr(*, "label")= chr "3-letter ISO country code"
## ..- attr(*, "format.stata")= chr "%9s"
## $ country : chr "Aruba" "Aruba" "Aruba" "Aruba" ...
## ..- attr(*, "label")= chr "Country name"
## ..- attr(*, "format.stata")= chr "%34s"
## $ currency_unit: chr "Aruban Guilder" "Aruban Guilder" "Aruban Guilder" "Aruban Guilder" ...
## ..- attr(*, "label")= chr "Currency unit"
## ..- attr(*, "format.stata")= chr "%29s"
## $ year : num 1950 1951 1952 1953 1954 ...
## ..- attr(*, "label")= chr "Year"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ rgdpe : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Expenditure-side real GDP at chained PPPs (in mil. 2011US$)"
## ..- attr(*, "format.stata")= chr "%14.3g"
## $ rgdpo : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Output-side real GDP at chained PPPs (in mil. 2011US$)"
## ..- attr(*, "format.stata")= chr "%14.3g"
## $ pop : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Population (in millions)"
## ..- attr(*, "format.stata")= chr "%10.0g"
## $ emp : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Number of persons engaged (in millions)"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ avh : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Average annual hours worked by persons engaged (source: The Conference Board)"
## ..- attr(*, "format.stata")= chr "%10.0gc"
## $ hc : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Human capital index, see note hc"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ ccon : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Real consumption of households and government, at current PPPs (in mil. 2011US$)"
## ..- attr(*, "format.stata")= chr "%14.3g"
## $ cda : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Real domestic absorption, see note cda"
## ..- attr(*, "format.stata")= chr "%14.3g"
## $ cgdpe : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Expenditure-side real GDP at current PPPs (in mil. 2011US$)"
## ..- attr(*, "format.stata")= chr "%14.3g"
## $ cgdpo : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Output-side real GDP at current PPPs (in mil. 2011US$)"
## ..- attr(*, "format.stata")= chr "%14.3g"
## $ ck : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Capital stock at current PPPs (in mil. 2011US$)"
## ..- attr(*, "format.stata")= chr "%14.3g"
## $ ctfp : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "TFP level at current PPPs (USA=1)"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ cwtfp : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Welfare-relevant TFP levels at current PPPs (USA=1)"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ rgdpna : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Real GDP at constant 2011 national prices (in mil. 2011US$)"
## ..- attr(*, "format.stata")= chr "%14.3g"
## $ rconna : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Real consumption at constant 2011 national prices (in mil. 2011US$)"
## ..- attr(*, "format.stata")= chr "%14.3g"
## $ rdana : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Real domestic absorption at constant 2011 national prices (in mil. 2011US$)"
## ..- attr(*, "format.stata")= chr "%14.3g"
## $ rkna : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Capital stock at constant 2011 national prices (in mil. 2011US$)"
## ..- attr(*, "format.stata")= chr "%14.3g"
## $ rtfpna : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "TFP at constant national prices (2011=1)"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ rwtfpna : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Welfare-relevant TFP at constant national prices (2011=1)"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ labsh : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Share of labour compensation in GDP at current national prices"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ delta : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Average depreciation rate of the capital stock"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ xr : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Exchange rate, national currency/USD (market+estimated)"
## ..- attr(*, "format.stata")= chr "%10.0g"
## $ pl_con : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Price level of CCON (PPP/XR), price level of USA GDPo in 2011=1"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ pl_da : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Price level of CDA (PPP/XR), price level of USA GDPo in 2011=1"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ pl_gdpo : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Price level of CGDPo (PPP/XR), price level of USA GDPo in 2011=1"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ i_cig : 'labelled' num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "0/1/2, see note i_cig"
## ..- attr(*, "format.stata")= chr "%12.0g"
## ..- attr(*, "labels")= Named num 0 1 2
## .. ..- attr(*, "names")= chr "Extrapolated" "Benchmark" "Interpolated"
## $ i_xm : 'labelled' num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "0/1/2, see note i_xm"
## ..- attr(*, "format.stata")= chr "%12.0g"
## ..- attr(*, "labels")= Named num 0 1 2
## .. ..- attr(*, "names")= chr "Extrapolated" "Benchmark" "Interpolated"
## $ i_xr : 'labelled' num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "0/1: the exchange rate is market-based (0) or estimated (1)"
## ..- attr(*, "format.stata")= chr "%9.0g"
## ..- attr(*, "labels")= Named num 0 1
## .. ..- attr(*, "names")= chr "Market" "Estimated"
## $ i_outlier : 'labelled' num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "0/1, see note i_outlier"
## ..- attr(*, "format.stata")= chr "%9.0g"
## ..- attr(*, "labels")= Named num 0 1
## .. ..- attr(*, "names")= chr "Regular" "Outlier"
## $ cor_exp : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Correlation between expenditure shares, see note cor_exp"
## ..- attr(*, "format.stata")= chr "%10.0g"
## $ statcap : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Statistical capacity indicator (source: World Bank, developing countries only)"
## ..- attr(*, "format.stata")= chr "%10.0g"
## $ csh_c : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Share of household consumption at current PPPs"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ csh_i : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Share of gross capital formation at current PPPs"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ csh_g : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Share of government consumption at current PPPs"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ csh_x : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Share of merchandise exports at current PPPs"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ csh_m : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Share of merchandise imports at current PPPs"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ csh_r : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Share of residual trade and GDP statistical discrepancy at current PPPs"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ pl_c : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Price level of household consumption, price level of USA GDPo in 2011=1"
## ..- attr(*, "format.stata")= chr "%10.0g"
## $ pl_i : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Price level of capital formation, price level of USA GDPo in 2011=1"
## ..- attr(*, "format.stata")= chr "%10.0g"
## $ pl_g : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Price level of government consumption, price level of USA GDPo in 2011=1"
## ..- attr(*, "format.stata")= chr "%10.0g"
## $ pl_x : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Price level of exports, price level of USA GDPo in 2011=1"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ pl_m : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Price level of imports, price level of USA GDPo in 2011=1"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ pl_k : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Price level of the capital stock, price level of USA 2011=1"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ gdpc : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Output-side real GDP at chained PPPs (in mil. 2011US$)"
## ..- attr(*, "format.stata")= chr "%14.3g"
## - attr(*, ".internal.selfref")=<externalptr>
# Drop Bermuda
pwt <- pwt[countrycode != "BMU"]
str(pwt)
## Classes 'data.table' and 'data.frame': 11765 obs. of 48 variables:
## $ countrycode : chr "ABW" "ABW" "ABW" "ABW" ...
## ..- attr(*, "label")= chr "3-letter ISO country code"
## ..- attr(*, "format.stata")= chr "%9s"
## $ country : chr "Aruba" "Aruba" "Aruba" "Aruba" ...
## ..- attr(*, "label")= chr "Country name"
## ..- attr(*, "format.stata")= chr "%34s"
## $ currency_unit: chr "Aruban Guilder" "Aruban Guilder" "Aruban Guilder" "Aruban Guilder" ...
## ..- attr(*, "label")= chr "Currency unit"
## ..- attr(*, "format.stata")= chr "%29s"
## $ year : num 1950 1951 1952 1953 1954 ...
## ..- attr(*, "label")= chr "Year"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ rgdpe : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Expenditure-side real GDP at chained PPPs (in mil. 2011US$)"
## ..- attr(*, "format.stata")= chr "%14.3g"
## $ rgdpo : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Output-side real GDP at chained PPPs (in mil. 2011US$)"
## ..- attr(*, "format.stata")= chr "%14.3g"
## $ pop : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Population (in millions)"
## ..- attr(*, "format.stata")= chr "%10.0g"
## $ emp : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Number of persons engaged (in millions)"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ avh : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Average annual hours worked by persons engaged (source: The Conference Board)"
## ..- attr(*, "format.stata")= chr "%10.0gc"
## $ hc : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Human capital index, see note hc"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ ccon : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Real consumption of households and government, at current PPPs (in mil. 2011US$)"
## ..- attr(*, "format.stata")= chr "%14.3g"
## $ cda : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Real domestic absorption, see note cda"
## ..- attr(*, "format.stata")= chr "%14.3g"
## $ cgdpe : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Expenditure-side real GDP at current PPPs (in mil. 2011US$)"
## ..- attr(*, "format.stata")= chr "%14.3g"
## $ cgdpo : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Output-side real GDP at current PPPs (in mil. 2011US$)"
## ..- attr(*, "format.stata")= chr "%14.3g"
## $ ck : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Capital stock at current PPPs (in mil. 2011US$)"
## ..- attr(*, "format.stata")= chr "%14.3g"
## $ ctfp : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "TFP level at current PPPs (USA=1)"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ cwtfp : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Welfare-relevant TFP levels at current PPPs (USA=1)"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ rgdpna : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Real GDP at constant 2011 national prices (in mil. 2011US$)"
## ..- attr(*, "format.stata")= chr "%14.3g"
## $ rconna : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Real consumption at constant 2011 national prices (in mil. 2011US$)"
## ..- attr(*, "format.stata")= chr "%14.3g"
## $ rdana : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Real domestic absorption at constant 2011 national prices (in mil. 2011US$)"
## ..- attr(*, "format.stata")= chr "%14.3g"
## $ rkna : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Capital stock at constant 2011 national prices (in mil. 2011US$)"
## ..- attr(*, "format.stata")= chr "%14.3g"
## $ rtfpna : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "TFP at constant national prices (2011=1)"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ rwtfpna : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Welfare-relevant TFP at constant national prices (2011=1)"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ labsh : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Share of labour compensation in GDP at current national prices"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ delta : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Average depreciation rate of the capital stock"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ xr : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Exchange rate, national currency/USD (market+estimated)"
## ..- attr(*, "format.stata")= chr "%10.0g"
## $ pl_con : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Price level of CCON (PPP/XR), price level of USA GDPo in 2011=1"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ pl_da : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Price level of CDA (PPP/XR), price level of USA GDPo in 2011=1"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ pl_gdpo : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Price level of CGDPo (PPP/XR), price level of USA GDPo in 2011=1"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ i_cig : 'labelled' num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "0/1/2, see note i_cig"
## ..- attr(*, "format.stata")= chr "%12.0g"
## ..- attr(*, "labels")= Named num 0 1 2
## .. ..- attr(*, "names")= chr "Extrapolated" "Benchmark" "Interpolated"
## $ i_xm : 'labelled' num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "0/1/2, see note i_xm"
## ..- attr(*, "format.stata")= chr "%12.0g"
## ..- attr(*, "labels")= Named num 0 1 2
## .. ..- attr(*, "names")= chr "Extrapolated" "Benchmark" "Interpolated"
## $ i_xr : 'labelled' num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "0/1: the exchange rate is market-based (0) or estimated (1)"
## ..- attr(*, "format.stata")= chr "%9.0g"
## ..- attr(*, "labels")= Named num 0 1
## .. ..- attr(*, "names")= chr "Market" "Estimated"
## $ i_outlier : 'labelled' num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "0/1, see note i_outlier"
## ..- attr(*, "format.stata")= chr "%9.0g"
## ..- attr(*, "labels")= Named num 0 1
## .. ..- attr(*, "names")= chr "Regular" "Outlier"
## $ cor_exp : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Correlation between expenditure shares, see note cor_exp"
## ..- attr(*, "format.stata")= chr "%10.0g"
## $ statcap : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Statistical capacity indicator (source: World Bank, developing countries only)"
## ..- attr(*, "format.stata")= chr "%10.0g"
## $ csh_c : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Share of household consumption at current PPPs"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ csh_i : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Share of gross capital formation at current PPPs"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ csh_g : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Share of government consumption at current PPPs"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ csh_x : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Share of merchandise exports at current PPPs"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ csh_m : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Share of merchandise imports at current PPPs"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ csh_r : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Share of residual trade and GDP statistical discrepancy at current PPPs"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ pl_c : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Price level of household consumption, price level of USA GDPo in 2011=1"
## ..- attr(*, "format.stata")= chr "%10.0g"
## $ pl_i : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Price level of capital formation, price level of USA GDPo in 2011=1"
## ..- attr(*, "format.stata")= chr "%10.0g"
## $ pl_g : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Price level of government consumption, price level of USA GDPo in 2011=1"
## ..- attr(*, "format.stata")= chr "%10.0g"
## $ pl_x : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Price level of exports, price level of USA GDPo in 2011=1"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ pl_m : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Price level of imports, price level of USA GDPo in 2011=1"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ pl_k : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Price level of the capital stock, price level of USA 2011=1"
## ..- attr(*, "format.stata")= chr "%9.0g"
## $ gdpc : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "Output-side real GDP at chained PPPs (in mil. 2011US$)"
## ..- attr(*, "format.stata")= chr "%14.3g"
## - attr(*, ".internal.selfref")=<externalptr>
# GDP per capita
pwt[year == 2010, mean(gdpc, na.rm = TRUE), by = list(year, pop > 1)]
## year pop V1
## 1: 2010 FALSE 21780.64
## 2: 2010 TRUE 16232.58
# GDP per capita
pwt[year > 2010, mean(gdpc, na.rm = TRUE), by = list(year, gdpc > 10000)]
## year gdpc > 10000 V1
## 1: 2011 TRUE 29432.780
## 2: 2012 TRUE 29638.135
## 3: 2013 TRUE 30292.655
## 4: 2014 TRUE 29189.817
## 5: 2011 FALSE 4367.621
## 6: 2012 FALSE 4413.643
## 7: 2013 FALSE 4377.382
## 8: 2014 FALSE 4096.922
# GDP per capita for selected countries
pwt4 <- pwt[countrycode %in% c("TUR", "USA", "CHN", "KOR"),
list(countrycode, gdpc, year)]
table(pwt4$countrycode)
##
## CHN KOR TUR USA
## 65 65 65 65
# Relative GDP per capita
pwt4[countrycode == "USA", usgdpc := gdpc]
pwt4[, relgdpc := gdpc / sum(usgdpc, na.rm = T), by = year]
ggplot(pwt4, aes(x = year, y = relgdpc, color = countrycode)) + geom_line(size = 1.5)
