Related
I have a panel dataset C_sf1 of 922 cities and 155 variables that looks like this:
structure(list(FID = c("AT001C1", "AT002C1"), E_pc_2000 = c(10.2,
10.2), E_pc_2001 = c(10.7, 10.7), E_pc_2002 = c(10.8, 10.8),
E_pc_2003 = c(11.5, 11.5), E_pc_2004 = c(11.4, 11.4), E_pc_2005 = c(11.5,
11.5), geometry = structure(list(structure(list(list(structure(c(16.4197,
16.4176, 16.4156, 16.4065, 16.4038, 16.4018, 16.4015, 16.4011,
16.398, 16.3951, 16.3941, 16.394, 16.3946, 16.3948, 16.3944,
16.3921, 16.3917, 16.3914, 16.3907, 16.3889, 16.3887, 16.3894,
16.3894, 16.3889, 16.3885, 16.3866, 16.3858, 16.3834, 16.3826,
16.3802, 16.3797, 16.3797, 16.3799, 16.3815, 16.3826, 16.3838,
16.3833, 16.3846, 16.3838, 16.3853, 16.3855, 16.3801, 16.3793,
16.3784, 16.378, 16.3773, 16.3763, 16.3755, 16.3751, 16.3752,
16.3783, 16.3782, 16.3775, 16.3763, 16.3723, 16.3694, 16.3691,
16.369, 16.3692, 16.3696, 16.3692, 16.3685, 16.3657, 16.3631,
16.3619, 16.3603, 16.3596, 16.3587, 16.3585, 16.3575, 16.3553,
16.3545, 16.3533, 16.3529, 16.3495, 16.3475, 16.3465, 16.3454,
16.3446, 16.3445, 16.3444, 16.3446, 16.3458, 16.348, 16.346,
16.3445, 16.3437, 16.3408, 16.3348, 16.3306, 16.329, 16.3236,
16.3212, 16.3193, 16.3179, 16.316, 16.3136, 16.3099, 16.3074,
16.3056, 16.3043, 16.3031, 16.3015, 16.2998, 16.2962, 16.2918,
16.2888, 16.2881, 16.2866, 16.2865, 16.2869, 16.2884, 16.2898,
16.2909, 16.2909, 16.291, 16.2902, 16.2856, 16.282, 16.2785,
16.2718, 16.2688, 16.2657, 16.2644, 16.2639, 16.2638, 16.2645,
16.2652, 16.268, 16.269, 16.2695, 16.2694, 16.2689, 16.2664,
16.2655, 16.2654, 16.2654, 16.265, 16.264, 16.2635, 16.2595,
16.2582, 16.2578, 16.2569, 16.2564, 16.2518, 16.2488, 16.2474,
16.2454, 16.2439, 16.2427, 16.2413, 16.2387, 16.2374, 16.2373,
16.2376, 16.2387, 16.2397, 16.2402, 16.2398, 16.2367, 16.2335,
16.2301, 16.2264, 16.2242, 16.2221, 16.2198, 16.2166, 16.2149,
16.2137, 16.2104, 16.2085, 16.2067, 16.2057, 16.2049, 16.2046,
16.2051, 16.2051, 16.2047, 16.2009, 16.1996, 16.1978, 16.1961,
16.1959, 16.196, 16.1957, 16.1951, 16.1951, 16.1958, 16.1969,
16.1982, 16.1984, 16.1978, 16.1967, 16.1955, 16.1938, 16.1901,
16.1848, 16.184, 16.184, 16.1839, 16.1841, 16.1844, 16.1847,
16.1854, 16.1856, 16.1869, 16.1886, 16.1932, 16.1951, 16.1965,
16.1975, 16.1978, 16.1965, 16.1961, 16.1962, 16.1968, 16.1971,
16.1977, 16.1994, 16.2008, 16.2035, 16.2057, 16.2067, 16.2077,
16.2085, 16.2073, 16.2065, 16.2063, 16.2054, 16.2042, 16.2027,
16.2037, 16.2086, 16.2088, 16.207, 16.2055, 16.2048, 16.2048,
16.2047, 16.2045, 16.204, 16.2035, 16.2016, 16.1997, 16.1986,
16.1975, 16.1961, 16.1947, 16.193, 16.1913, 16.1906, 16.19,
16.1894, 16.1885, 16.1875, 16.1867, 16.1858, 16.1851, 16.1844,
16.1842, 16.1845, 16.1848, 16.1856, 16.1864, 16.1888, 16.1911,
16.1936, 16.1946, 16.1956, 16.1961, 16.1967, 16.1969, 16.1971,
16.197, 16.1969, 16.1971, 16.198, 16.1988, 16.1992, 16.1996,
16.2006, 16.2016, 16.2032, 16.2048, 16.2105, 16.2114, 16.2181,
16.2238, 16.2239, 16.2226, 16.2197, 16.2193, 16.2179, 16.2154,
16.2137, 16.211, 16.2104, 16.2104, 16.2106, 16.2113, 16.2125,
16.2134, 16.216, 16.2184, 16.2194, 16.2199, 16.2205, 16.2208,
16.2215, 16.2241, 16.2314, 16.2348, 16.2369, 16.2376, 16.2364,
16.2385, 16.2383, 16.2312, 16.2265, 16.2186, 16.2169, 16.2154,
16.2151, 16.216, 16.2193, 16.2209, 16.2235, 16.2245, 16.2291,
16.2307, 16.2328, 16.2348, 16.2377, 16.2418, 16.2456, 16.2484,
16.2496, 16.25, 16.2512, 16.2522, 16.2527, 16.2534, 16.2559,
16.2569, 16.2594, 16.2599, 16.2668, 16.2717, 16.2729, 16.2818,
16.2834, 16.2846, 16.2915, 16.2975, 16.2985, 16.2985, 16.2986,
16.2987, 16.299, 16.2991, 16.2992, 16.2992, 16.2989, 16.3003,
16.3008, 16.3061, 16.3073, 16.3103, 16.3121, 16.3125, 16.3127,
16.314, 16.3146, 16.3178, 16.3183, 16.3187, 16.3197, 16.3206,
16.3222, 16.325, 16.3266, 16.3319, 16.3324, 16.3348, 16.3362,
16.3387, 16.3425, 16.3429, 16.3443, 16.3453, 16.3463, 16.349,
16.351, 16.3524, 16.3527, 16.3546, 16.3548, 16.3564, 16.3585,
16.3587, 16.3605, 16.3649, 16.3688, 16.3731, 16.3814, 16.3816,
16.3818, 16.3848, 16.3857, 16.3873, 16.3936, 16.3945, 16.3971,
16.398, 16.3994, 16.403, 16.405, 16.4081, 16.4088, 16.4099,
16.4109, 16.4115, 16.4119, 16.4119, 16.4121, 16.4123, 16.4139,
16.4186, 16.4192, 16.4208, 16.421, 16.4215, 16.4224, 16.4232,
16.424, 16.4251, 16.4268, 16.4272, 16.4296, 16.4333, 16.4351,
16.436, 16.4366, 16.4366, 16.4362, 16.4359, 16.4352, 16.4351,
16.4348, 16.4346, 16.4341, 16.4337, 16.4332, 16.4325, 16.4321,
16.4335, 16.4347, 16.437, 16.4437, 16.4457, 16.4492, 16.4517,
16.4536, 16.4551, 16.4645, 16.4658, 16.467, 16.4681, 16.4715,
16.4724, 16.4735, 16.4751, 16.4762, 16.477, 16.4774, 16.4776,
16.4795, 16.4811, 16.4845, 16.4884, 16.4926, 16.4943, 16.4976,
16.4992, 16.5036, 16.5087, 16.5105, 16.5114, 16.5118, 16.5142,
16.5196, 16.5284, 16.5309, 16.5331, 16.5378, 16.5401, 16.5404,
16.5426, 16.5473, 16.5506, 16.5529, 16.5623, 16.5684, 16.575,
16.576, 16.5762, 16.5767, 16.5769, 16.5771, 16.5775, 16.578,
16.578, 16.5782, 16.578, 16.578, 16.5772, 16.5767, 16.5757,
16.5747, 16.5742, 16.5743, 16.5743, 16.5744, 16.5742, 16.5737,
16.5733, 16.5725, 16.5714, 16.5706, 16.5695, 16.5689, 16.5674,
16.5666, 16.5655, 16.5647, 16.5643, 16.5636, 16.563, 16.5628,
16.5623, 16.5616, 16.5609, 16.5605, 16.5602, 16.5597, 16.559,
16.5581, 16.5575, 16.5567, 16.5561, 16.5553, 16.5548, 16.5544,
16.5543, 16.5545, 16.5549, 16.5551, 16.5553, 16.5552, 16.5548,
16.5542, 16.5535, 16.5529, 16.5524, 16.552, 16.5514, 16.5507,
16.5501, 16.5499, 16.5497, 16.5495, 16.5485, 16.5482, 16.5478,
16.5468, 16.5466, 16.5466, 16.5471, 16.5477, 16.5484, 16.5492,
16.5494, 16.5498, 16.5499, 16.5501, 16.5501, 16.55, 16.5498,
16.5497, 16.5495, 16.5492, 16.5486, 16.5478, 16.547, 16.5463,
16.5452, 16.5442, 16.543, 16.5424, 16.5416, 16.5411, 16.541,
16.5412, 16.5413, 16.5417, 16.5423, 16.5432, 16.5437, 16.5441,
16.5444, 16.5445, 16.5444, 16.5441, 16.5434, 16.5429, 16.5422,
16.5419, 16.5412, 16.5404, 16.5393, 16.5384, 16.5371, 16.5365,
16.5364, 16.5364, 16.5364, 16.5369, 16.5373, 16.538, 16.5387,
16.5392, 16.5396, 16.5398, 16.5396, 16.538, 16.5365, 16.5365,
16.5365, 16.5365, 16.5367, 16.5366, 16.5365, 16.5365, 16.5365,
16.5366, 16.5367, 16.5368, 16.5368, 16.5369, 16.537, 16.5371,
16.5372, 16.5373, 16.5373, 16.5374, 16.5376, 16.5377, 16.5378,
16.538, 16.5381, 16.5382, 16.5383, 16.5384, 16.5385, 16.5385,
16.5385, 16.5385, 16.5386, 16.5386, 16.5389, 16.5393, 16.5399,
16.54, 16.54, 16.5401, 16.5401, 16.5401, 16.5399, 16.5396,
16.5391, 16.5386, 16.5386, 16.5385, 16.5384, 16.5383, 16.5382,
16.538, 16.5378, 16.5375, 16.5373, 16.5371, 16.5369, 16.5367,
16.5364, 16.5362, 16.5363, 16.5364, 16.5365, 16.5375, 16.5377,
16.5385, 16.54, 16.5406, 16.5413, 16.5423, 16.543, 16.5439,
16.545, 16.5456, 16.5461, 16.5467, 16.5469, 16.5465, 16.5467,
16.548, 16.55, 16.5518, 16.5521, 16.552, 16.5502, 16.5488,
16.5484, 16.5482, 16.5479, 16.5474, 16.5463, 16.5445, 16.5423,
16.5417, 16.5414, 16.5424, 16.5419, 16.5416, 16.542, 16.5421,
16.5468, 16.5417, 16.5379, 16.5329, 16.5326, 16.5263, 16.5213,
16.5192, 16.5128, 16.514, 16.5087, 16.5074, 16.5079, 16.5083,
16.5088, 16.5095, 16.5098, 16.5116, 16.5129, 16.5137, 16.5145,
16.5142, 16.5125, 16.5128, 16.5094, 16.5041, 16.5021, 16.4983,
16.4951, 16.494, 16.4898, 16.484, 16.4813, 16.4809, 16.4809,
16.4809, 16.4819, 16.4824, 16.4833, 16.4835, 16.4835, 16.4834,
16.4836, 16.4829, 16.4813, 16.4809, 16.4816, 16.4772, 16.4748,
16.472, 16.4715, 16.4647, 16.4606, 16.4596, 16.452, 16.4483,
16.4402, 16.4401, 16.4401, 16.4402, 16.438, 16.4379, 16.438,
16.4363, 16.4339, 16.4336, 16.4299, 16.4271, 16.4247, 16.4226,
16.4197, 48.3222, 48.3221, 48.3218, 48.3196, 48.3191, 48.3185,
48.3173, 48.3171, 48.3177, 48.3184, 48.3188, 48.3191, 48.3201,
48.3215, 48.3217, 48.3219, 48.3216, 48.321, 48.3206, 48.3205,
48.3202, 48.3195, 48.3191, 48.3186, 48.3184, 48.3188, 48.3187,
48.3174, 48.3169, 48.3154, 48.3147, 48.3137, 48.3128, 48.31,
48.3085, 48.306, 48.3053, 48.3039, 48.3029, 48.302, 48.301,
48.3018, 48.3011, 48.2995, 48.2986, 48.2976, 48.2963, 48.2958,
48.2949, 48.2944, 48.2901, 48.2892, 48.2886, 48.2879, 48.2878,
48.2874, 48.2867, 48.2859, 48.285, 48.2839, 48.2835, 48.2833,
48.284, 48.2846, 48.2846, 48.2842, 48.2842, 48.2842, 48.2842,
48.2843, 48.2843, 48.2844, 48.2848, 48.2852, 48.2894, 48.2896,
48.2897, 48.2895, 48.2887, 48.2883, 48.288, 48.2878, 48.2867,
48.2845, 48.2838, 48.283, 48.2828, 48.2816, 48.2803, 48.2791,
48.2788, 48.2791, 48.2792, 48.279, 48.2785, 48.2776, 48.2754,
48.2737, 48.2728, 48.2725, 48.2725, 48.2726, 48.2732, 48.2734,
48.2728, 48.2717, 48.2706, 48.2701, 48.2679, 48.2668, 48.2663,
48.2658, 48.2654, 48.2646, 48.2641, 48.2634, 48.2628, 48.2607,
48.2595, 48.2586, 48.2564, 48.2559, 48.2556, 48.2552, 48.2548,
48.2544, 48.254, 48.2539, 48.2535, 48.2532, 48.2525, 48.252,
48.2514, 48.25, 48.2491, 48.2485, 48.2479, 48.2475, 48.2471,
48.2467, 48.2432, 48.2418, 48.2407, 48.2397, 48.2395, 48.2401,
48.24, 48.2398, 48.239, 48.2389, 48.2392, 48.2398, 48.2405,
48.2411, 48.242, 48.243, 48.2446, 48.2461, 48.2471, 48.2476,
48.2486, 48.2498, 48.2514, 48.2536, 48.255, 48.2567, 48.2578,
48.2585, 48.2587, 48.2593, 48.2619, 48.2626, 48.2626, 48.2623,
48.2613, 48.2593, 48.2563, 48.2556, 48.255, 48.2533, 48.252,
48.2488, 48.2451, 48.2436, 48.2418, 48.2378, 48.2357, 48.2347,
48.2339, 48.2332, 48.2327, 48.2324, 48.2316, 48.2307, 48.2291,
48.228, 48.2258, 48.2243, 48.2238, 48.2237, 48.2234, 48.2233,
48.2227, 48.2223, 48.222, 48.222, 48.222, 48.2224, 48.2226,
48.2225, 48.2222, 48.2218, 48.2214, 48.2203, 48.2196, 48.2185,
48.2165, 48.2161, 48.2151, 48.2137, 48.2131, 48.2125, 48.2117,
48.211, 48.2095, 48.2087, 48.208, 48.2073, 48.2072, 48.2067,
48.2062, 48.2061, 48.2045, 48.2031, 48.203, 48.1993, 48.1968,
48.1942, 48.1934, 48.1924, 48.1914, 48.1905, 48.1897, 48.1879,
48.1861, 48.1852, 48.1844, 48.183, 48.1816, 48.1802, 48.1789,
48.1783, 48.1777, 48.1762, 48.1746, 48.1731, 48.1727, 48.1724,
48.172, 48.1716, 48.1712, 48.1709, 48.1705, 48.1701, 48.1697,
48.1686, 48.1676, 48.1664, 48.1656, 48.1648, 48.164, 48.1632,
48.1621, 48.1609, 48.1601, 48.1593, 48.1581, 48.1569, 48.1556,
48.1556, 48.1555, 48.1558, 48.156, 48.1561, 48.1562, 48.1556,
48.1554, 48.1547, 48.1536, 48.1524, 48.1529, 48.1534, 48.1526,
48.1521, 48.1515, 48.151, 48.1501, 48.1494, 48.1485, 48.1479,
48.1472, 48.1463, 48.1453, 48.1417, 48.1406, 48.14, 48.1392,
48.1375, 48.1364, 48.1352, 48.1352, 48.1361, 48.1364, 48.1364,
48.1358, 48.1347, 48.1343, 48.1325, 48.1315, 48.1297, 48.1282,
48.1273, 48.1265, 48.1258, 48.1253, 48.1243, 48.1238, 48.1245,
48.1247, 48.1269, 48.1274, 48.1283, 48.129, 48.1297, 48.1301,
48.1308, 48.131, 48.1308, 48.1308, 48.1306, 48.1304, 48.1303,
48.1302, 48.1298, 48.1298, 48.1302, 48.1304, 48.1335, 48.1333,
48.1332, 48.1308, 48.1307, 48.1305, 48.1296, 48.1286, 48.1287,
48.1283, 48.128, 48.1279, 48.1275, 48.1273, 48.1269, 48.1266,
48.1262, 48.1251, 48.1245, 48.1217, 48.121, 48.1197, 48.1197,
48.1217, 48.1231, 48.1256, 48.1269, 48.1261, 48.1269, 48.1276,
48.1292, 48.1307, 48.1333, 48.138, 48.1376, 48.1368, 48.1367,
48.1362, 48.1359, 48.1353, 48.1345, 48.1344, 48.134, 48.1337,
48.1336, 48.1332, 48.1328, 48.1322, 48.1321, 48.1309, 48.1309,
48.1304, 48.1296, 48.1296, 48.1293, 48.1287, 48.1278, 48.1272,
48.1265, 48.1264, 48.1264, 48.126, 48.1259, 48.1257, 48.125,
48.1249, 48.1247, 48.1244, 48.1235, 48.122, 48.1209, 48.1189,
48.1187, 48.1186, 48.1187, 48.1189, 48.1194, 48.1202, 48.1204,
48.1206, 48.1207, 48.1209, 48.1209, 48.1216, 48.1217, 48.122,
48.1226, 48.1229, 48.1229, 48.1225, 48.122, 48.1219, 48.1212,
48.1204, 48.12, 48.1199, 48.1202, 48.1204, 48.1213, 48.1219,
48.123, 48.1232, 48.1237, 48.124, 48.1256, 48.1278, 48.1303,
48.131, 48.1314, 48.1379, 48.138, 48.1381, 48.1384, 48.1384,
48.1387, 48.1389, 48.1393, 48.14, 48.1462, 48.147, 48.1478,
48.1484, 48.1507, 48.1514, 48.1522, 48.1533, 48.1541, 48.1548,
48.1556, 48.1562, 48.157, 48.1568, 48.1561, 48.1554, 48.1586,
48.1581, 48.1572, 48.1568, 48.157, 48.1577, 48.1581, 48.1583,
48.1586, 48.1591, 48.1572, 48.1538, 48.1525, 48.151, 48.147,
48.1452, 48.145, 48.1438, 48.1422, 48.1413, 48.1407, 48.139,
48.1377, 48.1361, 48.1363, 48.1369, 48.1386, 48.1393, 48.1402,
48.1418, 48.1428, 48.1433, 48.144, 48.146, 48.1476, 48.1506,
48.152, 48.1538, 48.1554, 48.1563, 48.1573, 48.1582, 48.159,
48.1597, 48.1604, 48.161, 48.1616, 48.1623, 48.1627, 48.163,
48.1631, 48.1631, 48.1631, 48.1629, 48.1628, 48.1625, 48.1621,
48.1616, 48.1612, 48.161, 48.1609, 48.1609, 48.1613, 48.1616,
48.1622, 48.1629, 48.1638, 48.1641, 48.1643, 48.1644, 48.1644,
48.1642, 48.1639, 48.1636, 48.1633, 48.163, 48.1627, 48.1624,
48.1622, 48.162, 48.1621, 48.1621, 48.1623, 48.1624, 48.1626,
48.163, 48.1635, 48.1642, 48.165, 48.166, 48.1665, 48.1675,
48.1681, 48.1689, 48.1704, 48.1709, 48.1713, 48.1714, 48.1716,
48.1718, 48.172, 48.1723, 48.1727, 48.1732, 48.1737, 48.1742,
48.1749, 48.1756, 48.1761, 48.1767, 48.177, 48.1773, 48.1775,
48.1775, 48.1775, 48.1773, 48.1771, 48.1771, 48.1771, 48.1771,
48.1768, 48.1765, 48.1762, 48.176, 48.1758, 48.1756, 48.1752,
48.1748, 48.1745, 48.1739, 48.1735, 48.1733, 48.1731, 48.1732,
48.1733, 48.1735, 48.1734, 48.1732, 48.173, 48.1727, 48.1727,
48.173, 48.1732, 48.1734, 48.1736, 48.1739, 48.174, 48.174,
48.174, 48.174, 48.1738, 48.1738, 48.1741, 48.1745, 48.1757,
48.1771, 48.1774, 48.1778, 48.179, 48.1803, 48.1814, 48.182,
48.1826, 48.1828, 48.1831, 48.1833, 48.1835, 48.1837, 48.1839,
48.1842, 48.1846, 48.1851, 48.1855, 48.1856, 48.186, 48.1864,
48.1867, 48.1869, 48.1872, 48.1874, 48.1876, 48.1879, 48.188,
48.1882, 48.1884, 48.1887, 48.1892, 48.1895, 48.1897, 48.191,
48.1915, 48.1923, 48.1924, 48.1927, 48.193, 48.1935, 48.1938,
48.1942, 48.1947, 48.1954, 48.1961, 48.1963, 48.1965, 48.1967,
48.1972, 48.1976, 48.1981, 48.1984, 48.1988, 48.1991, 48.1993,
48.1995, 48.1997, 48.2001, 48.2003, 48.201, 48.2014, 48.202,
48.2039, 48.2043, 48.2055, 48.2075, 48.2082, 48.2089, 48.2099,
48.2109, 48.2122, 48.2137, 48.2148, 48.2161, 48.2184, 48.2196,
48.2215, 48.2233, 48.2283, 48.2333, 48.2367, 48.2377, 48.2388,
48.2395, 48.2401, 48.2404, 48.2416, 48.2423, 48.2426, 48.2428,
48.2429, 48.2429, 48.243, 48.2437, 48.2482, 48.2487, 48.2497,
48.2511, 48.2513, 48.2631, 48.2639, 48.2642, 48.2629, 48.262,
48.2632, 48.2645, 48.2669, 48.2679, 48.2712, 48.2734, 48.2737,
48.2752, 48.2761, 48.277, 48.2783, 48.2787, 48.2812, 48.2821,
48.2839, 48.2855, 48.2859, 48.2866, 48.2873, 48.288, 48.2896,
48.2904, 48.2914, 48.2924, 48.2926, 48.2931, 48.294, 48.2944,
48.2938, 48.2926, 48.2922, 48.2898, 48.2876, 48.286, 48.2849,
48.2848, 48.2794, 48.2789, 48.2776, 48.2757, 48.2749, 48.273,
48.2737, 48.2751, 48.2769, 48.2772, 48.2808, 48.2829, 48.283,
48.2864, 48.2882, 48.2919, 48.296, 48.2989, 48.3031, 48.3099,
48.3126, 48.3164, 48.317, 48.3174, 48.3176, 48.3193, 48.3204,
48.3212, 48.3217, 48.3222), .Dim = c(807L, 2L)))), class = c("XY",
"MULTIPOLYGON", "sfg")), structure(list(list(structure(c(15.4008,
15.3968, 15.3961, 15.3958, 15.3971, 15.3978, 15.3973, 15.3957,
15.3919, 15.3911, 15.3904, 15.3823, 15.3817, 15.3814, 15.3781,
15.3778, 15.3758, 15.3736, 15.3723, 15.3696, 15.367, 15.3668,
15.3661, 15.3619, 15.3585, 15.3583, 15.3582, 15.358, 15.3584,
15.3619, 15.364, 15.3638, 15.3634, 15.3626, 15.3522, 15.3515,
15.3513, 15.3512, 15.3517, 15.3526, 15.3533, 15.3581, 15.3594,
15.3602, 15.3607, 15.3642, 15.3651, 15.3652, 15.3652, 15.3652,
15.3675, 15.3718, 15.3728, 15.3739, 15.3723, 15.3682, 15.3682,
15.3684, 15.374, 15.3772, 15.3795, 15.3804, 15.3808, 15.3805,
15.3794, 15.3761, 15.3761, 15.3753, 15.3743, 15.3714, 15.3705,
15.3693, 15.3682, 15.3674, 15.3664, 15.3654, 15.364, 15.3618,
15.3612, 15.3598, 15.3592, 15.3592, 15.3609, 15.3613, 15.366,
15.3662, 15.3684, 15.3687, 15.3687, 15.3673, 15.3674, 15.3675,
15.3681, 15.3693, 15.3718, 15.3764, 15.378, 15.3858, 15.3897,
15.3912, 15.3924, 15.3941, 15.3956, 15.3986, 15.4033, 15.4076,
15.4096, 15.4102, 15.422, 15.4228, 15.434, 15.4373, 15.4396,
15.4404, 15.4409, 15.442, 15.4422, 15.4429, 15.4476, 15.4466,
15.4484, 15.4494, 15.4509, 15.4531, 15.4534, 15.457, 15.4575,
15.4577, 15.4585, 15.4594, 15.4605, 15.4616, 15.4639, 15.4647,
15.4688, 15.4734, 15.4777, 15.48, 15.4846, 15.4862, 15.4882,
15.4888, 15.4875, 15.4884, 15.489, 15.4885, 15.487, 15.4863,
15.486, 15.4836, 15.4833, 15.4829, 15.4834, 15.4849, 15.4867,
15.4865, 15.4846, 15.4901, 15.4934, 15.4982, 15.4982, 15.4981,
15.4979, 15.4995, 15.5034, 15.5001, 15.4988, 15.4981, 15.498,
15.4989, 15.5002, 15.5015, 15.5034, 15.5049, 15.5057, 15.5069,
15.5077, 15.5081, 15.5083, 15.5082, 15.507, 15.5065, 15.5062,
15.506, 15.506, 15.5062, 15.5068, 15.5093, 15.5108, 15.5119,
15.5123, 15.513, 15.5184, 15.5206, 15.522, 15.5255, 15.5264,
15.5271, 15.5308, 15.5336, 15.5327, 15.5307, 15.5287, 15.5277,
15.5264, 15.526, 15.526, 15.5283, 15.5302, 15.5302, 15.5303,
15.5303, 15.5266, 15.5201, 15.5147, 15.5136, 15.5095, 15.5102,
15.5118, 15.5132, 15.5145, 15.5157, 15.5176, 15.5182, 15.5184,
15.5196, 15.5199, 15.5194, 15.5186, 15.5159, 15.5148, 15.5118,
15.5099, 15.5054, 15.5023, 15.4991, 15.4953, 15.4922, 15.4911,
15.4838, 15.4822, 15.4813, 15.48, 15.4792, 15.4783, 15.4771,
15.4682, 15.4651, 15.4637, 15.4624, 15.4596, 15.4459, 15.4423,
15.4369, 15.4298, 15.4268, 15.4248, 15.4223, 15.4189, 15.4173,
15.4137, 15.4108, 15.4089, 15.4073, 15.4044, 15.4008, 47.1349,
47.1327, 47.1323, 47.1322, 47.1287, 47.1254, 47.123, 47.1199,
47.1151, 47.1147, 47.1145, 47.1141, 47.1144, 47.1147, 47.1212,
47.1213, 47.1206, 47.1198, 47.1194, 47.1194, 47.1195, 47.1195,
47.1195, 47.1199, 47.1202, 47.1197, 47.1195, 47.1187, 47.1171,
47.1155, 47.1113, 47.1109, 47.1106, 47.1104, 47.108, 47.1076,
47.1071, 47.1032, 47.1024, 47.1021, 47.102, 47.1028, 47.1028,
47.1026, 47.1022, 47.0977, 47.0962, 47.0948, 47.0946, 47.0914,
47.0877, 47.0878, 47.0876, 47.0872, 47.085, 47.0806, 47.0804,
47.0803, 47.0782, 47.0788, 47.0764, 47.0754, 47.0717, 47.0704,
47.0706, 47.0712, 47.0681, 47.0657, 47.0642, 47.0628, 47.0619,
47.0582, 47.054, 47.0534, 47.0529, 47.0528, 47.0529, 47.0534,
47.0533, 47.0525, 47.0515, 47.0513, 47.0492, 47.0486, 47.0422,
47.0421, 47.0411, 47.0402, 47.039, 47.0326, 47.0322, 47.031,
47.0294, 47.028, 47.0257, 47.0231, 47.0221, 47.0198, 47.0201,
47.0187, 47.0176, 47.0173, 47.0164, 47.0163, 47.0163, 47.0163,
47.0166, 47.0167, 47.0183, 47.0185, 47.0211, 47.0219, 47.0183,
47.0169, 47.0161, 47.0163, 47.0163, 47.0153, 47.0166, 47.0177,
47.0185, 47.0176, 47.0177, 47.0189, 47.019, 47.0182, 47.0171,
47.0166, 47.015, 47.013, 47.0125, 47.012, 47.0129, 47.0132,
47.0136, 47.0153, 47.0167, 47.0172, 47.0193, 47.0197, 47.0203,
47.0205, 47.0213, 47.0217, 47.0219, 47.0227, 47.0219, 47.0217,
47.0223, 47.0246, 47.0249, 47.0255, 47.0256, 47.026, 47.0265,
47.0268, 47.0308, 47.0328, 47.0339, 47.0363, 47.0366, 47.0374,
47.0378, 47.0402, 47.0412, 47.0456, 47.0473, 47.0491, 47.0525,
47.0533, 47.0543, 47.057, 47.0589, 47.0608, 47.0622, 47.066,
47.0673, 47.068, 47.0689, 47.07, 47.074, 47.075, 47.0756,
47.0767, 47.0778, 47.0786, 47.079, 47.0807, 47.082, 47.0837,
47.0863, 47.0903, 47.092, 47.0931, 47.0941, 47.0969, 47.0973,
47.0974, 47.0973, 47.0984, 47.0998, 47.1008, 47.1018, 47.1023,
47.1034, 47.1043, 47.1052, 47.1063, 47.1072, 47.1075, 47.1103,
47.1107, 47.1095, 47.1075, 47.1069, 47.1068, 47.1074, 47.1081,
47.1095, 47.1118, 47.1144, 47.1157, 47.1171, 47.1176, 47.1177,
47.1189, 47.1201, 47.1213, 47.1234, 47.1263, 47.127, 47.1289,
47.1292, 47.1289, 47.1288, 47.1286, 47.1245, 47.1271, 47.1272,
47.1295, 47.1295, 47.1292, 47.1285, 47.1276, 47.1269, 47.1269,
47.1285, 47.1305, 47.1309, 47.1309, 47.1305, 47.1278, 47.127,
47.1276, 47.1279, 47.1281, 47.128, 47.128, 47.1277, 47.128,
47.1309, 47.1332, 47.1338, 47.1339, 47.1336, 47.1349), .Dim = c(266L,
2L)))), class = c("XY", "MULTIPOLYGON", "sfg"))), class = c("sfc_MULTIPOLYGON",
"sfc"), precision = 0, bbox = structure(c(xmin = 15.3512,
ymin = 47.012, xmax = 16.5782, ymax = 48.3222), class = "bbox"), crs = structure(list(
input = "WGS 84", wkt = "GEOGCRS[\"WGS 84\",\n DATUM[\"World Geodetic System 1984\",\n ELLIPSOID[\"WGS 84\",6378137,298.257223563,\n LENGTHUNIT[\"metre\",1]]],\n PRIMEM[\"Greenwich\",0,\n ANGLEUNIT[\"degree\",0.0174532925199433]],\n CS[ellipsoidal,2],\n AXIS[\"latitude\",north,\n ORDER[1],\n ANGLEUNIT[\"degree\",0.0174532925199433]],\n AXIS[\"longitude\",east,\n ORDER[2],\n ANGLEUNIT[\"degree\",0.0174532925199433]],\n ID[\"EPSG\",4326]]"), class = "crs"), n_empty = 0L)), row.names = 1:2, class = c("sf",
"data.frame"), sf_column = "geometry", agr = structure(c(FID = NA_integer_,
E_pc_2000 = NA_integer_, E_pc_2001 = NA_integer_, E_pc_2002 = NA_integer_,
E_pc_2003 = NA_integer_, E_pc_2004 = NA_integer_, E_pc_2005 = NA_integer_
), .Label = c("constant", "aggregate", "identity"), class = "factor"))
My task is to calculate the spatial weights matrix and the spatial lags associated with the future dependent variable of my model: log(E).
I built the spatial weights matrix for the wide dataframe.
library(spdep)
coord <- st_coordinates(st_centroid(C_sf1))
nb2.knn1 <- knn2nb(knearneigh(coord,k=1,longlat=T))
all.linkedT <- max(unlist(nbdists(nb2.knn1, coord,longlat=T)))
all.linkedT
nb3 <- dnearneigh(coord, 0, all.linkedT,row.names=C_sf$FDI,longlat=TRUE)
n.nb3 <- nbdists(nb3, coord,longlat=TRUE);n.nb3
n.nb3.gl <- lapply(n.nb3, function(x) 1/x); n.nb3.gl
nb3.listw.gl <- nb2listw(nb3, glist=n.nb3.gl,zero.policy=F)
My problem is that I want to shift to a long format (17k+ observations), keeping the spatial weight matrix. The long format would allow me to have a variable called E (and not E_2001, E_2002, etc.) from which I can calculate the logarithm.
I tried to use the same code for the dataframe C_sf_long (obtained with reshape), but R returns me a spatial weight matrix of 0 distances. The reason is that in the panel dataset, the same coordinates are repeated for each time unit, so distance and neighbour functions only recognize a city as the neighbour of itself.
The code for the spatial lags is the following:
lag.E<- lag.listw(nb2listw(nb3),C_sf1$E)
but, as you can see, C_sf1 does not have the variable E, because it is still in the wide format.
So, I want C_sf1 to be in long format, in order to compute the spatial lags, but still I want to keep the spatial weight matrix computed with C_sf1 in wide format.
I cant run my MCa with a ind.sup, I dont know why. Here's my code :
library(FactoMineR)
data_ACM <- data[80000:81000,]
res.AMC <- MCA(data_ACM, ncp = 6, graph = F, ind.sup = 700)
Error in eigen(crossprod(t(X), t(X)), symmetric = TRUE) :
infinite or missing values in 'x'
dpu(data[1:10,]) returns this :
"SOHDZET", "SOHDZF", "SOHDZFT", "SOKFXD", "SOKFXF", "SOM31BL4B1D",
"SONFZD", "SONFZF", "SOVJZC", "SOVJZE", "SOVJZET", "SOVJZF",
"SOVJZFT", "SOVJZJ", "SOZAAZ", "SP01", "SP21", "SP308", "SP41",
"SP4DWC", "SP61", "SP81", "SR11", "SR1F0H", "SR1G0H", "SR1H06",
"SR1J0H", "SR1J22", "SR1J32", "SR1L22", "SR1L32", "SR1M22",
"SR270E", "SR31", "SR4DWC", "SR51", "SR81", "SR8J32", "SR8L32",
"SREU0H", "SRH1ELC", "SRH1ELCD", "SS11", "SS31", "SS51",
"SS55", "SS71", "SS91", "ST12", "ST72", "ST785E", "STA",
"STADRX01AG", "STADRX01SG", "STAFNX01AG", "STAFNX01SG", "STAGZX01SG",
"STAHLX01AG", "STAHLX01SG", "STAHUX01AG", "STAHUX01SG", "STAV",
"STB", "STBV", "STC", "STCDCA", "STCV", "SU31", "SU35", "SU55",
"SU55A", "SV10", "SV10EM", "SV1V12", "SV31", "SVA", "SVAV",
"SW11", "SW31", "SW51", "SW52X400", "SW52X700", "SW53K2",
"SW53K200", "SW53K9", "SW53K900", "SW5441", "SW5449", "SW58K2",
"SW58K200", "SW58K3", "SW58K7", "SW58K9", "SW58K900", "SW5952",
"SW5953", "SW5957", "SW5959", "SW59G200", "SW59G300", "SW59G700",
"SW59G900", "SW6142", "SW6149", "SW61P2", "SW61P200", "SW61P9",
"SW61P900", "SW61S200", "SW61S900", "SW6542", "SW6549", "SW65P2",
"SW65P200", "SW65P9", "SW65P900", "SW65S200", "SW65S900",
"SW6942", "SW71", "SW7140", "SW7141", "SW7144", "SW7146",
"SW7362", "SW743200", "SW743900", "SW7462", "SW7469", "SW793200",
"SW793700", "SW793900", "SW7962", "SW7964", "SW7967", "SW7969",
"SW8042", "SW8142", "SX11", "SX31", "SX71", "SXBHW6", "SXBHY6",
"SXHMP6", "SXHMZ6", "SXHNZ6", "SXHNZT", "SY11", "SY31", "SY33",
"SY34", "SY53", "SY54", "SY71", "SYN1E", "SZ58K3", "SZ58K7",
"SZ5953", "SZ5957", "SZ59G300", "SZ59G700", "SZ7144", "SZ7146",
"SZ793700", "SZ7964", "SZ7967", "SZ799400", "T11", "T12",
"T1X000", "T1X004", "T1X300", "T1X305", "T1X805", "T1XE05",
"T1XF05", "T27WWT270", "T27WWT271", "T27ZRT270", "T27ZRT271",
"T2DG1", "T2DH1", "T2DH2", "T2DJ1", "T2DJ2", "T2W300", "T2W305",
"T2W800", "T2W805", "T2WE05", "T2WF05", "T31DD03", "T31EE01",
"T32AA", "T32AA01", "T32AA02", "T32AA03", "T32AA04", "T32BB",
"T32BB01", "T32BB02", "T32BB05", "T32BB06", "T32CC", "T32CC01",
"T32CC02", "T32CC03", "T32CC04", "T32DD", "T32EE", "T3408",
"T3419", "T3435", "T3448", "T3467", "T3DB1", "T3DE1", "T3DE2",
"T3LA1", "T3LA2", "T3X000", "T3X100", "T3X300", "T3X305",
"T3X800", "T3X805", "T3XE05", "T3XF05", "T4BBD", "T4BBG",
"T4BCD", "T4BCE", "T4BCF", "T4DCA", "T4DCD", "T4DCT", "T4DCX",
"T4DCY", "T4X000", "T4X100", "T4X305", "T4X405", "T4X800",
"T4X805", "T4XE05", "T4XF05", "T4XG05", "T5408", "T5419",
"T5435", "T5448", "T5808", "T5819", "T5835", "T5848", "T5X000",
"T5X200", "T5X205", "T5X300", "T5X305", "T5X405", "T5X800",
"T5X805", "T5XF05", "T5XG05", "T6EB1", "T6EB4", "T6EB6",
"T6EB7", "T6ED5", "T6EE2", "T6MB1", "T6MB6", "T6MB7", "T6MC1",
"T6MD5", "T72", "T7JA1", "T7W200", "T7W205", "T7W300", "T7W405",
"T7WG05", "T824833", "T82W433", "T82W434", "T8344", "T8358A",
"T8378A", "T8380", "T8390A", "T8394", "T83T4", "T844433",
"T844434", "T844533", "T844544", "T845034", "T845043", "T845044",
"T8458", "T847433", "T847443", "T847833", "T8478A", "T8480",
"T8488", "T849033", "T849043", "T8490A", "T8494", "T8494B43",
"T8494B44", "T84T4", "T84T443", "T84T444", "T84T543", "T84T544",
"T84T833", "T84T843", "T855", "T857", "T85T4", "T85T8", "T9408",
"T9419", "T9435", "T9448", "T9508", "T9519", "T9535", "T9548",
"T9808", "T9819", "T9835", "T9848", "T9PG2", "TA", "TA12",
"TA12F2", "TA12J5", "TA12L2", "TA12L5", "TA13008", "TA13035",
"TA13068", "TA1308", "TA1335", "TA1368", "TA14", "TA16008",
"TA16035", "TA16068", "TA1608", "TA1635", "TA1668", "TA2",
"TA22", "TA23H", "TA23MB", "TA23MQG", "TA2422", "TA2422A",
"TA28MB", "TA32408", "TA32419", "TA32468", "TA40", "TA40AKS",
"TA40AMS", "TA43035", "TA43068", "TA4308", "TA4319", "TA4335",
"TA4368", "TA4408", "TA4419", "TA4435", "TA4468", "TA60AMS",
"TA60BLMB", "TA72L2", "TA72L5", "TA8408", "TA8419", "TA8435",
"TA8468", "TA9408", "TA9419", "TA9435", "TA9468", "TA94E08",
"TA94E19", "TA94E35", "TA94E68", "TAA408", "TAA419", "TAA435",
"TAA468", "TAA8408", "TAA9408", "TAA94E08", "TAAA408", "TAAP11",
"TAAW10", "TABCAS", "TABCDS", "TAC408", "TAC419", "TAC468",
"TAD308", "TAD3G08", "TAF408", "TAF419", "TAF435", "TAF468",
"TAJ3G008", "TAJ3G08", "TAK1367", "TAK4367", "TAKD467", "TAKU467",
"TAKV467", "TAM3019", "TAM3035", "TAM3068", "TAM319", "TAM335",
"TAM368", "TAM6019", "TAM6035", "TAM6068", "TAM619", "TAM635",
"TAM668", "TAN308", "TAN419", "TAN435", "TAN468", "TAN619",
"TAN635", "TAN668", "TANT30", "TAS335", "TAS3G08", "TATB",
"TATF", "TATHA", "TATHB", "TATJA", "TATJB", "TAU8470", "TAUA470",
"TAW4035", "TAW4068", "TAW6035", "TAW6068", "TAX100", "TAX300",
"TAZ100", "TAZ300", "TB", "TB2", "TB2412", "TB367", "TBAABS",
"TBABBS", "TBAP12", "TBAV10", "TBBCBS", "TBC1NPLY", "TBC1NRLY",
"TBCAAS", "TBCADS", "TBCAES", "TBCBAS", "TBCBDS", "TBCBES",
"TBDCAS", "TBDCDS", "TBDCES", "TBE1OVM", "TBE1OVN", "TBE2PZN",
"TBE2UZN", "TBE4TWN", "TBE4TYN", "TBE4TYNC", "TBE5IWN", "TBE5IWNC",
"TBE5TWM", "TBE5TWN", "TBE5TWNC", "TBE5TYN", "TBGMSLV", "TBGMTLVI",
"TBGNSLV", "TBGNSLVC", "TBGNSLY", "TBGNTLVI", "TBIF5P11M51AZ1",
"TBIF5P21M52CZ1", "TBNT30", "TBUR20", "TBX200", "TBX205",
"TBX300", "TBX305", "TBX405", "TBXE05", "TBXG05", "TBZ200",
"TBZ205", "TBZ300", "TBZ305", "TBZ405", "TBZ800", "TBZE05",
"TBZG05", "TC2", "TCAP11", "TCAP12", "TCAV10", "TCTG05",
"TCY405", "TCYG05", "TD2", "TD367", "TD4HSJ", "TD4HTH", "TDABKS",
"TDABTS", "TDABUS", "TDAP11", "TDAV10", "TDAW10", "TDAW10M",
"TDAW10U", "TDFBAS", "TDGCKS", "TDGCTS", "TDGCUS", "TDHCAS",
"TDSY61", "TDUHZJ", "TDX8ZA", "TDXFVJ", "TE2", "TE51", "TEABKL",
"TEABTL", "TEAP12", "TEBBAL", "TECCKL", "TECCTL", "TEDCAL",
"TENT30", "TESY61", "TF08II51", "TF08J551", "TF08M351", "TF35II51",
"TF35J551", "TF35M351", "TF408", "TF419", "TF435", "TF448",
"TF48II51", "TF48J551", "TF48M351", "TFABAL", "TFABKL", "TFABSL",
"TFABTL", "TFABUL", "TFAP12", "TFAY10", "TFBCAL", "TFBCKL",
"TFBCSL", "TFBCTL", "TFBCUL", "TG08GCA1", "TG08GCS1", "TG08ICA1",
"TG08ICS1", "TG08M351", "TG08M3A1", "TG2", "TG35GCA1", "TG35GCS1",
"TG35ICA1", "TG35ICS1", "TG35M351", "TG35M3A1", "TG408",
"TG419", "TG435", "TG448", "TG48GCA1", "TG48GCS1", "TG48ICA1",
"TG48ICS1", "TG48M351", "TG48M3A1", "TGABAL", "TGABKL", "TGABSL",
"TGABTE", "TGABTL", "TGBN", "TGCCAL", "TGCCKL", "TGCCSL",
"TGCCTE", "TGCCTL", "TGFCAL", "TGFCTL", "TGFCUL", "TH308",
"TH348", "TJAC9S", "TJAN91", "TJAN9S", "TJBBAL", "TJBBTE",
"TJBBTL", "TJBBUL", "TJBCAL", "TJBCTE", "TJBCTL", "TJBCUL",
"TK308", "TK348", "TKACAL", "TKACTE", "TKACTL", "TKACUL",
"TLB310", "TLEF5D14A63ZZ1", "TLEF5D24M62ZZ1", "TLEF5D31M62ZZ1",
"TLEF5D51M66ZZ1", "TLEF5D61D7", "TLEF5P31M64ZZ1", "TLEF5P41D71ZZ1",
"TLEF5P44D71ZZ1", "TM4BUL", "TM4DBC", "TP260", "TPAD", "TPADJ",
"TPAE", "TPAEJ", "TQF8D62M53AZ1", "TQF8D82A51AZ1", "TR160",
"TR160T", "TR467", "TR567", "TR7969", "TR9145", "TR9245",
"TS08ICS1", "TS08MCS1", "TS31", "TS48ICS1", "TS6142", "TS6149",
"TS61E2", "TS61P2", "TS61P200", "TS61P9", "TS61P900", "TS61S200",
"TS61S900", "TS65E2", "TS65P2", "TS65P200", "TS65P9", "TS65P900",
"TS65S200", "TS65S900", "TS7202", "TS7206", "TS7432", "TS743200",
"TS7439", "TS743900", "TS7462", "TS7469", "TS7932", "TS793200",
"TS793900", "TS7962", "TS7969", "TS90C5", "TS90K5", "TS9145",
"TS91J500", "TS91K5", "TS91K500", "TS9245", "TS92K5", "TS92K500",
"TS94K2", "TS94K5", "TS97C2", "TS97C5", "TT132", "TT408",
"TT419", "TT435", "TT448", "TT508", "TT519", "TT535", "TT548",
"TT808", "TT819", "TT835", "TT848", "TTFN44", "TU308", "TU319",
"TU348", "TU408", "TU419", "TU435", "TU448", "TU467", "TU508",
"TU519", "TU535", "TU548", "TU567", "TU808", "TU819", "TU835",
"TU848", "TVUF", "TVUR20", "TVUR20U", "TW308", "TW319", "TW348",
"TW3T08", "TW3T19", "TW3T35", "TW3T48", "TWUR20", "TX31",
"TX71", "TY260", "TY260T", "U11", "U11T", "U2", "U2MAC",
"U2NAC", "U6UA", "U6UB", "U6UC", "U6UE", "U6UF", "U6UG",
"U6UJ", "U6UJT", "U6UK", "U6UR", "U6UR9", "U6UT", "U6UT9",
"U6UU", "U6UU9", "U6UW", "U6UWT", "U857", "U858", "U859533",
"U85T", "U9C1G6", "U9C1GH", "U9C1K6", "U9C1KH", "U9C2G6",
"U9C2K6", "U9CYG6", "U9CYK6", "U9VCK5", "UA5FV81P", "UA5FWC",
"UA5FWC1", "UA5FXH1P", "UA6FYC", "UA71", "UA9HR8", "UA9HR82PS",
"UA9HR8P", "UA9HR8PS", "UA9HZC", "UA9HZC1", "UA9HZH1P", "UA9HZHP",
"UARFJHP", "UARHB8P", "UARHE8", "UARHHA", "UARHJH1P", "UARHJHP",
"UARHRJ", "UB11", "UB31", "UB51", "UB53", "UB71", "UB91",
"UC11", "UC31", "UC71", "UC91", "UD11", "UD1Y11", "UD51",
"UD5FS0", "UD5FV81P", "UD5FWC", "UD5FXH1P", "UD5FXHP", "UD6FYC",
"UD71", "UD91", "UD9HR8", "UD9HR82PS", "UD9HR8P", "UD9HR8PS",
"UD9HZC", "UD9HZC1", "UD9HZC1CU1", "UD9HZH1P", "UD9HZH1PCU1",
"UD9HZHP", "UDC1G6", "UDC1K6", "UDC2G6", "UDC2K6", "UDC3G6",
"UDC3K6", "UDCAG5", "UDCAK5", "UDCCG5", "UDCCK5", "UDCEG5",
"UDCEK5", "UDCGG5", "UDCGK5", "UDCJG5", "UDCJK5", "UDCMG5",
"UDCMK5", "UDCNG5", "UDCNK5", "UDCSG6", "UDCSK6", "UDCUG5",
"UDCUG6", "UDCUK5", "UDCUK6", "UDCVG5", "UDCVK5", "UDCYG6",
"UDCYK6", "UDRFJHP", "UDRHB8P", "UDRHE8", "UDRHHA", "UDRHJH1P",
"UDRHJHP", "UDRHRJ", "UE11", "UE31", "UE51", "UE9HR8", "UE9HR8P",
"UE9HZC1", "UE9HZH1P", "UF11", "UF1Y32", "UF31", "UF51",
"UF8U52", "UF8UL2", "UF91", "UG31", "UG51", "UH11", "UH31",
"UH51", "UK11", "UK31", "UK51", "UL91", "UM11", "UM51", "UM71",
"UM91", "UMC7D2", "UN10", "UN13", "UN1A22", "UN71", "UN8B42",
"UN8D32", "UN8F42", "UN91", "UP11", "UR11", "UR31", "UR51",
"UR91", "URMD21", "US31", "US91", "USD135", "USD145", "USD1F5",
"USD1K5", "USD1W5", "USDBL5", "USDUK5", "USDUW5", "UT31",
"UT71", "UU11", "UU51", "UU6M", "UU71", "UU91", "UV51", "UV71",
"UV91", "UW71", "UW91", "UX11", "UX31", "UX51", "UX71", "UX91",
"UXC1", "UY11", "UY31", "UY51", "UY71", "UY91", "UZ10BC",
"UZ68BC", "UZA2BC", "UZA8BC", "UZBABD", "UZJ100LGNAEKW",
"V1DKS", "V1DVS", "V1JKS", "V1JVS", "V23CGRHE", "V23WGNXE",
"V23WGNXES", "V24WGNXF", "V24WGNXFS", "V24WNDF", "V24WNDFS",
"V24WNHF", "V24WNHFS", "V25WGNX", "V25WGRX", "V26WGNX", "V26WGNXS",
"V2DAC", "V2JAC", "V2W200", "V2W300", "V36MM01", "V36NN01",
"V36PP01", "V36RR01", "V36SS01", "V36TT01", "V37CC03", "V37DD03",
"V37EE01", "V37FF01", "V37GG01", "V3X208", "V3X308", "V3XG08",
"V43WGRXE", "V44WGNXF", "V44WGNXFS", "V44WNHF", "V45WGRX",
"V46WGNX", "V46WGRX", "V46WNH", "V46WNHS", "V4X208", "V4X300",
"V4X308", "V4X408", "V4XE08", "V4XG08", "V4Y208", "V4Y408",
"V4YG08", "V521L22HCR167", "V521SLLDA5865", "V521SLLDAR165",
"V5W208", "V5W408", "V5WG08", "V6419", "V6435", "V64WMN",
"V7419", "V7435", "V9X300", "VA", "VA51", "VA71", "VA91",
"VABHSH", "VABHVB", "VABHXH", "VAJ", "VAJI", "VAJIA", "VAL",
"VAM", "VAN", "VANA", "VB11", "VB31", "VB71", "VB84035",
"VB84035M", "VB84069", "VB84069M", "VB8435", "VB8469", "VB86035",
"VB86035M", "VB86069", "VB86069M", "VB8635", "VB8669", "VB94035",
"VB94069", "VB9435", "VB9469", "VB96035", "VB96069", "VB9635",
"VB9669", "VBBHSH", "VBBHVB", "VBBHXH", "VBW4035", "VBW4069",
"VBW435", "VBW469", "VBW6035", "VBW6069", "VBW635", "VBW669",
"VBY4035", "VBY4069", "VBY6035", "VBY6069", "VBZ4035", "VBZ4069",
"VBZ6035", "VBZ6069", "VBZ6E035", "VBZ6E069", "VC11", "VC31",
"VC419", "VC435", "VC51", "VC619", "VC635", "VC91", "VCA419",
"VCA435", "VCA619", "VCA635", "VCDZ", "VCM36V", "VD0419",
"VD0619", "VD11", "VD21SDDAAN735", "VD21SEEAAN735", "VD31",
"VD419", "VD435", "VD51", "VD619", "VD635", "VD71", "VD91",
"VDGCS", "VDGDS", "VDL", "VDM", "VDNS", "VDPA", "VDPAF",
"XWYWMX", "XWYWTP", "XWYWTX", "XWYWXT", "XWZLUR", "XX11",
"XX31", "XX51", "XY11", "XY31", "XY51", "Y3AA", "Y3AB", "Y3AC",
"Y3ACA", "Y3AE", "Y3AF", "Y3AG", "Y3AH", "Y3AHA", "Y3AL",
"Y3AP", "Y3AR", "Y3AS", "Y3AT", "Y3AW", "Y3AX", "Y3AY", "Y3AZ",
"Y3CN", "Y4CW", "Y4CZ", "Y4GB", "Y4GG", "Y4GM", "Y4GN", "Y4GR",
"Y4GU", "Y4GV", "Y4GW", "Y4GX", "Y4GY", "Y4GZ", "Y4MZ", "Y4NW",
"Y4NX", "Y4NY", "Y4NZ", "Y4RM", "Y4RN", "Y4TD", "Y4TN", "Y4TR",
"Y4TS", "Y4TT", "Y4TU", "Y4TV", "Y4TX", "Y4WC", "Y4WG", "Y4WH",
"Y4WK", "Y51AA01", "Y51BB0", "Y51BB01", "Y51CC01", "Y51DD01",
"Y51FF0", "Y51HEE0", "Y51HEE01", "Y910", "YA01", "YA1MFA",
"YA1MFB", "YA1MRA", "YA2MFA", "YA2MFB", "YA2MRA", "YA3MFA",
"YA3MFB", "YA3MRA", "YA41", "YA61", "YA81", "YA9S", "YAAMFA",
"YAAMFAAX", "YAAMFB", "YAAMFBBX", "YAAMPA", "YAAMRA", "YABMFA",
"YABMFAAX", "YABMFB", "YABMFBBX", "YABMPA", "YABMRA", "YASMFA",
"YASMFB", "YASMPA", "YASMRA", "YATMFA", "YATMFB", "YATMPA",
"YATMRA", "YAUMPA", "YAUMRA", "YB01", "YB1MFA", "YB1MFB",
"YB1MFC", "YB1MRB", "YB2MFA", "YB2MFB", "YB2MFC", "YB2MRB",
"YB2R", "YB3MFA", "YB3MFB", "YB3MFC", "YB3MRB", "YB4R", "YB61",
"YB81", "YBAMAA", "YBAMAAAX", "YBAMAB", "YBAMABAX", "YBAMAC",
"YBAMACAX", "YBAMADAX", "YBAMCB", "YBAMCBAX", "YBAMCC", "YBAMCCAX",
"YBAMDBAX", "YBAMDCAX", "YBAMFA", "YBAMFAAX", "YBAMFB", "YBAMFBAX",
"YBAMFBBX", "YBAMFC", "YBAMFCBX", "YBAMFCCX", "YBAMGBAX",
"YBAMHBAX", "YBAMHCAX", "YBAMRB", "YBB5D11M6", "YBB5P11M6",
"YBB5P21M5", "YBB5P31M5", "YBB5P51A4", "YBBMAA", "YBBMAB",
"YBBMABAX", "YBBMAC", "YBBMACAX", "YBBMADAX", "YBBMCB", "YBBMCBAX",
"YBBMCC", "YBBMCCAX", "YBBMDBAX", "YBBMDC", "YBBMDCAX", "YBBMFA",
"YBBMFAAX", "YBBMFB", "YBBMFBAX", "YBBMFBBX", "YBBMFC", "YBBMFCBX",
"YBBMFCCX", "YBBMGBAX", "YBBMHBAX", "YBBMHCAX", "YBBMPB",
"YBBMRB", "YBC5D31M6", "YBCP11M6", "YBCP41M6", "YBDMAB",
"YBDMAC", "YBDMADAX", "YBDMDC", "YBDMDCAX", "YBDMFA", "YBDMFB",
"YBDMFBAX", "YBDMFBBX", "YBDMFC", "YBDMFCBX", "YBDMFCCX",
"YBDMHCAX", "YBDMPB", "YBPMFB", "YBPMFC", "YBPMPB", "YBSMFA",
"YBSMFB", "YBSMFC", "YBSMRB", "YBTMFA", "YBTMFB", "YBTMFC",
"YBTMPB", "YBTMRB", "YBUMFA", "YBUMFB", "YBUMFC", "YBUMPB",
"YBUMRB", "YC1MFA", "YC1MFB", "YC1MFC", "YC21", "YC2MFA",
"YC2MFB", "YC2MFC", "YC3MFB", "YC3MFC", "YC41", "YC81", "YCAMAB",
"YCAMABAX", "YCAMAC", "YCAMACAX", "YCAMACAXL", "YCAMAD",
"YCAMADAX", "YCAMCB", "YCAMCBAX", "YCAMCC", "YCAMCCAX", "YCAMCCAXL",
"YCAMDB", "YCAMDBAX", "YCAMDC", "YCAMDCAX", "YCAMFB", "YCAMFBAX",
"YCAMFBBX", "YCAMFC", "YCAMFCBX", "YCAMFCCX", "YCAMGCAX",
"YCAMHC", "YCAMHCAX", "YCBMAB", "YCBMABAX", "YCBMAC", "YCBMACAX",
"YCBMACAXL", "YCBMAD", "YCBMADAX", "YCBMCB", "YCBMCBAX",
"YCBMCC", "YCBMCCAX", "YCBMCCAXL", "YCBMDB", "YCBMDBAX",
"YCBMDC", "YCBMDCAX", "YCBMFB", "YCBMFBAX", "YCBMFBBX", "YCBMFC",
"YCBMFCBX", "YCBMFCCX", "YCBMGC", "YCBMGCAX", "YCBMHC", "YCBMHCAX",
"YCDMAB", "YCDMABAX", "YCDMAC", "YCDMACAX", "YCDMACAXL",
"YCDMAD", "YCDMADAX", "YCDMCB", "YCDMCC", "YCDMCCAX", "YCDMCCAXL",
"YCDMDB", "YCDMDC", "YCDMDCAX", "YCDMFB", "YCDMFBAX", "YCDMFBBX",
"YCDMFC", "YCDMFCBX", "YCDMFCCX", "YCPMFB", "YCPMFC", "YCSMFB",
"YCSMFC", "YCTMFA", "YCTMFB", "YCTMFC", "YCUMFB", "YCUMFC",
"YD2MFC", "YD3MFC", "YDBMACAX", "YDBMACAXL", "YDBMBCAX",
"YDBMCCAXL", "YDBMDCAXL", "YDBMFCBX", "YDBMFCCX", "YDBMFCCXL",
"YDBMGCAXL", "YDBMHCAXL", "YDDMACAX", "YDDMACAXL", "YDDMCCAXL",
"YDDMDCAXL", "YDDMFCBX", "YDDMFCBXL", "YDDMFCCX", "YDDMFCCXL",
"YDDMHCAXL", "YDPMFC", "YDTMFC", "YDUMFC", "YE01", "YE41",
"YE81", "YF61", "YF81", "YG01", "YG41", "YM", "YM11", "YM31",
"YM91", "YN", "YN1P", "YN2B", "YN2BCF", "YN2BDG", "YN2C",
"YN2M", "YN2MC", "YN2MD", "YN4M", "YN5B", "YN5M", "YN9P",
"YN9S", "YNA", "YNSF5DH1M6", "YNSF5P71M5", "YNSF5P91M6",
"YP", "YP11", "YP31", "YP51", "YP91", "YPA", "YS", "YS31",
"YS3B55D", "YS3F55E", "YS3F55L", "YS3F59E", "YS3GG4JA6A01",
"YS3GG4LA6F01", "YS3GG4LM6F02", "YS3GG4ZA6A01", "YS3GG4ZA6F01",
"YS3GG4ZM6A01", "YS3GG4ZM6F01", "YS51", "YS91", "YT11", "YT54Y",
"YT71", "YT91", "YU11", "YU31", "YU51", "YU71", "YV1PW10BD",
"YV1PW68BC", "YV1PW79B0", "YV1PW79B1", "YV1PWA2BC", "YV1PWA8B1",
"YV1PWA8BC", "YV1PWA8BD", "YV1PWARBC", "YV1PZ68BC", "YV1PZA2BC",
"YV1PZA8B4", "YV1PZA8BC", "YV71", "YV91", "YY23P", "YYCAD2M",
"Z10", "Z12", "Z12AA01", "Z12AA02", "Z12BB01", "Z16AMJ",
"Z16AMN", "Z34AA01", "Z34AA02", "Z34AA03", "Z34AA04", "Z34AA05",
"Z34BB01", "Z34BB02", "Z34BB03", "Z51AA01", "Z51BB01", "Z51CC01",
"ZA31", "ZAAMFA/BX", "ZAAMFAAX", "ZAAMFABX", "ZAAMFBAX",
"ZAAMNABX", "ZAAMPA", "ZAAMRA", "ZADA", "ZADAB", "ZADAC",
"ZADB", "ZADBB", "ZADBC", "ZADD", "ZADDB", "ZADE", "ZADEB",
"ZADFT", "ZADH", "ZADJ", "ZADM", "ZADN", "ZAKE", "ZAKET",
"ZAKF", "ZAKH", "ZAKJ", "ZALAT", "ZALB", "ZALC", "ZALF",
"ZALG", "ZALH", "ZALHT", "ZALJ", "ZAPAPA", "ZAPMFAAX", "ZAPMFAAY",
"ZAPMNAAX", "ZAPMPA", "ZARMFAAX", "ZARMFABX", "ZARMFBAX",
"ZARMNAAX", "ZARMPA", "ZARMRA", "ZATMFAAX", "ZATMFABX", "ZATMFBAX",
"ZATMNAAX", "ZAXZ", "ZAZA", "ZAZAT", "ZAZB", "ZAZBT", "ZAZC",
"ZAZD", "ZAZE", "ZAZF", "ZAZH", "ZAZJ", "ZAZK", "ZAZL", "ZAZT",
"ZAZTT", "ZAZU", "ZAZY", "ZAZZ", "ZB31", "ZB51", "ZB71",
"ZB8S", "ZB8Y", "ZB8Z", "ZBAMAAAX", "ZBAMABAX", "ZBAMDAAX",
"ZBAMDBAX", "ZBAMDCAX", "ZBAMFAAX", "ZBAMFABX", "ZBAMFBAX",
"ZBAMNAAX", "ZBAMNBAX", "ZBAMRB", "ZBC", "ZBD", "ZBPANBAX",
"ZBPANCAX", "ZBPARB", "ZBPMABAY", "ZBPMACAY", "ZBPMBBAY",
"ZBPMBCAY", "ZBPMCBAX", "ZBPMDCAY", "ZBPMFAAY", "ZBPMFABY",
"ZBPMFBAY", "ZBPMGBAY", "ZBPMGCAY", "ZBPMHCAY", "ZBPMNAAX",
"ZBPMNABY", "ZBPMNBAY", "ZBPMNCAX", "ZBPMNCAY", "ZBPMRB",
"ZBPWFAAY", "ZBPWFABY", "ZBPWNBAY", "ZBRMAAAX", "ZBRMABAX",
"ZBRMACAX", "ZBRMCBAX", "ZBRMDAAX", "ZBRMDBAX", "ZBRMDCAX",
"ZBRMFAAX", "ZBRMFABX", "ZBRMFBAX", "ZBRMGBAX", "ZBRMGCAX",
"ZBRMHCAX", "ZBRMNAAX", "ZBRMNBAX", "ZBRMNCAX", "ZBRMRB",
"ZBTMFAAX", "ZBTMFABX", "ZBTMFBAX", "ZBTMNAAX", "ZBTMNBAX",
"ZBTMNCAX", "ZBZMACAY", "ZBZMDCAY", "ZBZMFAAY", "ZBZMFABY",
"ZBZMFBAY", "ZBZMGCAY", "ZBZMHCAY", "ZBZMNBAY", "ZBZMNCAY",
"ZC31", "ZC51", "ZCBHZH", "ZCBHZW", "ZCPANBAX", "ZCPANBBX",
"ZCPANCAX", "ZCPANCBX", "ZCPMABAY", "ZCPMACAY", "ZCPMAGAY",
"ZCPMBBAY", "ZCPMBCAY", "ZCPMCBAY", "ZCPMCCAY", "ZCPMDBAY",
"ZCPMDCAY", "ZCPMFAAY", "ZCPMFABY", "ZCPMGCAY", "ZCPMHCAX",
"ZCPMHCAY", "ZCPMNBAY", "ZCPMNBBY", "ZCPMNCAY", "ZCPMNCBY",
"ZCPWABAY", "ZCPWACAY", "ZCPWCBAY", "ZCPWCCAY", "ZCPWDBAY",
"ZCPWDCAY", "ZCPWNBAY", "ZCPWNCAY", "ZCRMABAX", "ZCRMACAX",
"ZCRMAGAX", "ZCRMCBAX", "ZCRMCCAX", "ZCRMDBAX", "ZCRMDCAX",
"ZCRMFAAX", "ZCRMFABX", "ZCRMGCAX", "ZCRMHCAX", "ZCRMNBAX",
"ZCRMNBBX", "ZCRMNCAX", "ZCRMNCBX", "ZCTMCBAX", "ZCTMCCAX",
"ZCTMNBAX", "ZCTMNBBX", "ZCTMNCAX", "ZCTMNCBX", "ZCZMABAY",
"ZCZMACAY", "ZCZMAGAY", "ZCZMCBAY", "ZCZMCCAY", "ZCZMDBAY",
"ZCZMDCAY", "ZCZMDGAY", "ZCZMFAAY", "ZCZMFABY", "ZCZMGCAY",
"ZCZMGGAY", "ZCZMHCAY", "ZCZMHGAY", "ZCZMNBAY", "ZCZMNBBY",
"ZCZMNCAY", "ZCZMNCBY", "ZD11", "ZD31", "ZD51", "ZD71", "ZE2",
"ZE31", "ZF111", "ZF31", "ZFANF6EK6", "ZFANF6EKA", "ZG31",
"ZH51", "ZH71", "ZHNC9L", "ZHNE9L", "ZM31", "ZN31", "ZN57",
"ZN67", "ZN8M", "ZP31", "ZR31", "ZR71", "ZRHNYH", "ZRHNYW",
"ZT14", "ZT24", "ZT35", "ZT37", "ZT38", "ZT45", "ZT47", "ZT48",
"ZT55", "ZT57", "ZT58", "ZT65", "ZT67", "ZT68", "ZV18", "ZV30",
"ZV38", "ZV41", "ZV81", "ZW41", "ZW61", "ZW81", "ZZT220LAEMNKW",
"ZZT220LALMNKW", "ZZT221LAEMEKW", "ZZT221LAEPEKW", "ZZT221LALMEKW",
"ZZT221LALPEKW", "ZZT230LBLFGHW", "ZZT231LBLFVFW", "ZZT251LALMEKW",
"ZZT251LALMNKW", "ZZT251LALPEKW", "ZZW30LAKMQHW"), class = "factor")), .Names
= c("Lib_Marque__MRQ_",
"Lib_Modele__MOD_", "Lib_Carrosserie__CAR_", "Lib_Energie__ENE_",
"Puiss_Fiscale", "Nb_Places", "Nb_Cylindres", "Cylindree", "Vitesse_Max",
"Puiss_Reelle_kW", "Puiss_Reelle_DIN", "Regime_Puiss", "Couple_Max",
"Regime_Couple_Max", "Nb_Rapports", "Longueur", "Largeur", "Hauteur",
"Empattement", "Voie_AV", "Voie_AR", "Poids_Vide", "PTAC", "Charge_Utile",
"Valeur_Neuf_Orig_EUR", "Code_SRA"), row.names = c(NA, 10L), class =
"data.frame")
Thank you for your kind help
Any help would be greatly appreciated!!
I'm trying to create a choropleth map in R that shows the counties of texas, color-coded by their population ranges.
My problem is that the range of populations is too large. The highest population is over 4 million, but most of the counties have a population under 50,000. The criteria for the fill is: (0-1mil), (1-2mil), (2-3mil), (3-4mil), (4-5mil) but almost all fall under 0-1mil.
How can I change the legend to account for different ranges of numbers? For example, maybe:
(0-1,000), (1,000-10,000), (10,000-100,000), (100,000-1mil), (1mil-5mil)
Here's the code I wrote to plot the data:
txplot <- ggplot(txczpop, aes(fill=pop2014)) + geom_map(txmap)
tm_shape(txmap) +
tm_fill("pop2014", title="TX County Population", palette = "PRGn") +
tm_borders(alpha=.5) +
tm_style_beaver()
Here's the result:
[![enter image description here][1]][1]
I'm using a census county shapefile and population also retrieved from a census file.
Here's the output of my population data:
txczpop <- structure(list(county_fips = c(48001L, 48003L, 48005L, 48007L,
48009L, 48011L, 48013L, 48015L, 48017L, 48019L, 48021L, 48023L,
48025L, 48027L, 48029L, 48031L, 48033L, 48035L, 48037L, 48039L,
48041L, 48043L, 48045L, 48047L, 48049L, 48051L, 48053L, 48055L,
48057L, 48059L, 48061L, 48063L, 48065L, 48067L, 48069L, 48071L,
48073L, 48075L, 48077L, 48079L, 48081L, 48083L, 48085L, 48087L,
48089L, 48091L, 48093L, 48095L, 48097L, 48099L, 48101L, 48103L,
48105L, 48107L, 48109L, 48111L, 48113L, 48115L, 48117L, 48119L,
48121L, 48123L, 48125L, 48127L, 48129L, 48131L, 48133L, 48135L,
48137L, 48141L, 48139L, 48143L, 48145L, 48147L, 48149L, 48151L,
48153L, 48155L, 48157L, 48159L, 48161L, 48163L, 48165L, 48167L,
48169L, 48171L, 48173L, 48175L, 48177L, 48179L, 48181L, 48183L,
48185L, 48187L, 48189L, 48191L, 48193L, 48195L, 48197L, 48199L,
48201L, 48203L, 48205L, 48207L, 48209L, 48211L, 48213L, 48215L,
48217L, 48219L, 48221L, 48223L, 48225L, 48227L, 48229L, 48231L,
48233L, 48235L, 48237L, 48239L, 48241L, 48243L, 48245L, 48247L,
48249L, 48251L, 48253L, 48255L, 48257L, 48259L, 48261L, 48263L,
48265L, 48267L, 48269L, 48271L, 48273L, 48275L, 48283L, 48277L,
48279L, 48281L, 48285L, 48287L, 48289L, 48291L, 48293L, 48295L,
48297L, 48299L, 48301L, 48303L, 48305L, 48313L, 48315L, 48317L,
48319L, 48321L, 48323L, 48307L, 48309L, 48311L, 48325L, 48327L,
48329L, 48331L, 48333L, 48335L, 48337L, 48339L, 48341L, 48343L,
48345L, 48347L, 48349L, 48351L, 48353L, 48355L, 48357L, 48359L,
48361L, 48363L, 48365L, 48367L, 48369L, 48371L, 48373L, 48375L,
48377L, 48379L, 48381L, 48383L, 48385L, 48387L, 48389L, 48391L,
48393L, 48395L, 48397L, 48399L, 48401L, 48403L, 48405L, 48407L,
48409L, 48411L, 48413L, 48415L, 48417L, 48419L, 48421L, 48423L,
48425L, 48427L, 48429L, 48431L, 48433L, 48435L, 48437L, 48439L,
48441L, 48443L, 48445L, 48447L, 48449L, 48451L, 48453L, 48455L,
48457L, 48459L, 48461L, 48463L, 48465L, 48467L, 48469L, 48471L,
48473L, 48475L, 48477L, 48479L, 48481L, 48483L, 48485L, 48487L,
48489L, 48491L, 48493L, 48495L, 48497L, 48499L, 48501L, 48503L,
48505L, 48507L), county_name = c("Anderson", "Andrews", "Angelina",
"Aransas", "Archer", "Armstrong", "Atascosa", "Austin", "Bailey",
"Bandera", "Bastrop", "Baylor", "Bee", "Bell", "Bexar", "Blanco",
"Borden", "Bosque", "Bowie", "Brazoria", "Brazos", "Brewster",
"Briscoe", "Brooks", "Brown", "Burleson", "Burnet", "Caldwell",
"Calhoun", "Callahan", "Cameron", "Camp", "Carson", "Cass", "Castro",
"Chambers", "Cherokee", "Childress", "Clay", "Cochran", "Coke",
"Coleman", "Collin", "Collingsworth", "Colorado", "Comal", "Comanche",
"Concho", "Cooke", "Coryell", "Cottle", "Crane", "Crockett",
"Crosby", "Culberson", "Dallam", "Dallas", "Dawson", "Deaf Smith",
"Delta", "Denton", "DeWitt", "Dickens", "Dimmit", "Donley", "Duval",
"Eastland", "Ector", "Edwards", "El Paso", "Ellis", "Erath",
"Falls", "Fannin", "Fayette", "Fisher", "Floyd", "Foard", "Fort Bend",
"Franklin", "Freestone", "Frio", "Gaines", "Galveston", "Garza",
"Gillespie", "Glasscock", "Goliad", "Gonzales", "Gray", "Grayson",
"Gregg", "Grimes", "Guadalupe", "Hale", "Hall", "Hamilton", "Hansford",
"Hardeman", "Hardin", "Harris", "Harrison", "Hartley", "Haskell",
"Hays", "Hemphill", "Henderson", "Hidalgo", "Hill", "Hockley",
"Hood", "Hopkins", "Houston", "Howard", "Hudspeth", "Hunt", "Hutchinson",
"Irion", "Jack", "Jackson", "Jasper", "Jeff Davis", "Jefferson",
"Jim Hogg", "Jim Wells", "Johnson", "Jones", "Karnes", "Kaufman",
"Kendall", "Kenedy", "Kent", "Kerr", "Kimble", "King", "Kinney",
"Kleberg", "Knox", "La Salle", "Lamar", "Lamb", "Lampasas", "Lavaca",
"Lee", "Leon", "Liberty", "Limestone", "Lipscomb", "Live Oak",
"Llano", "Loving", "Lubbock", "Lynn", "Madison", "Marion", "Martin",
"Mason", "Matagorda", "Maverick", "McCulloch", "McLennan", "McMullen",
"Medina", "Menard", "Midland", "Milam", "Mills", "Mitchell",
"Montague", "Montgomery", "Moore", "Morris", "Motley", "Nacogdoches",
"Navarro", "Newton", "Nolan", "Nueces", "Ochiltree", "Oldham",
"Orange", "Palo Pinto", "Panola", "Parker", "Parmer", "Pecos",
"Polk", "Potter", "Presidio", "Rains", "Randall", "Reagan", "Real",
"Red River", "Reeves", "Refugio", "Roberts", "Robertson", "Rockwall",
"Runnels", "Rusk", "Sabine", "San Augustine", "San Jacinto",
"San Patricio", "San Saba", "Schleicher", "Scurry", "Shackelford",
"Shelby", "Sherman", "Smith", "Somervell", "Starr", "Stephens",
"Sterling", "Stonewall", "Sutton", "Swisher", "Tarrant", "Taylor",
"Terrell", "Terry", "Throckmorton", "Titus", "Tom Green", "Travis",
"Trinity", "Tyler", "Upshur", "Upton", "Uvalde", "Val Verde",
"Van Zandt", "Victoria", "Walker", "Waller", "Ward", "Washington",
"Webb", "Wharton", "Wheeler", "Wichita", "Wilbarger", "Willacy",
"Williamson", "Wilson", "Winkler", "Wise", "Wood", "Yoakum",
"Young", "Zapata", "Zavala"), pop2014 = c(57627L, 17477L, 87750L,
24972L, 8811L, 1955L, 47774L, 29114L, 6910L, 20892L, 78069L,
3592L, 32863L, 329140L, 1855866L, 10812L, 652L, 17780L, 93275L,
338124L, 209152L, 9173L, 1536L, 7194L, 37653L, 17253L, 44943L,
39810L, 21797L, 13513L, 420392L, 12621L, 6013L, 30261L, 7781L,
38145L, 50902L, 7089L, 10370L, 2935L, 3254L, 8430L, 885241L,
3017L, 20719L, 123694L, 13550L, 4050L, 38761L, 75562L, 1415L,
4950L, 3812L, 5899L, 2266L, 7135L, 2518638L, 13372L, 19195L,
5238L, 753363L, 20684L, 2218L, 11089L, 3543L, 11533L, 18176L,
153904L, 1879L, 833487L, 159317L, 40147L, 16989L, 33752L, 24833L,
3831L, 5949L, 1275L, 685345L, 10600L, 19762L, 18531L, 19425L,
314198L, 6435L, 25520L, 1291L, 7549L, 20462L, 23044L, 123534L,
123204L, 27172L, 147250L, 34720L, 3147L, 8199L, 5509L, 3928L,
55621L, 4441370L, 67336L, 6089L, 5769L, 185025L, 4180L, 79290L,
831073L, 34848L, 23577L, 53921L, 35921L, 22741L, 36651L, 3211L,
88493L, 21773L, 1574L, 8855L, 14739L, 35552L, 2204L, 252235L,
5255L, 41353L, 157456L, 19936L, 14906L, 111236L, 38880L, 400L,
785L, 50562L, 4438L, 262L, 3526L, 32190L, 3858L, 7474L, 49523L,
13574L, 20156L, 19721L, 16742L, 16861L, 78117L, 23524L, 3553L,
12091L, 19510L, 86L, 293974L, 5771L, 13861L, 10149L, 5460L, 4071L,
36519L, 57023L, 8199L, 243441L, 805L, 47894L, 2147L, 155830L,
24256L, 4870L, 9076L, 19416L, 518947L, 22148L, 12743L, 1153L,
65301L, 48195L, 14138L, 15093L, 356221L, 10758L, 2070L, 83433L,
28096L, 23769L, 123164L, 9908L, 15893L, 46079L, 121627L, 6976L,
11032L, 128220L, 3755L, 3371L, 12446L, 14349L, 7302L, 928L, 16500L,
87809L, 10416L, 53923L, 10350L, 8610L, 27099L, 66915L, 5622L,
3162L, 17328L, 3343L, 25515L, 3084L, 218842L, 8694L, 62955L,
9405L, 1339L, 1403L, 3972L, 7581L, 1945360L, 135143L, 927L, 12739L,
1608L, 32506L, 116608L, 1151145L, 14224L, 21418L, 40354L, 3454L,
27117L, 48974L, 52910L, 91081L, 69789L, 46820L, 11625L, 34438L,
266673L, 41168L, 5714L, 132355L, 12973L, 21903L, 489250L, 46402L,
7821L, 61638L, 42852L, 8286L, 18350L, 14319L, 12267L)), .Names = c("county_fips",
"county_name", "pop2014"), row.names = c(5100L, 5101L, 5103L,
5106L, 5107L, 5109L, 5112L, 5114L, 5116L, 5118L, 5120L, 5121L,
5124L, 5126L, 5128L, 5129L, 5131L, 5133L, 5136L, 5137L, 5140L,
5141L, 5143L, 5146L, 5147L, 5150L, 5152L, 5153L, 5156L, 5158L,
5159L, 5161L, 5163L, 5166L, 5168L, 5170L, 5171L, 5174L, 5176L,
5178L, 5179L, 5182L, 5183L, 5185L, 5188L, 5190L, 5192L, 5194L,
5195L, 5198L, 5200L, 5201L, 5203L, 5205L, 5208L, 5209L, 5212L,
5214L, 5215L, 5218L, 5219L, 5221L, 5224L, 5226L, 5228L, 5230L,
5232L, 5233L, 5235L, 5239L, 5237L, 5242L, 5244L, 5245L, 5248L,
5249L, 5251L, 5254L, 5256L, 5257L, 5260L, 5261L, 5264L, 5265L,
5268L, 5270L, 5272L, 5274L, 5276L, 5278L, 5280L, 5281L, 5284L,
5286L, 5288L, 5290L, 5292L, 5293L, 5296L, 5298L, 5300L, 5301L,
5303L, 5306L, 5308L, 5309L, 5312L, 5314L, 5316L, 5317L, 5319L,
5321L, 5323L, 5326L, 5327L, 5330L, 5332L, 5334L, 5335L, 5337L,
5339L, 5341L, 5343L, 5346L, 5348L, 5349L, 5352L, 5354L, 5356L,
5357L, 5360L, 5362L, 5364L, 5365L, 5368L, 5369L, 5372L, 5374L,
5382L, 5376L, 5378L, 5379L, 5383L, 5385L, 5388L, 5390L, 5392L,
5394L, 5396L, 5398L, 5400L, 5401L, 5404L, 5412L, 5413L, 5416L,
5418L, 5419L, 5421L, 5406L, 5407L, 5409L, 5423L, 5425L, 5427L,
5429L, 5432L, 5434L, 5435L, 5438L, 5440L, 5442L, 5443L, 5446L,
5448L, 5449L, 5451L, 5453L, 5456L, 5457L, 5460L, 5461L, 5464L,
5465L, 5468L, 5470L, 5472L, 5474L, 5476L, 5477L, 5480L, 5482L,
5484L, 5486L, 5488L, 5489L, 5491L, 5494L, 5496L, 5498L, 5499L,
5501L, 5504L, 5505L, 5508L, 5510L, 5511L, 5514L, 5516L, 5518L,
5520L, 5522L, 5524L, 5526L, 5527L, 5530L, 5531L, 5533L, 5536L,
5537L, 5540L, 5542L, 5544L, 5546L, 5547L, 5550L, 5552L, 5554L,
5555L, 5558L, 5559L, 5562L, 5563L, 5566L, 5568L, 5569L, 5571L,
5574L, 5575L, 5578L, 5579L, 5582L, 5584L, 5585L, 5587L, 5590L,
5592L, 5594L, 5595L, 5598L, 5600L, 5602L, 5604L, 5606L), class = "data.frame")
I just created a new column in the population dataframe that summarizes the population based on the ranges that I want to use, and then use that as the criteria for the fill:
txczpop$poprange[txczpop$pop2014 >= 0 & txczpop < 1000] <- "0-1,000"
txczpop$poprange[txczpop$pop2014 >= 1000 & txczpop < 10000] <- "1-10,000"
txczpop$poprange[txczpop$pop2014 >= 10000 & txczpop$pop2014 < 100000] <- "10,000-100,000"
txczpop$poprange[txczpop$pop2014 >= 100000 & txczpop$pop2014 < 1000000] <- "100,000 - 1,000,000"
txczpop$poprange[txczpop$pop2014 >= 1000000 & txczpop$pop2014 <= 5000000] <- "1,000,000 - 5,000,000"
I'm trying to fit some curves from the data below with a mono-exponential "decay". Graphical display is not as important as is pulling out the time constant. the y-axis is pA and the x is time in seconds.
dput(stackover_data)
structure(list(Time = c(0.09990001, 0.19990001, 0.29990001, 0.39990001,
0.49990001, 0.59990001, 0.69990001, 0.79990001, 0.89990001, 0.99990001,
1.09990001, 1.19990001, 1.29990001, 1.39990001, 1.49990001, 1.59990001,
1.69990001, 1.79990001, 1.89990001, 1.99990001, 2.09990001, 2.19990001,
2.29990001, 2.39990001, 2.49990001, 2.59990001, 2.69990001, 2.79990001,
2.89990001, 2.99990001, 3.09990001, 3.19990001, 3.29990001, 3.39990001,
3.49990001, 3.59990001, 3.69990001, 3.79990001, 3.89990001, 3.99990001,
4.09990001, 4.19990001, 4.29990001, 4.39990001, 4.49990001, 4.59990001,
4.69990001, 4.79990001, 4.89990001, 4.99990001, 5.09990001, 5.19990001,
5.29990001, 5.39990001, 5.49990001, 5.59990001, 5.69990001, 5.79990001,
5.89990001, 5.99990001, 6.09990001, 6.19990001, 6.29990001, 6.39990001,
6.49990001, 6.59990001, 6.69990001, 6.79990001, 6.89990001, 6.99990001,
7.09990001, 7.19990001, 7.29990001, 7.39990001, 7.49990001, 7.59990001,
7.69990001, 7.79990001, 7.89990001, 7.99990001, 8.09990001, 8.19990001,
8.29990001, 8.39990001, 8.49990001, 8.59990001, 8.69990001, 8.79990001,
8.89990001, 8.99990001, 9.09990001, 9.19990001, 9.29990001, 9.39990001,
9.49990001, 9.59990001, 9.69990001, 9.79990001, 9.89990001, 9.99990001,
10.09990001, 10.19990001, 10.29990001, 10.39990001, 10.49990001,
10.59990001, 10.69990001, 10.79990001, 10.89990001, 10.99990001,
11.09990001, 11.19990001, 11.29990001, 11.39990001, 11.49990001,
11.59990001, 11.69990001, 11.79990001, 11.89990001, 11.99990001,
12.09990001, 12.19990001, 12.29990001, 12.39990001, 12.49990001,
12.59990001, 12.69990001, 12.79990001, 12.89990001, 12.99990001,
13.09990001, 13.19990001, 13.29990001, 13.39990001, 13.49990001,
13.59990001, 13.69990001, 13.79990001, 13.89990001, 13.99990001,
14.09990001, 14.19990001, 14.29990001, 14.39990001, 14.49990001,
14.59990001, 14.69990001, 14.79990001, 14.89990001, 14.99990001,
15.09990001, 15.19990001, 15.29990001, 15.39990001, 15.49990001,
15.59990001, 15.69990001, 15.79990001, 15.89990001, 15.99990001,
16.09990001, 16.19990001, 16.29990001, 16.39990001, 16.49990001,
16.59990001, 16.69990001, 16.79990001, 16.89990001, 16.99990001,
17.09990001, 17.19990001, 17.29990001, 17.39990001, 17.49990001,
17.59990001, 17.69990001, 17.79990001, 17.89990001, 17.99990001,
18.09990001, 18.19990001, 18.29990001, 18.39990001, 18.49990001,
18.59990001, 18.69990001, 18.79990001, 18.89990001, 18.99990001,
19.09990001, 19.19990001, 19.29990001, 19.39990001, 19.49990001,
19.59990001, 19.69990001, 19.79990001, 19.89990001, 19.99990001
), `Trace 1` = c(-3.08656892325052, 9.36821982641837, 8.13806079083122,
10.7039590839898, 7.25670468903547, 4.31122291688919, 1.77905971163193,
-6.27606834721828, -8.65955381985049, -10.1445673910916, -10.6649772153892,
-6.52301948183154, -11.9757817854835, -0.976254254762154, 0.237467076202677,
-11.8114896779541, -11.0022757370468, -16.1845427042923, 2.70565927469852,
-18.9048281652216, -22.153682283437, -4363.32044948884, -3470.59111611883,
-2877.51064886248, -2253.51505229908, -1876.62974792002, -1541.14546478629,
-1288.26617158403, -1232.52313999354, -1042.98549163259, -938.795146277054,
-810.913567086442, -736.390770574588, -682.044521632168, -638.575324886466,
-596.587948389699, -629.282103146111, -569.200610245336, -587.387817942122,
-521.939307772762, -533.776693538631, -514.782451411908, -445.949395199026,
-451.026618716539, -461.669600651513, -442.372477296489, -426.132547857502,
-452.471101919398, -449.377497412324, -436.817659066873, -487.464805660851,
-442.915455035179, -458.666741489705, -454.990437793055, -455.173154690614,
-426.885702219019, -456.985443835707, -408.869318611773, -442.418742303429,
-407.42593099033, -414.538253423006, -316.036755461507, -248.395118743017,
-198.708015370115, -183.88480278352, -160.041754268138, -139.020785864805,
-123.701641615743, -137.253323547789, -124.240619047461, -121.512859107349,
-113.719386521321, -102.98429740535, -118.721098087137, -109.789471870234,
-122.145399997255, -109.542467905009, -96.6725326170008, -108.16233877188,
-94.1092968366083, -88.0987495250118, -89.2245425850472, -96.3495439993499,
-110.340434956898, -98.8777265479938, -88.0674013452629, -102.755317774957,
-91.0752842494157, -107.959830375198, -78.8424385901398, -93.7743443479161,
-98.3591889757604, -72.3214090579355, -85.6296201608712, -112.596656084088,
-115.518068650615, -110.973655632476, -78.8021307215932, -85.6443070182152,
-93.1744356389988, -100.483987323044, -87.6672104421484, -83.6481077757535,
-82.1465876740715, -86.7838666454595, -75.9066755520263, -86.6416980998645,
-88.9405806921788, -74.0592581080291, -86.3433338300531, -93.5114839187431,
-91.1875849041866, -82.7018083540351, -102.859075734953, -82.1494206590809,
-79.2197323780198, -78.1558787387238, -86.2649418863144, -93.5271994290692,
-84.6678528566242, -71.4828270654073, -69.5618263581887, -67.2920558863641,
-58.1490330439793, -59.8163238740351, -73.6957797622946, -61.7673947702343,
-63.6492255747164, -64.043638367468, -73.8301991524909, -69.3055992018769,
-65.9342860783478, -74.9891521715357, -74.5779808619617, -69.9029875902787,
-69.746935396023, -71.35455530782, -73.5279471991205, -77.8000626250279,
-74.6065254864801, -64.0834786591292, -45.1346212136811, -43.1615385011179,
-42.554841323715, -50.276566542849, -56.1940469314277, -49.6368473019083,
-53.9842269565738, -54.1480156577708, -48.2160751112714, -43.2902743874793,
-46.4358385791469, -58.4459814967147, -49.0443619038303, -56.8007684056415,
-56.3345802020277, -40.1704617527471, -38.2624410086947, -39.9339279963857,
-27.7637250414188, -40.8686520798649, -30.9246406597275, -32.70021580322,
-28.7536071219405, -32.4403271637661, -30.3389870650563, -23.2866770834185,
-40.349370835336, -39.7478049410975, -38.4648592612063, -37.4191914226093,
-52.1259539643142, -38.0993644763616, -37.5493511061199, -25.0794144286873,
-32.6165749053184, -17.2616237291667, -32.1515419380766, -27.0814392288745,
-26.2440551871993, -17.5882979851251, -19.852012057918, -22.1408312947152,
-21.3395092026716, -28.0789344732281, -19.6443281019371, -26.6585778398899,
-24.2994817465722, -27.847437388783, -22.0565059455712, -24.8747905836351,
-30.7544250721296, -29.9924061105416, -22.5728855329732, -24.9767037068316,
-18.8208568928653, -10.7306087216159, -16.4281210876173, -19.3057174287183,
-15.9745523586581), `Trace 2` = c(5.94927992143286, 1.42121402161905,
6.78788136514507, 3.33970424403748, -3.73956433922802, -7.3097330836793,
-9.18242380095097, 3.29952017048882, 5.17208246028, 1.53238537592179,
6.90832098860733, 3.16748380079213, 5.49988319742749, 3.86758484926656,
11.981378748128, -0.208159377585758, -7.2781503569274, 1.71389537416221,
-4.77614396689646, -0.561115871778583, -14.5385068323377, -4944.63006397872,
-3891.23141079918, -3065.66300899921, -2527.61853791531, -2060.05223709386,
-1763.27814995793, -1451.77884853301, -1368.42409410389, -1115.3385743493,
-990.942590410893, -852.951182088994, -920.516421149141, -860.486324532349,
-751.257935697359, -768.757880984424, -713.935725464811, -715.011716140783,
-616.355633698422, -672.133827760044, -629.418857335266, -607.342562427532,
-599.657392949772, -553.265193759536, -572.878306401819, -583.153736972572,
-589.661926862833, -527.552306430313, -553.672277977987, -539.906254466276,
-525.491319968336, -564.876755897633, -505.363603776286, -575.649358123809,
-497.897048596953, -544.360501913781, -483.501835830955, -533.085516954168,
-478.007776808595, -453.333986756674, -447.125338584774, -380.268661132557,
-253.127091546528, -227.877427383931, -191.488369613846, -166.799480062655,
-171.446527550207, -156.775073739292, -148.804014388301, -151.79670303645,
-152.599645619669, -113.40778617347, -133.161498039771, -129.929086626311,
-156.827022418459, -146.618934465797, -94.20997521855, -124.514491421933,
-89.8538131075507, -91.6500841191965, -113.887119052802, -109.971138752404,
-122.532101557384, -107.714865018275, -120.258046104427, -133.296582343699,
-111.726093380577, -121.533575165211, -93.2177144160483, -97.3820468959891,
-111.294540122284, -100.913019601926, -90.8456358331201, -94.143801546019,
-89.9179990101121, -79.3197988311624, -87.0020882554621, -103.315348613102,
-92.6091563297179, -87.3900461813473, -106.61173718685, -117.18534688083,
-106.417237368758, -114.977458547258, -101.759575332149, -93.2156781534656,
-79.4452885332048, -83.8035016923155, -85.7770898942133, -80.151862850598,
-72.4787832080729, -75.4887322569144, -70.8074563945233, -76.9714115362765,
-74.2305756099668, -72.7994824773122, -86.4890002605112, -68.5241155035978,
-76.723559709657, -60.5721619300473, -64.185161607482, -81.0777473166996,
-59.7023719535658, -48.0615168793405, -63.7668096308392, -69.0412888697281,
-70.9894443873153, -48.8366672284658, -43.1003899300401, -50.5591327356918,
-65.184353265464, -52.2850955061746, -66.0074402595226, -69.3081002474696,
-69.3876079446205, -53.0919105504865, -69.3289906081927, -63.5198462609014,
-70.9316181916525, -68.753726760587, -57.1360889764902, -45.7297495037924,
-45.0225501894778, -45.7487568762369, -46.486850874694, -46.7437348684622,
-41.1396209831975, -38.5056250833188, -49.8076229939374, -46.9868446480308,
-38.4639795698766, -40.8336417843587, -49.9810721301961, -58.7820616810459,
-51.4348042554539, -56.9132093579117, -53.527642322397, -49.7620292214747,
-54.6821412367706, -66.7962460110546, -52.0791798888695, -43.3022792658504,
-47.5584813904168, -36.3352478868053, -39.8699383685185, -43.9973042247199,
-49.3191665577514, -42.0818922168797, -37.7538461067954, -29.1591142899914,
-31.6815992654018, -35.3664789321372, -29.6435427630652, -45.5528322348158,
-37.0568327299924, -29.6739001502498, -37.0952936191428, -44.1786405275305,
-50.929131558954, -31.7842698176079, -31.4300485665082, -33.960404707688,
-34.0566040705914, -41.2299273666286, -30.1094672132779, -28.6802701495622,
-37.1973277887207, -31.7114680610898, -32.5450840790921, -28.5823680785893,
-27.5268500574853, -15.2389866754792, -19.2079199191329, -28.9156878395334,
-16.5863741697772, -31.6638747515138, -34.679018047114, -32.6358481677933,
-26.3463300811612, -18.542719903123), `Trace 3` = c(4.74984904312473,
-6.83237755613888, -4.15774846489683, 0.0116753456567835, 4.3974208810568,
8.56175743226069, -0.271281722577697, -10.7982854290066, -18.31704442037,
-15.7577103721081, -20.5775487342614, -17.7307614062074, -14.2844531791322,
-5.43267691733291, 2.52755064598145, -0.79977817636393, -8.07074631773528,
-14.6631228094247, -9.70657043283029, -11.4521183354252, -10.8401698276907,
-5495.45466682125, -4324.53856178012, -3525.03574657791, -2774.37269062517,
-2175.8494069729, -1790.54797875842, -1574.71991002551, -1342.90235937345,
-1212.59715782746, -1042.56659527622, -1006.18885759924, -962.27473179716,
-870.459562514915, -847.89784200018, -770.205392585522, -709.279500982242,
-740.484104629636, -739.054108223511, -645.039948927494, -712.347389757369,
-609.734365947307, -613.120930696929, -630.805982834916, -581.828460576409,
-608.843958294059, -602.432567890874, -547.023947351376, -571.475614300591,
-561.707380041648, -535.094375334464, -578.783728464712, -535.290908514374,
-489.922390143035, -522.824368605682, -520.322259788742, -556.271028784585,
-536.474251127913, -546.811709881227, -541.304029428261, -563.702125503968,
-460.326259578172, -348.718859066443, -261.740363110507, -244.084411205067,
-182.533320181224, -205.356578064485, -157.759947570376, -162.214202634405,
-180.011654146528, -176.904445787135, -167.116794266109, -163.416234638164,
-155.821745831568, -134.760994564946, -131.702791057298, -138.598664148971,
-139.818370538306, -134.700400696581, -137.275258643317, -119.921721377143,
-118.187525264885, -118.490007264421, -130.391338033262, -138.05114734155,
-127.505366899025, -132.636233878671, -123.490111075107, -113.42942794286,
-122.40925776411, -114.48254277334, -122.78054774327, -90.9528981496313,
-96.9508488566491, -107.842237007564, -91.0521871208113, -91.0347011279681,
-85.7054147126896, -89.7363815316723, -77.6588521980595, -76.6733051688234,
-85.821987540449, -90.3216600489126, -66.1137192454051, -77.7157754909198,
-89.8637060331546, -80.9149315468257, -71.7434481475478, -97.6600776026528,
-103.405900069071, -83.7502802699966, -77.4600968679799, -77.8705250264031,
-78.9012879762878, -87.4698127883497, -92.9587198547098, -96.1222986385502,
-91.1990314226718, -83.8744116889207, -61.7720026226848, -83.6616594930646,
-91.3899042204887, -74.2552104995002, -94.3117384183545, -81.130091800586,
-87.6638514042855, -79.061763018934, -78.7472359255489, -72.9642062301352,
-83.5759889935287, -74.066505526926, -62.6191749912397, -75.7150349466975,
-71.8704211546101, -68.4798935028027, -77.0470926055507, -76.6024078920054,
-66.5760393943856, -77.8848130194484, -77.5151449198185, -65.5482413040489,
-47.6087884294611, -35.6865291924552, -47.9648700291328, -62.8176093048552,
-48.969044800301, -46.853228499419, -52.562083189922, -58.7677963803543,
-63.6615259259388, -51.7351340481719, -54.5967510836755, -54.1294816023731,
-47.7621427567758, -40.9244974392914, -50.5286026229961, -45.9961132210891,
-43.143208905187, -43.8395623054062, -49.3326618452772, -48.308625153823,
-44.9833219970479, -32.4206723427767, -45.6628898472981, -41.4748918817457,
-43.6438737242104, -42.9844450147366, -30.3729810016378, -44.9519080631137,
-45.8241111308495, -40.4766896430523, -40.3007872752484, -46.4770136239361,
-49.3759491415156, -46.0844075639024, -34.3090862431636, -26.6047381016158,
-28.9119815377167, -33.1619464253006, -37.3813739690468, -46.5001186141503,
-54.5235621985407, -44.5233400090119, -55.2272273265567, -48.088333647706,
-55.3522589332341, -52.8146474401922, -60.1877313269088, -48.2512741329643,
-34.8540879147507, -55.2019332852645, -50.8488894007021, -49.9600753927381,
-48.5654335180739, -47.8918651630979, -45.1405419376454, -40.8504490169926,
-38.5815843253789, -45.001677748311, -43.0547862406721)), .Names = c("Time",
"Trace 1", "Trace 2", "Trace 3"), row.names = c(1000L, 2000L,
3000L, 4000L, 5000L, 6000L, 7000L, 8000L, 9000L, 10000L, 11000L,
12000L, 13000L, 14000L, 15000L, 16000L, 17000L, 18000L, 19000L,
20000L, 21000L, 22000L, 23000L, 24000L, 25000L, 26000L, 27000L,
28000L, 29000L, 30000L, 31000L, 32000L, 33000L, 34000L, 35000L,
36000L, 37000L, 38000L, 39000L, 40000L, 41000L, 42000L, 43000L,
44000L, 45000L, 46000L, 47000L, 48000L, 49000L, 50000L, 51000L,
52000L, 53000L, 54000L, 55000L, 56000L, 57000L, 58000L, 59000L,
60000L, 61000L, 62000L, 63000L, 64000L, 65000L, 66000L, 67000L,
68000L, 69000L, 70000L, 71000L, 72000L, 73000L, 74000L, 75000L,
76000L, 77000L, 78000L, 79000L, 80000L, 81000L, 82000L, 83000L,
84000L, 85000L, 86000L, 87000L, 88000L, 89000L, 90000L, 91000L,
92000L, 93000L, 94000L, 95000L, 96000L, 97000L, 98000L, 99000L,
100000L, 101000L, 102000L, 103000L, 104000L, 105000L, 106000L,
107000L, 108000L, 109000L, 110000L, 111000L, 112000L, 113000L,
114000L, 115000L, 116000L, 117000L, 118000L, 119000L, 120000L,
121000L, 122000L, 123000L, 124000L, 125000L, 126000L, 127000L,
128000L, 129000L, 130000L, 131000L, 132000L, 133000L, 134000L,
135000L, 136000L, 137000L, 138000L, 139000L, 140000L, 141000L,
142000L, 143000L, 144000L, 145000L, 146000L, 147000L, 148000L,
149000L, 150000L, 151000L, 152000L, 153000L, 154000L, 155000L,
156000L, 157000L, 158000L, 159000L, 160000L, 161000L, 162000L,
163000L, 164000L, 165000L, 166000L, 167000L, 168000L, 169000L,
170000L, 171000L, 172000L, 173000L, 174000L, 175000L, 176000L,
177000L, 178000L, 179000L, 180000L, 181000L, 182000L, 183000L,
184000L, 185000L, 186000L, 187000L, 188000L, 189000L, 190000L,
191000L, 192000L, 193000L, 194000L, 195000L, 196000L, 197000L,
198000L, 199000L, 200000L), class = "data.frame")
I've tried doing lm(y~x)but it doesn't seem to get the right answer (verified the right answer in Igor) and obviously this is because its a linear model and not a exponential. Any all suggestions are welcomed. I'm struggling on this!
Thanks all!
names(DF) <- make.names(names(DF))
plot(Trace.1 ~ Time, data = DF)
#remove the initial values that clearly don't follow the model
DF1 <- DF[-seq_len(which((diff(DF$Trace.1) < -1e3))),]
plot(Trace.1 ~ Time, data = DF1)
fit <- nls(Trace.1 ~ SSasymp(Time, Asym, R0, lrc), data = DF1)
summary(fit)
coef(fit)
help("SSasymp") #for an explanation of the parameters
lines(DF1$Time, predict(fit))
I apologize in advance for the specificity of this question, I'm hoping my neuroscience jargon doesn't confuse things.
I have current traces from single cell recordings and I need to fit a tau from the peak to 4 seconds after the peak out on each trace. In reality, this is just an exponential decay.
So what I'm looking for is a ggplot overlay of the fit (so really a fit line over the period requested.) something along the lines of stat_smooth()
My data looks like follow:
dput(stackover_data)
structure(list(Time = c(0.09990001, 0.19990001, 0.29990001, 0.39990001,
0.49990001, 0.59990001, 0.69990001, 0.79990001, 0.89990001, 0.99990001,
1.09990001, 1.19990001, 1.29990001, 1.39990001, 1.49990001, 1.59990001,
1.69990001, 1.79990001, 1.89990001, 1.99990001, 2.09990001, 2.19990001,
2.29990001, 2.39990001, 2.49990001, 2.59990001, 2.69990001, 2.79990001,
2.89990001, 2.99990001, 3.09990001, 3.19990001, 3.29990001, 3.39990001,
3.49990001, 3.59990001, 3.69990001, 3.79990001, 3.89990001, 3.99990001,
4.09990001, 4.19990001, 4.29990001, 4.39990001, 4.49990001, 4.59990001,
4.69990001, 4.79990001, 4.89990001, 4.99990001, 5.09990001, 5.19990001,
5.29990001, 5.39990001, 5.49990001, 5.59990001, 5.69990001, 5.79990001,
5.89990001, 5.99990001, 6.09990001, 6.19990001, 6.29990001, 6.39990001,
6.49990001, 6.59990001, 6.69990001, 6.79990001, 6.89990001, 6.99990001,
7.09990001, 7.19990001, 7.29990001, 7.39990001, 7.49990001, 7.59990001,
7.69990001, 7.79990001, 7.89990001, 7.99990001, 8.09990001, 8.19990001,
8.29990001, 8.39990001, 8.49990001, 8.59990001, 8.69990001, 8.79990001,
8.89990001, 8.99990001, 9.09990001, 9.19990001, 9.29990001, 9.39990001,
9.49990001, 9.59990001, 9.69990001, 9.79990001, 9.89990001, 9.99990001,
10.09990001, 10.19990001, 10.29990001, 10.39990001, 10.49990001,
10.59990001, 10.69990001, 10.79990001, 10.89990001, 10.99990001,
11.09990001, 11.19990001, 11.29990001, 11.39990001, 11.49990001,
11.59990001, 11.69990001, 11.79990001, 11.89990001, 11.99990001,
12.09990001, 12.19990001, 12.29990001, 12.39990001, 12.49990001,
12.59990001, 12.69990001, 12.79990001, 12.89990001, 12.99990001,
13.09990001, 13.19990001, 13.29990001, 13.39990001, 13.49990001,
13.59990001, 13.69990001, 13.79990001, 13.89990001, 13.99990001,
14.09990001, 14.19990001, 14.29990001, 14.39990001, 14.49990001,
14.59990001, 14.69990001, 14.79990001, 14.89990001, 14.99990001,
15.09990001, 15.19990001, 15.29990001, 15.39990001, 15.49990001,
15.59990001, 15.69990001, 15.79990001, 15.89990001, 15.99990001,
16.09990001, 16.19990001, 16.29990001, 16.39990001, 16.49990001,
16.59990001, 16.69990001, 16.79990001, 16.89990001, 16.99990001,
17.09990001, 17.19990001, 17.29990001, 17.39990001, 17.49990001,
17.59990001, 17.69990001, 17.79990001, 17.89990001, 17.99990001,
18.09990001, 18.19990001, 18.29990001, 18.39990001, 18.49990001,
18.59990001, 18.69990001, 18.79990001, 18.89990001, 18.99990001,
19.09990001, 19.19990001, 19.29990001, 19.39990001, 19.49990001,
19.59990001, 19.69990001, 19.79990001, 19.89990001, 19.99990001
), `Trace 1` = c(-3.08656892325052, 9.36821982641837, 8.13806079083122,
10.7039590839898, 7.25670468903547, 4.31122291688919, 1.77905971163193,
-6.27606834721828, -8.65955381985049, -10.1445673910916, -10.6649772153892,
-6.52301948183154, -11.9757817854835, -0.976254254762154, 0.237467076202677,
-11.8114896779541, -11.0022757370468, -16.1845427042923, 2.70565927469852,
-18.9048281652216, -22.153682283437, -4363.32044948884, -3470.59111611883,
-2877.51064886248, -2253.51505229908, -1876.62974792002, -1541.14546478629,
-1288.26617158403, -1232.52313999354, -1042.98549163259, -938.795146277054,
-810.913567086442, -736.390770574588, -682.044521632168, -638.575324886466,
-596.587948389699, -629.282103146111, -569.200610245336, -587.387817942122,
-521.939307772762, -533.776693538631, -514.782451411908, -445.949395199026,
-451.026618716539, -461.669600651513, -442.372477296489, -426.132547857502,
-452.471101919398, -449.377497412324, -436.817659066873, -487.464805660851,
-442.915455035179, -458.666741489705, -454.990437793055, -455.173154690614,
-426.885702219019, -456.985443835707, -408.869318611773, -442.418742303429,
-407.42593099033, -414.538253423006, -316.036755461507, -248.395118743017,
-198.708015370115, -183.88480278352, -160.041754268138, -139.020785864805,
-123.701641615743, -137.253323547789, -124.240619047461, -121.512859107349,
-113.719386521321, -102.98429740535, -118.721098087137, -109.789471870234,
-122.145399997255, -109.542467905009, -96.6725326170008, -108.16233877188,
-94.1092968366083, -88.0987495250118, -89.2245425850472, -96.3495439993499,
-110.340434956898, -98.8777265479938, -88.0674013452629, -102.755317774957,
-91.0752842494157, -107.959830375198, -78.8424385901398, -93.7743443479161,
-98.3591889757604, -72.3214090579355, -85.6296201608712, -112.596656084088,
-115.518068650615, -110.973655632476, -78.8021307215932, -85.6443070182152,
-93.1744356389988, -100.483987323044, -87.6672104421484, -83.6481077757535,
-82.1465876740715, -86.7838666454595, -75.9066755520263, -86.6416980998645,
-88.9405806921788, -74.0592581080291, -86.3433338300531, -93.5114839187431,
-91.1875849041866, -82.7018083540351, -102.859075734953, -82.1494206590809,
-79.2197323780198, -78.1558787387238, -86.2649418863144, -93.5271994290692,
-84.6678528566242, -71.4828270654073, -69.5618263581887, -67.2920558863641,
-58.1490330439793, -59.8163238740351, -73.6957797622946, -61.7673947702343,
-63.6492255747164, -64.043638367468, -73.8301991524909, -69.3055992018769,
-65.9342860783478, -74.9891521715357, -74.5779808619617, -69.9029875902787,
-69.746935396023, -71.35455530782, -73.5279471991205, -77.8000626250279,
-74.6065254864801, -64.0834786591292, -45.1346212136811, -43.1615385011179,
-42.554841323715, -50.276566542849, -56.1940469314277, -49.6368473019083,
-53.9842269565738, -54.1480156577708, -48.2160751112714, -43.2902743874793,
-46.4358385791469, -58.4459814967147, -49.0443619038303, -56.8007684056415,
-56.3345802020277, -40.1704617527471, -38.2624410086947, -39.9339279963857,
-27.7637250414188, -40.8686520798649, -30.9246406597275, -32.70021580322,
-28.7536071219405, -32.4403271637661, -30.3389870650563, -23.2866770834185,
-40.349370835336, -39.7478049410975, -38.4648592612063, -37.4191914226093,
-52.1259539643142, -38.0993644763616, -37.5493511061199, -25.0794144286873,
-32.6165749053184, -17.2616237291667, -32.1515419380766, -27.0814392288745,
-26.2440551871993, -17.5882979851251, -19.852012057918, -22.1408312947152,
-21.3395092026716, -28.0789344732281, -19.6443281019371, -26.6585778398899,
-24.2994817465722, -27.847437388783, -22.0565059455712, -24.8747905836351,
-30.7544250721296, -29.9924061105416, -22.5728855329732, -24.9767037068316,
-18.8208568928653, -10.7306087216159, -16.4281210876173, -19.3057174287183,
-15.9745523586581), `Trace 2` = c(5.94927992143286, 1.42121402161905,
6.78788136514507, 3.33970424403748, -3.73956433922802, -7.3097330836793,
-9.18242380095097, 3.29952017048882, 5.17208246028, 1.53238537592179,
6.90832098860733, 3.16748380079213, 5.49988319742749, 3.86758484926656,
11.981378748128, -0.208159377585758, -7.2781503569274, 1.71389537416221,
-4.77614396689646, -0.561115871778583, -14.5385068323377, -4944.63006397872,
-3891.23141079918, -3065.66300899921, -2527.61853791531, -2060.05223709386,
-1763.27814995793, -1451.77884853301, -1368.42409410389, -1115.3385743493,
-990.942590410893, -852.951182088994, -920.516421149141, -860.486324532349,
-751.257935697359, -768.757880984424, -713.935725464811, -715.011716140783,
-616.355633698422, -672.133827760044, -629.418857335266, -607.342562427532,
-599.657392949772, -553.265193759536, -572.878306401819, -583.153736972572,
-589.661926862833, -527.552306430313, -553.672277977987, -539.906254466276,
-525.491319968336, -564.876755897633, -505.363603776286, -575.649358123809,
-497.897048596953, -544.360501913781, -483.501835830955, -533.085516954168,
-478.007776808595, -453.333986756674, -447.125338584774, -380.268661132557,
-253.127091546528, -227.877427383931, -191.488369613846, -166.799480062655,
-171.446527550207, -156.775073739292, -148.804014388301, -151.79670303645,
-152.599645619669, -113.40778617347, -133.161498039771, -129.929086626311,
-156.827022418459, -146.618934465797, -94.20997521855, -124.514491421933,
-89.8538131075507, -91.6500841191965, -113.887119052802, -109.971138752404,
-122.532101557384, -107.714865018275, -120.258046104427, -133.296582343699,
-111.726093380577, -121.533575165211, -93.2177144160483, -97.3820468959891,
-111.294540122284, -100.913019601926, -90.8456358331201, -94.143801546019,
-89.9179990101121, -79.3197988311624, -87.0020882554621, -103.315348613102,
-92.6091563297179, -87.3900461813473, -106.61173718685, -117.18534688083,
-106.417237368758, -114.977458547258, -101.759575332149, -93.2156781534656,
-79.4452885332048, -83.8035016923155, -85.7770898942133, -80.151862850598,
-72.4787832080729, -75.4887322569144, -70.8074563945233, -76.9714115362765,
-74.2305756099668, -72.7994824773122, -86.4890002605112, -68.5241155035978,
-76.723559709657, -60.5721619300473, -64.185161607482, -81.0777473166996,
-59.7023719535658, -48.0615168793405, -63.7668096308392, -69.0412888697281,
-70.9894443873153, -48.8366672284658, -43.1003899300401, -50.5591327356918,
-65.184353265464, -52.2850955061746, -66.0074402595226, -69.3081002474696,
-69.3876079446205, -53.0919105504865, -69.3289906081927, -63.5198462609014,
-70.9316181916525, -68.753726760587, -57.1360889764902, -45.7297495037924,
-45.0225501894778, -45.7487568762369, -46.486850874694, -46.7437348684622,
-41.1396209831975, -38.5056250833188, -49.8076229939374, -46.9868446480308,
-38.4639795698766, -40.8336417843587, -49.9810721301961, -58.7820616810459,
-51.4348042554539, -56.9132093579117, -53.527642322397, -49.7620292214747,
-54.6821412367706, -66.7962460110546, -52.0791798888695, -43.3022792658504,
-47.5584813904168, -36.3352478868053, -39.8699383685185, -43.9973042247199,
-49.3191665577514, -42.0818922168797, -37.7538461067954, -29.1591142899914,
-31.6815992654018, -35.3664789321372, -29.6435427630652, -45.5528322348158,
-37.0568327299924, -29.6739001502498, -37.0952936191428, -44.1786405275305,
-50.929131558954, -31.7842698176079, -31.4300485665082, -33.960404707688,
-34.0566040705914, -41.2299273666286, -30.1094672132779, -28.6802701495622,
-37.1973277887207, -31.7114680610898, -32.5450840790921, -28.5823680785893,
-27.5268500574853, -15.2389866754792, -19.2079199191329, -28.9156878395334,
-16.5863741697772, -31.6638747515138, -34.679018047114, -32.6358481677933,
-26.3463300811612, -18.542719903123), `Trace 3` = c(4.74984904312473,
-6.83237755613888, -4.15774846489683, 0.0116753456567835, 4.3974208810568,
8.56175743226069, -0.271281722577697, -10.7982854290066, -18.31704442037,
-15.7577103721081, -20.5775487342614, -17.7307614062074, -14.2844531791322,
-5.43267691733291, 2.52755064598145, -0.79977817636393, -8.07074631773528,
-14.6631228094247, -9.70657043283029, -11.4521183354252, -10.8401698276907,
-5495.45466682125, -4324.53856178012, -3525.03574657791, -2774.37269062517,
-2175.8494069729, -1790.54797875842, -1574.71991002551, -1342.90235937345,
-1212.59715782746, -1042.56659527622, -1006.18885759924, -962.27473179716,
-870.459562514915, -847.89784200018, -770.205392585522, -709.279500982242,
-740.484104629636, -739.054108223511, -645.039948927494, -712.347389757369,
-609.734365947307, -613.120930696929, -630.805982834916, -581.828460576409,
-608.843958294059, -602.432567890874, -547.023947351376, -571.475614300591,
-561.707380041648, -535.094375334464, -578.783728464712, -535.290908514374,
-489.922390143035, -522.824368605682, -520.322259788742, -556.271028784585,
-536.474251127913, -546.811709881227, -541.304029428261, -563.702125503968,
-460.326259578172, -348.718859066443, -261.740363110507, -244.084411205067,
-182.533320181224, -205.356578064485, -157.759947570376, -162.214202634405,
-180.011654146528, -176.904445787135, -167.116794266109, -163.416234638164,
-155.821745831568, -134.760994564946, -131.702791057298, -138.598664148971,
-139.818370538306, -134.700400696581, -137.275258643317, -119.921721377143,
-118.187525264885, -118.490007264421, -130.391338033262, -138.05114734155,
-127.505366899025, -132.636233878671, -123.490111075107, -113.42942794286,
-122.40925776411, -114.48254277334, -122.78054774327, -90.9528981496313,
-96.9508488566491, -107.842237007564, -91.0521871208113, -91.0347011279681,
-85.7054147126896, -89.7363815316723, -77.6588521980595, -76.6733051688234,
-85.821987540449, -90.3216600489126, -66.1137192454051, -77.7157754909198,
-89.8637060331546, -80.9149315468257, -71.7434481475478, -97.6600776026528,
-103.405900069071, -83.7502802699966, -77.4600968679799, -77.8705250264031,
-78.9012879762878, -87.4698127883497, -92.9587198547098, -96.1222986385502,
-91.1990314226718, -83.8744116889207, -61.7720026226848, -83.6616594930646,
-91.3899042204887, -74.2552104995002, -94.3117384183545, -81.130091800586,
-87.6638514042855, -79.061763018934, -78.7472359255489, -72.9642062301352,
-83.5759889935287, -74.066505526926, -62.6191749912397, -75.7150349466975,
-71.8704211546101, -68.4798935028027, -77.0470926055507, -76.6024078920054,
-66.5760393943856, -77.8848130194484, -77.5151449198185, -65.5482413040489,
-47.6087884294611, -35.6865291924552, -47.9648700291328, -62.8176093048552,
-48.969044800301, -46.853228499419, -52.562083189922, -58.7677963803543,
-63.6615259259388, -51.7351340481719, -54.5967510836755, -54.1294816023731,
-47.7621427567758, -40.9244974392914, -50.5286026229961, -45.9961132210891,
-43.143208905187, -43.8395623054062, -49.3326618452772, -48.308625153823,
-44.9833219970479, -32.4206723427767, -45.6628898472981, -41.4748918817457,
-43.6438737242104, -42.9844450147366, -30.3729810016378, -44.9519080631137,
-45.8241111308495, -40.4766896430523, -40.3007872752484, -46.4770136239361,
-49.3759491415156, -46.0844075639024, -34.3090862431636, -26.6047381016158,
-28.9119815377167, -33.1619464253006, -37.3813739690468, -46.5001186141503,
-54.5235621985407, -44.5233400090119, -55.2272273265567, -48.088333647706,
-55.3522589332341, -52.8146474401922, -60.1877313269088, -48.2512741329643,
-34.8540879147507, -55.2019332852645, -50.8488894007021, -49.9600753927381,
-48.5654335180739, -47.8918651630979, -45.1405419376454, -40.8504490169926,
-38.5815843253789, -45.001677748311, -43.0547862406721)), .Names = c("Time",
"Trace 1", "Trace 2", "Trace 3"), row.names = c(1000L, 2000L,
3000L, 4000L, 5000L, 6000L, 7000L, 8000L, 9000L, 10000L, 11000L,
12000L, 13000L, 14000L, 15000L, 16000L, 17000L, 18000L, 19000L,
20000L, 21000L, 22000L, 23000L, 24000L, 25000L, 26000L, 27000L,
28000L, 29000L, 30000L, 31000L, 32000L, 33000L, 34000L, 35000L,
36000L, 37000L, 38000L, 39000L, 40000L, 41000L, 42000L, 43000L,
44000L, 45000L, 46000L, 47000L, 48000L, 49000L, 50000L, 51000L,
52000L, 53000L, 54000L, 55000L, 56000L, 57000L, 58000L, 59000L,
60000L, 61000L, 62000L, 63000L, 64000L, 65000L, 66000L, 67000L,
68000L, 69000L, 70000L, 71000L, 72000L, 73000L, 74000L, 75000L,
76000L, 77000L, 78000L, 79000L, 80000L, 81000L, 82000L, 83000L,
84000L, 85000L, 86000L, 87000L, 88000L, 89000L, 90000L, 91000L,
92000L, 93000L, 94000L, 95000L, 96000L, 97000L, 98000L, 99000L,
100000L, 101000L, 102000L, 103000L, 104000L, 105000L, 106000L,
107000L, 108000L, 109000L, 110000L, 111000L, 112000L, 113000L,
114000L, 115000L, 116000L, 117000L, 118000L, 119000L, 120000L,
121000L, 122000L, 123000L, 124000L, 125000L, 126000L, 127000L,
128000L, 129000L, 130000L, 131000L, 132000L, 133000L, 134000L,
135000L, 136000L, 137000L, 138000L, 139000L, 140000L, 141000L,
142000L, 143000L, 144000L, 145000L, 146000L, 147000L, 148000L,
149000L, 150000L, 151000L, 152000L, 153000L, 154000L, 155000L,
156000L, 157000L, 158000L, 159000L, 160000L, 161000L, 162000L,
163000L, 164000L, 165000L, 166000L, 167000L, 168000L, 169000L,
170000L, 171000L, 172000L, 173000L, 174000L, 175000L, 176000L,
177000L, 178000L, 179000L, 180000L, 181000L, 182000L, 183000L,
184000L, 185000L, 186000L, 187000L, 188000L, 189000L, 190000L,
191000L, 192000L, 193000L, 194000L, 195000L, 196000L, 197000L,
198000L, 199000L, 200000L), class = "data.frame")
If you
plot(stackover_data$Time,stackover_data$'Trace 1',type='l'
you can see from time = 2 is where the peak is and I'm looking to fit it from 2 seconds to 6 seconds in ggplot. ( yes I know I used base graphics for the example, its just easier to plot one graph in that vs ggplot2.
Thanks all!
library(ggplot)
# rename the variables, so there's no space in them (is easier in aes())
names(neurons)[2:4] <- paste0('Trace', 1:3)
Now we make a plot, where we use a selection of the data based on the lowest point in the trace:
ggplot(neurons, aes(x = Time, y = Trace1)) +
geom_line() +
geom_smooth(data = subset(neurons, Time >= Time[which.min(Trace1)] &
Time < Time[which.min(Trace1)] + 4)) +
theme_classic() +
xlim(0, 10)
Results in:
Note that the ggplot equivalent of your simple base plot is this:
qplot(neurons$Time, neurons$Trace1, geom = 'line')
Basically, just as easy ;)