search

eliocamp / shceof

5 stars 2 forks source link

Hacer todas los promedios pesados por la cantidad de días #23

Closed eliocamp closed 3 years ago

eliocamp commented 3 years ago

Esto:

 .[, .(hgt = mean(hgt)), by = .(lev, lon, lat, time = seasonally(time))] %>%

tiene que ser 

 .[, .(hgt = weighted.mean(hgt, days_in_month)), by = .(lev, lon, lat, time = seasonally(time))] %>%
eliocamp commented 3 years ago

Esto se complica un poco porque es weighted.mean muchísimo más lento que mean (probablemente porque data.table optimiza un montón la primera expresión).

La alternativa es calcular los pesos a parte, multiplicar y luego tomar la media, pero es un cambio más rompebolas.