Advertisement
Guest User

Team Two

a guest
May 30th, 2017
95
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 5.43 KB | None | 0 0
  1. > library(dplyr)
  2. >
  3. > bety_src <- src_postgres(dbname = "bety",
  4. + password = "bety",
  5. + host = "terra-bety.default",
  6. + user = "bety",
  7. + port = 5432)
  8. >
  9. > traits <- tbl(bety_src, "traits")
  10. >
  11. > variables <- tbl(bety_src, "variables") %>%
  12. + mutate(variable_id = id)
  13. >
  14. > traits_variables <- traits %>%
  15. + left_join(variables, by = "variable_id") %>%
  16. + select(name, mean)
  17. >
  18. > head(traits_variables)
  19. > library(dplyr)
  20. >
  21. > bety_src <- src_postgres(dbname = "bety",
  22. + password = "bety",
  23. + host = "terra-bety.default",
  24. + user = "bety",
  25. + port = 5432)
  26. >
  27. > traits <- tbl(bety_src, "traits")
  28. >
  29. > variables <- tbl(bety_src, "variables")
  30. >
  31. > traits_variables <- traits %>%
  32. + left_join(variables, by = "id") %>%
  33. + select(name, mean)
  34. >
  35. > head(traits_variables)
  36. > library(dplyr)
  37. >
  38. > bety_src <- src_postgres(dbname = "bety",
  39. + password = "bety",
  40. + host = "terra-bety.default",
  41. + user = "bety",
  42. + port = 5432)
  43. >
  44. > traits <- tbl(bety_src, "traits")
  45. >
  46. > variables <- tbl(bety_src, "variables") %>%
  47. + mutate(variable_id = id)
  48. >
  49. > traits_variables <- traits %>%
  50. + left_join(variables, by = "variable_id") %>%
  51. + select(name, mean)
  52. >
  53. > head(traits_variables)
  54. > library(traits)
  55. > library(dplyr)
  56. >
  57. > bety_src <- src_postgres(dbname = "bety",
  58. + password = 'bety',
  59. + host = 'terra-bety.default',
  60. + user = 'bety',
  61. + port = 5432)
  62. >
  63. > species <- betydb_query(table = 'species', limit = 'none') %>%
  64. + mutate(specie_id = id)
  65.  
  66. Error in curl::curl_fetch_memory(url, handle = handle) :
  67. Operation was aborted by an application callback
  68. > library(traits)
  69. > library(dplyr)
  70. >
  71. > bety <- src_postgres(dbname = "bety",
  72. + password = 'bety',
  73. + host = 'terra-bety.default',
  74. + user = 'bety',
  75. + port = 5432)
  76. >
  77. > species <- betydb_query(table = 'species', limit = 'none') %>%
  78. + mutate(specie_id = id)
  79. Auto-disconnecting postgres connection (231, 4)
  80. Auto-disconnecting postgres connection (231, 3)
  81. Auto-disconnecting postgres connection (231, 2)
  82. >
  83. > sites <- betydb_query(table = 'sites', limit = 'none') %>%
  84. + mutate(site_id = id)
  85. Error in mutate_impl(.data, dots) :
  86. Unsupported type CLOSXP for column "site_id"
  87. >
  88. > yields <- betydb_query(table = 'yields', limit = 'none') %>%
  89. + select(id, date, mean, n, statname, stat, site_id, specie_id, treatment_id, citation_id, cultivar_id) %>%
  90. + left_join(species, by = 'specie_id') %>%
  91. + left_join(sites, by = 'site_id')
  92. Error in is.data.frame(y) : object 'sites' not found
  93. >
  94. > managements_treatments <- tbl(bety, 'managements_treatments') %>%
  95. + select(treatment_id, management_id)
  96. >
  97. > treatments <- tbl(bety, 'treatments') %>%
  98. + dplyr::mutate(treatment_id = id) %>%
  99. + dplyr::select(treatment_id, name, definition, control)
  100. >
  101. > managements <- tbl(bety, 'managements') %>%
  102. + filter(mgmttype %in% c('fertilizer_N', 'fertilizer_N_rate', 'planting', 'irrigation')) %>%
  103. + dplyr::mutate(management_id = id) %>%
  104. + dplyr::select(management_id, date, mgmttype, level, units) %>%
  105. + left_join(managements_treatments, by = 'management_id') %>%
  106. + left_join(treatments, by = 'treatment_id')
  107. >
  108. >
  109. > nitrogen <- managements %>%
  110. + filter(mgmttype == "fertilizer_N_rate") %>%
  111. + select(treatment_id, nrate = level)
  112. >
  113. > planting <- managements %>% filter(mgmttype == "planting") %>%
  114. + select(treatment_id, planting_date = date)
  115. >
  116. > planting_rate <- managements %>% filter(mgmttype == "planting") %>%
  117. + select(treatment_id, planting_date = date, planting_density = level)
  118. >
  119. > irrigation <- managements %>%
  120. + filter(mgmttype == 'irrigation')
  121. >
  122. > irrigation_rate <- irrigation %>%
  123. + filter(units == 'mm', !is.na(treatment_id)) %>%
  124. + group_by(treatment_id, year = sql("extract(year from date)"), units) %>%
  125. + summarise(irrig.mm = sum(level)) %>%
  126. + group_by(treatment_id) %>%
  127. + summarise(irrig.mm.y = mean(irrig.mm))
  128. >
  129. > irrigation_boolean <- irrigation %>%
  130. + collect %>%
  131. + group_by(treatment_id) %>%
  132. + mutate(irrig = as.logical(mean(level))) %>%
  133. + select(treatment_id, irrig = irrig)
  134. Error in resolve_vars(new_groups, tbl_vars(.data)) :
  135. unknown variable to group by : treatment_id
  136. >
  137. > irrigation_all <- irrigation_boolean %>%
  138. + full_join(irrigation_rate, copy = TRUE, by = 'treatment_id')
  139. Error in eval(lhs, parent, parent) :
  140. object 'irrigation_boolean' not found
  141. >
  142. > grass_yields <- yields %>%
  143. + filter(genus %in% c('Miscanthus', 'Panicum')) %>%
  144. + left_join(nitrogen, by = 'treatment_id') %>%
  145. + #left_join(planting, by = 'treatment_id') %>%
  146. + left_join(planting_rate, by = 'treatment_id') %>%
  147. + left_join(irrigation_all, by = 'treatment_id', copy = TRUE) %>%
  148. + collect %>%
  149. + mutate(age = year(date)- year(planting_date),
  150. + nrate = ifelse(is.na(nrate), 0, nrate),
  151. + SE = ifelse(statname == "SE", stat, ifelse(statname == 'SD', stat / sqrt(n), NA)),
  152. + continent = ifelse(lon < -30, 'united_states', ifelse(lon < 75, 'europe', 'asia'))) %>%
  153. + filter(!duplicated(.))
  154. Error in eval(lhs, parent, parent) : object 'yields' not found
  155. >
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement