How to Calculate Confidence Intervals for Crosstabs with Categorical Variable Sets

// Choose which stats to replace based on those available for categorical data
var upper_statistic_to_replace = "Index";
var lower_statistic_to_replace = "Missing n";

// Confirm the table is a crosstab using categorical data
if (table.availableStatistics.indexOf("Column %") == -1) {
table.suppressOutput('Confidence intervals not computed - it only applies to tables containing Column % in the Statistics - Cells');
form.ruleNotApplicable('it only applies to tables containing Column % in the Statistics - Cells');
}

// Confirm the categorical data used is either pick one or pick any -- other types not supported
if (["Pick One", "Pick Any"].indexOf(table.blueQuestion.questionType) == -1 || ["Pick One", "Pick Any"].indexOf(table.brownQuestion.questionType) == -1){
table.suppressOutput('Confidence intervals not computed - the table should have a Nominal or Binary-Multi variable set selected in both the Rows and Columns of the table.');
form.ruleNotApplicable('the table should have a Nominal or Binary-Multi variable set selected in both the Rows and Columns of the table.');
}

//add titles to rule
form.setSummary('Calculate confidence intervals on crosstabs with percentages');
form.setHeading('Calculate confidence intervals on crosstabs with percentages');

// Set up controls for user input to select the significance level
var label = form.newLabel('Significance level:');
var p_level = form.newNumericUpDown('Lower'); // create numeric control
p_level.setIncrement(0.01);
p_level.setDefault(0.05);
p_level.setMaximum(1);
p_level.setMinimum(0);
p_level.lineBreakAfter = true;

let t_box = form.newCheckBox('uset', 'Use Survey Reporter / Quantum formula');
t_box.setDefault(false);

form.setInputControls([label, p_level, t_box]); // add controls to settings shown


// get the value from the input form above and set it as the cutoff variable
var p_cutoff = p_level.getValue();
let use_t_dist = t_box.getValue();

// Get relevant table stats
var n = table.get("n");
var col_n = table.get("Column n");
var col_percent = table.get("Column %");
var col_se = table.get("Column Standard Error");
var upper_values = table.get(upper_statistic_to_replace);
var lower_values = table.get(lower_statistic_to_replace);
var isWeight = table.weight != null;

// Do the calculations on each cell in the table
for (var row = 0; row < table.numberRows; row++) {
for (var col = 0; col < table.numberColumns; col++) {
if (use_t_dist) {
let df = col_n[row][col] - 1;
let t = jStat.studentt.inv(p_cutoff/2, df);
let increment = t * col_se[row][col] * 100;
lower_values[row][col] = col_percent[row][col] + increment;
upper_values[row][col] = col_percent[row][col] - increment;

} else {
var z = jStat.normal.inv(p_cutoff/2,0,1);

if (isWeight) {
var increment = z*col_se[row][col]*100;
lower_values[row][col] = col_percent[row][col] + increment;
upper_values[row][col] = col_percent[row][col] - increment;
} else {
var calc1 = col_n[row][col] + Math.pow(z,2);
var calc2 = (n[row][col] + Math.pow(z,2)/2)/calc1;
var calc3 = Math.sqrt(calc2*(1-calc2)/calc1);
var increment = z*calc3;
lower_values[row][col] = (calc2 + increment)*100;
upper_values[row][col] = (calc2 - increment)*100;
}
}

}
}

//set the values of the statistics we're overwriting to the new values on the table
table.set(lower_statistic_to_replace, lower_values);
table.set(upper_statistic_to_replace, upper_values);

//see if statistics for upper and lower CI are showing on table and if not, add them
var stats = table.statistics;
if(stats.indexOf(upper_statistic_to_replace) == -1) stats.push(upper_statistic_to_replace);
if(stats.indexOf(lower_statistic_to_replace) == -1) stats.push(lower_statistic_to_replace);
table.statistics = stats;

//rename the statistics we're overwriting to upper and lower CI
form.setTranslation(upper_statistic_to_replace, "Upper Confidence Interval");
form.setTranslation(lower_statistic_to_replace, "Lower Confidence Interval");