JustPaste.it

sensWth<-
structure(list(W = c(5665987.3454, 5872688.2009, 6090384.3529,
6304430.6732, 6506105.9085, 6697331.4024, 6880658.1781, 7054469.5193,
7224339.2914, 7388265.7955, 7547861.6658, 7703290.8471, 7853705.3315,
8001215.3499, 8144108.3417, 5868344.0375, 6082427.0642, 6307898.0787,
6529588.9092, 6738466.8349, 6936521.8118, 7126395.9718, 7306414.8578,
7482351.4095, 7652132.433, 7817428.1539, 7978408.3776, 8134194.8097,
8286973.0438, 8434969.3576, 6070700.728, 6292165.9286, 6525411.8078,
6754747.1496, 6970827.7608, 7175712.2165, 7372133.7647, 7558360.1983,
7740363.5297, 7915999.0664, 8086994.6407, 8253525.9088, 8414684.2843,
8572730.7345, 8725830.3682, 6273057.4189, 6501904.7916, 6742925.5327,
6979905.3854, 7203188.6862, 7414902.6277, 7617871.5595, 7810305.5402,
7998375.6481, 8179865.7051, 8356561.1316, 8528643.4385, 8695173.7599,
8858488.4261, 9016691.3802, 6475414.1093, 6711643.6574, 6960439.2596,
7205063.6259, 7435549.6119, 7654093.0331, 7863609.349, 8062250.8783,
8256387.763, 8443732.3373, 8626127.6195, 8803760.9676, 8975663.2369,
9144246.1189, 9307552.3946, 6677770.7993, 6921382.5218, 7177952.9843,
7430221.8645, 7667910.5355, 7893283.4418, 8109347.1418, 8314196.2205,
8514399.8819, 8707598.9732, 8895694.1055, 9078878.5007, 9256152.7126,
9430003.8076, 9598413.4057, 6880127.4911, 7131121.3839, 7395466.7129,
7655380.1029, 7900271.4625, 8132473.8454, 8355084.9333, 8566141.558,
8772412.0012, 8971465.6087, 9165260.5943, 9353996.0279, 9536642.1866,
9715761.4983, 9889274.4159, 7082484.1829, 7340860.2524, 7612980.4413,
7880538.3386, 8132632.3881, 8371664.2532, 8600822.7237, 8818086.9011,
9030424.117, 9235332.2415, 9434827.0857, 9629113.5598, 9817131.6631,
10001519.1915, 10180135.4302, 7284840.8715, 7550599.1146, 7830494.1691,
8105696.5763, 8364993.3112, 8610854.6608, 8846560.5191, 9070032.238,
9288436.2322, 9499198.8814, 9704393.5688, 9904231.0891, 10097621.1421,
10287276.881, 10470996.4438, 7487197.5637, 7760337.9768, 8048007.8942,
8330854.8131, 8597354.2388, 8850045.0703, 9092298.3092, 9321977.5759,
9546448.35, 9763065.5142, 9973960.0565, 10179348.6194, 10378110.6197,
10573034.5731, 10761857.4552, 7689554.2554, 7970076.8434, 8265521.6201,
8556013.05, 8829715.1641, 9089235.4752, 9338036.1035, 9573922.9204,
9804460.4684, 10026932.1525, 10243526.5475, 10454466.1529, 10658600.0927,
10858792.2647, 11052718.468, 7891910.9467, 8179815.7083, 8483035.3492,
8781171.2923, 9062076.0875, 9328425.8821, 9583773.8971, 9825868.2549,
10062472.5889, 10290798.7875, 10513093.0339, 10729583.6786, 10939089.5695,
11144549.957, 11343579.4827, 8013324.9598, 8305659.029, 8613543.5823,
8916266.2359, 9201492.6458, 9471940.1286, 9731216.5693, 9977035.4639,
10217279.8573, 10449118.7667, 10674832.9283, 10894654.1968, 11107383.2569,
11316004.57, 11518096.0894, 8094267.6373, 8389554.5691, 8700549.0751,
9006329.5283, 9294437.0128, 9567616.2919, 9829511.6861, 10077813.5965,
10320484.7081, 10554665.4269, 10782659.5238, 11004701.2105, 11219579.042,
11430307.6475, 11634440.4931, 8175210.3114, 8473450.1176, 8787554.5639,
9096392.8265, 9387381.3836, 9663292.4523, 9927806.8018, 10178591.7345,
10423689.5526, 10660212.0768, 10890486.1181, 11114748.2252, 11331774.8373,
11544610.7229, 11750784.897, 8296624.3288, 8599293.4337, 8918062.7988,
9231487.7689, 9526797.9398, 9806806.6986, 10075249.4773, 10329758.938,
10578496.8221, 10818532.0551, 11052226.0114, 11279818.7392, 11500068.521,
11716065.3363, 11925301.5068, 8498981.0166, 8809032.3015, 9135576.5301,
9456646.0082, 9759158.8669, 10045997.1047, 10320987.2713, 10581704.278,
10836508.9401, 11082398.6907, 11321792.4962, 11554936.2703, 11780557.9978,
12001823.0299, 12216162.5167, 8701337.7086, 9018771.16, 9353090.2567,
9681804.2452, 9991519.7891, 10285187.5142, 10566725.0652, 10833649.6176,
11094521.0545, 11346265.3293, 11591358.9874, 11830053.7977, 12061047.4748,
12287580.721, 12507023.5283, 8903694.3987, 9228510.0277, 9570603.9818,
9906962.4846, 10223880.7154, 10524377.9188, 10812462.8545, 11085594.9565,
11352533.1744, 11610131.9644, 11860925.4727, 12105171.3314, 12341536.9491,
12573338.4101, 12797884.5431, 9106051.0917, 9438248.8941, 9788117.7081,
10132120.7214, 10456241.6422, 10763568.3264, 11058200.6467, 11337540.2972,
11610545.2911, 11873998.5996, 12130491.9645, 12380288.8639, 12622026.4284,
12859096.1009, 13088745.5527, 9308407.7803, 9647987.7565, 10005631.4349,
10357278.9613, 10688602.5676, 11002758.7355, 11303938.4391, 11589485.636,
11868557.4106, 12137865.2368, 12400058.4525, 12655406.3911, 12902515.9034,
13144853.794, 13379606.5653, 9510764.4725, 9857726.6209, 10223145.1599,
10582437.1964, 10920963.4899, 11241949.1401, 11549676.2347, 11841430.9779,
12126569.5278, 12401731.8694, 12669624.9388, 12930523.923, 13183005.3779,
13430611.4858, 13670467.5806, 9713121.1659, 10067465.4835, 10440658.8884,
10807595.4348, 11153324.4177, 11481139.5513, 11795414.0241, 12093376.3166,
12384581.6477, 12665598.5055, 12939191.4293, 13205641.4523, 13463494.8554,
13716369.1739, 13961328.5909, 9915477.8537, 10277204.3495, 10658172.6177,
11032753.6754, 11385685.3423, 11720329.9559, 12041151.8168, 12345321.6608,
12642593.7646, 12929465.1432, 13208757.9169, 13480758.9817, 13743984.3321,
14002126.8646, 14252189.6036, 10117834.5457, 10486943.2147, 10875686.3432,
11257911.9133, 11618046.2671, 11959520.364, 12286889.6089, 12597266.9975,
12900605.8814, 13193331.78, 13478324.4022, 13755876.5114, 14024473.809,
14287884.5577, 14543050.6153, 10320191.2366, 10696682.0791, 11093200.0679,
11483070.1503, 11850407.1921, 12198710.7711, 12532627.4011, 12849212.3386,
13158617.9995, 13457198.4139, 13747890.8924, 14030994.0442, 14304963.2859,
14573642.2515, 14833911.6273, 10522547.9279, 10906420.9451, 11310713.798,
11708228.3888, 12082768.1196, 12437901.1781, 12778365.1912, 13101157.678,
13416630.1166, 13721065.0495, 14017457.3796, 14306111.5775, 14585452.7598,
14859399.9412, 15124772.6396), t = structure(c(1L, 2L, 3L, 4L,
5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 1L, 2L, 3L,
4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 1L, 2L,
3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 1L,
2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L,
1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L,
15L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L,
14L, 15L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L,
13L, 14L, 15L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L,
12L, 13L, 14L, 15L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L,
11L, 12L, 13L, 14L, 15L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L,
10L, 11L, 12L, 13L, 14L, 15L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L,
9L, 10L, 11L, 12L, 13L, 14L, 15L, 1L, 2L, 3L, 4L, 5L, 6L, 7L,
8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 1L, 2L, 3L, 4L, 5L, 6L,
7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 1L, 2L, 3L, 4L, 5L,
6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 1L, 2L, 3L, 4L,
5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 1L, 2L, 3L,
4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 1L, 2L,
3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 1L,
2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L,
1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L,
15L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L,
14L, 15L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L,
13L, 14L, 15L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L,
12L, 13L, 14L, 15L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L,
11L, 12L, 13L, 14L, 15L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L,
10L, 11L, 12L, 13L, 14L, 15L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L,
9L, 10L, 11L, 12L, 13L, 14L, 15L, 1L, 2L, 3L, 4L, 5L, 6L, 7L,
8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 1L, 2L, 3L, 4L, 5L, 6L,
7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L), .Label = c("t_m1_tc",
"t_m0.5_tc", "t_p0_tc", "t_p0.5_tc", "t_p1_tc", "t_p1.5_tc",
"t_p2_tc", "t_p2.5_tc", "t_p3_tc", "t_p3.5_tc", "t_p4_tc", "t_p4.5_tc",
"t_p5_tc", "t_p5.5_tc", "t_p6_tc"), class = "factor"), p = structure(c(1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 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, 5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 6L, 6L,
6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 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, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L,
9L, 9L, 9L, 9L, 9L, 9L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L,
10L, 10L, 10L, 10L, 10L, 10L, 10L, 11L, 11L, 11L, 11L, 11L, 11L,
11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 12L, 12L, 12L, 12L,
12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 13L, 13L,
13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L,
14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L,
14L, 14L, 15L, 15L, 15L, 15L, 15L, 15L, 15L, 15L, 15L, 15L, 15L,
15L, 15L, 15L, 15L, 16L, 16L, 16L, 16L, 16L, 16L, 16L, 16L, 16L,
16L, 16L, 16L, 16L, 16L, 16L, 17L, 17L, 17L, 17L, 17L, 17L, 17L,
17L, 17L, 17L, 17L, 17L, 17L, 17L, 17L, 18L, 18L, 18L, 18L, 18L,
18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 19L, 19L, 19L,
19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 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, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L,
22L, 22L, 22L, 22L, 22L, 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, 25L, 25L, 25L, 25L,
25L, 25L, 25L, 25L, 25L, 25L, 25L, 25L, 25L, 25L, 25L, 26L, 26L,
26L, 26L, 26L, 26L, 26L, 26L, 26L, 26L, 26L, 26L, 26L, 26L, 26L,
27L, 27L, 27L, 27L, 27L, 27L, 27L, 27L, 27L, 27L, 27L, 27L, 27L,
27L, 27L), .Label = c("h_30", "h_27.5", "h_25", "h_22.5", "h_20",
"h_17.5", "h_15", "h_12.5", "h_10", "h_7.5", "h_5", "h_2.5",
"h_1", "h0", "h1", "h2.5", "h5", "h7.5", "h10", "h12.5", "h15",
"h17.5", "h20", "h22.5", "h25", "h27.5", "h30"), class = "factor"),
tt = c(-1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6, -1, -0.5, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5,
5, 5.5, 6), hh = c(-30, -30, -30, -30, -30, -30, -30, -30,
-30, -30, -30, -30, -30, -30, -30, -27.5, -27.5, -27.5, -27.5,
-27.5, -27.5, -27.5, -27.5, -27.5, -27.5, -27.5, -27.5, -27.5,
-27.5, -27.5, -25, -25, -25, -25, -25, -25, -25, -25, -25,
-25, -25, -25, -25, -25, -25, -22.5, -22.5, -22.5, -22.5,
-22.5, -22.5, -22.5, -22.5, -22.5, -22.5, -22.5, -22.5, -22.5,
-22.5, -22.5, -20, -20, -20, -20, -20, -20, -20, -20, -20,
-20, -20, -20, -20, -20, -20, -17.5, -17.5, -17.5, -17.5,
-17.5, -17.5, -17.5, -17.5, -17.5, -17.5, -17.5, -17.5, -17.5,
-17.5, -17.5, -15, -15, -15, -15, -15, -15, -15, -15, -15,
-15, -15, -15, -15, -15, -15, -12.5, -12.5, -12.5, -12.5,
-12.5, -12.5, -12.5, -12.5, -12.5, -12.5, -12.5, -12.5, -12.5,
-12.5, -12.5, -10, -10, -10, -10, -10, -10, -10, -10, -10,
-10, -10, -10, -10, -10, -10, -7.5, -7.5, -7.5, -7.5, -7.5,
-7.5, -7.5, -7.5, -7.5, -7.5, -7.5, -7.5, -7.5, -7.5, -7.5,
-5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5,
-2.5, -2.5, -2.5, -2.5, -2.5, -2.5, -2.5, -2.5, -2.5, -2.5,
-2.5, -2.5, -2.5, -2.5, -2.5, -1, -1, -1, -1, -1, -1, -1,
-1, -1, -1, -1, -1, -1, -1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5,
2.5, 2.5, 2.5, 2.5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
5, 5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
7.5, 7.5, 7.5, 7.5, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
10, 10, 10, 10, 10, 12.5, 12.5, 12.5, 12.5, 12.5, 12.5, 12.5,
12.5, 12.5, 12.5, 12.5, 12.5, 12.5, 12.5, 12.5, 15, 15, 15,
15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 17.5, 17.5,
17.5, 17.5, 17.5, 17.5, 17.5, 17.5, 17.5, 17.5, 17.5, 17.5,
17.5, 17.5, 17.5, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20,
20, 20, 20, 20, 20, 22.5, 22.5, 22.5, 22.5, 22.5, 22.5, 22.5,
22.5, 22.5, 22.5, 22.5, 22.5, 22.5, 22.5, 22.5, 25, 25, 25,
25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 27.5, 27.5,
27.5, 27.5, 27.5, 27.5, 27.5, 27.5, 27.5, 27.5, 27.5, 27.5,
27.5, 27.5, 27.5, 30, 30, 30, 30, 30, 30, 30, 30, 30, 30,
30, 30, 30, 30, 30), pChange = c(-34.877818669911, -32.502096704368,
-29.9999999962071, -27.5398527290355, -25.2218928674312,
-23.024037395904, -20.9169660591689, -18.9192606304687, -16.9668577345854,
-15.0827639528591, -13.2484444297758, -11.4620148612694,
-9.73322184945282, -8.03781139746013, -6.39546686763104,
-32.5520264658406, -30.0914573126514, -27.5000000085914,
-24.9519903532645, -22.5512461715233, -20.2748958493717,
-18.0925719711765, -16.0235199556526, -14.001388361642, -12.0500055002327,
-10.1501745875716, -8.29994396062525, -6.50940832022707,
-4.75344748624667, -3.05244778470429, -30.2262342801598,
-27.6808179082918, -24.999999983047, -22.3641279269221, -19.8805994813623,
-17.5257543568591, -15.2681778923789, -13.1277792578496,
-11.0359189645622, -9.01724709472985, -7.05190476030902,
-5.13787305193564, -3.28559483237802, -1.46908361181251,
0.290571237306776, -27.9004420898815, -25.2701785200232,
-22.5000000057755, -19.7762655534498, -17.2099527969479,
-14.7766127896385, -12.4437837917437, -10.2320385439555,
-8.07044958817072, -5.98448862831125, -3.95363488592294,
-1.97580216048633, -0.0617813330354367, 1.81528027296583,
3.63359027540877, -25.5746499053501, -22.8595390995727, -20.0000000055169,
-17.188403125958, -14.5393061090855, -12.0274712890804, -9.61938975202412,
-7.33629787373693, -5.10498025200661, -2.95173023660058,
-0.85536504601744, 1.18626872406686, 3.16203218239808, 5.09964417153639,
6.97660934109523, -23.248857725416, -20.4488996952132, -17.500000030544,
-14.6005407203039, -11.8686594453595, -9.27832975059359,
-6.79499567437592, -4.44055715639485, -2.13951086986841,
0.0810281976361441, 2.24290477205035, 4.34833965459419, 6.38584568289002,
8.38400802298351, 10.319628368853, -20.9230655247936, -18.0382603172888,
-15.0000000107464, -12.0126783169485, -9.19801274255557,
-6.52918827072382, -3.97060161166931, -1.54481649307235,
0.825958516867222, 3.11378662727544, 5.34117462230003, 7.51041051730968,
9.60965916384296, 11.6683718974177, 13.6626473862667, -18.5972733241712,
-15.6276208658058, -12.4999999932475, -9.42481594462559,
-6.52736605584254, -3.78004674258123, -1.14620756160556,
1.3509242346139, 3.79142786337551, 6.14654502588222, 8.43944450243287,
10.6724814340448, 12.8334726735297, 14.9527358005857, 17.0056664508038,
-16.2714811603281, -13.216981486732, -9.99999998264477, -6.83695354931565,
-3.85671939786331, -1.03090521673734, 1.6781865459258, 4.24666489104027,
6.75689720298767, 9.17930350609306, 11.5377142871694, 13.8345523208967,
16.0572862119503, 18.2370996612276, 20.3486855072954, -13.9456889551083,
-10.8063421076582, -7.50000000307451, -4.24909116434988,
-1.18607268816322, 1.71823633094423, 4.50258059254149, 7.14240555896018,
9.72236657248299, 12.2120619046998, 14.6359841247762, 16.9966232192421,
19.2810997342799, 21.5214635517527, 23.6917045385013, -11.6198967556353,
-8.39570267801293, -5.00000001430943, -1.66122877823475,
1.48457399510176, 4.46737782575559, 7.32697468742998, 10.0381463027374,
12.6878359488744, 15.244820366521, 17.7342540003117, 20.1586941543668,
22.5049132037394, 24.805827436531, 27.0347235857981, -9.29410456075963,
-5.98506326790662, -2.49999998876508, 0.926633669945423,
4.15522065652904, 7.216519343554, 10.151368774273, 12.9338869315793,
15.6533053494023, 18.2775787904136, 20.8325238229769, 23.320764999842,
25.7287267168741, 28.0901913293548, 30.3777426549326, -7.89862926314338,
-4.53867960161423, -1.00000002355023, 2.47935111839503, 5.75760870234783,
8.86600428135778, 11.8460051808645, 14.6713314042807, 17.4325869448943,
20.0972338240607, 22.6914857459994, 25.2180075390792, 27.6630148399265,
30.0608096376945, 32.3835540720471, -6.96831237392947, -3.57442390492378,
0, 3.51449604571637, 6.82586733979402, 9.96566089468365,
12.9757628082458, 15.8296276420252, 18.6187747349886, 21.3103372648776,
23.9307936858694, 26.4828359165771, 28.9525401805868, 31.3745552014903,
33.720761674645, -6.03799552379355, -2.6101681116877, 0.999999977576133,
4.54964103970015, 7.89412602091561, 11.0653174746783, 14.1055204229843,
16.9879239418348, 19.8049624526737, 22.5234405873112, 25.1701016119472,
27.7476643055686, 30.2420656384811, 32.6883007411496, 35.0579692795416,
-4.64252017675512, -1.16378449826552, 2.49999996347934, 6.10235847435753,
9.49651404259799, 12.7148024101834, 15.8001568675044, 18.7253683513218,
21.5842440608085, 24.3430956106142, 27.0290635223269, 29.6449067965329,
32.1763537190074, 34.6589190540867, 37.063780732286, -2.31672802210684,
1.246854945172, 5.00000001430944, 8.69022088805711, 12.1671607465513,
15.4639439187869, 18.6245509589448, 21.6211090433782, 24.5497134326025,
27.3758540414028, 30.1273333266024, 32.8069777040732, 35.4001672321422,
37.9432829618521, 40.4067997462516, 0.00906418081425968,
3.65749428171973, 7.50000001112, 11.2780832753216, 14.8378073941864,
18.2130854664685, 21.4489450492359, 24.5168497308371, 27.5151827630199,
30.4086125066721, 33.2256032044366, 35.9690485690874, 38.6239807475757,
41.2276468408837, 43.7498187797562, 2.33485636189767, 6.0681337240079,
9.99999999069025, 13.8659456901705, 17.5084540889449, 20.9622269578318,
24.2733390866567, 27.4125904102505, 30.4806521566516, 33.441370931714,
36.3238730144589, 39.1311195065108, 41.8477942319767, 44.5120106969282,
47.09283785004, 4.66064857631231, 8.47877315135449, 12.4999999840527,
16.4538080751363, 20.1791007894501, 23.7113684836757, 27.0977331574089,
30.3083311103523, 33.4461215135041, 36.4741293579052, 39.4221428991891,
42.293190430142, 45.0716077738454, 47.7963745725117, 50.4358568605576,
6.98644074015542, 10.889412532727, 14.999999983162, 19.041670495732,
22.8497474738645, 26.4605100267599, 29.9221272304596, 33.2040717886164,
36.4115909025384, 39.5068878070835, 42.520412740244, 45.4552612928575,
48.2954212662918, 51.0807384745304, 53.7788759055557, 9.31223294537522,
13.3000519370865, 17.4999999615829, 21.6295328611588, 25.5203941226489,
29.2096515181232, 32.7465213402897, 36.0998125025104, 39.3770602651376,
42.5396462033916, 45.6186825617598, 48.6173322095926, 51.5192347529915,
54.3651023616074, 57.1218949815863, 11.6380251643873, 15.7106913207577,
19.9999999802311, 24.2173952645142, 28.1910408346477, 31.9587930853438,
35.5709153788599, 38.9955531796251, 42.3425296587694, 45.572404639927,
48.7169524315485, 51.7794030964445, 54.7430482741718, 57.6494662061584,
60.4649140001493, 13.9638173190356, 18.1213307435069, 22.5000000080742,
26.8052576931554, 30.8616875098672, 34.7079345767071, 38.3953094553588,
41.8912939199543, 45.3079990179205, 48.6051630948521, 51.815222268006,
54.9414739844457, 57.9668617861573, 60.9338300805926, 63.8079330462968,
16.2896095219567, 20.5319701570613, 24.9999999922419, 29.393120090764,
33.5323341873854, 37.4570761082977, 41.2197035249615, 44.7870345740819,
48.2734683759223, 51.637921539433, 54.9134920780283, 58.1035448758951,
61.1906753004414, 64.2181939826112, 67.1509520809507, 18.6154017122349,
22.9426095614208, 27.4999999672147, 31.9809824780285, 36.2029808672023,
40.2062176283948, 44.0440975957136, 47.6827752787811, 51.2389377488657,
54.6706799506827, 58.0117619443689, 61.2656158029743, 64.4144888147256,
67.5025578926753, 70.4939711190527, 20.9411939071105, 25.35324898417,
30.0000000042526, 34.5688448825333, 38.8736275757531, 42.9553591473426,
46.8684916423293, 50.5785159639413, 54.2044071103156, 57.7034383814713,
61.1100317762289, 64.4276867358003, 67.6383022945292, 70.7869217556159,
73.8369901606028)), class = "data.frame", row.names = c(NA,
-405L))

 

 

delta <- structure(list(rcp = c("rcp45", "rcp45", "rcp45", "rcp45", "rcp45",
"rcp45", "rcp45", "rcp45", "rcp45", "rcp45", "rcp45", "rcp45",
"rcp45", "rcp45", "rcp45", "rcp45", "rcp45", "rcp45", "rcp45",
"rcp45", "rcp45", "rcp45", "rcp45", "rcp45", "rcp45", "rcp45",
"rcp45", "rcp45", "rcp45", "rcp45", "rcp45", "rcp45", "rcp45",
"rcp45", "rcp45", "rcp45", "rcp45", "rcp45", "rcp45", "rcp45",
"rcp45", "rcp45", "rcp45", "rcp45", "rcp45", "rcp45", "rcp45",
"rcp45", "rcp45", "rcp45", "rcp45", "rcp45", "rcp45", "rcp45",
"rcp45", "rcp45", "rcp45", "rcp45", "rcp45", "rcp45", "rcp45",
"rcp45", "rcp45", "rcp45", "rcp45", "rcp45", "rcp45", "rcp45",
"rcp85", "rcp85", "rcp85", "rcp85", "rcp85", "rcp85", "rcp85",
"rcp85", "rcp85", "rcp85", "rcp85", "rcp85", "rcp85", "rcp85",
"rcp85", "rcp85", "rcp85", "rcp85", "rcp85", "rcp85", "rcp85",
"rcp85", "rcp85", "rcp85", "rcp85", "rcp85", "rcp85", "rcp85",
"rcp85", "rcp85", "rcp85", "rcp85", "rcp85", "rcp85", "rcp85",
"rcp85", "rcp85", "rcp85", "rcp85", "rcp85", "rcp85", "rcp85",
"rcp85", "rcp85", "rcp85", "rcp85", "rcp85", "rcp85", "rcp85",
"rcp85", "rcp85", "rcp85", "rcp85", "rcp85", "rcp85", "rcp85",
"rcp85", "rcp85", "rcp85", "rcp85", "rcp85", "rcp85", "rcp85",
"rcp85", "rcp85", "rcp85", "rcp85", "rcp85"), gcm = c("ACCESS1-0",
"ACCESS1-3", "bcc-csm1-1", "bcc-csm1-1-m", "CanESM2", "CCSM4",
"CESM1-BGC", "CESM1-CAM5", "CMCC-CM", "CMCC-CMS", "CNRM-CM5",
"CSIRO-Mk3-6-0", "FGOALS-g2", "FIO-ESM", "GFDL-CM3", "GFDL-ESM2G",
"GFDL-ESM2M", "GISS-E2-H", "GISS-E2-R", "HadGEM2-AO", "HadGEM2-CC",
"HadGEM2-ES", "inmcm4", "IPSL-CM5A-LR", "IPSL-CM5A-MR", "IPSL-CM5B-LR",
"MIROC-ESM", "MIROC-ESM-CHEM", "MIROC5", "MPI-ESM-LR", "MPI-ESM-MR",
"MRI-CGCM3", "NorESM1-M", "NorESM1-ME", "ACCESS1-0", "ACCESS1-3",
"bcc-csm1-1", "bcc-csm1-1-m", "CanESM2", "CCSM4", "CESM1-BGC",
"CESM1-CAM5", "CMCC-CM", "CMCC-CMS", "CNRM-CM5", "CSIRO-Mk3-6-0",
"FGOALS-g2", "FIO-ESM", "GFDL-CM3", "GFDL-ESM2G", "GFDL-ESM2M",
"GISS-E2-H", "GISS-E2-R", "HadGEM2-AO", "HadGEM2-CC", "HadGEM2-ES",
"inmcm4", "IPSL-CM5A-LR", "IPSL-CM5A-MR", "IPSL-CM5B-LR", "MIROC-ESM",
"MIROC-ESM-CHEM", "MIROC5", "MPI-ESM-LR", "MPI-ESM-MR", "MRI-CGCM3",
"NorESM1-M", "NorESM1-ME", "ACCESS1-0", "ACCESS1-3", "bcc-csm1-1",
"bcc-csm1-1-m", "CanESM2", "CCSM4", "CESM1-BGC", "CESM1-CAM5",
"CMCC-CM", "CMCC-CMS", "CNRM-CM5", "CSIRO-Mk3-6-0", "FGOALS-g2",
"FIO-ESM", "GFDL-CM3", "GFDL-ESM2G", "GFDL-ESM2M", "GISS-E2-H",
"GISS-E2-R", "HadGEM2-AO", "HadGEM2-CC", "HadGEM2-ES", "inmcm4",
"IPSL-CM5A-LR", "IPSL-CM5A-MR", "IPSL-CM5B-LR", "MIROC-ESM",
"MIROC-ESM-CHEM", "MIROC5", "MPI-ESM-LR", "MPI-ESM-MR", "MRI-CGCM3",
"NorESM1-M", "NorESM1-ME", "ACCESS1-0", "ACCESS1-3", "bcc-csm1-1",
"bcc-csm1-1-m", "CanESM2", "CCSM4", "CESM1-BGC", "CESM1-CAM5",
"CMCC-CM", "CMCC-CMS", "CNRM-CM5", "CSIRO-Mk3-6-0", "FGOALS-g2",
"FIO-ESM", "GFDL-CM3", "GFDL-ESM2G", "GFDL-ESM2M", "GISS-E2-H",
"GISS-E2-R", "HadGEM2-AO", "HadGEM2-CC", "HadGEM2-ES", "inmcm4",
"IPSL-CM5A-LR", "IPSL-CM5A-MR", "IPSL-CM5B-LR", "MIROC-ESM",
"MIROC-ESM-CHEM", "MIROC5", "MPI-ESM-LR", "MPI-ESM-MR", "MRI-CGCM3",
"NorESM1-M", "NorESM1-ME"), future = c("Near", "Near", "Near",
"Near", "Near", "Near", "Near", "Near", "Near", "Near", "Near",
"Near", "Near", "Near", "Near", "Near", "Near", "Near", "Near",
"Near", "Near", "Near", "Near", "Near", "Near", "Near", "Near",
"Near", "Near", "Near", "Near", "Near", "Near", "Near", "Far",
"Far", "Far", "Far", "Far", "Far", "Far", "Far", "Far", "Far",
"Far", "Far", "Far", "Far", "Far", "Far", "Far", "Far", "Far",
"Far", "Far", "Far", "Far", "Far", "Far", "Far", "Far", "Far",
"Far", "Far", "Far", "Far", "Far", "Far", "Near", "Near", "Near",
"Near", "Near", "Near", "Near", "Near", "Near", "Near", "Near",
"Near", "Near", "Near", "Near", "Near", "Near", "Near", "Near",
"Near", "Near", "Near", "Near", "Near", "Near", "Near", "Near",
"Near", "Near", "Near", "Near", "Near", "Near", "Near", "Far",
"Far", "Far", "Far", "Far", "Far", "Far", "Far", "Far", "Far",
"Far", "Far", "Far", "Far", "Far", "Far", "Far", "Far", "Far",
"Far", "Far", "Far", "Far", "Far", "Far", "Far", "Far", "Far",
"Far", "Far", "Far", "Far", "Far", "Far"), dP = c(-3.25601208225227,
4.25422575313546, 3.47856932609794, 6.07325308463769, 6.13007803837502,
3.7621513993393, 13.713339640492, 10.2070499979329, 0.577718791474743,
-1.09176058823822, 8.96908501620543, 7.59563902682608, -1.09812004690998,
-3.71905411491341, 5.24520341268675, 9.16074352957066, 6.15524950613622,
0.807207207207211, 1.34566862910008, -4.51506437950621, -4.08101204233061,
-5.39728674911793, -2.69215440031884, 1.15328520828994, 8.40475894717955,
-4.07115502025864, 6.08013359849063, 7.057408574626, 2.01259854997822,
-1.90455419455918, -0.902599845908517, 8.59364585647634, -0.164483012875046,
2.82756313256128, 0.271076172404446, 9.65568966235642, 15.8954542366997,
3.56600434296903, 20.642681394342, 9.85687780515799, 14.330867807689,
23.8083068053264, -6.07905899497749, -4.29246854228087, 21.4528795811518,
12.3292738184407, 5.01202364452658, -9.56785087232639, 38.1477228727452,
12.0792602972542, 12.590308517961, 14.0396396396396, 12.6072329688814,
-0.455401787611376, -1.24072497263106, -2.20760639134636, 2.4530188977207,
3.8655276800706, 18.7349040067354, 8.45867791530222, 12.7983019733823,
14.9721511171711, 12.5, -3.69382922660589, -0.862733152962687,
4.44493816188995, 7.86531299037679, 8.55302678904106, -6.42992680943345,
3.39032486562313, 8.2779232737345, -1.03490692741217, 9.54729294486627,
4.00816195557822, 9.2332501935151, 11.8442069973804, 2.4744506280926,
-2.61853507206316, 12.8153826975816, 14.145421255583, 3.61420093109208,
-6.14250934985085, 5.13571972175119, 4.99119267500688, 9.6092180443554,
-3.35135135135136, 5.4247266610597, 7.4904032749267, -5.6934138858479,
-8.46111781142083, -1.11671767722722, 4.7687914404422, 6.38413868300829,
-2.77330669182809, 8.93503645913007, 4.87194659188845, 2.79604611544708,
-1.08736291126619, -2.8690580709896, 7.28578834333112, 2.46270773070161,
1.86983400375641, -6.96898114084344, 19.4502701426409, 24.4669690634505,
7.83554456564033, 28.658774829192, 8.87077863575217, 17.918637653737,
23.3422657861999, -4.51877214309878, -7.51356000457387, 31.4454001495886,
22.5381677624583, 14.9164267860739, -14.6229709355072, 45.9891584442634,
32.6658900463158, 25.2656791256512, 6.81801801801802, 14.6173254835997,
7.30838420896072, -6.61179364500128, -9.24692680515084, 6.47169853886704,
-1.97933646268675, 26.8794942928519, 6.50959279199276, 18.1148335118046,
20.2061975019616, 19.3187670852977, -3.67290984190286, -4.62677608356775,
13.2044730800563, 11.812905299378, 10.2593738831563), dT = c(1.74,
1.116, 1.286, 1.314, 1.952, 1.418, 1.376, 1.826, 1.657, 1.99,
1.769, 1.703, 1.547, 1.183, 3.343, 1.46, 1.105, 1.274, 1.295,
1.797, 1.149, 1.427, 0.785, 1.956, 2.267, 1.011, 1.972, 1.808,
1.713, 1.769, 1.864, 1.202, 1.486, 1.273, 3.366, 2.504, 2.273,
2.284, 3.671, 2.461, 2.344, 3.256, 3.594, 4.136, 3.901, 3.174,
2.833, 1.701, 5.998, 2.184, 2.082, 2.407, 2.315, 3.354, 2.656,
3.135, 1.881, 3.931, 4.02, 2.117, 4.002, 4.151, 3.641, 3.283,
3.018, 2.318, 2.447, 2.742, 1.713, 1.268, 1.574, 1.499, 1.879,
1.66, 1.658, 1.928, 1.781, 2.569, 2.036, 1.49, 1.75, 1.272, 3.008,
1.635, 1.481, 1.511, 1.447, 1.143, 1.343, 1.427, 1.069, 2.352,
2.337, 1.366, 2.156, 2.417, 1.759, 2.158, 2.205, 1.319, 1.446,
1.522, 4.822, 3.837, 3.997, 3.815, 5.742, 4.306, 4.197, 5.109,
6.036, 6.457, 5.376, 4.614, 4.244, 3.638, 8.075, 3.981, 3.69,
4.115, 3.78, 4.29, 4.737, 4.604, 3.605, 6.276, 6.486, 3.897,
6.876, 7.334, 4.911, 5.451, 5.488, 4.125, 4.141, 4.202)), row.names = c(NA,
-136L), class = "data.frame")

library(directlabels)
model <- loess(W ~ tt + hh, data = sensWth)
new.data <- expand.grid(tt = seq(from = min(sensWth$tt), to = max(sensWth$tt), length.out = 500),
hh = seq(from = min(sensWth$hh), to = max(sensWth$hh), length.out = 500))

gg <- predict(model, newdata = new.data)
new.data = data.frame(W = as.vector(gg),
new.data)

 

plt<-ggplot(new.data, aes(x = hh/100, y = tt, z = floor(W))) +
geom_tile(aes(fill = W)) +
stat_contour(bins = 10, aes(color = ..level..), show.legend = FALSE) +
#geom_contour(aes(z = W, colour = ..level..), breaks = 9000000, color = "red") +
ylab("Change in temperature in degree Celsius") +
xlab("percentage change in precipitation") +
scale_fill_gradientn(name = "W (in m3/year)",
values = scales::rescale(quantile(new.data$W)),
limits = c(min(new.data$W),max(new.data$W)),
# breaks = unname(quantile(new.data$W)),
breaks = seq(round(min(new.data$W)/1000000)*1000000,
round(max(new.data$W)/1000000)*1000000,
(round(max(new.data$W)/1000000)*1000000-round(min(new.data$W)/1000000)*1000000)/3),
colors = rainbow(7), guide = "colorbar") +
scale_x_continuous(breaks = seq(-0.3,0.3, 0.1), label = scales::percent)+
# scale_y_continuous(breaks = seq(7.24, 13.25, 0.53), label = scales::percent )+
scale_y_continuous(breaks = seq(-1, 6, 1))+
ggtitle("Variation of average annual sediment production with \n temperature and precipitation")+
guides(fill = guide_colorbar(barwidth = 0.5, barheight = 10))

plt