Advertisement
Guest User

Untitled

a guest
Jun 23rd, 2017
58
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.16 KB | None | 0 0
  1. Cumulative elevation trends
  2. ---------------------------
  3.  
  4. ### SH3 deep - including all data points (pre & post Wilma)
  5.  
  6. ``` r
  7. DCSET_all <- lm(Site_avg~Days, data=subset(SH3, Type=="Deep_CU"))
  8. ```
  9.  
  10. ##
  11. ## Call:
  12. ## lm(formula = Site_avg ~ Days, data = subset(SH3, Type == "Deep_CU"))
  13. ##
  14. ## Residuals:
  15. ## Min 1Q Median 3Q Max
  16. ## -20.3496 -3.4127 -0.7529 2.3280 26.7322
  17. ##
  18. ## Coefficients:
  19. ## Estimate Std. Error t value Pr(>|t|)
  20. ## (Intercept) 9.6417345 2.0161167 4.782 2.15e-05 ***
  21. ## Days 0.0082418 0.0008241 10.001 1.12e-12 ***
  22. ## ---
  23. ## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
  24. ##
  25. ## Residual standard error: 8.376 on 42 degrees of freedom
  26. ## Multiple R-squared: 0.7043, Adjusted R-squared: 0.6972
  27. ## F-statistic: 100 on 1 and 42 DF, p-value: 1.118e-12
  28.  
  29. ##### rate of elevation change for deep = 3.01 mm/yr
  30.  
  31. ``` r
  32. coef(DCSET_all)["Days"] * 365.25
  33. ```
  34.  
  35. ## Days
  36. ## 3.010314
  37.  
  38. ------------------------------------------------------------------------
  39.  
  40. ### SH3 shallow - including all data points (pre & post Wilma)
  41.  
  42. ``` r
  43. SCSET_all <- lm(Site_avg~Days, data=subset(SH3, Type="Shallow_CU"))
  44. ```
  45.  
  46. ##
  47. ## Call:
  48. ## lm(formula = Site_avg ~ Days, data = subset(SH3, Type = "Shallow_CU"))
  49. ##
  50. ## Residuals:
  51. ## Min 1Q Median 3Q Max
  52. ## -42.221 -13.955 -1.929 9.733 65.137
  53. ##
  54. ## Coefficients:
  55. ## Estimate Std. Error t value Pr(>|t|)
  56. ## (Intercept) 3.1016718 1.5934502 1.947 0.0524 .
  57. ## Days 0.0049936 0.0006328 7.892 3.79e-14 ***
  58. ## ---
  59. ## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
  60. ##
  61. ## Residual standard error: 18.82 on 353 degrees of freedom
  62. ## Multiple R-squared: 0.15, Adjusted R-squared: 0.1476
  63. ## F-statistic: 62.28 on 1 and 353 DF, p-value: 3.79e-14
  64.  
  65. ##### rate of elevation change for shallow = 1.8 mm/yr
  66.  
  67. ``` r
  68. coef(SCSET_all)["Days"] * 365.25
  69. ```
  70.  
  71. ## Days
  72. ## 1.823916
  73.  
  74. #### Cumulative Elevation - all dates
  75.  
  76. ![](SH3_files/figure-markdown_github/unnamed-chunk-11-1.png)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement