data(data_HPI)
data_HPI = data_HPI[complete.cases(data_HPI),]
data_HPI_2019 = data_HPI[data_HPI$year==2019,]
Indic_name = c("Life_Expectancy","Ladder_of_life","Ecological_Footprint")
Indic_norm = data.frame("ISO"=data_HPI_2019$ISO,
normalise_ci(data_HPI_2019[, Indic_name],
c(1:3),
c("POS","POS","NEG"),
method=2)$ci_norm)
Indic_norm = Indic_norm[Indic_norm$Life_Expectancy>0 &
Indic_norm$Ladder_of_life>0 &
Indic_norm$Ecological_Footprint >0 ,]
atleast = 2
CI_owa_n = ci_owa(Indic_norm, id="ISO",
indic_col=c(2:4),
atleastjp=atleast)$CI_OWA_n
CI_owa_p = ci_owa(Indic_norm, id="ISO",
indic_col=c(2:4),
atleastjp=atleast)$CI_OWA_p
Run the code above in your browser using DataLab