/
missingdatastatsALQ_9_20_2021.m
117 lines (95 loc) · 4.81 KB
/
missingdatastatsALQ_9_20_2021.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
%% Load Data
clear all;close all
clc;
load('USALQ2020_2019_calcresp.mat')
%% Replace missing placeholder values with NaNs
USALQ = [USALQHH2019; USALQHH2020];
for i= 3:width(USALQ)
USALQ.(i)(USALQ.(i)==-9999) = NaN;
end
%% Calculate Average Monthly FCH4 and NEE
x = month(USALQ.TIMESTAMP_END);
monthlyNEE = accumarray(x,USALQ.NEE_F,[],@nanmean);
monthlyFCH4 = accumarray(x,USALQ.FCH4_RF_filled,[],@nanmean);
monthlyNEEnums = accumarray(x,~isnan(USALQ.NEE_F),[]);
monthlyFCH4nums = accumarray(x,~isnan(USALQ.FCH4_RF_filled),[]);
monthlyNEEstd = accumarray(x,USALQ.NEE_F,[],@nanstd);
monthlyFCH4std = accumarray(x,USALQ.FCH4_RF_filled,[],@nanstd);
err = monthlyNEEstd./sqrt(monthlyNEEnums);
err2 = monthlyFCH4std./sqrt(monthlyFCH4nums);
FCH4variance = accumarray(x,USALQ.FCH4_RF_filled,[],@nanvar);
NEEvariance = accumarray(x,USALQ.NEE_F,[],@nanvar);
nanmean(USALQ.FCH4_RF_filled)
%% Create Bar Plots with Error Bars
months = 1:12;
figure()
hBar = bar(months,monthlyNEE,'hist');
hold on
er = errorbar(months,monthlyNEE,err);
er.Color = [0 0 0];
er.LineStyle = 'none';
hold off
title('Bar Plot of Monthly NEE ALQ 2019-20')
xlabel('Month')
ylabel('CO_2 flux (\mumol CO_2 m^-^2s^-^1)')
% Add a colormap to represent variance
CMap = colormap(spring);
set(hBar,'CData',NEEvariance,'CDataMapping','scaled');
c = colorbar;
c.Label.String = 'Variance';
caxis([min(NEEvariance) max(NEEvariance)]);
figure()
hBar2 = bar(months,monthlyFCH4,'hist');
hold on
er = errorbar(months,monthlyFCH4,err2');
er.Color = [0 0 0];
er.LineStyle = 'none';
hold off
title('Bar Plot of Monthly FCH_4 ALQ 2019-20')
xlabel('Month')
ylabel('CH_4 flux (\etamol m^-^2s^-^1)')
% Add a colormap to represent variance
CMap = colormap(winter);
set(hBar2,'CData',FCH4variance,'CDataMapping','scaled');
c = colorbar;
c.Label.String = 'Variance';
caxis([min(FCH4variance) max(FCH4variance)]);
%% Calculate Missing Data
missingNEE19 = sum(isnan(USALQ.NEE(year(USALQ.TIMESTAMP_END)==2019)));
NEE19 = sum(~isnan(USALQ.NEE_F(year(USALQ.TIMESTAMP_END)==2019)));
missingtotalNEE19 = (missingNEE19./length(USALQ.NEE_F(year(USALQ.TIMESTAMP_END)==2019)))*100;
totalNEE19 = (NEE19./length(USALQ.NEE_F(year(USALQ.TIMESTAMP_END)==2019)))*100;
missingNEE20 = sum(isnan(USALQ.NEE(year(USALQ.TIMESTAMP_END)==2020)));
NEE20 = sum(~isnan(USALQ.NEE_F(year(USALQ.TIMESTAMP_END)==2020)));
missingtotalNEE20 = (missingNEE20./length(USALQ.NEE_F(year(USALQ.TIMESTAMP_END)==2020)))*100;
totalNEE20 = (NEE20./length(USALQ.NEE_F(year(USALQ.TIMESTAMP_END)==2020)))*100;
mo = month(USALQ.TIMESTAMP_END);
missingNEE = sum(isnan(USALQ.NEE(ismember(mo,4:10))));
missingtotalNEEseason = (missingNEE./length(USALQ.NEE_F(ismember(mo,4:10))))*100;
missingNEE = sum(isnan(USALQ.NEE(~ismember(mo,4:10))));
missingtotalNEEoffseason = (missingNEE./length(USALQ.NEE_F(~ismember(mo,4:10))))*100;
%% Calculate Missing Data Each Month
TF = 100*sum(isnan(USALQ.FCH4)/length(USALQ.FCH4))
TF = 100*sum(isnan(USALQ.FCH4(ismember(x,4:10))))/length(USALQ.FCH4(ismember(x,4:10)))
TF = 100*sum(isnan(USALQ.FCH4(~ismember(x,4:10))))/length(USALQ.FCH4(~ismember(x,4:10)))
TF = 100*sum(isnan(USALQ.FCH4_RF_filled)/length(USALQ.FCH4_RF_filled))
x = month(USALQ.TIMESTAMP_END(year(USALQ.TIMESTAMP_END)==2019));
missingNEE19 = accumarray(x,isnan(USALQ.NEE(year(USALQ.TIMESTAMP_END)==2019)),[]);
monthlyNEE19 = accumarray(x,~isnan(USALQ.NEE(year(USALQ.TIMESTAMP_END)==2019)),[]);
missingtotalNEE19 = (missingNEE19./monthlyNEE19)*100;
missingFCH419 = accumarray(x,isnan(USALQ.FCH4(year(USALQ.TIMESTAMP_END)==2019)),[]);
monthlyFCH419 = accumarray(x,~isnan(USALQ.FCH4(year(USALQ.TIMESTAMP_END)==2019)),[]);
missingtotalFCH419 = (missingFCH419./monthlyFCH419)*100;
TF = 100*sum(isnan(USALQ.FCH4(year(USALQ.TIMESTAMP_END)==2019))/length(USALQ.FCH4(year(USALQ.TIMESTAMP_END)==2019)))
yr = 2019*ones(12,1);
x2 = month(USALQ.TIMESTAMP_END(year(USALQ.TIMESTAMP_END)==2020));
missingNEE20 = accumarray(x2,isnan(USALQ.NEE(year(USALQ.TIMESTAMP_END)==2020)),[]);
monthlyNEE20 = accumarray(x2,~isnan(USALQ.NEE(year(USALQ.TIMESTAMP_END)==2020)),[]);
missingtotalNEE20 = (missingNEE20./monthlyNEE20)*100;
missingFCH420 = accumarray(x2,isnan(USALQ.FCH4(year(USALQ.TIMESTAMP_END)==2020)),[]);
monthlyFCH420 = accumarray(x2,~isnan(USALQ.FCH4(year(USALQ.TIMESTAMP_END)==2020)),[]);
missingtotalFCH420 = (missingFCH420./monthlyFCH420)*100;
TF = 100*sum(isnan(USALQ.FCH4(year(USALQ.TIMESTAMP_END)==2020))/length(USALQ.FCH4(year(USALQ.TIMESTAMP_END)==2020)))
yr2 = 2020*ones(length(unique(x2)),1);
T = table([yr;yr2], [unique(x);unique(x2)] ,[missingtotalNEE19;missingtotalNEE20] ,[missingtotalFCH419;missingtotalFCH420]);
T.Properties.VariableNames = {'Year','Month','Missing NEE','Missing FCH4'};