[R] Systemfit Question
Spencer Graves
@pencer@gr@ve@ @end|ng |rom e||ect|vede|en@e@org
Thu May 17 08:24:39 CEST 2018
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 using 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.
More information about the R-help
mailing list