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@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.