Continued updates to diagnostics.

derivatives/lossless/code/tools/diagnostics/main_diagnostics.m

@@ -8,7 +8,7 @@ function main_diagnostics(filePath)
     fOut = fopen(fName,'w');
     
     % Update this eventually...
-    header = 'fname,nbchan,pnts,srate,manual,ch_s_sd-m,ch_s_sd_b-m,ch_s_sd_b-m-ch_sd_b,ch_s_sd_b-m-ic_sd1_b,ch_sd-m,ch_sd_b-m,ch_sd_b-m-ic_sd1_b,low_r-m,ic_sd1-m,ic_sd1_b-m,ic_sd2,ch_s_sd,ch_sd,low_r,bridge,LtModelMean,LtModelStd,Lt_quant_0.05,Lt_quant_0.15,Lt_quant_0.25,Lt_quant_0.5,Lt_quant_0.75,Lt_quant_0.85,Lt_quant_0.95,manual_comp_1,ic_rt_1,linkpvaluesMean,linkpvaluesStd,lpv_quant_0.05,lpv_quant_0.15,lpv_quant_0.25,lpv_quant_0.5,lpv_quant_0.75,lpv_quant_0.85,lpv_quant_0.95';
+    header = 'fname,nbchan,pnts,srate,manual,ch_s_sd-m,ch_s_sd_b-m,ch_s_sd_b-m-ch_sd_b,ch_s_sd_b-m-ic_sd1_b,ch_sd-m,ch_sd_b-m,ch_sd_b-m-ic_sd1_b,low_r-m,ic_sd1-m,ic_sd1_b-m,ic_sd2,ch_s_sd,ch_sd,low_r,bridge,LtModelMean,LtModelStd,Lt_quant_0.05,Lt_quant_0.15,Lt_quant_0.25,Lt_quant_0.5,Lt_quant_0.75,Lt_quant_0.85,Lt_quant_0.95,manual_comp_1,ic_rt_1,linkpvaluesMean,linkpvaluesStd,lpv_quant_0.05,lpv_quant_0.15,lpv_quant_0.25,lpv_quant_0.5,lpv_quant_0.75,lpv_quant_0.85,lpv_quant_0.95,tvar_mean, tvar_std, tvar_quant_0.05,tvar_quant_0.15,tvar_quant_0.25,tvar_quant_0.5,tvar_quant_0.75,tvar_quant_0.85,tvar_quant_0.95, pvar_mean, pvar_std,pvar_quant_0.05,pvar_quant_0.15,pvar_quant_0.25,pvar_quant_0.5,pvar_quant_0.75,pvar_quant_0.85,pvar_quant_0.95';
     fprintf(fOut,'%s\n',header);
     
     % Read line by line of the specified diag list

derivatives/lossless/code/tools/diagnostics/single_diagnostic.m

@@ -74,6 +74,28 @@ function single_diagnostic(fOutID, singleSetFile)
         outString = [outString, num2str(quants(i)), ','];
     end
     
+    % Varianc measures:
+    % Total data:
+    EEG = pop_marks_select_data(EEG,'channel marks',[],'labels',{'manual'},'remove','on');
+    EEG = pop_marks_select_data(EEG,'time marks',[],'labels',{'manual'},'remove','on');
+    tvar=std(EEG.data,[],1);
+    outString = [outString, num2str(mean(tvar)), ','];
+    outString = [outString, num2str(std(tvar)), ','];
+    quants = quantile(tvar,quantBreakdown);
+    for i=1:length(quants)
+        outString = [outString, num2str(quants(i)), ','];
+    end
+
+    % Projected data:
+    EEG = pop_marks_select_data(EEG,'component marks',[],'labels',{'manual'},'remove','on');
+    pvar=std(EEG.data,[],1);
+    outString = [outString, num2str(mean(pvar)), ','];
+    outString = [outString, num2str(std(pvar)), ','];
+    quants = quantile(pvar,quantBreakdown);
+    for i=1:length(quants)
+        outString = [outString, num2str(quants(i)), ','];
+    end
+    
     % disp('hold');
     
     % Final write - includes newline