Ef_emfac

[Joaobazzo]

• [x] Improve efficiency (too slow)

• [x] Produce a different type of data(emfac) so it can be faster

• [x] Check EF data from 2005, which seems to be wrong

• [x] Export in 'units' class

• [ ] if aggregate == FALSE, no need for veh input

[Joaobazzo]

Now it's looking like this @rafapereirabr @pedro-andrade-inpe . I think it's looking a bit better

library(gtfs2emis)
gtfs2emis::ef_emfac(pol = c("CO","PM10"),calendar_year = "2019",
                    model_year = c(2005,2006,2007),aggregate = FALSE,
                    veh = rep(0.33,1), # sum = 1
                    speed = units::set_units(c(10,20,30),'km/h')) 
#>      CO_2005   CO_2006   CO_2007   PM10_2005   PM10_2006   PM10_2007
#> 1: 0.6421423 0.9269853 0.6439274 0.004818207 0.005758708 0.008073102
#> 2: 0.3725280 0.5374732 0.4534922 0.004397409 0.005254205 0.006937188
#> 3: 0.2218897 0.3199493 0.3280440 0.003991091 0.004767238 0.006188870
gtfs2emis::ef_emfac(pol = c("CO","PM10"),calendar_year = "2019",
                    model_year = c(2005,2006,2007),aggregate = TRUE,
                    veh = rep(0.33,1), # sum = 1
                    speed = units::set_units(c(10,20,30),'km/h'))
#>       CO_avg    PM10_avg
#> 1: 0.2119069 0.001590008
#> 2: 0.1229342 0.001451145
#> 3: 0.0732236 0.001317060

^{Created on 2020-05-18 by the reprex package (v0.3.0)}

[Joaobazzo]

check cut function