Systemfit Question

Previous Topic Next Topic
 
classic Classic list List threaded Threaded
2 messages Options
Reply | Threaded
Open this post in threaded view
|

Systemfit Question

Ngwinui Azenui
I can't get my simultaneous equations to work using system fit. Please help.

#Reproducible script
Empdata<- read.csv("/Users/ngwinuiazenui/Documents/UPLOADemp.csv")
View(Empdata)
str(Empdata)
Empdata$gnipc<-as.numeric(Empdata$gnipc)
install.packages("systemfit")
library("systemfit")
pdata <- plm.data(Empdata, index=c("country","year"))
install.packages("plm")
library("plm")

#Simultaneous equation 1
eqprd<- prd ~ in_emp +gnipc
eqemp<- in_emp ~ prd
inst <- ~ gnipc
system<- list(productivity = eqprd, industrialemp = eqemp)
fit2sls<- systemfit(system, "2SLS", inst = inst, data = Empdata)
#I get an error code when I run the last line

#Simultaneous equation 2
eqprd2<- prd ~  in_emp + lag(prd, -1) + gnipc
eqemp2<- in_emp ~ prd + lag(prd, -1) + conflict
system2<- list(productivity = eqprd2, industrialemp = eqemp2)
inst <- ~ lag(prd, -1) + gnipc + conflict
fit2sls<- systemfit(system2, "2SLS", inst = inst, data = Empdata)
#I also wanted to include year and country fix effects(string
variables) to the model in "Simultaneous equation 2)

EmpSur<- systemfit(prd ~ in_emp + ag_emp, method = "SUR", data = pdata)
#When I run the code below, I get a table with reults by country
print(EmpSur)
#An except of the output from running <print(EmpSur)
#caf_(Intercept)      caf_in_emp      caf_ag_emp cod_(Intercept)
#     -0.13243213     14.25003985     -8.57346398     -0.11581997
#     cod_in_emp      cod_ag_emp eri_(Intercept)      eri_in_emp
#     11.58115777    -13.29628753     -0.55615848     23.69284487
#     eri_ag_emp eth_(Intercept)      eth_in_emp      eth_ag_emp
#     -6.26339787      0.17436006     34.87782981     -4.80423119
#gin_(Intercept)      gin_in_emp      gin_ag_emp gmb_(Intercept)
#     0.01122243    -11.22002009     -0.64218211     -0.07104175
#     gmb_in_emp      gmb_ag_emp gnb_(Intercept)      gnb_in_emp
#     9.74177468      3.29511466     -0.07681684     17.00346867
#     gnb_ag_emp gnq_(Intercept)      gnq_in_emp      gnq_ag_emp
#     -3.40077876     -0.25706504     -1.09274682     -7.39122680

#Below is the data. Sorry that it's a long panel dataset
dput(Empdata)
Empdata <- structure(list(country = structure(c(1L, 1L, 1L, 1L, 1L, 1L,
                                                1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
                                                1L, 1L, 1L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
                                                3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L,
                                                4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
                                                4L, 4L, 4L, 4L, 4L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
                                                2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 5L, 5L,
                                                5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
                                                5L, 5L, 5L, 5L, 5L,
5L, 5L, 26L, 26L, 26L, 26L, 26L, 26L, 26L,
                                                26L, 26L, 26L, 26L,
26L, 26L, 26L, 26L, 26L, 26L, 26L, 26L, 26L,
                                                26L, 26L, 26L, 26L,
26L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L,
                                                6L, 6L, 6L, 6L, 6L,
6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L,
                                                12L, 12L, 12L, 12L,
12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L,
                                                12L, 12L, 12L, 12L,
12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 7L,
                                                7L, 7L, 7L, 7L, 7L,
7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L,
                                                7L, 7L, 7L, 7L, 7L,
7L, 7L, 7L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L,
                                                8L, 8L, 8L, 8L, 8L,
8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L,
                                                8L, 10L, 10L, 10L,
10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L,
                                                10L, 10L, 10L, 10L,
10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L,
                                                9L, 9L, 9L, 9L, 9L,
9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L,
                                                9L, 9L, 9L, 9L, 9L,
9L, 9L, 9L, 9L, 11L, 11L, 11L, 11L, 11L,
                                                11L, 11L, 11L, 11L,
11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L,
                                                11L, 11L, 11L, 11L,
11L, 11L, 11L, 14L, 14L, 14L, 14L, 14L, 14L,
                                                14L, 14L, 14L, 14L,
14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L,
                                                14L, 14L, 14L, 14L,
14L, 14L, 13L, 13L, 13L, 13L, 13L, 13L, 13L,
                                                13L, 13L, 13L, 13L,
13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L,
                                                13L, 13L, 13L, 13L,
13L, 15L, 15L, 15L, 15L, 15L, 15L, 15L, 15L,
                                                15L, 15L, 15L, 15L,
15L, 15L, 15L, 15L, 15L, 15L, 15L, 15L, 15L,
                                                15L, 15L, 15L, 15L,
19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L,
                                                19L, 19L, 19L, 19L,
19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L,
                                                19L, 19L, 19L, 16L,
16L, 16L, 16L, 16L, 16L, 16L, 16L, 16L, 16L,
                                                16L, 16L, 16L, 16L,
16L, 16L, 16L, 16L, 16L, 16L, 16L, 16L, 16L,
                                                16L, 16L, 18L, 18L,
18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L,
                                                18L, 18L, 18L, 18L,
18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L,
                                                18L, 17L, 17L, 17L,
17L, 17L, 17L, 17L, 17L, 17L, 17L, 17L, 17L,
                                                17L, 17L, 17L, 17L,
17L, 17L, 17L, 17L, 17L, 17L, 17L, 17L, 17L,
                                                20L, 20L, 20L, 20L,
20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L,
                                                20L, 20L, 20L, 20L,
20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 21L,
                                                21L, 21L, 21L, 21L,
21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L,
                                                21L, 21L, 21L, 21L,
21L, 21L, 21L, 21L, 21L, 21L, 21L, 23L, 23L,
                                                23L, 23L, 23L, 23L,
23L, 23L, 23L, 23L, 23L, 23L, 23L, 23L, 23L,
                                                23L, 23L, 23L, 23L,
23L, 23L, 23L, 23L, 23L, 23L, 24L, 24L, 24L,
                                                24L, 24L, 24L, 24L,
24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L,
                                                24L, 24L, 24L, 24L,
24L, 24L, 24L, 24L, 24L, 25L, 25L, 25L, 25L,
                                                25L, 25L, 25L, 25L,
25L, 25L, 25L, 25L, 25L, 25L, 25L, 25L, 25L,
                                                25L, 25L, 25L, 25L,
25L, 25L, 25L, 25L, 22L, 22L, 22L, 22L, 22L,
                                                22L, 22L, 22L, 22L,
22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L,
                                                22L, 22L, 22L, 22L,
22L, 22L, 22L, 28L, 28L, 28L, 28L, 28L, 28L,
                                                28L, 28L, 28L, 28L,
28L, 28L, 28L, 28L, 28L, 28L, 28L, 28L, 28L,
                                                28L, 28L, 28L, 28L,
28L, 28L, 27L, 27L, 27L, 27L, 27L, 27L, 27L,
                                                27L, 27L, 27L, 27L,
27L, 27L, 27L, 27L, 27L, 27L, 27L, 27L, 27L,
                                                27L, 27L, 27L, 27L,
27L, 29L, 29L, 29L, 29L, 29L, 29L, 29L, 29L,
                                                29L, 29L, 29L, 29L,
29L, 29L, 29L, 29L, 29L, 29L, 29L, 29L, 29L,
                                                29L, 29L, 29L, 29L,
30L, 30L, 30L, 30L, 30L, 30L, 30L, 30L, 30L,
                                                30L, 30L, 30L, 30L,
30L, 30L, 30L, 30L, 30L, 30L, 30L, 30L, 30L,
                                                30L, 30L, 30L), .Label
= c("ago", "bdi", "ben", "bfa", "caf",

    "cod", "eri", "eth", "gin", "gmb", "gnb", "gnq", "lbr", "lso",

    "mdg", "mli", "moz", "mrt", "mwi", "ner", "rwa", "sdn", "sen",

    "sle", "som", "tcd", "tgo", "tza", "uga", "zmb"), class =
"factor"),
                          year = c(1991L, 1992L, 1993L, 1994L, 1995L,
1996L, 1997L,
                                   1998L, 1999L, 2000L, 2001L, 2002L,
2003L, 2004L, 2005L, 2006L,
                                   2007L, 2008L, 2009L, 2010L, 2011L,
2012L, 2013L, 2014L, 2015L,
                                   1991L, 1992L, 1993L, 1994L, 1995L,
1996L, 1997L, 1998L, 1999L,
                                   2000L, 2001L, 2002L, 2003L, 2004L,
2005L, 2006L, 2007L, 2008L,
                                   2009L, 2010L, 2011L, 2012L, 2013L,
2014L, 2015L, 1991L, 1992L,
                                   1993L, 1994L, 1995L, 1996L, 1997L,
1998L, 1999L, 2000L, 2001L,
                                   2002L, 2003L, 2004L, 2005L, 2006L,
2007L, 2008L, 2009L, 2010L,
                                   2011L, 2012L, 2013L, 2014L, 2015L,
1991L, 1992L, 1993L, 1994L,
                                   1995L, 1996L, 1997L, 1998L, 1999L,
2000L, 2001L, 2002L, 2003L,
                                   2004L, 2005L, 2006L, 2007L, 2008L,
2009L, 2010L, 2011L, 2012L,
                                   2013L, 2014L, 2015L, 1991L, 1992L,
1993L, 1994L, 1995L, 1996L,
                                   1997L, 1998L, 1999L, 2000L, 2001L,
2002L, 2003L, 2004L, 2005L,
                                   2006L, 2007L, 2008L, 2009L, 2010L,
2011L, 2012L, 2013L, 2014L,
                                   2015L, 1991L, 1992L, 1993L, 1994L,
1995L, 1996L, 1997L, 1998L,
                                   1999L, 2000L, 2001L, 2002L, 2003L,
2004L, 2005L, 2006L, 2007L,
                                   2008L, 2009L, 2010L, 2011L, 2012L,
2013L, 2014L, 2015L, 1991L,
                                   1992L, 1993L, 1994L, 1995L, 1996L,
1997L, 1998L, 1999L, 2000L,
                                   2001L, 2002L, 2003L, 2004L, 2005L,
2006L, 2007L, 2008L, 2009L,
                                   2010L, 2011L, 2012L, 2013L, 2014L,
2015L, 1991L, 1992L, 1993L,
                                   1994L, 1995L, 1996L, 1997L, 1998L,
1999L, 2000L, 2001L, 2002L,
                                   2003L, 2004L, 2005L, 2006L, 2007L,
2008L, 2009L, 2010L, 2011L,
                                   2012L, 2013L, 2014L, 2015L, 1991L,
1992L, 1993L, 1994L, 1995L,
                                   1996L, 1997L, 1998L, 1999L, 2000L,
2001L, 2002L, 2003L, 2004L,
                                   2005L, 2006L, 2007L, 2008L, 2009L,
2010L, 2011L, 2012L, 2013L,
                                   2014L, 2015L, 1991L, 1992L, 1993L,
1994L, 1995L, 1996L, 1997L,
                                   1998L, 1999L, 2000L, 2001L, 2002L,
2003L, 2004L, 2005L, 2006L,
                                   2007L, 2008L, 2009L, 2010L, 2011L,
2012L, 2013L, 2014L, 2015L,
                                   1991L, 1992L, 1993L, 1994L, 1995L,
1996L, 1997L, 1998L, 1999L,
                                   2000L, 2001L, 2002L, 2003L, 2004L,
2005L, 2006L, 2007L, 2008L,
                                   2009L, 2010L, 2011L, 2012L, 2013L,
2014L, 2015L, 1991L, 1992L,
                                   1993L, 1994L, 1995L, 1996L, 1997L,
1998L, 1999L, 2000L, 2001L,
                                   2002L, 2003L, 2004L, 2005L, 2006L,
2007L, 2008L, 2009L, 2010L,
                                   2011L, 2012L, 2013L, 2014L, 2015L,
1991L, 1992L, 1993L, 1994L,
                                   1995L, 1996L, 1997L, 1998L, 1999L,
2000L, 2001L, 2002L, 2003L,
                                   2004L, 2005L, 2006L, 2007L, 2008L,
2009L, 2010L, 2011L, 2012L,
                                   2013L, 2014L, 2015L, 1991L, 1992L,
1993L, 1994L, 1995L, 1996L,
                                   1997L, 1998L, 1999L, 2000L, 2001L,
2002L, 2003L, 2004L, 2005L,
                                   2006L, 2007L, 2008L, 2009L, 2010L,
2011L, 2012L, 2013L, 2014L,
                                   2015L, 1991L, 1992L, 1993L, 1994L,
1995L, 1996L, 1997L, 1998L,
                                   1999L, 2000L, 2001L, 2002L, 2003L,
2004L, 2005L, 2006L, 2007L,
                                   2008L, 2009L, 2010L, 2011L, 2012L,
2013L, 2014L, 2015L, 1991L,
                                   1992L, 1993L, 1994L, 1995L, 1996L,
1997L, 1998L, 1999L, 2000L,
                                   2001L, 2002L, 2003L, 2004L, 2005L,
2006L, 2007L, 2008L, 2009L,
                                   2010L, 2011L, 2012L, 2013L, 2014L,
2015L, 1991L, 1992L, 1993L,
                                   1994L, 1995L, 1996L, 1997L, 1998L,
1999L, 2000L, 2001L, 2002L,
                                   2003L, 2004L, 2005L, 2006L, 2007L,
2008L, 2009L, 2010L, 2011L,
                                   2012L, 2013L, 2014L, 2015L, 1991L,
1992L, 1993L, 1994L, 1995L,
                                   1996L, 1997L, 1998L, 1999L, 2000L,
2001L, 2002L, 2003L, 2004L,
                                   2005L, 2006L, 2007L, 2008L, 2009L,
2010L, 2011L, 2012L, 2013L,
                                   2014L, 2015L, 1991L, 1992L, 1993L,
1994L, 1995L, 1996L, 1997L,
                                   1998L, 1999L, 2000L, 2001L, 2002L,
2003L, 2004L, 2005L, 2006L,
                                   2007L, 2008L, 2009L, 2010L, 2011L,
2012L, 2013L, 2014L, 2015L,
                                   1991L, 1992L, 1993L, 1994L, 1995L,
1996L, 1997L, 1998L, 1999L,
                                   2000L, 2001L, 2002L, 2003L, 2004L,
2005L, 2006L, 2007L, 2008L,
                                   2009L, 2010L, 2011L, 2012L, 2013L,
2014L, 2015L, 1991L, 1992L,
                                   1993L, 1994L, 1995L, 1996L, 1997L,
1998L, 1999L, 2000L, 2001L,
                                   2002L, 2003L, 2004L, 2005L, 2006L,
2007L, 2008L, 2009L, 2010L,
                                   2011L, 2012L, 2013L, 2014L, 2015L,
1991L, 1992L, 1993L, 1994L,
                                   1995L, 1996L, 1997L, 1998L, 1999L,
2000L, 2001L, 2002L, 2003L,
                                   2004L, 2005L, 2006L, 2007L, 2008L,
2009L, 2010L, 2011L, 2012L,
                                   2013L, 2014L, 2015L, 1991L, 1992L,
1993L, 1994L, 1995L, 1996L,
                                   1997L, 1998L, 1999L, 2000L, 2001L,
2002L, 2003L, 2004L, 2005L,
                                   2006L, 2007L, 2008L, 2009L, 2010L,
2011L, 2012L, 2013L, 2014L,
                                   2015L, 1991L, 1992L, 1993L, 1994L,
1995L, 1996L, 1997L, 1998L,
                                   1999L, 2000L, 2001L, 2002L, 2003L,
2004L, 2005L, 2006L, 2007L,
                                   2008L, 2009L, 2010L, 2011L, 2012L,
2013L, 2014L, 2015L, 1991L,
                                   1992L, 1993L, 1994L, 1995L, 1996L,
1997L, 1998L, 1999L, 2000L,
                                   2001L, 2002L, 2003L, 2004L, 2005L,
2006L, 2007L, 2008L, 2009L,
                                   2010L, 2011L, 2012L, 2013L, 2014L,
2015L, 1991L, 1992L, 1993L,
                                   1994L, 1995L, 1996L, 1997L, 1998L,
1999L, 2000L, 2001L, 2002L,
                                   2003L, 2004L, 2005L, 2006L, 2007L,
2008L, 2009L, 2010L, 2011L,
                                   2012L, 2013L, 2014L, 2015L, 1991L,
1992L, 1993L, 1994L, 1995L,
                                   1996L, 1997L, 1998L, 1999L, 2000L,
2001L, 2002L, 2003L, 2004L,
                                   2005L, 2006L, 2007L, 2008L, 2009L,
2010L, 2011L, 2012L, 2013L,
                                   2014L, 2015L, 1991L, 1992L, 1993L,
1994L, 1995L, 1996L, 1997L,
                                   1998L, 1999L, 2000L, 2001L, 2002L,
2003L, 2004L, 2005L, 2006L,
                                   2007L, 2008L, 2009L, 2010L, 2011L,
2012L, 2013L, 2014L, 2015L,
                                   1991L, 1992L, 1993L, 1994L, 1995L,
1996L, 1997L, 1998L, 1999L,
                                   2000L, 2001L, 2002L, 2003L, 2004L,
2005L, 2006L, 2007L, 2008L,
                                   2009L, 2010L, 2011L, 2012L, 2013L,
2014L, 2015L, 1991L, 1992L,
                                   1993L, 1994L, 1995L, 1996L, 1997L,
1998L, 1999L, 2000L, 2001L,
                                   2002L, 2003L, 2004L, 2005L, 2006L,
2007L, 2008L, 2009L, 2010L,
                                   2011L, 2012L, 2013L, 2014L, 2015L),
ag_emp = c(NA, 0.002,

           0.019, 0.002, -0.006, -0.014, -0.012, -0.01, -0.006,
-0.005,

           -0.008, -0.013, -0.01, -0.011, -0.017, -0.021, -0.022,
-0.017,

           -0.007, -0.004, -0.007, -0.007, -0.007, -0.005, -0.003, NA,

           0.002, 0.002, 0.002, 0.002, -0.002, -0.001, -0.001, -0.001,

           -0.001, 0, 0.001, 0.006, 0.004, 0.004, 0.003, 0.003, 0.003,

           0.003, 0.004, 0.007, 0.004, 0.004, 0.004, -0.002, NA,
0.001,

           -0.008, 0.003, 0, -0.005, -0.002, -0.004, -0.003, -0.005,

           -0.001, -0.01, -0.003, 0, 0, -0.038, -0.138, -0.055,
-0.024,

           -0.018, -0.016, -0.014, -0.012, -0.009, -0.007, NA, -0.014,

           -0.01, -0.008, -0.007, -0.008, -0.008, -0.008, -0.001, 0,

           0.009, 0.009, 0.009, 0.008, 0.007, 0.004, 0.004, 0.004,
0.004,

           0.004, -0.002, -0.001, -0.001, -0.002, 0, NA, 0.006, 0.003,

           0.002, -0.003, 0.001, 0.001, -0.003, -0.001, 0.003, -0.003,

           -0.002, 0.004, 0, 0.001, 0, -0.001, 0.001, 0.001, 0.001,

           0.001, 0, 0.027, 0.011, 0.001, NA, -0.01, 0.009, -0.003,

           -0.003, -0.004, -0.003, -0.003, -0.002, 0, -0.003, -0.004,

           -0.004, -0.006, -0.01, -0.002, 0.001, 0.001, 0.002, 0.001,

           -0.001, 0.002, 0.001, 0.002, 0.002, NA, 0.006, 0.008,
0.002,

           -0.001, -0.002, 0.003, -0.003, 0.001, 0.005, -0.002,
-0.005,

           -0.006, -0.005, -0.006, -0.007, -0.006, -0.006, -0.004,
-0.005,

           -0.007, -0.006, -0.008, -0.007, -0.002, NA, -0.01, -0.015,

           -0.014, -0.007, -0.004, -0.006, -0.014, 0, 0.001, -0.011,

           -0.005, 0, -0.003, -0.005, 0.003, -0.002, -0.002, 0.004,

           0.002, -0.001, 0, 0.001, -0.002, 0.038, NA, -0.009, -0.015,

           -0.022, -0.009, -0.003, 0.006, 0.008, 0.009, 0.012, 0.004,

           0.017, 0.012, 0.008, 0.009, 0.001, 0.001, 0.004, 0.001,
-0.003,

           -0.005, -0.006, -0.004, -0.005, -0.004, NA, -0.002, 0.001,

           -0.001, -0.001, -0.004, -0.009, -0.008, -0.01, 0.034,
-0.012,

           0.01, 0.01, 0.009, -0.001, -0.006, -0.006, -0.005, -0.004,

           -0.005, -0.005, -0.008, 0.002, -0.007, 0, NA, -0.008,
-0.013,

           -0.016, -0.009, -0.015, -0.015, -0.011, -0.014, -0.014,
-0.016,

           -0.018, -0.012, -0.016, -0.019, -0.016, -0.016, -0.016,
-0.016,

           -0.016, -0.019, -0.014, -0.016, -0.017, -0.003, NA, 0.001,

           0.002, 0.002, 0.005, 0.001, 0.002, 0.003, 0.004, 0.006,
0.001,

           0.003, 0.005, 0.004, 0.004, 0.004, 0.004, 0.003, 0.006,
0.007,

           0.016, 0.002, 0.087, 0.002, 0, NA, -0.009, -0.007, -0.006,

           -0.006, -0.013, -0.008, 0.019, 0.006, -0.01, -0.004,
-0.002,

           -0.002, -0.004, -0.005, -0.006, -0.006, -0.005, -0.005,
-0.006,

           -0.01, -0.005, -0.003, -0.005, -0.004, NA, -0.006, 0.011,

           -0.006, -0.083, 0.084, -0.125, 0.066, 0.064, -0.005,
-0.115,

           0.101, 0.006, -0.064, -0.054, 0.098, -0.013, -0.007, 0.005,

           -0.007, 0.007, -0.004, 0.006, -0.016, -0.008, NA, -0.027,

           0.055, 0.044, 0.029, 0.023, -0.045, -0.014, -0.015, -0.012,

           -0.003, -0.003, 0.02, 0.003, -0.009, -0.001, -0.013,
-0.003,

           0, 0.001, -0.006, -0.006, -0.006, -0.002, 0, NA, -0.002,

           0.012, -0.01, 0.018, 0.006, 0.004, 0.002, 0, 0.009, -0.005,

           0.006, 0.017, 0.036, 0.029, -0.006, -0.001, -0.002, 0.001,

           -0.002, -0.007, -0.008, -0.008, -0.006, -0.003, NA, -0.004,

           -0.003, -0.002, -0.128, -0.009, -0.006, -0.007, 0.005,
0.005,

           0.008, 0.007, 0.006, -0.001, 0.006, 0.007, -0.005, 0.007,

           -0.004, -0.003, -0.002, 0, -0.004, -0.001, -0.001, NA,
-0.012,

           -0.003, 0.007, 0.004, 0.001, -0.02, -0.005, 0.013, -0.012,

           -0.017, -0.016, 0.026, -0.011, 0.02, 0.013, 0.011, 0.028,

           0.026, 0.031, -0.005, -0.004, -0.007, -0.013, -0.007, NA,

           -0.001, -0.006, -0.001, -0.008, -0.007, 0.001, -0.003,
-0.005,

           -0.001, -0.002, -0.004, -0.007, -0.007, -0.004, -0.011,
-0.003,

           -0.002, -0.002, -0.002, -0.005, -0.005, -0.005, -0.006,
-0.003,

           NA, 0.005, 0.016, 0.004, 0.004, 0, -0.007, 0.003, 0.002,

           0, 0.003, -0.005, -0.004, -0.007, -0.008, -0.009, -0.011,

           -0.011, -0.01, -0.012, -0.014, -0.016, -0.003, 0, 0.001,

           NA, -0.009, 0.021, 0.006, -0.013, 0.02, 0.002, -0.003,
0.007,

           -0.001, 0.001, 0.003, -0.001, 0, 0.001, -0.004, -0.001,
-0.005,

           -0.002, -0.004, -0.003, -0.006, -0.004, -0.006, -0.005, NA,

           -0.006, 0.029, 0.029, 0.02, 0.001, -0.014, -0.033, -0.002,

           -0.016, -0.015, 0.014, -0.015, -0.034, -0.023, -0.001,
-0.003,

           0.013, -0.015, -0.003, -0.008, -0.008, -0.002, 0.007, 0,

           NA, 0, -0.015, -0.015, -0.015, -0.024, -0.015, -0.009,
-0.016,

           -0.026, -0.019, -0.007, -0.037, -0.027, -0.026, -0.033,
0.042,

           0.034, 0.033, 0.041, 0.028, 0.016, 0.006, 0.009, 0.001, NA,

           0.009, 0, -0.002, 0.003, -0.005, 0, 0.001, -0.001, 0.003,

           0.006, -0.018, -0.008, 0.013, 0.001, -0.002, -0.003,
-0.003,

           0, 0, -0.003, -0.006, -0.009, -0.005, 0.016, NA, -0.007,

           -0.005, -0.004, -0.006, -0.01, -0.01, -0.01, -0.008,
-0.006,

           -0.009, -0.006, -0.003, -0.005, -0.004, -0.006, -0.006,
-0.005,

           -0.003, -0.003, -0.005, -0.001, -0.001, -0.001, -0.001, NA,

           -0.019, -0.14, 0.006, 0.009, 0.007, 0.044, -0.035, -0.004,

           -0.002, -0.017, -0.005, -0.001, -0.004, -0.008, 0.004,
-0.003,

           -0.011, 0.006, -0.001, -0.004, -0.002, 0.007, -0.003,
-0.001,

           NA, 0.001, -0.002, -0.001, -0.005, -0.007, -0.002, 0.004,

           0.001, 0.002, -0.012, 0.004, -0.009, -0.014, -0.017,
-0.029,

           -0.044, -0.014, -0.013, -0.017, -0.02, -0.017, -0.025,
-0.006,

           -0.008, NA, 0.004, 0.009, 0.001, -0.001, -0.008, -0.004,

           0.003, 0.003, 0.003, 0.003, 0.003, 0.001, 0, 0.001, -0.002,

           -0.002, -0.002, -0.002, -0.003, -0.004, -0.004, -0.004,
-0.004,

           -0.003, NA, -0.006, -0.005, 0.002, -0.001, -0.004, 0,
-0.012,

           -0.016, 0.008, -0.017, -0.011, -0.014, 0.015, 0.008, 0.003,

           0.009, 0.014, 0.014, 0.012, 0.01, 0.009, 0.009, 0.006,
0.001,

           NA, -0.001, 0, 0.015, 0.014, 0.025, 0.022, 0.027, 0.007,

           0.007, -0.006, 0.013, 0.016, -0.01, -0.031, -0.003, 0.001,

           0.016, -0.008, -0.043, 0.024, 0.025, 0, -0.015, -0.008),
                          in_emp = c(NA, -0.001, -0.004, 0.005, 0.004,
0.001, 0, 0.001,
                                     0.001, 0.001, 0.001, 0.003, 0,
0.003, 0.004, 0.003, 0.003,
                                     0.001, -0.001, 0.002, 0.001,
0.002, 0.002, 0.002, 0.002,
                                     NA, 0, 0.001, 0.001, 0.001, 0, 0,
0, 0, 0, 0, 0, 0.001, 0.002,
                                     0.002, 0.002, 0.002, 0.002,
0.002, 0.002, 0.002, 0.001, 0.001,
                                     0.001, 0.002, NA, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, -0.002,
                                     -0.002, 0.024, 0.046, 0.021,
0.01, 0.008, 0.007, 0.006, 0.006,
                                     0.006, 0.004, NA, -0.001, -0.001,
0, -0.001, 0, 0, 0, 0,
                                     0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, NA, -0.001,
                                     0, 0.002, 0, -0.002, 0.003, 0, 0,
-0.001, 0, 0, -0.001, 0.002,
                                     0, 0.001, 0.001, 0, 0.001, 0.001,
0.001, 0.001, -0.008, 0.007,
                                     0.001, NA, 0, -0.001, 0, 0, 0, 0,
0, 0, 0, 0.001, 0, 0.001,
                                     0.002, -0.001, 0, 0, 0, 0, 0.001,
-0.001, 0.001, 0, 0, 0.001,
                                     NA, 0, -0.001, 0.001, 0.001,
0.001, 0.002, -0.002, 0, 0.002,
                                     0, 0, 0.001, 0.002, 0, 0, 0.001,
0.001, 0, 0.002, 0.001,
                                     0.001, 0, 0.002, 0.002, NA,
0.001, -0.001, 0.002, 0, 0, 0.002,
                                     -0.003, -0.001, 0.001, 0.002, 0,
0.001, 0.004, 0, 0, 0.002,
                                     0.001, -0.001, 0.001, 0, 0.001,
-0.001, 0, -0.028, NA, -0.001,
                                     0.002, 0.001, 0.001, -0.001,
0.002, 0.001, 0, 0, 0, 0.004,
                                     0, 0, 0.002, 0, 0.001, -0.002,
0.002, 0, 0.002, 0.001, 0,
                                     0, 0.001, NA, 0, 0, 0, 0, 0,
-0.001, -0.001, -0.001, 0.004,
                                     -0.001, 0.001, 0.001, 0.001,
0.001, 0.001, 0.001, 0.001,
                                     0.001, 0.001, 0.001, 0.001,
0.002, 0.001, 0.003, NA, 0.001,
                                     0.003, 0.002, 0.003, 0.003,
0.003, 0.003, 0.003, 0.003, 0.003,
                                     0.003, 0.004, 0.004, 0.003,
0.004, 0.004, 0.004, 0.004, 0.004,
                                     0.003, 0.005, 0.004, 0.004,
0.001, NA, 0, 0, 0, 0, 0, 0,
                                     0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, -0.002, 0, 0.004,
                                     0, 0.001, NA, 0, 0.001, 0.001,
0.002, 0.002, 0.002, -0.005,
                                     0.008, -0.001, 0.002, 0.001,
0.002, 0.002, 0.002, 0.001,
                                     0.002, 0.002, 0.002, 0.002,
0.003, -0.001, 0.002, 0.002,
                                     0.003, NA, 0.001, 0.003, 0.001,
-0.009, 0.014, -0.016, 0.01,
                                     0.01, 0.002, -0.015, 0.015,
0.004, -0.009, -0.006, 0.018,
                                     0, 0.001, 0.002, 0.001, 0.003,
0.001, 0.003, -0.001, 0.001,
                                     NA, -0.017, -0.008, -0.007,
-0.006, -0.004, 0.006, 0.002,
                                     0.001, 0.003, 0, 0.001, -0.007,
0.006, 0, 0.003, 0, 0.001,
                                     0.002, 0.003, 0.001, 0.001,
0.001, -0.001, 0.001, NA, -0.004,
                                     -0.001, -0.006, -0.001, -0.003,
-0.004, -0.004, -0.005, -0.005,
                                     -0.007, -0.008, -0.021, -0.02,
-0.012, 0.004, 0.005, 0.005,
                                     0.007, 0.008, 0.009, 0.01, 0.008,
0.003, 0.002, NA, -0.002,
                                     0.003, -0.002, 0.013, 0, 0, 0,
0.001, 0, 0, 0.002, 0.002,
                                     0, 0.001, 0, 0, 0.001, 0.001, 0,
0, 0, 0.001, 0.001, 0.001,
                                     NA, -0.001, 0, 0.001, 0.001,
0.001, 0, 0.001, 0.002, 0, 0.002,
                                     -0.001, 0.004, 0, 0.003, 0.003,
0.002, 0.004, 0.004, 0.004,
                                     0.001, 0, 0.001, 0.001, 0.002,
NA, 0.003, 0.003, -0.001,
                                     0.005, 0.001, -0.001, 0.003,
0.004, 0.001, 0.003, 0.001,
                                     0.003, 0.002, 0.005, 0.006, 0,
0.002, 0.001, 0.004, 0.002,
                                     0.003, 0.002, 0.002, 0.002, NA,
-0.001, 0.002, 0, 0, 0.001,
                                     0, 0.001, 0.001, 0, 0.001, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0,
                                     -0.001, 0.001, 0.001, 0.001, NA,
-0.003, 0.006, 0.003, -0.003,
                                     0.007, 0.001, 0.004, 0, 0, 0.004,
0, 0, -0.002, 0.002, 0,
                                     0, 0.001, -0.002, 0.002, -0.001,
0.002, -0.001, 0.001, 0.001,
                                     NA, 0, 0.001, 0.001, 0.001, 0, 0,
0, 0.001, 0, 0.001, 0.003,
                                     0.004, 0.003, 0.004, 0.004,
0.004, 0.005, 0.004, 0.005, 0.005,
                                     0.005, 0.006, 0.003, 0.004, NA,
0, -0.004, -0.004, -0.003,
                                     -0.005, -0.003, -0.001, -0.002,
-0.004, 0.003, 0.009, 0.001,
                                     0.004, 0.005, 0.002, 0.012,
0.009, 0.01, 0.013, 0.01, 0.008,
                                     0.003, 0.004, 0.003, NA, -0.003,
0.002, 0, -0.001, 0, -0.001,
                                     0.002, -0.001, 0.004, -0.001,
0.001, -0.001, 0.002, 0, 0,
                                     0.001, 0.001, 0, 0.001, 0.001,
0.002, 0.003, 0, -0.004, NA,
                                     0, 0.001, -0.001, 0.002, 0, 0, 0,
0, 0, 0, 0, 0, 0.001, 0.001,
                                     0.001, 0.001, 0, 0.001, 0.002,
0.001, 0.001, 0.002, 0.002,
                                     0.002, NA, -0.001, 0.066, 0.001,
0.004, 0.003, 0.027, -0.018,
                                     -0.001, 0.001, -0.006, 0, 0.003,
0, -0.001, 0.006, 0.003,
                                     -0.003, 0.004, 0, 0.001, 0.003,
0.004, 0, 0.001, NA, -0.002,
                                     -0.002, -0.002, -0.002, -0.002,
-0.001, -0.001, -0.001, -0.001,
                                     -0.002, 0.014, 0.003, 0.003,
0.003, 0.002, -0.007, 0.002,
                                     0.002, 0.002, 0.002, 0.002,
0.002, 0.003, 0.003, NA, 0, -0.002,
                                     0.005, 0.002, -0.002, 0.002, 0,
0.002, 0, 0.001, 0.001, 0.002,
                                     0.001, 0.001, 0.001, 0.001,
0.001, 0.001, 0.001, 0.001, 0.001,
                                     0.001, 0.001, 0.002, NA, 0,
0.001, 0.001, 0.001, 0.001, 0.002,
                                     0.002, 0.006, 0.007, 0.006,
0.008, 0.008, -0.019, -0.015,
                                     0.005, 0.006, 0.007, 0.008,
-0.006, -0.006, -0.005, -0.005,
                                     -0.003, 0.002, NA, -0.002,
-0.002, 0, -0.001, -0.001, -0.004,
                                     -0.003, -0.005, -0.004, 0.002,
0.004, 0.004, 0.002, 0, 0.001,
                                     0.002, 0.003, 0.002, -0.003,
0.004, 0.004, 0.001, 0,  prd = c(NA, -0.095, -0.276, -0.032, 0.06,

                            0.083, 0.053, 0.028, -0.006, -0.003,
-0.046, 0.099, -0.031,

                            0.105, 0.143, 0.123, 0.127, 0.048, 0.075,
0.114, 0.101, 0.164,

                            0.081, 0.012, -0.033, NA, -0.017, -0.019,
-0.019, 0.033,

                            0.007, 0.023, -0.005, -0.006, 0.039,
0.029, -0.011, -0.01,

                            0.008, -0.016, 0.008, 0.009, 0.008,
-0.002, -0.018, -0.004,

                            -0.013, 0.038, 0.031, 0.01, NA, -0.023,
0.023, -0.059, 0.024,

                            0.066, 0.032, 0.051, 0.042, 0.018, 0.031,
0.023, 0.097, 0.007,

                            0.033, 0.089, -0.002, 0.038, -0.039,
0.051, 0.021, 0.031,

                            0.036, 0.023, 0.018, NA, 0.025, -0.074,
-0.098, -0.088, -0.041,

                            0.014, 0.027, -0.013, -0.044, -0.018, 0,
-0.037, -0.023,

                            -0.01, -0.003, -0.013, 0.005, -0.009,
0.004, 0.006, 0.011,

                            0.007, 0.015, -0.07, NA, -0.05, -0.037,
0.015, 0.013, -0.091,

                            0.007, 0.008, 0.076, 0.005, -0.018, -0.03,
-0.097, -0.003,

                            0.005, 0.02, 0.016, 0.001, -0.006, -0.002,
0.008, 0.017,

                            -0.36, -0.012, 0.007, NA, 0.007, -0.088,
0.051, -0.024, -0.011,

                            0.029, 0.038, -0.043, -0.04, 0.08, 0.042,
0.191, 0.264, 0.112,

                            -0.046, 0.133, -0.018, 0.058, 0.076,
-0.069, 0.071, 0.038,

                            0.194, 0.026, NA, -0.132, -0.182, -0.075,
-0.026, -0.017,

                            -0.1, -0.02, -0.05, -0.101, -0.062, 0.001,
0.024, 0.033,

                            0.029, 0.019, 0.041, 0.032, -0.003, 0.037,
0.036, 0.036,

                            0.049, 0.062, 0.034, NA, 0.149, 0.064,
-0.018, 0.062, 0.307,

                            0.605, 0.18, 0.187, 0.119, 0.575, 0.171,
0.093, 0.165, 0.061,

                            0.006, 0.111, 0.15, -0.014, -0.12, 0.043,
0.051, -0.073,

                            -0.031, -0.093, NA, 0.137, 0.11, 0.19,
0.213, 0.009, 0.053,

                            0.063, -0.017, -0.028, -0.175, -0.004,
0.01, -0.071, -0.015,

                            -0.029, 0.003, -0.108, 0.022, -0.006,
0.067, 0.05, 0.008,

                            0.031, 0.025, NA, -0.12, 0.091, -0.009,
0.032, 0.087, 0.008,

                            -0.062, 0.033, -0.012, 0.056, -0.03,
-0.06, 0.067, 0.085,

                            0.084, 0.083, 0.081, 0.065, 0.071, 0.104,
0.061, 0.058, 0.074,

                            0.071, NA, 0.002, 0.029, -0.025, -0.004,
0.016, -0.02, 0,

                            0.037, 0.021, 0.026, -0.018, -0.015,
-0.037, -0.032, -0.018,

                            -0.001, 0.023, 0.038, 0.019, -0.059,
0.017, 0.019, -0.002,

                            0.011, NA, -0.029, -0.014, -0.018, -0.01,
0.011, 0.014, 0.025,

                            0.021, -0.002, 0.016, 0.014, -0.016,
-0.006, -0.007, -0.015,

                            0.008, 0.001, -0.039, 0.014, 0.019, 0.009,
-0.1, -0.023,

                            0.002, NA, 0.003, -0.004, 0.01, 0.007,
0.022, 0.034, -0.189,

                            0.089, -0.026, -0.005, -0.048, -0.008,
-0.01, 0.021, 0.005,

                            0.004, 0.011, 0.008, 0.018, 0.062, -0.039,
0.001, -0.026,

                            -0.017, NA, 0.058, -0.003, 0.034, 0.128,
-0.075, 0.245, -0.073,

                            -0.066, 0.041, 0.249, -0.12, 0.022, 0.117,
0.089, -0.103,

                            0.051, 0.047, 0.009, 0.072, 0.002, 0.034,
0.011, 0.029, 0.033,

                            NA, -0.339, -0.325, -0.224, -0.073, 0.047,
0.903, 0.19, 0.155,

                            0.176, 0.239, 0.287, -0.338, -0.083,
0.063, 0.041, 0.074,

                            0.122, 0.046, 0.062, 0.022, 0.061, 0.051,
-0.011, -0.049,

                            NA, -0.013, -0.017, -0.007, -0.032,
-0.014, 0.005, 0.013,

                            0.023, 0.008, 0.043, -0.145, 0.057, 0.005,
0, 0.021, 0.027,

                            0.022, -0.154, 0.055, -0.027, 0.001,
-0.005, 0.009, 0.001,

                            NA, -0.092, 0.122, -0.134, 0.095, 0.088,
0.039, 0.015, -0.017,

                            -0.023, -0.077, 0.001, 0.017, 0.057,
-0.03, -0.026, 0.074,

                            0.029, 0.05, 0.046, 0.011, -0.013, 0.023,
0.026, -0.01, NA,

                            0.063, 0.01, 0, -0.028, 0.019, 0.022,
0.044, -0.003, 0.027,

                            0.134, 0.071, 0.036, -0.098, 0.03, 0.066,
0.025, 0.029, -0.005,

                            0.038, 0.033, 0.109, 0.035, 0.066, 0.041,
NA, -0.028, 0.092,

                            -0.005, 0.029, -0.028, -0.098, -0.012,
0.054, -0.049, -0.013,

                            -0.031, 0.022, 0.014, 0.054, 0.173,
-0.024, -0.032, 0.023,

                            -0.04, 0.006, 0.016, 0.04, 0.034, -0.016,
NA, -0.073, 0.029,

                            0.024, -0.03, 0.056, 0.09, 0.055, 0.047,
-0.018, 0.097, 0.049,

                            0.038, 0.066, 0.056, 0.071, 0.058, 0.061,
0.048, 0.061, 0.06,

                            0.064, 0.035, 0.039, 0.042, NA, -0.006,
-0.057, -0.026, 0.006,

                            -0.03, -0.034, 0.066, -0.038, -0.067,
0.024, 0.002, 0.001,

                            -0.053, 0.036, 0.026, -0.006, 0.059,
-0.048, 0.049, -0.021,

                            0.102, 0.003, 0.033, -0.017, NA, 0.119,
-0.115, -0.469, 0.33,

                            0.067, 0.096, 0.061, 0.005, 0.048, 0.062,
0.066, -0.007,

                            0.061, 0.069, 0.051, 0.055, 0.053, 0.039,
0.034, 0.051, 0.073,

                            0.022, 0.035, 0.032, NA, -0.02, -0.042,
0.011, 0.047, 0.053,

                            0.021, 0.03, 0.055, 0.044, 0.014, -0.047,
0.091, 0.055, 0.049,

                            0.032, -0.014, 0.002, -0.007, -0.015,
-0.027, 0.016, 0.007,

                            0.001, 0.025, NA, -0.092, 0, 0.061,
-0.113, -0.227, -0.154,

                            -0.039, -0.091, -0.012, 0.11, 0.213,
0.044, -0.001, -0.002,

                            0.006, 0.047, 0.024, 0.005, 0.025, 0.036,
0.135, 0.194, 0.023,

                            -0.206, NA, -0.099, 0.007, -0.186, -0.005,
0.026, -0.047,

                            -0.006, 0.006, -0.004, 0.017, 0.008,
0.003, 0.005, 0.004,

                            0.007, -0.001, 0.004, 0.002, -0.002,
0.002, -0.026, -0.004,

                            0.014, -0.008, NA, 0.005, -0.019, -0.002,
0.007, 0.069, 0.059,

                            0.067, 0.037, 0.068, 0.065, 0.034, 0.043,
0.026, 0.045, 0.056,

                            0.07, -0.179, -0.009, 0.033, -0.006,
-0.061, 0.024, -0.004,

                            0.019, NA, 0.018, 0.008, -0.017, -0.008,
0.011, 0.012, 0.02,

                            0.018, 0.018, 0.046, 0.025, 0.033, 0.049,
0.045, 0.038, 0.023,

                            0.042, 0.031, 0.055, 0.073, 0.042, 0.068,
0.04, 0.038, NA,

                            -0.065, -0.163, 0.098, 0.04, 0.061, 0.101,
-0.052, -0.008,

                            -0.05, -0.061, -0.065, -0.006, -0.005,
-0.026, 0.04, -0.02,

                            -0.011, 0.012, 0.001, 0.052, -0.033,
0.028, 0.024, 0.02,

                            NA, 0.023, 0.035, 0.061, 0.051, 0.033,
0.023, 0.069, 0.055,

                            0.002, 0.076, 0.035, 0.024, 0.02, 0.05,
0.026, 0.018, 0.035,

                            0.019, 0.042, 0.011, -0.01, 0.006, 0.004,
0.014, NA, -0.082,

                            0.087, -0.157, 0, 0.004, -0.017, -0.044,
0.025, 0.011, 0.034,

                            0.006, 0.021, 0.059, 0.09, 0.055, 0.05,
0.024, 0.067, 0.113,

                            -0.006, 0.009, 0.012, 0.034, 0.003), gnipc
= structure(c(15L,


              487L, 397L, 323L, 233L, 315L, 370L, 309L, 268L, 391L,
381L,


              456L, 488L, 14L, 46L, 122L, 124L, 197L, 125L, 196L,
272L,


              273L, 274L, 335L, 199L, 264L, 286L, 285L, 183L, 265L,
284L,


              263L, 282L, 277L, 254L, 257L, 288L, 343L, 379L, 386L,
402L,


              441L, 472L, 461L, 453L, 473L, 474L, 496L, 502L, 465L,
240L,


              250L, 230L, 116L, 151L, 164L, 146L, 167L, 172L, 143L,
149L,


              168L, 224L, 252L, 279L, 293L, 321L, 372L, 360L, 367L,
411L,


              418L, 429L, 425L, 368L, 132L, 112L, 96L, 103L, 95L, 80L,


              90L, 81L, 68L, 62L, 56L, 500L, 478L, 500L, 83L, 94L,
97L,


              118L, 130L, 145L, 157L, 161L, 165L, 179L, 180L, 318L,
318L,


              287L, 173L, 233L, 200L, 179L, 181L, 184L, 162L, 163L,
171L,


              195L, 220L, 232L, 248L, 281L, 316L, 314L, 314L, 342L,
339L,


              239L, 266L, 246L, 194L, 195L, 157L, 119L, 148L, 163L,
146L,


              167L, 137L, 109L, 142L, 156L, 181L, 302L, 373L, 450L,
457L,


              497L, 467L, 486L, 499L, 9L, 1L, 4L, 470L, 247L, 216L,
284L,


              127L, 108L, 107L, 129L, 134L, 105L, 99L, 95L, 97L, 97L,
111L,


              138L, 160L, 180L, 211L, 187L, 216L, 253L, 283L, 291L,
323L,


              327L, 295L, 296L, 292L, 210L, 270L, 374L, 3L, 484L, 42L,


              47L, 51L, 121L, 123L, 198L, 337L, 336L, 385L, 55L, 60L,
54L,


              71L, 78L, 70L, 64L, 436L, 67L, 79L, 93L, 129L, 133L,
152L,


              158L, 163L, 154L, 131L, 136L, 126L, 145L, 185L, 179L,
194L,


              214L, 219L, 296L, 322L, 378L, 422L, 452L, 480L, 508L,
166L,


              119L, 76L, 68L, 72L, 77L, 79L, 67L, 62L, 65L, 63L, 59L,
63L,


              74L, 92L, 112L, 149L, 206L, 225L, 195L, 227L, 308L,
329L,


              363L, 405L, 444L, 446L, 448L, 440L, 447L, 455L, 432L,
439L,


              417L, 393L, 344L, 289L, 234L, 267L, 283L, 287L, 330L,
383L,


              358L, 365L, 340L, 333L, 320L, 297L, 313L, 394L, 406L,
388L,


              387L, 416L, 418L, 403L, 365L, 356L, 319L, 304L, 306L,
346L,


              354L, 294L, 295L, 400L, 419L, 401L, 404L, 398L, 410L,
425L,


              437L, 424L, 382L, 391L, 397L, 407L, 415L, 298L, 225L,
191L,


              209L, 184L, 188L, 207L, 241L, 264L, 290L, 288L, 323L,
377L,


              357L, 359L, 430L, 389L, 401L, 384L, 373L, 459L, 491L,
482L,


              477L, 494L, 479L, 489L, 507L, 498L, 471L, 438L, 413L,
506L,


              16L, 18L, 17L, 26L, 25L, 21L, 39L, 45L, 38L, 33L, 30L,
22L,


              146L, 86L, 61L, 495L, 493L, 56L, 91L, 102L, 98L, 104L,
84L,


              90L, 58L, 73L, 88L, 104L, 136L, 156L, 172L, 184L, 237L,
241L,


              261L, 266L, 341L, 161L, 182L, 201L, 165L, 170L, 221L,
186L,


              191L, 186L, 186L, 220L, 203L, 251L, 186L, 217L, 231L,
306L,


              367L, 320L, 326L, 367L, 350L, 361L, 329L, 311L, 271L,
230L,


              251L, 144L, 148L, 128L, 195L, 193L, 187L, 179L, 169L,
91L,


              79L, 82L, 89L, 191L, 213L, 247L, 285L, 310L, 331L, 230L,


              200L, 232L, 245L, 234L, 236L, 229L, 150L, 189L, 212L,
183L,


              208L, 216L, 174L, 191L, 212L, 258L, 295L, 318L, 338L,
379L,


              424L, 428L, 428L, 469L, 452L, 460L, 475L, 447L, 470L,
468L,


              396L, 423L, 431L, 427L, 424L, 375L, 358L, 324L, 309L,
341L,


              362L, 421L, 445L, 501L, 512L, 19L, 6L, 20L, 32L, 31L,
36L,


              27L, 12L, 146L, 89L, 91L, 89L, 89L, 112L, 141L, 144L,
160L,


              141L, 129L, 127L, 155L, 186L, 215L, 219L, 246L, 300L,
293L,


              280L, 351L, 384L, 390L, 397L, 347L, 200L, 195L, 175L,
91L,


              106L, 114L, 97L, 110L, 102L, 85L, 91L, 100L, 130L, 140L,


              159L, 167L, 191L, 243L, 231L, 236L, 255L, 262L, 281L,
294L,


              244L, 166L, 186L, 191L, 111L, 131L, 137L, 178L, 176L,
152L,


              136L, 120L, 116L, 132L, 147L, 186L, 231L, 278L, 334L,
355L,


              366L, 394L, 420L, 426L, 432L, 434L, 451L, 464L, 435L,
306L,


              384L, 392L, 352L, 374L, 364L, 317L, 329L, 349L, 409L,
449L,


              458L, 462L, 503L, 10L, 2L, 511L, 7L, 513L, 5L, 8L, 492L,


              151L, 126L, 153L, 174L, 179L, 203L, 182L, 132L, 131L,
115L,


              149L, 169L, 177L, 171L, 193L, 223L, 243L, 284L, 271L,
282L,


              303L, 380L, 442L, 433L, 379L, 463L, 454L, 69L, 80L, 77L,


              87L, 101L, 139L, 132L, 138L, 75L, 66L, 83L, 110L, 135L,
136L,


              137L, 139L, 94L, 481L, 476L, 509L, 57L, 509L, 500L,
249L,


              307L, 279L, 262L, 281L, 181L, 223L, 232L, 218L, 252L,
301L,


              328L, 363L, 422L, 485L, 11L, 37L, 24L, 23L, 35L, 50L,
40L,


              41L, 44L, 53L, 190L, 158L, 141L, 142L, 164L, 205L, 240L,


              256L, 257L, 263L, 262L, 265L, 282L, 299L, 314L, 312L,
348L,


              412L, 414L, 428L, 443L, 466L, 490L, 504L, 483L, 344L,
353L,


              277L, 189L, 260L, 281L, 276L, 252L, 244L, 192L, 185L,
217L,


              264L, 275L, 267L, 269L, 305L, 355L, 352L, 349L, 332L,
323L,


              345L, 371L, 342L, 113L, 117L, 120L, 202L, 235L, 231L,
245L,


              227L, 204L, 182L, 180L, 184L, 188L, 226L, 259L, 279L,
325L,


              369L, 375L, 376L, 395L, 418L, 423L, 441L, 399L, 248L,
242L,


              238L, 262L, 266L, 246L, 292L, 228L, 214L, 222L, 248L,
249L,


              290L, 333L, 408L, 505L, 510L, 28L, 13L, 34L, 43L, 49L,
52L,


              48L, 29L), .Label = c("1,003", "1,005", "1,011",
"1,035",


                                    "1,051", "1,059", "1,063",
"1,068", "1,076", "1,088", "1,091",


                                    "1,106", "1,108", "1,117",
"1,128", "1,149", "1,170", "1,175",


                                    "1,176", "1,183", "1,197",
"1,202", "1,245", "1,246", "1,247",


                                    "1,253", "1,258", "1,262",
"1,298", "1,305", "1,316", "1,347",


                                    "1,362", "1,365", "1,377",
"1,378", "1,390", "1,480", "1,486",


                                    "1,494", "1,508", "1,512",
"1,564", "1,597", "1,599", "1,684",


                                    "1,686", "1,705", "1,712",
"1,718", "1,726", "1,775", "1,789",


                                    "10,061", "10,693", "100", "102",
"107", "109", "11,437",


                                    "110", "113", "116", "12,711",
"120", "122", "123", "127",


                                    "128", "13,546", "13,880", "131",
"132", "133", "135", "136",


                                    "139", "14,293", "140", "142",
"143", "146", "148", "149",


                                    "151", "152", "153", "155", "156",
"157", "158", "159", "163",


                                    "165", "166", "167", "170", "171",
"172", "174", "175", "179",


                                    "180", "181", "183", "184", "185",
"186", "187", "188", "189",


                                    "190", "191", "192", "193", "194",
"195", "196", "197", "198",


                                    "2,205", "2,280", "2,510",
"2,692", "2,761", "200", "201",


                                    "204", "205", "207", "208", "209",
"210", "211", "213", "214",


                                    "216", "218", "219", "220", "221",
"222", "225", "228", "230",


                                    "232", "233", "234", "235", "236",
"237", "238", "241", "242",


                                    "243", "244", "245", "246", "247",
"248", "249", "251", "252",


                                    "253", "255", "256", "258", "259",
"260", "262", "263", "265",


                                    "266", "267", "268", "270", "271",
"272", "274", "275", "276",


                                    "277", "278", "282", "283", "284",
"286", "288", "290", "291",


                                    "292", "293", "294", "296", "299",
"3,143", "3,400", "3,420",


                                    "3,974", "300", "301", "302",
"304", "305", "307", "308",


                                    "309", "310", "311", "312", "313",
"314", "315", "316", "317",


                                    "318", "319", "321", "323", "324",
"325", "327", "328", "329",


                                    "330", "331", "332", "333", "336",
"338", "339", "340", "341",


                                    "343", "344", "345", "347", "348",
"349", "351", "352", "353",


                                    "354", "355", "356", "361", "362",
"363", "364", "365", "368",


                                    "369", "371", "372", "373", "376",
"377", "378", "379", "380",


                                    "382", "384", "385", "386", "389",
"390", "394", "395", "397",


                                    "398", "399", "4,191", "4,660",
"4,847", "400", "401", "402",


                                    "403", "404", "405", "406", "408",
"410", "411", "412", "413",


                                    "414", "415", "416", "417", "418",
"419", "420", "422", "425",


                                    "427", "428", "430", "431", "433",
"434", "436", "439", "440",


                                    "441", "442", "443", "449", "451",
"452", "453", "458", "460",


                                    "462", "464", "465", "469", "472",
"473", "474", "475", "478",


                                    "479", "482", "483", "486", "487",
"488", "489", "490", "491",


                                    "493", "496", "497", "5,036",
"5,551", "5,713", "502", "504",


                                    "509", "511", "512", "514", "516",
"517", "518", "520", "521",


                                    "523", "526", "527", "528", "533",
"534", "535", "537", "538",


                                    "544", "545", "548", "550", "553",
"555", "556", "558", "559",


                                    "560", "562", "564", "566", "569",
"570", "571", "572", "573",


                                    "574", "575", "578", "579", "580",
"581", "585", "586", "593",


                                    "6,671", "600", "601", "602",
"603", "604", "606", "608",


                                    "611", "612", "613", "615", "616",
"617", "618", "619", "621",


                                    "622", "626", "627", "629", "630",
"637", "639", "640", "641",


                                    "642", "644", "649", "652", "653",
"654", "657", "661", "668",


                                    "669", "670", "672", "673", "674",
"676", "677", "680", "681",


                                    "685", "688", "690", "691", "692",
"693", "694", "7,183",


                                    "700", "701", "702", "703", "704",
"707", "712", "716", "718",


                                    "719", "723", "730", "732", "735",
"737", "746", "749", "75",


                                    "759", "764", "765", "769", "775",
"777", "785", "786", "79",


                                    "790", "797", "800", "801", "802",
"805", "809", "813", "815",


                                    "820", "824", "825", "83", "831",
"84", "840", "846", "85",


                                    "852", "853", "863", "867", "879",
"882", "883", "885", "886",


                                    "887", "888", "89", "891", "90",
"904", "910", "914", "918",


                                    "92", "921", "925", "933", "935",
"936", "939", "941", "978",


                                    "98", "983", "991", "995", "997"),





             ))

______________________________________________
[hidden email] mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Reply | Threaded
Open this post in threaded view
|

Re: Systemfit Question

Spencer Graves-4
       I got an error when I tried to source the "structure" below
beginning with "Empdata <- structure(list(country = structure(c(1L, 1L,
1L, 1L, 1L, 1L,":  The source ended with " "98", "983", "991", "995",
"997"),
+
+
+
+
+
+                                      ))
+"


       I supplied ")" and got another "+".  When I supplied a second
")", I got, "Error in structure(c(15L, 487L, 397L, 323L, 233L, 315L,
370L, 309L, 268L,  :
   argument is missing, with no default".


       Please excuse if I missed something obvious that I might have
seen if I had read this entire thread, but it looked to me like this
example was intended to be self contained and does not seem to be so.


       Spencer


On 2018-05-16 17:21, Ngwinui Azenui wrote:

> I can't get my simultaneous equations to work using system fit. Please help.
>
> #Reproducible script
> Empdata<- read.csv("/Users/ngwinuiazenui/Documents/UPLOADemp.csv")
> View(Empdata)
> str(Empdata)
> Empdata$gnipc<-as.numeric(Empdata$gnipc)
> install.packages("systemfit")
> library("systemfit")
> pdata <- plm.data(Empdata, index=c("country","year"))
> install.packages("plm")
> library("plm")
>
> #Simultaneous equation 1
> eqprd<- prd ~ in_emp +gnipc
> eqemp<- in_emp ~ prd
> inst <- ~ gnipc
> system<- list(productivity = eqprd, industrialemp = eqemp)
> fit2sls<- systemfit(system, "2SLS", inst = inst, data = Empdata)
> #I get an error code when I run the last line
>
> #Simultaneous equation 2
> eqprd2<- prd ~  in_emp + lag(prd, -1) + gnipc
> eqemp2<- in_emp ~ prd + lag(prd, -1) + conflict
> system2<- list(productivity = eqprd2, industrialemp = eqemp2)
> inst <- ~ lag(prd, -1) + gnipc + conflict
> fit2sls<- systemfit(system2, "2SLS", inst = inst, data = Empdata)
> #I also wanted to include year and country fix effects(string
> variables) to the model in "Simultaneous equation 2)
>
> EmpSur<- systemfit(prd ~ in_emp + ag_emp, method = "SUR", data = pdata)
> #When I run the code below, I get a table with reults by country
> print(EmpSur)
> #An except of the output from running <print(EmpSur)
> #caf_(Intercept)      caf_in_emp      caf_ag_emp cod_(Intercept)
> #     -0.13243213     14.25003985     -8.57346398     -0.11581997
> #     cod_in_emp      cod_ag_emp eri_(Intercept)      eri_in_emp
> #     11.58115777    -13.29628753     -0.55615848     23.69284487
> #     eri_ag_emp eth_(Intercept)      eth_in_emp      eth_ag_emp
> #     -6.26339787      0.17436006     34.87782981     -4.80423119
> #gin_(Intercept)      gin_in_emp      gin_ag_emp gmb_(Intercept)
> #     0.01122243    -11.22002009     -0.64218211     -0.07104175
> #     gmb_in_emp      gmb_ag_emp gnb_(Intercept)      gnb_in_emp
> #     9.74177468      3.29511466     -0.07681684     17.00346867
> #     gnb_ag_emp gnq_(Intercept)      gnq_in_emp      gnq_ag_emp
> #     -3.40077876     -0.25706504     -1.09274682     -7.39122680
>
> #Below is the data. Sorry that it's a long panel dataset
> dput(Empdata)
> Empdata <- structure(list(country = structure(c(1L, 1L, 1L, 1L, 1L, 1L,
>                                                  1L, 1L, 1L, 1L, 1L,
> 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
>                                                  1L, 1L, 1L, 3L, 3L,
> 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
>                                                  3L, 3L, 3L, 3L, 3L,
> 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L,
>                                                  4L, 4L, 4L, 4L, 4L,
> 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
>                                                  4L, 4L, 4L, 4L, 4L,
> 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
>                                                  2L, 2L, 2L, 2L, 2L,
> 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 5L, 5L,
>                                                  5L, 5L, 5L, 5L, 5L,
> 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
>                                                  5L, 5L, 5L, 5L, 5L,
> 5L, 5L, 26L, 26L, 26L, 26L, 26L, 26L, 26L,
>                                                  26L, 26L, 26L, 26L,
> 26L, 26L, 26L, 26L, 26L, 26L, 26L, 26L, 26L,
>                                                  26L, 26L, 26L, 26L,
> 26L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L,
>                                                  6L, 6L, 6L, 6L, 6L,
> 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L,
>                                                  12L, 12L, 12L, 12L,
> 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L,
>                                                  12L, 12L, 12L, 12L,
> 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 7L,
>                                                  7L, 7L, 7L, 7L, 7L,
> 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L,
>                                                  7L, 7L, 7L, 7L, 7L,
> 7L, 7L, 7L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L,
>                                                  8L, 8L, 8L, 8L, 8L,
> 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L,
>                                                  8L, 10L, 10L, 10L,
> 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L,
>                                                  10L, 10L, 10L, 10L,
> 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L,
>                                                  9L, 9L, 9L, 9L, 9L,
> 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L,
>                                                  9L, 9L, 9L, 9L, 9L,
> 9L, 9L, 9L, 9L, 11L, 11L, 11L, 11L, 11L,
>                                                  11L, 11L, 11L, 11L,
> 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L,
>                                                  11L, 11L, 11L, 11L,
> 11L, 11L, 11L, 14L, 14L, 14L, 14L, 14L, 14L,
>                                                  14L, 14L, 14L, 14L,
> 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L,
>                                                  14L, 14L, 14L, 14L,
> 14L, 14L, 13L, 13L, 13L, 13L, 13L, 13L, 13L,
>                                                  13L, 13L, 13L, 13L,
> 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L,
>                                                  13L, 13L, 13L, 13L,
> 13L, 15L, 15L, 15L, 15L, 15L, 15L, 15L, 15L,
>                                                  15L, 15L, 15L, 15L,
> 15L, 15L, 15L, 15L, 15L, 15L, 15L, 15L, 15L,
>                                                  15L, 15L, 15L, 15L,
> 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L,
>                                                  19L, 19L, 19L, 19L,
> 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L,
>                                                  19L, 19L, 19L, 16L,
> 16L, 16L, 16L, 16L, 16L, 16L, 16L, 16L, 16L,
>                                                  16L, 16L, 16L, 16L,
> 16L, 16L, 16L, 16L, 16L, 16L, 16L, 16L, 16L,
>                                                  16L, 16L, 18L, 18L,
> 18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L,
>                                                  18L, 18L, 18L, 18L,
> 18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L,
>                                                  18L, 17L, 17L, 17L,
> 17L, 17L, 17L, 17L, 17L, 17L, 17L, 17L, 17L,
>                                                  17L, 17L, 17L, 17L,
> 17L, 17L, 17L, 17L, 17L, 17L, 17L, 17L, 17L,
>                                                  20L, 20L, 20L, 20L,
> 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L,
>                                                  20L, 20L, 20L, 20L,
> 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 21L,
>                                                  21L, 21L, 21L, 21L,
> 21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L,
>                                                  21L, 21L, 21L, 21L,
> 21L, 21L, 21L, 21L, 21L, 21L, 21L, 23L, 23L,
>                                                  23L, 23L, 23L, 23L,
> 23L, 23L, 23L, 23L, 23L, 23L, 23L, 23L, 23L,
>                                                  23L, 23L, 23L, 23L,
> 23L, 23L, 23L, 23L, 23L, 23L, 24L, 24L, 24L,
>                                                  24L, 24L, 24L, 24L,
> 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 24L,
>                                                  24L, 24L, 24L, 24L,
> 24L, 24L, 24L, 24L, 24L, 25L, 25L, 25L, 25L,
>                                                  25L, 25L, 25L, 25L,
> 25L, 25L, 25L, 25L, 25L, 25L, 25L, 25L, 25L,
>                                                  25L, 25L, 25L, 25L,
> 25L, 25L, 25L, 25L, 22L, 22L, 22L, 22L, 22L,
>                                                  22L, 22L, 22L, 22L,
> 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L,
>                                                  22L, 22L, 22L, 22L,
> 22L, 22L, 22L, 28L, 28L, 28L, 28L, 28L, 28L,
>                                                  28L, 28L, 28L, 28L,
> 28L, 28L, 28L, 28L, 28L, 28L, 28L, 28L, 28L,
>                                                  28L, 28L, 28L, 28L,
> 28L, 28L, 27L, 27L, 27L, 27L, 27L, 27L, 27L,
>                                                  27L, 27L, 27L, 27L,
> 27L, 27L, 27L, 27L, 27L, 27L, 27L, 27L, 27L,
>                                                  27L, 27L, 27L, 27L,
> 27L, 29L, 29L, 29L, 29L, 29L, 29L, 29L, 29L,
>                                                  29L, 29L, 29L, 29L,
> 29L, 29L, 29L, 29L, 29L, 29L, 29L, 29L, 29L,
>                                                  29L, 29L, 29L, 29L,
> 30L, 30L, 30L, 30L, 30L, 30L, 30L, 30L, 30L,
>                                                  30L, 30L, 30L, 30L,
> 30L, 30L, 30L, 30L, 30L, 30L, 30L, 30L, 30L,
>                                                  30L, 30L, 30L), .Label
> = c("ago", "bdi", "ben", "bfa", "caf",
>
>      "cod", "eri", "eth", "gin", "gmb", "gnb", "gnq", "lbr", "lso",
>
>      "mdg", "mli", "moz", "mrt", "mwi", "ner", "rwa", "sdn", "sen",
>
>      "sle", "som", "tcd", "tgo", "tza", "uga", "zmb"), class =
> "factor"),
>                            year = c(1991L, 1992L, 1993L, 1994L, 1995L,
> 1996L, 1997L,
>                                     1998L, 1999L, 2000L, 2001L, 2002L,
> 2003L, 2004L, 2005L, 2006L,
>                                     2007L, 2008L, 2009L, 2010L, 2011L,
> 2012L, 2013L, 2014L, 2015L,
>                                     1991L, 1992L, 1993L, 1994L, 1995L,
> 1996L, 1997L, 1998L, 1999L,
>                                     2000L, 2001L, 2002L, 2003L, 2004L,
> 2005L, 2006L, 2007L, 2008L,
>                                     2009L, 2010L, 2011L, 2012L, 2013L,
> 2014L, 2015L, 1991L, 1992L,
>                                     1993L, 1994L, 1995L, 1996L, 1997L,
> 1998L, 1999L, 2000L, 2001L,
>                                     2002L, 2003L, 2004L, 2005L, 2006L,
> 2007L, 2008L, 2009L, 2010L,
>                                     2011L, 2012L, 2013L, 2014L, 2015L,
> 1991L, 1992L, 1993L, 1994L,
>                                     1995L, 1996L, 1997L, 1998L, 1999L,
> 2000L, 2001L, 2002L, 2003L,
>                                     2004L, 2005L, 2006L, 2007L, 2008L,
> 2009L, 2010L, 2011L, 2012L,
>                                     2013L, 2014L, 2015L, 1991L, 1992L,
> 1993L, 1994L, 1995L, 1996L,
>                                     1997L, 1998L, 1999L, 2000L, 2001L,
> 2002L, 2003L, 2004L, 2005L,
>                                     2006L, 2007L, 2008L, 2009L, 2010L,
> 2011L, 2012L, 2013L, 2014L,
>                                     2015L, 1991L, 1992L, 1993L, 1994L,
> 1995L, 1996L, 1997L, 1998L,
>                                     1999L, 2000L, 2001L, 2002L, 2003L,
> 2004L, 2005L, 2006L, 2007L,
>                                     2008L, 2009L, 2010L, 2011L, 2012L,
> 2013L, 2014L, 2015L, 1991L,
>                                     1992L, 1993L, 1994L, 1995L, 1996L,
> 1997L, 1998L, 1999L, 2000L,
>                                     2001L, 2002L, 2003L, 2004L, 2005L,
> 2006L, 2007L, 2008L, 2009L,
>                                     2010L, 2011L, 2012L, 2013L, 2014L,
> 2015L, 1991L, 1992L, 1993L,
>                                     1994L, 1995L, 1996L, 1997L, 1998L,
> 1999L, 2000L, 2001L, 2002L,
>                                     2003L, 2004L, 2005L, 2006L, 2007L,
> 2008L, 2009L, 2010L, 2011L,
>                                     2012L, 2013L, 2014L, 2015L, 1991L,
> 1992L, 1993L, 1994L, 1995L,
>                                     1996L, 1997L, 1998L, 1999L, 2000L,
> 2001L, 2002L, 2003L, 2004L,
>                                     2005L, 2006L, 2007L, 2008L, 2009L,
> 2010L, 2011L, 2012L, 2013L,
>                                     2014L, 2015L, 1991L, 1992L, 1993L,
> 1994L, 1995L, 1996L, 1997L,
>                                     1998L, 1999L, 2000L, 2001L, 2002L,
> 2003L, 2004L, 2005L, 2006L,
>                                     2007L, 2008L, 2009L, 2010L, 2011L,
> 2012L, 2013L, 2014L, 2015L,
>                                     1991L, 1992L, 1993L, 1994L, 1995L,
> 1996L, 1997L, 1998L, 1999L,
>                                     2000L, 2001L, 2002L, 2003L, 2004L,
> 2005L, 2006L, 2007L, 2008L,
>                                     2009L, 2010L, 2011L, 2012L, 2013L,
> 2014L, 2015L, 1991L, 1992L,
>                                     1993L, 1994L, 1995L, 1996L, 1997L,
> 1998L, 1999L, 2000L, 2001L,
>                                     2002L, 2003L, 2004L, 2005L, 2006L,
> 2007L, 2008L, 2009L, 2010L,
>                                     2011L, 2012L, 2013L, 2014L, 2015L,
> 1991L, 1992L, 1993L, 1994L,
>                                     1995L, 1996L, 1997L, 1998L, 1999L,
> 2000L, 2001L, 2002L, 2003L,
>                                     2004L, 2005L, 2006L, 2007L, 2008L,
> 2009L, 2010L, 2011L, 2012L,
>                                     2013L, 2014L, 2015L, 1991L, 1992L,
> 1993L, 1994L, 1995L, 1996L,
>                                     1997L, 1998L, 1999L, 2000L, 2001L,
> 2002L, 2003L, 2004L, 2005L,
>                                     2006L, 2007L, 2008L, 2009L, 2010L,
> 2011L, 2012L, 2013L, 2014L,
>                                     2015L, 1991L, 1992L, 1993L, 1994L,
> 1995L, 1996L, 1997L, 1998L,
>                                     1999L, 2000L, 2001L, 2002L, 2003L,
> 2004L, 2005L, 2006L, 2007L,
>                                     2008L, 2009L, 2010L, 2011L, 2012L,
> 2013L, 2014L, 2015L, 1991L,
>                                     1992L, 1993L, 1994L, 1995L, 1996L,
> 1997L, 1998L, 1999L, 2000L,
>                                     2001L, 2002L, 2003L, 2004L, 2005L,
> 2006L, 2007L, 2008L, 2009L,
>                                     2010L, 2011L, 2012L, 2013L, 2014L,
> 2015L, 1991L, 1992L, 1993L,
>                                     1994L, 1995L, 1996L, 1997L, 1998L,
> 1999L, 2000L, 2001L, 2002L,
>                                     2003L, 2004L, 2005L, 2006L, 2007L,
> 2008L, 2009L, 2010L, 2011L,
>                                     2012L, 2013L, 2014L, 2015L, 1991L,
> 1992L, 1993L, 1994L, 1995L,
>                                     1996L, 1997L, 1998L, 1999L, 2000L,
> 2001L, 2002L, 2003L, 2004L,
>                                     2005L, 2006L, 2007L, 2008L, 2009L,
> 2010L, 2011L, 2012L, 2013L,
>                                     2014L, 2015L, 1991L, 1992L, 1993L,
> 1994L, 1995L, 1996L, 1997L,
>                                     1998L, 1999L, 2000L, 2001L, 2002L,
> 2003L, 2004L, 2005L, 2006L,
>                                     2007L, 2008L, 2009L, 2010L, 2011L,
> 2012L, 2013L, 2014L, 2015L,
>                                     1991L, 1992L, 1993L, 1994L, 1995L,
> 1996L, 1997L, 1998L, 1999L,
>                                     2000L, 2001L, 2002L, 2003L, 2004L,
> 2005L, 2006L, 2007L, 2008L,
>                                     2009L, 2010L, 2011L, 2012L, 2013L,
> 2014L, 2015L, 1991L, 1992L,
>                                     1993L, 1994L, 1995L, 1996L, 1997L,
> 1998L, 1999L, 2000L, 2001L,
>                                     2002L, 2003L, 2004L, 2005L, 2006L,
> 2007L, 2008L, 2009L, 2010L,
>                                     2011L, 2012L, 2013L, 2014L, 2015L,
> 1991L, 1992L, 1993L, 1994L,
>                                     1995L, 1996L, 1997L, 1998L, 1999L,
> 2000L, 2001L, 2002L, 2003L,
>                                     2004L, 2005L, 2006L, 2007L, 2008L,
> 2009L, 2010L, 2011L, 2012L,
>                                     2013L, 2014L, 2015L, 1991L, 1992L,
> 1993L, 1994L, 1995L, 1996L,
>                                     1997L, 1998L, 1999L, 2000L, 2001L,
> 2002L, 2003L, 2004L, 2005L,
>                                     2006L, 2007L, 2008L, 2009L, 2010L,
> 2011L, 2012L, 2013L, 2014L,
>                                     2015L, 1991L, 1992L, 1993L, 1994L,
> 1995L, 1996L, 1997L, 1998L,
>                                     1999L, 2000L, 2001L, 2002L, 2003L,
> 2004L, 2005L, 2006L, 2007L,
>                                     2008L, 2009L, 2010L, 2011L, 2012L,
> 2013L, 2014L, 2015L, 1991L,
>                                     1992L, 1993L, 1994L, 1995L, 1996L,
> 1997L, 1998L, 1999L, 2000L,
>                                     2001L, 2002L, 2003L, 2004L, 2005L,
> 2006L, 2007L, 2008L, 2009L,
>                                     2010L, 2011L, 2012L, 2013L, 2014L,
> 2015L, 1991L, 1992L, 1993L,
>                                     1994L, 1995L, 1996L, 1997L, 1998L,
> 1999L, 2000L, 2001L, 2002L,
>                                     2003L, 2004L, 2005L, 2006L, 2007L,
> 2008L, 2009L, 2010L, 2011L,
>                                     2012L, 2013L, 2014L, 2015L, 1991L,
> 1992L, 1993L, 1994L, 1995L,
>                                     1996L, 1997L, 1998L, 1999L, 2000L,
> 2001L, 2002L, 2003L, 2004L,
>                                     2005L, 2006L, 2007L, 2008L, 2009L,
> 2010L, 2011L, 2012L, 2013L,
>                                     2014L, 2015L, 1991L, 1992L, 1993L,
> 1994L, 1995L, 1996L, 1997L,
>                                     1998L, 1999L, 2000L, 2001L, 2002L,
> 2003L, 2004L, 2005L, 2006L,
>                                     2007L, 2008L, 2009L, 2010L, 2011L,
> 2012L, 2013L, 2014L, 2015L,
>                                     1991L, 1992L, 1993L, 1994L, 1995L,
> 1996L, 1997L, 1998L, 1999L,
>                                     2000L, 2001L, 2002L, 2003L, 2004L,
> 2005L, 2006L, 2007L, 2008L,
>                                     2009L, 2010L, 2011L, 2012L, 2013L,
> 2014L, 2015L, 1991L, 1992L,
>                                     1993L, 1994L, 1995L, 1996L, 1997L,
> 1998L, 1999L, 2000L, 2001L,
>                                     2002L, 2003L, 2004L, 2005L, 2006L,
> 2007L, 2008L, 2009L, 2010L,
>                                     2011L, 2012L, 2013L, 2014L, 2015L),
> ag_emp = c(NA, 0.002,
>
>             0.019, 0.002, -0.006, -0.014, -0.012, -0.01, -0.006,
> -0.005,
>
>             -0.008, -0.013, -0.01, -0.011, -0.017, -0.021, -0.022,
> -0.017,
>
>             -0.007, -0.004, -0.007, -0.007, -0.007, -0.005, -0.003, NA,
>
>             0.002, 0.002, 0.002, 0.002, -0.002, -0.001, -0.001, -0.001,
>
>             -0.001, 0, 0.001, 0.006, 0.004, 0.004, 0.003, 0.003, 0.003,
>
>             0.003, 0.004, 0.007, 0.004, 0.004, 0.004, -0.002, NA,
> 0.001,
>
>             -0.008, 0.003, 0, -0.005, -0.002, -0.004, -0.003, -0.005,
>
>             -0.001, -0.01, -0.003, 0, 0, -0.038, -0.138, -0.055,
> -0.024,
>
>             -0.018, -0.016, -0.014, -0.012, -0.009, -0.007, NA, -0.014,
>
>             -0.01, -0.008, -0.007, -0.008, -0.008, -0.008, -0.001, 0,
>
>             0.009, 0.009, 0.009, 0.008, 0.007, 0.004, 0.004, 0.004,
> 0.004,
>
>             0.004, -0.002, -0.001, -0.001, -0.002, 0, NA, 0.006, 0.003,
>
>             0.002, -0.003, 0.001, 0.001, -0.003, -0.001, 0.003, -0.003,
>
>             -0.002, 0.004, 0, 0.001, 0, -0.001, 0.001, 0.001, 0.001,
>
>             0.001, 0, 0.027, 0.011, 0.001, NA, -0.01, 0.009, -0.003,
>
>             -0.003, -0.004, -0.003, -0.003, -0.002, 0, -0.003, -0.004,
>
>             -0.004, -0.006, -0.01, -0.002, 0.001, 0.001, 0.002, 0.001,
>
>             -0.001, 0.002, 0.001, 0.002, 0.002, NA, 0.006, 0.008,
> 0.002,
>
>             -0.001, -0.002, 0.003, -0.003, 0.001, 0.005, -0.002,
> -0.005,
>
>             -0.006, -0.005, -0.006, -0.007, -0.006, -0.006, -0.004,
> -0.005,
>
>             -0.007, -0.006, -0.008, -0.007, -0.002, NA, -0.01, -0.015,
>
>             -0.014, -0.007, -0.004, -0.006, -0.014, 0, 0.001, -0.011,
>
>             -0.005, 0, -0.003, -0.005, 0.003, -0.002, -0.002, 0.004,
>
>             0.002, -0.001, 0, 0.001, -0.002, 0.038, NA, -0.009, -0.015,
>
>             -0.022, -0.009, -0.003, 0.006, 0.008, 0.009, 0.012, 0.004,
>
>             0.017, 0.012, 0.008, 0.009, 0.001, 0.001, 0.004, 0.001,
> -0.003,
>
>             -0.005, -0.006, -0.004, -0.005, -0.004, NA, -0.002, 0.001,
>
>             -0.001, -0.001, -0.004, -0.009, -0.008, -0.01, 0.034,
> -0.012,
>
>             0.01, 0.01, 0.009, -0.001, -0.006, -0.006, -0.005, -0.004,
>
>             -0.005, -0.005, -0.008, 0.002, -0.007, 0, NA, -0.008,
> -0.013,
>
>             -0.016, -0.009, -0.015, -0.015, -0.011, -0.014, -0.014,
> -0.016,
>
>             -0.018, -0.012, -0.016, -0.019, -0.016, -0.016, -0.016,
> -0.016,
>
>             -0.016, -0.019, -0.014, -0.016, -0.017, -0.003, NA, 0.001,
>
>             0.002, 0.002, 0.005, 0.001, 0.002, 0.003, 0.004, 0.006,
> 0.001,
>
>             0.003, 0.005, 0.004, 0.004, 0.004, 0.004, 0.003, 0.006,
> 0.007,
>
>             0.016, 0.002, 0.087, 0.002, 0, NA, -0.009, -0.007, -0.006,
>
>             -0.006, -0.013, -0.008, 0.019, 0.006, -0.01, -0.004,
> -0.002,
>
>             -0.002, -0.004, -0.005, -0.006, -0.006, -0.005, -0.005,
> -0.006,
>
>             -0.01, -0.005, -0.003, -0.005, -0.004, NA, -0.006, 0.011,
>
>             -0.006, -0.083, 0.084, -0.125, 0.066, 0.064, -0.005,
> -0.115,
>
>             0.101, 0.006, -0.064, -0.054, 0.098, -0.013, -0.007, 0.005,
>
>             -0.007, 0.007, -0.004, 0.006, -0.016, -0.008, NA, -0.027,
>
>             0.055, 0.044, 0.029, 0.023, -0.045, -0.014, -0.015, -0.012,
>
>             -0.003, -0.003, 0.02, 0.003, -0.009, -0.001, -0.013,
> -0.003,
>
>             0, 0.001, -0.006, -0.006, -0.006, -0.002, 0, NA, -0.002,
>
>             0.012, -0.01, 0.018, 0.006, 0.004, 0.002, 0, 0.009, -0.005,
>
>             0.006, 0.017, 0.036, 0.029, -0.006, -0.001, -0.002, 0.001,
>
>             -0.002, -0.007, -0.008, -0.008, -0.006, -0.003, NA, -0.004,
>
>             -0.003, -0.002, -0.128, -0.009, -0.006, -0.007, 0.005,
> 0.005,
>
>             0.008, 0.007, 0.006, -0.001, 0.006, 0.007, -0.005, 0.007,
>
>             -0.004, -0.003, -0.002, 0, -0.004, -0.001, -0.001, NA,
> -0.012,
>
>             -0.003, 0.007, 0.004, 0.001, -0.02, -0.005, 0.013, -0.012,
>
>             -0.017, -0.016, 0.026, -0.011, 0.02, 0.013, 0.011, 0.028,
>
>             0.026, 0.031, -0.005, -0.004, -0.007, -0.013, -0.007, NA,
>
>             -0.001, -0.006, -0.001, -0.008, -0.007, 0.001, -0.003,
> -0.005,
>
>             -0.001, -0.002, -0.004, -0.007, -0.007, -0.004, -0.011,
> -0.003,
>
>             -0.002, -0.002, -0.002, -0.005, -0.005, -0.005, -0.006,
> -0.003,
>
>             NA, 0.005, 0.016, 0.004, 0.004, 0, -0.007, 0.003, 0.002,
>
>             0, 0.003, -0.005, -0.004, -0.007, -0.008, -0.009, -0.011,
>
>             -0.011, -0.01, -0.012, -0.014, -0.016, -0.003, 0, 0.001,
>
>             NA, -0.009, 0.021, 0.006, -0.013, 0.02, 0.002, -0.003,
> 0.007,
>
>             -0.001, 0.001, 0.003, -0.001, 0, 0.001, -0.004, -0.001,
> -0.005,
>
>             -0.002, -0.004, -0.003, -0.006, -0.004, -0.006, -0.005, NA,
>
>             -0.006, 0.029, 0.029, 0.02, 0.001, -0.014, -0.033, -0.002,
>
>             -0.016, -0.015, 0.014, -0.015, -0.034, -0.023, -0.001,
> -0.003,
>
>             0.013, -0.015, -0.003, -0.008, -0.008, -0.002, 0.007, 0,
>
>             NA, 0, -0.015, -0.015, -0.015, -0.024, -0.015, -0.009,
> -0.016,
>
>             -0.026, -0.019, -0.007, -0.037, -0.027, -0.026, -0.033,
> 0.042,
>
>             0.034, 0.033, 0.041, 0.028, 0.016, 0.006, 0.009, 0.001, NA,
>
>             0.009, 0, -0.002, 0.003, -0.005, 0, 0.001, -0.001, 0.003,
>
>             0.006, -0.018, -0.008, 0.013, 0.001, -0.002, -0.003,
> -0.003,
>
>             0, 0, -0.003, -0.006, -0.009, -0.005, 0.016, NA, -0.007,
>
>             -0.005, -0.004, -0.006, -0.01, -0.01, -0.01, -0.008,
> -0.006,
>
>             -0.009, -0.006, -0.003, -0.005, -0.004, -0.006, -0.006,
> -0.005,
>
>             -0.003, -0.003, -0.005, -0.001, -0.001, -0.001, -0.001, NA,
>
>             -0.019, -0.14, 0.006, 0.009, 0.007, 0.044, -0.035, -0.004,
>
>             -0.002, -0.017, -0.005, -0.001, -0.004, -0.008, 0.004,
> -0.003,
>
>             -0.011, 0.006, -0.001, -0.004, -0.002, 0.007, -0.003,
> -0.001,
>
>             NA, 0.001, -0.002, -0.001, -0.005, -0.007, -0.002, 0.004,
>
>             0.001, 0.002, -0.012, 0.004, -0.009, -0.014, -0.017,
> -0.029,
>
>             -0.044, -0.014, -0.013, -0.017, -0.02, -0.017, -0.025,
> -0.006,
>
>             -0.008, NA, 0.004, 0.009, 0.001, -0.001, -0.008, -0.004,
>
>             0.003, 0.003, 0.003, 0.003, 0.003, 0.001, 0, 0.001, -0.002,
>
>             -0.002, -0.002, -0.002, -0.003, -0.004, -0.004, -0.004,
> -0.004,
>
>             -0.003, NA, -0.006, -0.005, 0.002, -0.001, -0.004, 0,
> -0.012,
>
>             -0.016, 0.008, -0.017, -0.011, -0.014, 0.015, 0.008, 0.003,
>
>             0.009, 0.014, 0.014, 0.012, 0.01, 0.009, 0.009, 0.006,
> 0.001,
>
>             NA, -0.001, 0, 0.015, 0.014, 0.025, 0.022, 0.027, 0.007,
>
>             0.007, -0.006, 0.013, 0.016, -0.01, -0.031, -0.003, 0.001,
>
>             0.016, -0.008, -0.043, 0.024, 0.025, 0, -0.015, -0.008),
>                            in_emp = c(NA, -0.001, -0.004, 0.005, 0.004,
> 0.001, 0, 0.001,
>                                       0.001, 0.001, 0.001, 0.003, 0,
> 0.003, 0.004, 0.003, 0.003,
>                                       0.001, -0.001, 0.002, 0.001,
> 0.002, 0.002, 0.002, 0.002,
>                                       NA, 0, 0.001, 0.001, 0.001, 0, 0,
> 0, 0, 0, 0, 0, 0.001, 0.002,
>                                       0.002, 0.002, 0.002, 0.002,
> 0.002, 0.002, 0.002, 0.001, 0.001,
>                                       0.001, 0.002, NA, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, -0.002,
>                                       -0.002, 0.024, 0.046, 0.021,
> 0.01, 0.008, 0.007, 0.006, 0.006,
>                                       0.006, 0.004, NA, -0.001, -0.001,
> 0, -0.001, 0, 0, 0, 0,
>                                       0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 0, 0, 0, 0, NA, -0.001,
>                                       0, 0.002, 0, -0.002, 0.003, 0, 0,
> -0.001, 0, 0, -0.001, 0.002,
>                                       0, 0.001, 0.001, 0, 0.001, 0.001,
> 0.001, 0.001, -0.008, 0.007,
>                                       0.001, NA, 0, -0.001, 0, 0, 0, 0,
> 0, 0, 0, 0.001, 0, 0.001,
>                                       0.002, -0.001, 0, 0, 0, 0, 0.001,
> -0.001, 0.001, 0, 0, 0.001,
>                                       NA, 0, -0.001, 0.001, 0.001,
> 0.001, 0.002, -0.002, 0, 0.002,
>                                       0, 0, 0.001, 0.002, 0, 0, 0.001,
> 0.001, 0, 0.002, 0.001,
>                                       0.001, 0, 0.002, 0.002, NA,
> 0.001, -0.001, 0.002, 0, 0, 0.002,
>                                       -0.003, -0.001, 0.001, 0.002, 0,
> 0.001, 0.004, 0, 0, 0.002,
>                                       0.001, -0.001, 0.001, 0, 0.001,
> -0.001, 0, -0.028, NA, -0.001,
>                                       0.002, 0.001, 0.001, -0.001,
> 0.002, 0.001, 0, 0, 0, 0.004,
>                                       0, 0, 0.002, 0, 0.001, -0.002,
> 0.002, 0, 0.002, 0.001, 0,
>                                       0, 0.001, NA, 0, 0, 0, 0, 0,
> -0.001, -0.001, -0.001, 0.004,
>                                       -0.001, 0.001, 0.001, 0.001,
> 0.001, 0.001, 0.001, 0.001,
>                                       0.001, 0.001, 0.001, 0.001,
> 0.002, 0.001, 0.003, NA, 0.001,
>                                       0.003, 0.002, 0.003, 0.003,
> 0.003, 0.003, 0.003, 0.003, 0.003,
>                                       0.003, 0.004, 0.004, 0.003,
> 0.004, 0.004, 0.004, 0.004, 0.004,
>                                       0.003, 0.005, 0.004, 0.004,
> 0.001, NA, 0, 0, 0, 0, 0, 0,
>                                       0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 0, -0.002, 0, 0.004,
>                                       0, 0.001, NA, 0, 0.001, 0.001,
> 0.002, 0.002, 0.002, -0.005,
>                                       0.008, -0.001, 0.002, 0.001,
> 0.002, 0.002, 0.002, 0.001,
>                                       0.002, 0.002, 0.002, 0.002,
> 0.003, -0.001, 0.002, 0.002,
>                                       0.003, NA, 0.001, 0.003, 0.001,
> -0.009, 0.014, -0.016, 0.01,
>                                       0.01, 0.002, -0.015, 0.015,
> 0.004, -0.009, -0.006, 0.018,
>                                       0, 0.001, 0.002, 0.001, 0.003,
> 0.001, 0.003, -0.001, 0.001,
>                                       NA, -0.017, -0.008, -0.007,
> -0.006, -0.004, 0.006, 0.002,
>                                       0.001, 0.003, 0, 0.001, -0.007,
> 0.006, 0, 0.003, 0, 0.001,
>                                       0.002, 0.003, 0.001, 0.001,
> 0.001, -0.001, 0.001, NA, -0.004,
>                                       -0.001, -0.006, -0.001, -0.003,
> -0.004, -0.004, -0.005, -0.005,
>                                       -0.007, -0.008, -0.021, -0.02,
> -0.012, 0.004, 0.005, 0.005,
>                                       0.007, 0.008, 0.009, 0.01, 0.008,
> 0.003, 0.002, NA, -0.002,
>                                       0.003, -0.002, 0.013, 0, 0, 0,
> 0.001, 0, 0, 0.002, 0.002,
>                                       0, 0.001, 0, 0, 0.001, 0.001, 0,
> 0, 0, 0.001, 0.001, 0.001,
>                                       NA, -0.001, 0, 0.001, 0.001,
> 0.001, 0, 0.001, 0.002, 0, 0.002,
>                                       -0.001, 0.004, 0, 0.003, 0.003,
> 0.002, 0.004, 0.004, 0.004,
>                                       0.001, 0, 0.001, 0.001, 0.002,
> NA, 0.003, 0.003, -0.001,
>                                       0.005, 0.001, -0.001, 0.003,
> 0.004, 0.001, 0.003, 0.001,
>                                       0.003, 0.002, 0.005, 0.006, 0,
> 0.002, 0.001, 0.004, 0.002,
>                                       0.003, 0.002, 0.002, 0.002, NA,
> -0.001, 0.002, 0, 0, 0.001,
>                                       0, 0.001, 0.001, 0, 0.001, 0, 0,
> 0, 0, 0, 0, 0, 0, 0, 0,
>                                       -0.001, 0.001, 0.001, 0.001, NA,
> -0.003, 0.006, 0.003, -0.003,
>                                       0.007, 0.001, 0.004, 0, 0, 0.004,
> 0, 0, -0.002, 0.002, 0,
>                                       0, 0.001, -0.002, 0.002, -0.001,
> 0.002, -0.001, 0.001, 0.001,
>                                       NA, 0, 0.001, 0.001, 0.001, 0, 0,
> 0, 0.001, 0, 0.001, 0.003,
>                                       0.004, 0.003, 0.004, 0.004,
> 0.004, 0.005, 0.004, 0.005, 0.005,
>                                       0.005, 0.006, 0.003, 0.004, NA,
> 0, -0.004, -0.004, -0.003,
>                                       -0.005, -0.003, -0.001, -0.002,
> -0.004, 0.003, 0.009, 0.001,
>                                       0.004, 0.005, 0.002, 0.012,
> 0.009, 0.01, 0.013, 0.01, 0.008,
>                                       0.003, 0.004, 0.003, NA, -0.003,
> 0.002, 0, -0.001, 0, -0.001,
>                                       0.002, -0.001, 0.004, -0.001,
> 0.001, -0.001, 0.002, 0, 0,
>                                       0.001, 0.001, 0, 0.001, 0.001,
> 0.002, 0.003, 0, -0.004, NA,
>                                       0, 0.001, -0.001, 0.002, 0, 0, 0,
> 0, 0, 0, 0, 0, 0.001, 0.001,
>                                       0.001, 0.001, 0, 0.001, 0.002,
> 0.001, 0.001, 0.002, 0.002,
>                                       0.002, NA, -0.001, 0.066, 0.001,
> 0.004, 0.003, 0.027, -0.018,
>                                       -0.001, 0.001, -0.006, 0, 0.003,
> 0, -0.001, 0.006, 0.003,
>                                       -0.003, 0.004, 0, 0.001, 0.003,
> 0.004, 0, 0.001, NA, -0.002,
>                                       -0.002, -0.002, -0.002, -0.002,
> -0.001, -0.001, -0.001, -0.001,
>                                       -0.002, 0.014, 0.003, 0.003,
> 0.003, 0.002, -0.007, 0.002,
>                                       0.002, 0.002, 0.002, 0.002,
> 0.002, 0.003, 0.003, NA, 0, -0.002,
>                                       0.005, 0.002, -0.002, 0.002, 0,
> 0.002, 0, 0.001, 0.001, 0.002,
>                                       0.001, 0.001, 0.001, 0.001,
> 0.001, 0.001, 0.001, 0.001, 0.001,
>                                       0.001, 0.001, 0.002, NA, 0,
> 0.001, 0.001, 0.001, 0.001, 0.002,
>                                       0.002, 0.006, 0.007, 0.006,
> 0.008, 0.008, -0.019, -0.015,
>                                       0.005, 0.006, 0.007, 0.008,
> -0.006, -0.006, -0.005, -0.005,
>                                       -0.003, 0.002, NA, -0.002,
> -0.002, 0, -0.001, -0.001, -0.004,
>                                       -0.003, -0.005, -0.004, 0.002,
> 0.004, 0.004, 0.002, 0, 0.001,
>                                       0.002, 0.003, 0.002, -0.003,
> 0.004, 0.004, 0.001, 0,  prd = c(NA, -0.095, -0.276, -0.032, 0.06,
>
>                              0.083, 0.053, 0.028, -0.006, -0.003,
> -0.046, 0.099, -0.031,
>
>                              0.105, 0.143, 0.123, 0.127, 0.048, 0.075,
> 0.114, 0.101, 0.164,
>
>                              0.081, 0.012, -0.033, NA, -0.017, -0.019,
> -0.019, 0.033,
>
>                              0.007, 0.023, -0.005, -0.006, 0.039,
> 0.029, -0.011, -0.01,
>
>                              0.008, -0.016, 0.008, 0.009, 0.008,
> -0.002, -0.018, -0.004,
>
>                              -0.013, 0.038, 0.031, 0.01, NA, -0.023,
> 0.023, -0.059, 0.024,
>
>                              0.066, 0.032, 0.051, 0.042, 0.018, 0.031,
> 0.023, 0.097, 0.007,
>
>                              0.033, 0.089, -0.002, 0.038, -0.039,
> 0.051, 0.021, 0.031,
>
>                              0.036, 0.023, 0.018, NA, 0.025, -0.074,
> -0.098, -0.088, -0.041,
>
>                              0.014, 0.027, -0.013, -0.044, -0.018, 0,
> -0.037, -0.023,
>
>                              -0.01, -0.003, -0.013, 0.005, -0.009,
> 0.004, 0.006, 0.011,
>
>                              0.007, 0.015, -0.07, NA, -0.05, -0.037,
> 0.015, 0.013, -0.091,
>
>                              0.007, 0.008, 0.076, 0.005, -0.018, -0.03,
> -0.097, -0.003,
>
>                              0.005, 0.02, 0.016, 0.001, -0.006, -0.002,
> 0.008, 0.017,
>
>                              -0.36, -0.012, 0.007, NA, 0.007, -0.088,
> 0.051, -0.024, -0.011,
>
>                              0.029, 0.038, -0.043, -0.04, 0.08, 0.042,
> 0.191, 0.264, 0.112,
>
>                              -0.046, 0.133, -0.018, 0.058, 0.076,
> -0.069, 0.071, 0.038,
>
>                              0.194, 0.026, NA, -0.132, -0.182, -0.075,
> -0.026, -0.017,
>
>                              -0.1, -0.02, -0.05, -0.101, -0.062, 0.001,
> 0.024, 0.033,
>
>                              0.029, 0.019, 0.041, 0.032, -0.003, 0.037,
> 0.036, 0.036,
>
>                              0.049, 0.062, 0.034, NA, 0.149, 0.064,
> -0.018, 0.062, 0.307,
>
>                              0.605, 0.18, 0.187, 0.119, 0.575, 0.171,
> 0.093, 0.165, 0.061,
>
>                              0.006, 0.111, 0.15, -0.014, -0.12, 0.043,
> 0.051, -0.073,
>
>                              -0.031, -0.093, NA, 0.137, 0.11, 0.19,
> 0.213, 0.009, 0.053,
>
>                              0.063, -0.017, -0.028, -0.175, -0.004,
> 0.01, -0.071, -0.015,
>
>                              -0.029, 0.003, -0.108, 0.022, -0.006,
> 0.067, 0.05, 0.008,
>
>                              0.031, 0.025, NA, -0.12, 0.091, -0.009,
> 0.032, 0.087, 0.008,
>
>                              -0.062, 0.033, -0.012, 0.056, -0.03,
> -0.06, 0.067, 0.085,
>
>                              0.084, 0.083, 0.081, 0.065, 0.071, 0.104,
> 0.061, 0.058, 0.074,
>
>                              0.071, NA, 0.002, 0.029, -0.025, -0.004,
> 0.016, -0.02, 0,
>
>                              0.037, 0.021, 0.026, -0.018, -0.015,
> -0.037, -0.032, -0.018,
>
>                              -0.001, 0.023, 0.038, 0.019, -0.059,
> 0.017, 0.019, -0.002,
>
>                              0.011, NA, -0.029, -0.014, -0.018, -0.01,
> 0.011, 0.014, 0.025,
>
>                              0.021, -0.002, 0.016, 0.014, -0.016,
> -0.006, -0.007, -0.015,
>
>                              0.008, 0.001, -0.039, 0.014, 0.019, 0.009,
> -0.1, -0.023,
>
>                              0.002, NA, 0.003, -0.004, 0.01, 0.007,
> 0.022, 0.034, -0.189,
>
>                              0.089, -0.026, -0.005, -0.048, -0.008,
> -0.01, 0.021, 0.005,
>
>                              0.004, 0.011, 0.008, 0.018, 0.062, -0.039,
> 0.001, -0.026,
>
>                              -0.017, NA, 0.058, -0.003, 0.034, 0.128,
> -0.075, 0.245, -0.073,
>
>                              -0.066, 0.041, 0.249, -0.12, 0.022, 0.117,
> 0.089, -0.103,
>
>                              0.051, 0.047, 0.009, 0.072, 0.002, 0.034,
> 0.011, 0.029, 0.033,
>
>                              NA, -0.339, -0.325, -0.224, -0.073, 0.047,
> 0.903, 0.19, 0.155,
>
>                              0.176, 0.239, 0.287, -0.338, -0.083,
> 0.063, 0.041, 0.074,
>
>                              0.122, 0.046, 0.062, 0.022, 0.061, 0.051,
> -0.011, -0.049,
>
>                              NA, -0.013, -0.017, -0.007, -0.032,
> -0.014, 0.005, 0.013,
>
>                              0.023, 0.008, 0.043, -0.145, 0.057, 0.005,
> 0, 0.021, 0.027,
>
>                              0.022, -0.154, 0.055, -0.027, 0.001,
> -0.005, 0.009, 0.001,
>
>                              NA, -0.092, 0.122, -0.134, 0.095, 0.088,
> 0.039, 0.015, -0.017,
>
>                              -0.023, -0.077, 0.001, 0.017, 0.057,
> -0.03, -0.026, 0.074,
>
>                              0.029, 0.05, 0.046, 0.011, -0.013, 0.023,
> 0.026, -0.01, NA,
>
>                              0.063, 0.01, 0, -0.028, 0.019, 0.022,
> 0.044, -0.003, 0.027,
>
>                              0.134, 0.071, 0.036, -0.098, 0.03, 0.066,
> 0.025, 0.029, -0.005,
>
>                              0.038, 0.033, 0.109, 0.035, 0.066, 0.041,
> NA, -0.028, 0.092,
>
>                              -0.005, 0.029, -0.028, -0.098, -0.012,
> 0.054, -0.049, -0.013,
>
>                              -0.031, 0.022, 0.014, 0.054, 0.173,
> -0.024, -0.032, 0.023,
>
>                              -0.04, 0.006, 0.016, 0.04, 0.034, -0.016,
> NA, -0.073, 0.029,
>
>                              0.024, -0.03, 0.056, 0.09, 0.055, 0.047,
> -0.018, 0.097, 0.049,
>
>                              0.038, 0.066, 0.056, 0.071, 0.058, 0.061,
> 0.048, 0.061, 0.06,
>
>                              0.064, 0.035, 0.039, 0.042, NA, -0.006,
> -0.057, -0.026, 0.006,
>
>                              -0.03, -0.034, 0.066, -0.038, -0.067,
> 0.024, 0.002, 0.001,
>
>                              -0.053, 0.036, 0.026, -0.006, 0.059,
> -0.048, 0.049, -0.021,
>
>                              0.102, 0.003, 0.033, -0.017, NA, 0.119,
> -0.115, -0.469, 0.33,
>
>                              0.067, 0.096, 0.061, 0.005, 0.048, 0.062,
> 0.066, -0.007,
>
>                              0.061, 0.069, 0.051, 0.055, 0.053, 0.039,
> 0.034, 0.051, 0.073,
>
>                              0.022, 0.035, 0.032, NA, -0.02, -0.042,
> 0.011, 0.047, 0.053,
>
>                              0.021, 0.03, 0.055, 0.044, 0.014, -0.047,
> 0.091, 0.055, 0.049,
>
>                              0.032, -0.014, 0.002, -0.007, -0.015,
> -0.027, 0.016, 0.007,
>
>                              0.001, 0.025, NA, -0.092, 0, 0.061,
> -0.113, -0.227, -0.154,
>
>                              -0.039, -0.091, -0.012, 0.11, 0.213,
> 0.044, -0.001, -0.002,
>
>                              0.006, 0.047, 0.024, 0.005, 0.025, 0.036,
> 0.135, 0.194, 0.023,
>
>                              -0.206, NA, -0.099, 0.007, -0.186, -0.005,
> 0.026, -0.047,
>
>                              -0.006, 0.006, -0.004, 0.017, 0.008,
> 0.003, 0.005, 0.004,
>
>                              0.007, -0.001, 0.004, 0.002, -0.002,
> 0.002, -0.026, -0.004,
>
>                              0.014, -0.008, NA, 0.005, -0.019, -0.002,
> 0.007, 0.069, 0.059,
>
>                              0.067, 0.037, 0.068, 0.065, 0.034, 0.043,
> 0.026, 0.045, 0.056,
>
>                              0.07, -0.179, -0.009, 0.033, -0.006,
> -0.061, 0.024, -0.004,
>
>                              0.019, NA, 0.018, 0.008, -0.017, -0.008,
> 0.011, 0.012, 0.02,
>
>                              0.018, 0.018, 0.046, 0.025, 0.033, 0.049,
> 0.045, 0.038, 0.023,
>
>                              0.042, 0.031, 0.055, 0.073, 0.042, 0.068,
> 0.04, 0.038, NA,
>
>                              -0.065, -0.163, 0.098, 0.04, 0.061, 0.101,
> -0.052, -0.008,
>
>                              -0.05, -0.061, -0.065, -0.006, -0.005,
> -0.026, 0.04, -0.02,
>
>                              -0.011, 0.012, 0.001, 0.052, -0.033,
> 0.028, 0.024, 0.02,
>
>                              NA, 0.023, 0.035, 0.061, 0.051, 0.033,
> 0.023, 0.069, 0.055,
>
>                              0.002, 0.076, 0.035, 0.024, 0.02, 0.05,
> 0.026, 0.018, 0.035,
>
>                              0.019, 0.042, 0.011, -0.01, 0.006, 0.004,
> 0.014, NA, -0.082,
>
>                              0.087, -0.157, 0, 0.004, -0.017, -0.044,
> 0.025, 0.011, 0.034,
>
>                              0.006, 0.021, 0.059, 0.09, 0.055, 0.05,
> 0.024, 0.067, 0.113,
>
>                              -0.006, 0.009, 0.012, 0.034, 0.003), gnipc
> = structure(c(15L,
>
>
>                487L, 397L, 323L, 233L, 315L, 370L, 309L, 268L, 391L,
> 381L,
>
>
>                456L, 488L, 14L, 46L, 122L, 124L, 197L, 125L, 196L,
> 272L,
>
>
>                273L, 274L, 335L, 199L, 264L, 286L, 285L, 183L, 265L,
> 284L,
>
>
>                263L, 282L, 277L, 254L, 257L, 288L, 343L, 379L, 386L,
> 402L,
>
>
>                441L, 472L, 461L, 453L, 473L, 474L, 496L, 502L, 465L,
> 240L,
>
>
>                250L, 230L, 116L, 151L, 164L, 146L, 167L, 172L, 143L,
> 149L,
>
>
>                168L, 224L, 252L, 279L, 293L, 321L, 372L, 360L, 367L,
> 411L,
>
>
>                418L, 429L, 425L, 368L, 132L, 112L, 96L, 103L, 95L, 80L,
>
>
>                90L, 81L, 68L, 62L, 56L, 500L, 478L, 500L, 83L, 94L,
> 97L,
>
>
>                118L, 130L, 145L, 157L, 161L, 165L, 179L, 180L, 318L,
> 318L,
>
>
>                287L, 173L, 233L, 200L, 179L, 181L, 184L, 162L, 163L,
> 171L,
>
>
>                195L, 220L, 232L, 248L, 281L, 316L, 314L, 314L, 342L,
> 339L,
>
>
>                239L, 266L, 246L, 194L, 195L, 157L, 119L, 148L, 163L,
> 146L,
>
>
>                167L, 137L, 109L, 142L, 156L, 181L, 302L, 373L, 450L,
> 457L,
>
>
>                497L, 467L, 486L, 499L, 9L, 1L, 4L, 470L, 247L, 216L,
> 284L,
>
>
>                127L, 108L, 107L, 129L, 134L, 105L, 99L, 95L, 97L, 97L,
> 111L,
>
>
>                138L, 160L, 180L, 211L, 187L, 216L, 253L, 283L, 291L,
> 323L,
>
>
>                327L, 295L, 296L, 292L, 210L, 270L, 374L, 3L, 484L, 42L,
>
>
>                47L, 51L, 121L, 123L, 198L, 337L, 336L, 385L, 55L, 60L,
> 54L,
>
>
>                71L, 78L, 70L, 64L, 436L, 67L, 79L, 93L, 129L, 133L,
> 152L,
>
>
>                158L, 163L, 154L, 131L, 136L, 126L, 145L, 185L, 179L,
> 194L,
>
>
>                214L, 219L, 296L, 322L, 378L, 422L, 452L, 480L, 508L,
> 166L,
>
>
>                119L, 76L, 68L, 72L, 77L, 79L, 67L, 62L, 65L, 63L, 59L,
> 63L,
>
>
>                74L, 92L, 112L, 149L, 206L, 225L, 195L, 227L, 308L,
> 329L,
>
>
>                363L, 405L, 444L, 446L, 448L, 440L, 447L, 455L, 432L,
> 439L,
>
>
>                417L, 393L, 344L, 289L, 234L, 267L, 283L, 287L, 330L,
> 383L,
>
>
>                358L, 365L, 340L, 333L, 320L, 297L, 313L, 394L, 406L,
> 388L,
>
>
>                387L, 416L, 418L, 403L, 365L, 356L, 319L, 304L, 306L,
> 346L,
>
>
>                354L, 294L, 295L, 400L, 419L, 401L, 404L, 398L, 410L,
> 425L,
>
>
>                437L, 424L, 382L, 391L, 397L, 407L, 415L, 298L, 225L,
> 191L,
>
>
>                209L, 184L, 188L, 207L, 241L, 264L, 290L, 288L, 323L,
> 377L,
>
>
>                357L, 359L, 430L, 389L, 401L, 384L, 373L, 459L, 491L,
> 482L,
>
>
>                477L, 494L, 479L, 489L, 507L, 498L, 471L, 438L, 413L,
> 506L,
>
>
>                16L, 18L, 17L, 26L, 25L, 21L, 39L, 45L, 38L, 33L, 30L,
> 22L,
>
>
>                146L, 86L, 61L, 495L, 493L, 56L, 91L, 102L, 98L, 104L,
> 84L,
>
>
>                90L, 58L, 73L, 88L, 104L, 136L, 156L, 172L, 184L, 237L,
> 241L,
>
>
>                261L, 266L, 341L, 161L, 182L, 201L, 165L, 170L, 221L,
> 186L,
>
>
>                191L, 186L, 186L, 220L, 203L, 251L, 186L, 217L, 231L,
> 306L,
>
>
>                367L, 320L, 326L, 367L, 350L, 361L, 329L, 311L, 271L,
> 230L,
>
>
>                251L, 144L, 148L, 128L, 195L, 193L, 187L, 179L, 169L,
> 91L,
>
>
>                79L, 82L, 89L, 191L, 213L, 247L, 285L, 310L, 331L, 230L,
>
>
>                200L, 232L, 245L, 234L, 236L, 229L, 150L, 189L, 212L,
> 183L,
>
>
>                208L, 216L, 174L, 191L, 212L, 258L, 295L, 318L, 338L,
> 379L,
>
>
>                424L, 428L, 428L, 469L, 452L, 460L, 475L, 447L, 470L,
> 468L,
>
>
>                396L, 423L, 431L, 427L, 424L, 375L, 358L, 324L, 309L,
> 341L,
>
>
>                362L, 421L, 445L, 501L, 512L, 19L, 6L, 20L, 32L, 31L,
> 36L,
>
>
>                27L, 12L, 146L, 89L, 91L, 89L, 89L, 112L, 141L, 144L,
> 160L,
>
>
>                141L, 129L, 127L, 155L, 186L, 215L, 219L, 246L, 300L,
> 293L,
>
>
>                280L, 351L, 384L, 390L, 397L, 347L, 200L, 195L, 175L,
> 91L,
>
>
>                106L, 114L, 97L, 110L, 102L, 85L, 91L, 100L, 130L, 140L,
>
>
>                159L, 167L, 191L, 243L, 231L, 236L, 255L, 262L, 281L,
> 294L,
>
>
>                244L, 166L, 186L, 191L, 111L, 131L, 137L, 178L, 176L,
> 152L,
>
>
>                136L, 120L, 116L, 132L, 147L, 186L, 231L, 278L, 334L,
> 355L,
>
>
>                366L, 394L, 420L, 426L, 432L, 434L, 451L, 464L, 435L,
> 306L,
>
>
>                384L, 392L, 352L, 374L, 364L, 317L, 329L, 349L, 409L,
> 449L,
>
>
>                458L, 462L, 503L, 10L, 2L, 511L, 7L, 513L, 5L, 8L, 492L,
>
>
>                151L, 126L, 153L, 174L, 179L, 203L, 182L, 132L, 131L,
> 115L,
>
>
>                149L, 169L, 177L, 171L, 193L, 223L, 243L, 284L, 271L,
> 282L,
>
>
>                303L, 380L, 442L, 433L, 379L, 463L, 454L, 69L, 80L, 77L,
>
>
>                87L, 101L, 139L, 132L, 138L, 75L, 66L, 83L, 110L, 135L,
> 136L,
>
>
>                137L, 139L, 94L, 481L, 476L, 509L, 57L, 509L, 500L,
> 249L,
>
>
>                307L, 279L, 262L, 281L, 181L, 223L, 232L, 218L, 252L,
> 301L,
>
>
>                328L, 363L, 422L, 485L, 11L, 37L, 24L, 23L, 35L, 50L,
> 40L,
>
>
>                41L, 44L, 53L, 190L, 158L, 141L, 142L, 164L, 205L, 240L,
>
>
>                256L, 257L, 263L, 262L, 265L, 282L, 299L, 314L, 312L,
> 348L,
>
>
>                412L, 414L, 428L, 443L, 466L, 490L, 504L, 483L, 344L,
> 353L,
>
>
>                277L, 189L, 260L, 281L, 276L, 252L, 244L, 192L, 185L,
> 217L,
>
>
>                264L, 275L, 267L, 269L, 305L, 355L, 352L, 349L, 332L,
> 323L,
>
>
>                345L, 371L, 342L, 113L, 117L, 120L, 202L, 235L, 231L,
> 245L,
>
>
>                227L, 204L, 182L, 180L, 184L, 188L, 226L, 259L, 279L,
> 325L,
>
>
>                369L, 375L, 376L, 395L, 418L, 423L, 441L, 399L, 248L,
> 242L,
>
>
>                238L, 262L, 266L, 246L, 292L, 228L, 214L, 222L, 248L,
> 249L,
>
>
>                290L, 333L, 408L, 505L, 510L, 28L, 13L, 34L, 43L, 49L,
> 52L,
>
>
>                48L, 29L), .Label = c("1,003", "1,005", "1,011",
> "1,035",
>
>
>                                      "1,051", "1,059", "1,063",
> "1,068", "1,076", "1,088", "1,091",
>
>
>                                      "1,106", "1,108", "1,117",
> "1,128", "1,149", "1,170", "1,175",
>
>
>                                      "1,176", "1,183", "1,197",
> "1,202", "1,245", "1,246", "1,247",
>
>
>                                      "1,253", "1,258", "1,262",
> "1,298", "1,305", "1,316", "1,347",
>
>
>                                      "1,362", "1,365", "1,377",
> "1,378", "1,390", "1,480", "1,486",
>
>
>                                      "1,494", "1,508", "1,512",
> "1,564", "1,597", "1,599", "1,684",
>
>
>                                      "1,686", "1,705", "1,712",
> "1,718", "1,726", "1,775", "1,789",
>
>
>                                      "10,061", "10,693", "100", "102",
> "107", "109", "11,437",
>
>
>                                      "110", "113", "116", "12,711",
> "120", "122", "123", "127",
>
>
>                                      "128", "13,546", "13,880", "131",
> "132", "133", "135", "136",
>
>
>                                      "139", "14,293", "140", "142",
> "143", "146", "148", "149",
>
>
>                                      "151", "152", "153", "155", "156",
> "157", "158", "159", "163",
>
>
>                                      "165", "166", "167", "170", "171",
> "172", "174", "175", "179",
>
>
>                                      "180", "181", "183", "184", "185",
> "186", "187", "188", "189",
>
>
>                                      "190", "191", "192", "193", "194",
> "195", "196", "197", "198",
>
>
>                                      "2,205", "2,280", "2,510",
> "2,692", "2,761", "200", "201",
>
>
>                                      "204", "205", "207", "208", "209",
> "210", "211", "213", "214",
>
>
>                                      "216", "218", "219", "220", "221",
> "222", "225", "228", "230",
>
>
>                                      "232", "233", "234", "235", "236",
> "237", "238", "241", "242",
>
>
>                                      "243", "244", "245", "246", "247",
> "248", "249", "251", "252",
>
>
>                                      "253", "255", "256", "258", "259",
> "260", "262", "263", "265",
>
>
>                                      "266", "267", "268", "270", "271",
> "272", "274", "275", "276",
>
>
>                                      "277", "278", "282", "283", "284",
> "286", "288", "290", "291",
>
>
>                                      "292", "293", "294", "296", "299",
> "3,143", "3,400", "3,420",
>
>
>                                      "3,974", "300", "301", "302",
> "304", "305", "307", "308",
>
>
>                                      "309", "310", "311", "312", "313",
> "314", "315", "316", "317",
>
>
>                                      "318", "319", "321", "323", "324",
> "325", "327", "328", "329",
>
>
>                                      "330", "331", "332", "333", "336",
> "338", "339", "340", "341",
>
>
>                                      "343", "344", "345", "347", "348",
> "349", "351", "352", "353",
>
>
>                                      "354", "355", "356", "361", "362",
> "363", "364", "365", "368",
>
>
>                                      "369", "371", "372", "373", "376",
> "377", "378", "379", "380",
>
>
>                                      "382", "384", "385", "386", "389",
> "390", "394", "395", "397",
>
>
>                                      "398", "399", "4,191", "4,660",
> "4,847", "400", "401", "402",
>
>
>                                      "403", "404", "405", "406", "408",
> "410", "411", "412", "413",
>
>
>                                      "414", "415", "416", "417", "418",
> "419", "420", "422", "425",
>
>
>                                      "427", "428", "430", "431", "433",
> "434", "436", "439", "440",
>
>
>                                      "441", "442", "443", "449", "451",
> "452", "453", "458", "460",
>
>
>                                      "462", "464", "465", "469", "472",
> "473", "474", "475", "478",
>
>
>                                      "479", "482", "483", "486", "487",
> "488", "489", "490", "491",
>
>
>                                      "493", "496", "497", "5,036",
> "5,551", "5,713", "502", "504",
>
>
>                                      "509", "511", "512", "514", "516",
> "517", "518", "520", "521",
>
>
>                                      "523", "526", "527", "528", "533",
> "534", "535", "537", "538",
>
>
>                                      "544", "545", "548", "550", "553",
> "555", "556", "558", "559",
>
>
>                                      "560", "562", "564", "566", "569",
> "570", "571", "572", "573",
>
>
>                                      "574", "575", "578", "579", "580",
> "581", "585", "586", "593",
>
>
>                                      "6,671", "600", "601", "602",
> "603", "604", "606", "608",
>
>
>                                      "611", "612", "613", "615", "616",
> "617", "618", "619", "621",
>
>
>                                      "622", "626", "627", "629", "630",
> "637", "639", "640", "641",
>
>
>                                      "642", "644", "649", "652", "653",
> "654", "657", "661", "668",
>
>
>                                      "669", "670", "672", "673", "674",
> "676", "677", "680", "681",
>
>
>                                      "685", "688", "690", "691", "692",
> "693", "694", "7,183",
>
>
>                                      "700", "701", "702", "703", "704",
> "707", "712", "716", "718",
>
>
>                                      "719", "723", "730", "732", "735",
> "737", "746", "749", "75",
>
>
>                                      "759", "764", "765", "769", "775",
> "777", "785", "786", "79",
>
>
>                                      "790", "797", "800", "801", "802",
> "805", "809", "813", "815",
>
>
>                                      "820", "824", "825", "83", "831",
> "84", "840", "846", "85",
>
>
>                                      "852", "853", "863", "867", "879",
> "882", "883", "885", "886",
>
>
>                                      "887", "888", "89", "891", "90",
> "904", "910", "914", "918",
>
>
>                                      "92", "921", "925", "933", "935",
> "936", "939", "941", "978",
>
>
>                                      "98", "983", "991", "995", "997"),
>
>
>
>
>
>               ))
>
> ______________________________________________
> [hidden email] mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

______________________________________________
[hidden email] mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.