In a Statistical Analysis Plan (SAP) the following statement is written:
To test the association between the pain intensity reported by the patient and the treatment the patient has been administered, a Cochran-Mantel-Haenszel test will be performed adjusting for the country of the patient's investigational site.
The following variables and levels are used:
- Pain intensity (PAIN): "No pain", "Mild pain", "Moderate pain", "Severe pain", "Most severe pain"
- Treatment group (TREAT): "Placebo", "Drug A", "Drug B"
- Country (COUNTRY): Eight different values.
The following program is written to perform this analysis:
proc freq data=pain;
by COUNTRY;
tables PAIN * TREAT / cmh;
run;
This code does not produce the analysis requested by the SAP.
What is wrong with this code?