dput(sample[1:50, ])
structure(list(year = c(2005L, 2006L, 2007L, 2008L, 2009L, 2010L,
2011L, 2005L, 2006L, 2007L, 2008L, 2009L, 2010L, 2011L, 2005L,
2006L, 2007L, 2008L, 2009L, 2010L, 2011L, 2005L, 2006L, 2007L,
2008L, 2009L, 2010L, 2011L, 2005L, 2006L, 2007L, 2008L, 2009L,
2010L, 2011L, 2005L, 2006L, 2007L, 2008L, 2009L, 2010L, 2011L,
2005L, 2006L, 2007L, 2008L, 2009L, 2010L, 2011L, 2005L), Y = c(4173.217045,
4277.996451, 4319.290767, 4549.83007, 4435.450694, 4280.088368,
4020.781806, 5877.274684, 5976.119041, 6014.478399, 6216.265093,
6236.914349, 6114.028861, 5889.205762, 6079.081436, 6205.04293,
6146.156668, 6404.103999, 6459.669141, 6366.698522, 6068.932465,
6046.077215, 6147.887524, 6109.346106, 6361.920033, 6372.34715,
6235.21232, 6079.461274, 5685.307721, 5647.604075, 5694.551862,
5985.826031, 6017.33036, 5964.862862, 5760.725342, 5704.078622,
5431.702292, 5582.624809, 5883.925628, 5832.585978, 5891.687208,
5638.702178, 5869.447414, 5945.162792, 5954.481229, 6159.511579,
6189.853063, 6019.20501, 5841.450154, 6081.946856), municipality = c("mu_1",
"mu_1", "mu_1", "mu_1", "mu_1", "mu_1", "mu_1", "mu_2", "mu_2",
"mu_2", "mu_2", "mu_2", "mu_2", "mu_2", "mu_3", "mu_3", "mu_3",
"mu_3", "mu_3", "mu_3", "mu_3", "mu_4", "mu_4", "mu_4", "mu_4",
"mu_4", "mu_4", "mu_4", "mu_5", "mu_5", "mu_5", "mu_5", "mu_5",
"mu_5", "mu_5", "mu_6", "mu_6", "mu_6", "mu_6", "mu_6", "mu_6",
"mu_6", "mu_7", "mu_7", "mu_7", "mu_7", "mu_7", "mu_7", "mu_7",
"mu_8"), region = c("re_re_1", "re_re_1", "re_re_1", "re_re_1",
"re_re_1", "re_re_1", "re_re_1", "re_re_1", "re_re_1", "re_re_1",
"re_re_1", "re_re_1", "re_re_1", "re_re_1", "re_re_1", "re_re_1",
"re_re_1", "re_re_1", "re_re_1", "re_re_1", "re_re_1", "re_re_1",
"re_re_1", "re_re_1", "re_re_1", "re_re_1", "re_re_1", "re_re_1",
"re_re_1", "re_re_1", "re_re_1", "re_re_1", "re_re_1", "re_re_1",
"re_re_1", "re_re_1", "re_re_1", "re_re_1", "re_re_1", "re_re_1",
"re_re_1", "re_re_1", "re_re_1", "re_re_1", "re_re_1", "re_re_1",
"re_re_1", "re_re_1", "re_re_1", "re_re_1"), treatment = c(0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L)), row.names = c(NA, 50L), class = "data.frame")