Upload New File

weighted_age_metrics.m

+%
+% Copyright 2019-2020 LIST (Luxembourg Institute of Science and Technology), all right reserved.
+%
+% Author: Nicolas Rodriguez (nicolas.bjorn.rodriguez@gmail.com)
+%
+function [FYW,MA,mA]=weighted_age_metrics(MD,frc,dT,weights)
+
+    % MD: master distribution
+    % MA: mean age
+    % mA: median age
+
+    temp=find(cumsum(MD*dT)<=0.5);
+    if ~isempty(temp)
+        mA=temp(end)*dT; % hours
+    else
+        mA=NaN;
+    end
+    temp=(dT:dT:dT*length(MD)); % vector of ages (hours)
+    MA=trapz(temp,temp.*MD); % hours
+    [r,~]=size(frc);
+    FYW=NaN(r,1);
+    for k=1:r
+        FYW(k,1)=sum(frc(k,:).*weights')/sum(weights);
+    end
+    
+end
\ No newline at end of file