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"),
))
????? 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"), > > > > > > )) > > ______________________________________________ > R-help at r-project.org 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.