Commit acbf64e4 authored by Grelot Frederic's avatar Grelot Frederic :swimmer_tone5:
Browse files

script/onrn

	- ajout de balance
	- réorganisation des variables

Refs #7
Showing with 10 additions and 6 deletions
+10 -6
...@@ -266,15 +266,19 @@ result = merge(result, temp, by = "commune", all.x = TRUE) ...@@ -266,15 +266,19 @@ result = merge(result, temp, by = "commune", all.x = TRUE)
### Bilan Sinistre - Prime : estimation ### Bilan Sinistre - Prime : estimation
# balance = ratio * population * premium_hab pop_france = 66992159 # INSEE (2018)
population = 66992159 # INSEE (2018) premium_france = 1670000000 # (CCR2019a pour 2018)
premium = 1670000000 # (CCR2019a pour 2018) premium_hab = premium_france / pop_france
premium_hab = premium / population
result[["balance"]] = (1 - result[["ratio"]]) * geau::so_ii_population[ , "2018"] * premium_hab result[["balance"]] = (1 - result[["ratio"]]) * geau::so_ii_population[ , "2018"] * premium_hab
result = result[c(
"commune", "n_catnat", "freq_sin", "cost", "cost_mean", "cost_hab", "ratio", "balance",
"ppri_year", "ppri_state", "ppri_state_sub", "ppri_state_age", "ppri_age_min", "ppri_age_max",
"freq_sin_min", "freq_sin_max", "cost_min", "cost_max", "cost_mean_min", "cost_mean_max",
"cost_hab_min", "cost_hab_max", "ratio_min", "ratio_max"
)]
write.csv2(result, sprintf("data-common/so-ii/onrn/onrn-%s.csv", today), row.names = FALSE) write.csv2(result, sprintf("data-common/so-ii/onrn/onrn-%s.csv", today), row.names = FALSE)
unlink(onrn_raw, recursive = TRUE, force = TRUE) unlink(onrn_raw, recursive = TRUE, force = TRUE)
metropole@data$loss.absolute = metropole@data$POPULATION * (metropole@data$loss.ratio - 1) / metropole@data$area
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment