about the "attribute_health.R" function

[enricopisoni]

Dears

I am checking the package...very well structured and clear, thanks for this !

I am also thinking how to apply it to my work

my first reaction is that probably (as you also said) I would mainly need the function "attribute_health.R" to be run for all the cells of my grid

in my case... what I would have are:

• data that are not spatially varying (as RR and cut-off) and
• others daa that are spatially varying (exposure, and i.e. baseline mortality rate)

it seems at the moment that, through the input to the "attribute_health.R" function( input "exp_central"), I can define a vector of exposure (in my case it would be an element for each cell) could I do the same for the "baseline mortality rate"? as now it seems as input I can only provide one value for all exposure elements (while in my case i.e. I could have a different "baseline mortality rate" by NUTS3

let me know if this is possible...otherwise I will simply loop on the "attribute_health.R" function, passing each time the 4 requested info, by cell:

• exposure (spatially varying)
• baseline health data (spatially varying)
• RR (fixed)
• cut-off (fixed)

let me know thanks Enrico

PS in any case, as soon as I have a code working on my sample dataset, I will share it with you, fyi...thanks!

[enricopisoni]

Dears @luytax @ungatoverde have you seen my previous message? it's not urgent...but I am curious to see your opinion on this ;-) thanks bests Enrico

[ungatoverde]

Hi @enricopisoni! Nice to see that you find our R pacakge useful for your work.

If you want to assign data to geographic units, you just need to insert the data as list elements. E.g. as.list(c(8,9,10)) means that you have the exposure in your three geographic units are 8, 9 and 10 ug/m3 respectively.

The geographic areas need an id to join all data refering to the same geographic unit. For this purpose you have to use the argument geo_id_raw with a vector e.g. c(01, 02, 03). In this way, the R package can interpret that the exposure in 01 is 8 and so on. Am I answering your question?

Sorry, we need to document the whole package much better (TBD next year).