Remove printing to permit bootci

mAGREE.m

@@ -45,17 +45,17 @@
 %% Check for valid data from multiple raters
 if n < 1
     P_O = NaN;
-    fprintf('ERROR: At least 1 valid item is required. \n')
+    disp('ERROR: At least 1 valid item is required.');
     return;
 end
 if r < 2
     P_O = NaN;
-    fprintf('ERROR: At least 2 raters are required. \n');
+    disp('ERROR: At least 2 raters are required.');
     return;
 end
 if any(ismember(x, CATEGORIES) == 0)
     P_O = NaN;
-    fprintf('ERROR: Unexpected category in CODES. \n');
+    disp('ERROR: Unexpected category in CODES.');
     return;
 end
 %% Get weights from mWEIGHTING function
@@ -73,7 +73,5 @@
 nprime = sum(r_i >= 2);
 possible = r_i .* (r_i - 1);
 P_O = sum(observed(r_i >= 2) ./ (possible(r_i >= 2))) ./ nprime;
-%% Output reliability and variance components
-fprintf('Percent observed agreement = %.3f \n', P_O);
 
 end

mALPHAK.m

@@ -87,9 +87,5 @@
 P_C = sum(sum(weights .* (pihat * transpose(pihat))));
 %% Calculate reliability point estimate
 ALPHAK = (P_O - P_C) / (1 - P_C);
-%% Output reliability and variance components
-fprintf('Percent observed agreement = %.3f \n',P_O);
-fprintf('Percent chance agreement = %.3f \n\n',P_C);
-fprintf('Krippendorff''s alpha coefficient = %.3f \n', ALPHAK);
 
 end

mGAMMA.m

@@ -82,9 +82,5 @@
 P_C = sum(weights(:)) .* sum(pihat .* (1 - pihat)) ./ (q * (q - 1));
 %% Calculate reliability point estimate
 GAMMA = (P_O - P_C) / (1 - P_C);
-%% Output reliability and variance components
-fprintf('Percent observed agreement = %.3f \n', P_O);
-fprintf('Percent chance agreement = %.3f \n\n', P_C);
-fprintf('Gwet''s gamma coefficient = %.3f \n', GAMMA);
 
 end

mKAPPA.m

@@ -37,9 +37,9 @@
 x(~isfinite(x)) = [];
 if nargin < 2
     CATEGORIES = x;
-    WEIGHTING = 'nominal';
+    WEIGHTING = 'identity';
 elseif nargin < 3
-    WEIGHTING = 'nominal';
+    WEIGHTING = 'identity';
 end
 if isempty(CATEGORIES)
     CATEGORIES = x;
@@ -91,9 +91,5 @@
 P_C = sum(sum(weights .* (pbar_k * transpose(pbar_k) - ssq_kl ./ r)));
 %% Calculate reliability point estimate
 KAPPA = (P_O - P_C) / (1 - P_C);
-%% Output reliability and variance components
-fprintf('Percent observed agreement = %.3f \n', P_O);
-fprintf('Percent chance agreement = %.3f \n\n', P_C);
-fprintf('Cohen''s kappa coefficient = %.3f \n', KAPPA);
 
 end

mSCOTTPI.m

@@ -83,9 +83,5 @@
 P_C = sum(sum(weights .* (pihat * transpose(pihat))));
 %% Calculate reliability point estimate
 SCOTTPI = (P_O - P_C) / (1 - P_C);
-%% Output reliability and variance components
-fprintf('Percent observed agreement = %.3f \n', P_O);
-fprintf('Percent chance agreement = %.3f \n\n', P_C);
-fprintf('Scott''s pi coefficient = %.3f \n', SCOTTPI);
 
 end

mSSCORE.m

@@ -84,9 +84,5 @@
 P_C = sum(weights(:)) / (q ^ 2);
 %% Calculate reliability point estimate
 SSCORE = (P_O - P_C) / (1 - P_C);
-%% Output reliability and variance components
-fprintf('Percent observed agreement = %.3f \n', P_O);
-fprintf('Percent chance agreement = %.3f \n\n', P_C);
-fprintf('Bennett et al.''s S score = %.3f \n', SSCORE);
 
 end