So I have a 22*9 matrix. with 22 regions and 9 years. I am doing a decomposition analysis where need to calculate different factors contribution to the overalla change in total carbon emissions, between each set of years, for each region.
For some reason however, my code returns complex answers (for both Ln(r,m) and DeltaCO2CI) when there should only be real numbers. Can't figure out what i have done wrong?
here's the problematic part of the code: