Skip to content

Instantly share code, notes, and snippets.

@mgrabovsky
Last active April 9, 2024 10:25
Show Gist options
  • Select an option

  • Save mgrabovsky/2104a47fbb7264071ac8992906074088 to your computer and use it in GitHub Desktop.

Select an option

Save mgrabovsky/2104a47fbb7264071ac8992906074088 to your computer and use it in GitHub Desktop.
Simulating observations in the Lighthouse problem in R
Display the source blob
Display the rendered blob
Raw
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Display the source blob
Display the rendered blob
Raw
<?xml version='1.0' encoding='UTF-8' ?>
<svg xmlns='http://www.w3.org/2000/svg' xmlns:xlink='http://www.w3.org/1999/xlink' viewBox='0 0 504.00 360.00'>
<defs>
<style type='text/css'><![CDATA[
line, polyline, path, rect, circle {
fill: none;
stroke: #000000;
stroke-linecap: round;
stroke-linejoin: round;
stroke-miterlimit: 10.00;
}
]]></style>
</defs>
<rect width='100%' height='100%' style='stroke: none; fill: #FFFFFF;'/>
<defs>
<clipPath id='cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ=='>
<rect x='40.68' y='5.48' width='457.84' height='322.15' />
</clipPath>
</defs>
<polyline points='40.68,270.62 498.52,270.62 ' style='stroke-width: 0.53; stroke: #EBEBEB; stroke-linecap: butt;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<polyline points='40.68,185.90 498.52,185.90 ' style='stroke-width: 0.53; stroke: #EBEBEB; stroke-linecap: butt;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<polyline points='40.68,101.17 498.52,101.17 ' style='stroke-width: 0.53; stroke: #EBEBEB; stroke-linecap: butt;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<polyline points='40.68,16.45 498.52,16.45 ' style='stroke-width: 0.53; stroke: #EBEBEB; stroke-linecap: butt;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<polyline points='73.05,327.63 73.05,5.48 ' style='stroke-width: 0.53; stroke: #EBEBEB; stroke-linecap: butt;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<polyline points='188.67,327.63 188.67,5.48 ' style='stroke-width: 0.53; stroke: #EBEBEB; stroke-linecap: butt;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<polyline points='304.28,327.63 304.28,5.48 ' style='stroke-width: 0.53; stroke: #EBEBEB; stroke-linecap: butt;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<polyline points='419.90,327.63 419.90,5.48 ' style='stroke-width: 0.53; stroke: #EBEBEB; stroke-linecap: butt;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<polyline points='40.68,312.98 498.52,312.98 ' style='stroke-width: 1.07; stroke: #EBEBEB; stroke-linecap: butt;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<polyline points='40.68,228.26 498.52,228.26 ' style='stroke-width: 1.07; stroke: #EBEBEB; stroke-linecap: butt;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<polyline points='40.68,143.54 498.52,143.54 ' style='stroke-width: 1.07; stroke: #EBEBEB; stroke-linecap: butt;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<polyline points='40.68,58.81 498.52,58.81 ' style='stroke-width: 1.07; stroke: #EBEBEB; stroke-linecap: butt;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<polyline points='130.86,327.63 130.86,5.48 ' style='stroke-width: 1.07; stroke: #EBEBEB; stroke-linecap: butt;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<polyline points='246.48,327.63 246.48,5.48 ' style='stroke-width: 1.07; stroke: #EBEBEB; stroke-linecap: butt;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<polyline points='362.09,327.63 362.09,5.48 ' style='stroke-width: 1.07; stroke: #EBEBEB; stroke-linecap: butt;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<polyline points='477.71,327.63 477.71,5.48 ' style='stroke-width: 1.07; stroke: #EBEBEB; stroke-linecap: butt;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<polyline points='61.49,310.17 62.30,309.85 63.12,309.58 63.93,309.35 64.75,309.17 65.56,309.03 66.38,308.93 67.19,308.85 68.00,308.79 68.82,308.73 69.63,308.68 70.45,308.63 71.26,308.58 72.08,308.54 72.89,308.51 73.71,308.50 74.52,308.51 75.34,308.54 76.15,308.59 76.96,308.65 77.78,308.73 78.59,308.81 79.41,308.89 80.22,308.96 81.04,309.03 81.85,309.08 82.67,309.12 83.48,309.15 84.30,309.15 85.11,309.14 85.92,309.11 86.74,309.07 87.55,309.01 88.37,308.94 89.18,308.87 90.00,308.79 90.81,308.71 91.63,308.63 92.44,308.56 93.25,308.49 94.07,308.44 94.88,308.39 95.70,308.36 96.51,308.33 97.33,308.31 98.14,308.29 98.96,308.26 99.77,308.22 100.59,308.15 101.40,308.06 102.21,307.94 103.03,307.79 103.84,307.60 104.66,307.39 105.47,307.16 106.29,306.93 107.10,306.70 107.92,306.48 108.73,306.29 109.55,306.13 110.36,305.99 111.17,305.89 111.99,305.82 112.80,305.78 113.62,305.77 114.43,305.77 115.25,305.78 116.06,305.81 116.88,305.83 117.69,305.85 118.51,305.86 119.32,305.86 120.13,305.85 120.95,305.83 121.76,305.79 122.58,305.73 123.39,305.67 124.21,305.59 125.02,305.50 125.84,305.40 126.65,305.28 127.46,305.15 128.28,305.01 129.09,304.86 129.91,304.69 130.72,304.51 131.54,304.31 132.35,304.11 133.17,303.91 133.98,303.70 134.80,303.48 135.61,303.27 136.42,303.06 137.24,302.86 138.05,302.67 138.87,302.50 139.68,302.36 140.50,302.25 141.31,302.19 142.13,302.16 142.94,302.19 143.76,302.26 144.57,302.38 145.38,302.52 146.20,302.66 147.01,302.80 147.83,302.90 148.64,302.94 149.46,302.92 150.27,302.82 151.09,302.63 151.90,302.35 152.72,302.01 153.53,301.61 154.34,301.17 155.16,300.71 155.97,300.25 156.79,299.82 157.60,299.40 158.42,299.02 159.23,298.65 160.05,298.31 160.86,297.97 161.67,297.64 162.49,297.31 163.30,297.00 164.12,296.70 164.93,296.42 165.75,296.19 166.56,296.00 167.38,295.87 168.19,295.79 169.01,295.75 169.82,295.74 170.63,295.76 171.45,295.80 172.26,295.85 173.08,295.91 173.89,295.98 174.71,296.05 175.52,296.13 176.34,296.20 177.15,296.25 177.97,296.27 178.78,296.23 179.59,296.13 180.41,295.96 181.22,295.70 182.04,295.34 182.85,294.88 183.67,294.32 184.48,293.67 185.30,292.93 186.11,292.11 186.92,291.22 187.74,290.28 188.55,289.32 189.37,288.35 190.18,287.41 191.00,286.51 191.81,285.66 192.63,284.87 193.44,284.13 194.26,283.45 195.07,282.82 195.88,282.24 196.70,281.69 197.51,281.17 198.33,280.68 199.14,280.21 199.96,279.76 200.77,279.31 201.59,278.86 202.40,278.39 203.22,277.90 204.03,277.37 204.84,276.79 205.66,276.16 206.47,275.47 207.29,274.72 208.10,273.90 208.92,273.02 209.73,272.07 210.55,271.05 211.36,269.97 212.18,268.81 212.99,267.58 213.80,266.29 214.62,264.94 215.43,263.53 216.25,262.07 217.06,260.55 217.88,258.97 218.69,257.31 219.51,255.57 220.32,253.74 221.13,251.80 221.95,249.74 222.76,247.56 223.58,245.23 224.39,242.75 225.21,240.12 226.02,237.33 226.84,234.39 227.65,231.31 228.47,228.09 229.28,224.74 230.09,221.31 230.91,217.82 231.72,214.30 232.54,210.78 233.35,207.29 234.17,203.83 234.98,200.42 235.80,197.03 236.61,193.64 237.43,190.20 238.24,186.64 239.05,182.88 239.87,178.86 240.68,174.52 241.50,169.81 242.31,164.73 243.13,159.25 243.94,153.42 244.76,147.28 245.57,140.89 246.39,134.33 247.20,127.69 248.01,121.04 248.83,114.46 249.64,108.04 250.46,101.79 251.27,95.73 252.09,89.88 252.90,84.22 253.72,78.74 254.53,73.43 255.34,68.30 256.16,63.37 256.97,58.68 257.79,54.26 258.60,50.19 259.42,46.54 260.23,43.43 261.05,40.87 261.86,38.90 262.68,37.56 263.49,36.83 264.30,36.71 265.12,37.15 265.93,38.10 266.75,39.52 267.56,41.35 268.38,43.54 269.19,46.04 270.01,48.82 270.82,51.86 271.64,55.11 272.45,58.57 273.26,62.22 274.08,66.07 274.89,70.11 275.71,74.36 276.52,78.84 277.34,83.57 278.15,88.56 278.97,93.84 279.78,99.41 280.60,105.29 281.41,111.47 282.22,117.89 283.04,124.52 283.85,131.30 284.67,138.16 285.48,145.02 286.30,151.80 287.11,158.42 287.93,164.80 288.74,170.86 289.55,176.55 290.37,181.82 291.18,186.60 292.00,190.94 292.81,194.86 293.63,198.41 294.44,201.64 295.26,204.62 296.07,207.42 296.89,210.09 297.70,212.71 298.51,215.29 299.33,217.88 300.14,220.49 300.96,223.12 301.77,225.75 302.59,228.37 303.40,230.95 304.22,233.47 305.03,235.89 305.85,238.20 306.66,240.37 307.47,242.40 308.29,244.27 309.10,245.98 309.92,247.54 310.73,248.95 311.55,250.22 312.36,251.36 313.18,252.41 313.99,253.39 314.80,254.33 315.62,255.29 316.43,256.29 317.25,257.38 318.06,258.60 318.88,259.95 319.69,261.45 320.51,263.09 321.32,264.84 322.14,266.66 322.95,268.50 323.76,270.31 324.58,272.07 325.39,273.72 326.21,275.27 327.02,276.71 327.84,278.02 328.65,279.23 329.47,280.35 330.28,281.38 331.10,282.34 331.91,283.25 332.72,284.11 333.54,284.92 334.35,285.69 335.17,286.42 335.98,287.10 336.80,287.72 337.61,288.26 338.43,288.73 339.24,289.12 340.06,289.44 340.87,289.70 341.68,289.91 342.50,290.12 343.31,290.34 344.13,290.61 344.94,290.93 345.76,291.33 346.57,291.79 347.39,292.30 348.20,292.85 349.01,293.41 349.83,293.95 350.64,294.45 351.46,294.91 352.27,295.30 353.09,295.63 353.90,295.89 354.72,296.11 355.53,296.29 356.35,296.47 357.16,296.65 357.97,296.86 358.79,297.10 359.60,297.39 360.42,297.73 361.23,298.10 362.05,298.50 362.86,298.90 363.68,299.29 364.49,299.63 365.31,299.93 366.12,300.16 366.93,300.33 367.75,300.44 368.56,300.49 369.38,300.49 370.19,300.47 371.01,300.42 371.82,300.36 372.64,300.30 373.45,300.25 374.27,300.22 375.08,300.22 375.89,300.26 376.71,300.33 377.52,300.45 378.34,300.61 379.15,300.83 379.97,301.09 380.78,301.39 381.60,301.73 382.41,302.09 383.22,302.46 384.04,302.83 384.85,303.18 385.67,303.49 386.48,303.76 387.30,303.96 388.11,304.10 388.93,304.17 389.74,304.17 390.56,304.10 391.37,303.98 392.18,303.82 393.00,303.63 393.81,303.43 394.63,303.23 395.44,303.04 396.26,302.89 397.07,302.78 397.89,302.72 398.70,302.73 399.52,302.80 400.33,302.94 401.14,303.15 401.96,303.42 402.77,303.74 403.59,304.11 404.40,304.49 405.22,304.87 406.03,305.23 406.85,305.55 407.66,305.82 408.48,306.05 409.29,306.22 410.10,306.35 410.92,306.46 411.73,306.55 412.55,306.65 413.36,306.76 414.18,306.90 414.99,307.06 415.81,307.25 416.62,307.45 417.43,307.64 418.25,307.82 419.06,307.96 419.88,308.06 420.69,308.11 421.51,308.12 422.32,308.07 423.14,307.99 423.95,307.88 424.77,307.76 425.58,307.65 426.39,307.55 427.21,307.47 428.02,307.43 428.84,307.41 429.65,307.42 430.47,307.45 431.28,307.50 432.10,307.54 432.91,307.57 433.73,307.59 434.54,307.58 435.35,307.55 436.17,307.51 436.98,307.45 437.80,307.39 438.61,307.34 439.43,307.31 440.24,307.29 441.06,307.30 441.87,307.34 442.69,307.40 443.50,307.47 444.31,307.56 445.13,307.65 445.94,307.74 446.76,307.83 447.57,307.92 448.39,307.99 449.20,308.07 450.02,308.13 450.83,308.20 451.64,308.26 452.46,308.31 453.27,308.37 454.09,308.42 454.90,308.47 455.72,308.51 456.53,308.54 457.35,308.56 458.16,308.56 458.98,308.54 459.79,308.52 460.60,308.48 461.42,308.45 462.23,308.43 463.05,308.42 463.86,308.43 464.68,308.47 465.49,308.54 466.31,308.64 467.12,308.77 467.94,308.93 468.75,309.10 469.56,309.28 470.38,309.47 471.19,309.68 472.01,309.88 472.82,310.10 473.64,310.32 474.45,310.54 475.27,310.78 476.08,311.02 476.89,311.26 477.71,311.50 477.71,312.98 476.89,312.98 476.08,312.98 475.27,312.98 474.45,312.98 473.64,312.98 472.82,312.98 472.01,312.98 471.19,312.98 470.38,312.98 469.56,312.98 468.75,312.98 467.94,312.98 467.12,312.98 466.31,312.98 465.49,312.98 464.68,312.98 463.86,312.98 463.05,312.98 462.23,312.98 461.42,312.98 460.60,312.98 459.79,312.98 458.98,312.98 458.16,312.98 457.35,312.98 456.53,312.98 455.72,312.98 454.90,312.98 454.09,312.98 453.27,312.98 452.46,312.98 451.64,312.98 450.83,312.98 450.02,312.98 449.20,312.98 448.39,312.98 447.57,312.98 446.76,312.98 445.94,312.98 445.13,312.98 444.31,312.98 443.50,312.98 442.69,312.98 441.87,312.98 441.06,312.98 440.24,312.98 439.43,312.98 438.61,312.98 437.80,312.98 436.98,312.98 436.17,312.98 435.35,312.98 434.54,312.98 433.73,312.98 432.91,312.98 432.10,312.98 431.28,312.98 430.47,312.98 429.65,312.98 428.84,312.98 428.02,312.98 427.21,312.98 426.39,312.98 425.58,312.98 424.77,312.98 423.95,312.98 423.14,312.98 422.32,312.98 421.51,312.98 420.69,312.98 419.88,312.98 419.06,312.98 418.25,312.98 417.43,312.98 416.62,312.98 415.81,312.98 414.99,312.98 414.18,312.98 413.36,312.98 412.55,312.98 411.73,312.98 410.92,312.98 410.10,312.98 409.29,312.98 408.48,312.98 407.66,312.98 406.85,312.98 406.03,312.98 405.22,312.98 404.40,312.98 403.59,312.98 402.77,312.98 401.96,312.98 401.14,312.98 400.33,312.98 399.52,312.98 398.70,312.98 397.89,312.98 397.07,312.98 396.26,312.98 395.44,312.98 394.63,312.98 393.81,312.98 393.00,312.98 392.18,312.98 391.37,312.98 390.56,312.98 389.74,312.98 388.93,312.98 388.11,312.98 387.30,312.98 386.48,312.98 385.67,312.98 384.85,312.98 384.04,312.98 383.22,312.98 382.41,312.98 381.60,312.98 380.78,312.98 379.97,312.98 379.15,312.98 378.34,312.98 377.52,312.98 376.71,312.98 375.89,312.98 375.08,312.98 374.27,312.98 373.45,312.98 372.64,312.98 371.82,312.98 371.01,312.98 370.19,312.98 369.38,312.98 368.56,312.98 367.75,312.98 366.93,312.98 366.12,312.98 365.31,312.98 364.49,312.98 363.68,312.98 362.86,312.98 362.05,312.98 361.23,312.98 360.42,312.98 359.60,312.98 358.79,312.98 357.97,312.98 357.16,312.98 356.35,312.98 355.53,312.98 354.72,312.98 353.90,312.98 353.09,312.98 352.27,312.98 351.46,312.98 350.64,312.98 349.83,312.98 349.01,312.98 348.20,312.98 347.39,312.98 346.57,312.98 345.76,312.98 344.94,312.98 344.13,312.98 343.31,312.98 342.50,312.98 341.68,312.98 340.87,312.98 340.06,312.98 339.24,312.98 338.43,312.98 337.61,312.98 336.80,312.98 335.98,312.98 335.17,312.98 334.35,312.98 333.54,312.98 332.72,312.98 331.91,312.98 331.10,312.98 330.28,312.98 329.47,312.98 328.65,312.98 327.84,312.98 327.02,312.98 326.21,312.98 325.39,312.98 324.58,312.98 323.76,312.98 322.95,312.98 322.14,312.98 321.32,312.98 320.51,312.98 319.69,312.98 318.88,312.98 318.06,312.98 317.25,312.98 316.43,312.98 315.62,312.98 314.80,312.98 313.99,312.98 313.18,312.98 312.36,312.98 311.55,312.98 310.73,312.98 309.92,312.98 309.10,312.98 308.29,312.98 307.47,312.98 306.66,312.98 305.85,312.98 305.03,312.98 304.22,312.98 303.40,312.98 302.59,312.98 301.77,312.98 300.96,312.98 300.14,312.98 299.33,312.98 298.51,312.98 297.70,312.98 296.89,312.98 296.07,312.98 295.26,312.98 294.44,312.98 293.63,312.98 292.81,312.98 292.00,312.98 291.18,312.98 290.37,312.98 289.55,312.98 288.74,312.98 287.93,312.98 287.11,312.98 286.30,312.98 285.48,312.98 284.67,312.98 283.85,312.98 283.04,312.98 282.22,312.98 281.41,312.98 280.60,312.98 279.78,312.98 278.97,312.98 278.15,312.98 277.34,312.98 276.52,312.98 275.71,312.98 274.89,312.98 274.08,312.98 273.26,312.98 272.45,312.98 271.64,312.98 270.82,312.98 270.01,312.98 269.19,312.98 268.38,312.98 267.56,312.98 266.75,312.98 265.93,312.98 265.12,312.98 264.30,312.98 263.49,312.98 262.68,312.98 261.86,312.98 261.05,312.98 260.23,312.98 259.42,312.98 258.60,312.98 257.79,312.98 256.97,312.98 256.16,312.98 255.34,312.98 254.53,312.98 253.72,312.98 252.90,312.98 252.09,312.98 251.27,312.98 250.46,312.98 249.64,312.98 248.83,312.98 248.01,312.98 247.20,312.98 246.39,312.98 245.57,312.98 244.76,312.98 243.94,312.98 243.13,312.98 242.31,312.98 241.50,312.98 240.68,312.98 239.87,312.98 239.05,312.98 238.24,312.98 237.43,312.98 236.61,312.98 235.80,312.98 234.98,312.98 234.17,312.98 233.35,312.98 232.54,312.98 231.72,312.98 230.91,312.98 230.09,312.98 229.28,312.98 228.47,312.98 227.65,312.98 226.84,312.98 226.02,312.98 225.21,312.98 224.39,312.98 223.58,312.98 222.76,312.98 221.95,312.98 221.13,312.98 220.32,312.98 219.51,312.98 218.69,312.98 217.88,312.98 217.06,312.98 216.25,312.98 215.43,312.98 214.62,312.98 213.80,312.98 212.99,312.98 212.18,312.98 211.36,312.98 210.55,312.98 209.73,312.98 208.92,312.98 208.10,312.98 207.29,312.98 206.47,312.98 205.66,312.98 204.84,312.98 204.03,312.98 203.22,312.98 202.40,312.98 201.59,312.98 200.77,312.98 199.96,312.98 199.14,312.98 198.33,312.98 197.51,312.98 196.70,312.98 195.88,312.98 195.07,312.98 194.26,312.98 193.44,312.98 192.63,312.98 191.81,312.98 191.00,312.98 190.18,312.98 189.37,312.98 188.55,312.98 187.74,312.98 186.92,312.98 186.11,312.98 185.30,312.98 184.48,312.98 183.67,312.98 182.85,312.98 182.04,312.98 181.22,312.98 180.41,312.98 179.59,312.98 178.78,312.98 177.97,312.98 177.15,312.98 176.34,312.98 175.52,312.98 174.71,312.98 173.89,312.98 173.08,312.98 172.26,312.98 171.45,312.98 170.63,312.98 169.82,312.98 169.01,312.98 168.19,312.98 167.38,312.98 166.56,312.98 165.75,312.98 164.93,312.98 164.12,312.98 163.30,312.98 162.49,312.98 161.67,312.98 160.86,312.98 160.05,312.98 159.23,312.98 158.42,312.98 157.60,312.98 156.79,312.98 155.97,312.98 155.16,312.98 154.34,312.98 153.53,312.98 152.72,312.98 151.90,312.98 151.09,312.98 150.27,312.98 149.46,312.98 148.64,312.98 147.83,312.98 147.01,312.98 146.20,312.98 145.38,312.98 144.57,312.98 143.76,312.98 142.94,312.98 142.13,312.98 141.31,312.98 140.50,312.98 139.68,312.98 138.87,312.98 138.05,312.98 137.24,312.98 136.42,312.98 135.61,312.98 134.80,312.98 133.98,312.98 133.17,312.98 132.35,312.98 131.54,312.98 130.72,312.98 129.91,312.98 129.09,312.98 128.28,312.98 127.46,312.98 126.65,312.98 125.84,312.98 125.02,312.98 124.21,312.98 123.39,312.98 122.58,312.98 121.76,312.98 120.95,312.98 120.13,312.98 119.32,312.98 118.51,312.98 117.69,312.98 116.88,312.98 116.06,312.98 115.25,312.98 114.43,312.98 113.62,312.98 112.80,312.98 111.99,312.98 111.17,312.98 110.36,312.98 109.55,312.98 108.73,312.98 107.92,312.98 107.10,312.98 106.29,312.98 105.47,312.98 104.66,312.98 103.84,312.98 103.03,312.98 102.21,312.98 101.40,312.98 100.59,312.98 99.77,312.98 98.96,312.98 98.14,312.98 97.33,312.98 96.51,312.98 95.70,312.98 94.88,312.98 94.07,312.98 93.25,312.98 92.44,312.98 91.63,312.98 90.81,312.98 90.00,312.98 89.18,312.98 88.37,312.98 87.55,312.98 86.74,312.98 85.92,312.98 85.11,312.98 84.30,312.98 83.48,312.98 82.67,312.98 81.85,312.98 81.04,312.98 80.22,312.98 79.41,312.98 78.59,312.98 77.78,312.98 76.96,312.98 76.15,312.98 75.34,312.98 74.52,312.98 73.71,312.98 72.89,312.98 72.08,312.98 71.26,312.98 70.45,312.98 69.63,312.98 68.82,312.98 68.00,312.98 67.19,312.98 66.38,312.98 65.56,312.98 64.75,312.98 63.93,312.98 63.12,312.98 62.30,312.98 61.49,312.98 61.49,310.17 ' style='stroke-width: 2.13; stroke: #CD6839;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='61.49' y='306.57' width='4.62' height='6.42' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='66.11' y='310.23' width='4.62' height='2.75' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='70.74' y='307.02' width='4.62' height='5.96' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='75.36' y='309.32' width='4.62' height='3.67' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='79.99' y='307.94' width='4.62' height='5.04' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='84.61' y='309.77' width='4.62' height='3.21' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='89.24' y='309.77' width='4.62' height='3.21' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='93.86' y='307.02' width='4.62' height='5.96' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='98.49' y='309.32' width='4.62' height='3.67' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='103.11' y='307.48' width='4.62' height='5.50' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='107.74' y='305.65' width='4.62' height='7.33' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='112.36' y='305.19' width='4.62' height='7.79' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='116.98' y='306.57' width='4.62' height='6.42' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='121.61' y='305.19' width='4.62' height='7.79' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='126.23' y='305.65' width='4.62' height='7.33' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='130.86' y='303.82' width='4.62' height='9.17' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='135.48' y='301.98' width='4.62' height='11.00' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='140.11' y='301.07' width='4.62' height='11.92' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='144.73' y='304.27' width='4.62' height='8.71' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='149.36' y='303.82' width='4.62' height='9.17' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='153.98' y='296.94' width='4.62' height='16.04' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='158.61' y='301.52' width='4.62' height='11.46' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='163.23' y='293.73' width='4.62' height='19.25' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='167.86' y='297.40' width='4.62' height='15.58' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='172.48' y='294.19' width='4.62' height='18.79' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='177.11' y='297.86' width='4.62' height='15.12' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='181.73' y='294.19' width='4.62' height='18.79' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='186.35' y='290.07' width='4.62' height='22.92' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='190.98' y='285.48' width='4.62' height='27.50' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='195.60' y='277.23' width='4.62' height='35.75' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='200.23' y='280.44' width='4.62' height='32.54' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='204.85' y='274.94' width='4.62' height='38.04' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='209.48' y='273.11' width='4.62' height='39.87' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='214.10' y='261.19' width='4.62' height='51.79' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='218.73' y='252.03' width='4.62' height='60.96' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='223.35' y='239.65' width='4.62' height='73.33' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='227.98' y='225.90' width='4.62' height='87.08' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='232.60' y='192.91' width='4.62' height='120.08' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='237.23' y='192.45' width='4.62' height='120.54' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='241.85' y='161.74' width='4.62' height='151.24' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='246.48' y='103.53' width='4.62' height='209.45' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='251.10' y='80.16' width='4.62' height='232.82' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='255.72' y='52.66' width='4.62' height='260.32' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='260.35' y='20.12' width='4.62' height='292.86' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='264.97' y='37.08' width='4.62' height='275.90' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='269.60' y='50.83' width='4.62' height='262.15' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='274.22' y='80.16' width='4.62' height='232.82' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='278.85' y='105.83' width='4.62' height='207.16' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='283.47' y='154.87' width='4.62' height='158.12' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='288.10' y='186.95' width='4.62' height='126.04' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='292.72' y='202.53' width='4.62' height='110.45' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='297.35' y='221.32' width='4.62' height='91.66' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='301.97' y='235.07' width='4.62' height='77.91' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='306.60' y='247.90' width='4.62' height='65.08' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='311.22' y='255.69' width='4.62' height='57.29' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='315.85' y='249.74' width='4.62' height='63.25' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='320.47' y='272.65' width='4.62' height='40.33' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='325.10' y='278.61' width='4.62' height='34.37' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='329.72' y='285.03' width='4.62' height='27.96' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='334.34' y='285.48' width='4.62' height='27.50' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='338.97' y='293.73' width='4.62' height='19.25' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='343.59' y='288.23' width='4.62' height='24.75' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='348.22' y='296.02' width='4.62' height='16.96' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='352.84' y='298.77' width='4.62' height='14.21' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='357.47' y='293.73' width='4.62' height='19.25' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='362.09' y='302.44' width='4.62' height='10.54' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='366.72' y='300.15' width='4.62' height='12.83' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='371.34' y='299.23' width='4.62' height='13.75' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='375.97' y='300.15' width='4.62' height='12.83' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='380.59' y='301.98' width='4.62' height='11.00' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='385.22' y='306.11' width='4.62' height='6.87' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='389.84' y='304.27' width='4.62' height='8.71' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='394.47' y='301.52' width='4.62' height='11.46' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='399.09' y='302.44' width='4.62' height='10.54' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='403.71' y='306.11' width='4.62' height='6.87' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='408.34' y='306.57' width='4.62' height='6.42' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='412.96' y='305.65' width='4.62' height='7.33' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='417.59' y='309.77' width='4.62' height='3.21' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='422.21' y='307.94' width='4.62' height='5.04' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='426.84' y='306.11' width='4.62' height='6.87' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='431.46' y='308.86' width='4.62' height='4.12' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='436.09' y='306.57' width='4.62' height='6.42' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='440.71' y='307.48' width='4.62' height='5.50' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='445.34' y='307.94' width='4.62' height='5.04' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='449.96' y='308.86' width='4.62' height='4.12' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='454.59' y='308.40' width='4.62' height='4.58' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='459.21' y='308.40' width='4.62' height='4.58' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='463.84' y='307.94' width='4.62' height='5.04' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='468.46' y='310.23' width='4.62' height='2.75' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<rect x='473.08' y='309.77' width='4.62' height='3.21' style='stroke-width: 1.07; stroke: #FFFFFF; stroke-linecap: butt; fill: #708090;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<line x1='264.97' y1='327.63' x2='264.97' y2='5.48' style='stroke-width: 1.71; stroke: #FFC125; stroke-dasharray: 9.10,9.10; stroke-linecap: butt;' clip-path='url(#cpNDAuNjc3NHw0OTguNTIxfDMyNy42MjV8NS40Nzk0NQ==)' />
<defs>
<clipPath id='cpMHw1MDR8MzYwfDA='>
<rect x='0.00' y='0.00' width='504.00' height='360.00' />
</clipPath>
</defs>
<g clip-path='url(#cpMHw1MDR8MzYwfDA=)'><text x='18.62' y='316.13' style='font-size: 8.80px; fill: #4D4D4D; font-family: Arial;' textLength='17.13px' lengthAdjust='spacingAndGlyphs'>0.00</text></g>
<g clip-path='url(#cpMHw1MDR8MzYwfDA=)'><text x='18.62' y='231.41' style='font-size: 8.80px; fill: #4D4D4D; font-family: Arial;' textLength='17.13px' lengthAdjust='spacingAndGlyphs'>0.02</text></g>
<g clip-path='url(#cpMHw1MDR8MzYwfDA=)'><text x='18.62' y='146.69' style='font-size: 8.80px; fill: #4D4D4D; font-family: Arial;' textLength='17.13px' lengthAdjust='spacingAndGlyphs'>0.04</text></g>
<g clip-path='url(#cpMHw1MDR8MzYwfDA=)'><text x='18.62' y='61.96' style='font-size: 8.80px; fill: #4D4D4D; font-family: Arial;' textLength='17.13px' lengthAdjust='spacingAndGlyphs'>0.06</text></g>
<g clip-path='url(#cpMHw1MDR8MzYwfDA=)'><text x='124.50' y='338.86' style='font-size: 8.80px; fill: #4D4D4D; font-family: Arial;' textLength='12.72px' lengthAdjust='spacingAndGlyphs'>-25</text></g>
<g clip-path='url(#cpMHw1MDR8MzYwfDA=)'><text x='244.03' y='338.86' style='font-size: 8.80px; fill: #4D4D4D; font-family: Arial;' textLength='4.89px' lengthAdjust='spacingAndGlyphs'>0</text></g>
<g clip-path='url(#cpMHw1MDR8MzYwfDA=)'><text x='357.20' y='338.86' style='font-size: 8.80px; fill: #4D4D4D; font-family: Arial;' textLength='9.79px' lengthAdjust='spacingAndGlyphs'>25</text></g>
<g clip-path='url(#cpMHw1MDR8MzYwfDA=)'><text x='472.82' y='338.86' style='font-size: 8.80px; fill: #4D4D4D; font-family: Arial;' textLength='9.79px' lengthAdjust='spacingAndGlyphs'>50</text></g>
<g clip-path='url(#cpMHw1MDR8MzYwfDA=)'><text x='266.85' y='352.00' style='font-size: 11.00px; font-family: Arial;' textLength='5.50px' lengthAdjust='spacingAndGlyphs'>x</text></g>
<g clip-path='url(#cpMHw1MDR8MzYwfDA=)'><text transform='translate(13.35,183.98) rotate(-90)' style='font-size: 11.00px; font-family: Arial;' textLength='34.85px' lengthAdjust='spacingAndGlyphs'>density</text></g>
</svg>
library(ggplot2)
# Metaparameters.
d <- 5 # Distance of the lighthouse from the coast.
x0 <- 4 # Perpendicular position of the lighthouse along the coast.
set.seed(1987)
# Simulate the data.
N <- 1e4
phi <- runif(N, -pi / 2, pi / 2)
x <- d * tan(phi) + x0
df_sims <- tibble(x = x)
# Plot the distribution.
ggplot(df_sims, aes(x, after_stat(density))) +
geom_density(colour = "sienna3", linewidth = 1) +
geom_histogram(binwidth = 1, center = .5, fill = "slategrey", colour = "white") +
geom_vline(xintercept = x0, colour = "goldenrod1", linetype = 2, linewidth = .8) +
xlim(-40, 50) +
theme_minimal()
# Distribution of absolute deviations from the median.
absolute_deviation <- abs(x - median(x))
mad_obs <- median(absolute_deviation)
df_absdev <- tibble(absolute_deviation = absolute_deviation)
ggplot(df_absdev, aes(absolute_deviation, after_stat(density))) +
geom_histogram(binwidth = 1, center = .5, fill = "sienna3", colour = "white") +
geom_vline(xintercept = mad_obs, colour = "slategrey", linetype = 2, linewidth = 1) +
annotate("text", x = mad_obs + 1.3, y = 0.114,
label = "Median absolute deviation",
colour = "slategrey", size = 3.5, hjust = 0) +
xlim(0, 60) +
theme_minimal()
ggplot(df_absdev, aes(absolute_deviation, after_stat(y))) +
stat_ecdf(colour = 'sienna3') +
geom_vline(xintercept = mad_obs, colour = "slategrey", linetype = 2, linewidth = 1) +
annotate("text", x = mad_obs + 1.3, y = 0.114,
label = "Median absolute deviation",
colour = "slategrey", size = 3.5, hjust = 0) +
scale_x_log10() +
ylab("Cumulative density") +
theme_minimal()
Display the source blob
Display the rendered blob
Raw
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment