{
data(net)
net <- net[1:100, ]
data(pc_profile)
x <- age_ldv(net$ldv)
pc_week <- temp_fact(net$ldv+net$hdv, pc_profile[[1]])
df <- netspeed(pc_week,
net$ps,
net$ffs,
net$capacity,
net$lkm,
alpha = 1)
s <- do.call("rbind",lapply(1:ncol(df), function(i) {
as.data.frame(replicate(ncol(x), df[, i]))
}))
ef <- ef_wear(wear = "tyre",
type = "PC",
pol = "PM10",
speed = as.data.frame(s))
e <- emis_long(x = x,
lkm = net$lkm,
ef = ef,
tfs = pc_profile[[1]],
speed = df)
ae <- emis_long(x = x,
lkm = net$lkm,
ef = ef,
tfs = pc_profile[[1]],
speed = df,
array = TRUE)
}
Run the code above in your browser using DataLab