0% found this document useful (0 votes)
37 views23 pages

Malaria Analysis with SAS Code

The document outlines a SAS program for analyzing data from a malaria survey, detailing the execution environment and the steps taken to import, clean, and analyze the data. Key steps include data import, variable creation, descriptive analysis, bivariate associations, and logistic regression modeling. The program aims to investigate the relationship between malaria positivity and various factors such as IRS spraying and bednet ownership.

Uploaded by

Pratik Pujari
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
37 views23 pages

Malaria Analysis with SAS Code

The document outlines a SAS program for analyzing data from a malaria survey, detailing the execution environment and the steps taken to import, clean, and analyze the data. Key steps include data import, variable creation, descriptive analysis, bivariate associations, and logistic regression modeling. The program aims to investigate the relationship between malaria positivity and various factors such as IRS spraying and bednet ownership.

Uploaded by

Pratik Pujari
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd

5/15/25, 12:13 PM Program Summary - code.

sas

Program Summary - [Link]

Execution Environment

Author: u64151009
File: /home/u64151009/sasuser.v94/ZWhite/PROJ/[Link]
SAS Platform: Linux LIN X64 5.14.0-284.30.1.el9_2.x86_64
SAS Host: [Link]
SAS Version: 9.04.01M7P08062020
SAS Locale: en_US
Submission Time: 5/15/2025, 12:13:25 PM
Browser Host: [Link]
User Agent: Mozilla/5.0 (Windows NT 10.0; Win64; x64) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/[Link] Safari/537.36
Application Server: [Link]

Code: [Link]

/* STEP 1: IMPORT DATA */


proc import
datafile="/home/u64151009/sasuser.v94/BioStat/Files/ngenirs-mopeia-cross-sectional-survey_2018.csv"
out=mopeia_raw
dbms=csv
replace;
getnames=yes;
guessingrows=max;
run;

/* STEP 2: CLEANING & VARIABLE CREATION */


data mopeia;
set mopeia_raw;

/* EXCLUSIONS: Remove missing RDT or spray status */


if missing('[Link]'n) or missing('PCD Reference File.2017 Spray St'n) then delete;

/* OUTCOME: Malaria (RDT positive) */


MalariaPos = '[Link]'n;

/* EXPOSURE: IRS in 2017 */


if 'PCD Reference File.2017 Spray St'n = 'Spray' then IRSsprayed = 1;
else if 'PCD Reference File.2017 Spray St'n = 'No_Spray' then IRSsprayed = 0;
else IRSsprayed = .;

/* BEDNET OWNERSHIP */
if 'MAIN_MALARIA_CONTROL_HAS_NETS'n = 1 then BedNet = 1;
else if 'MAIN_MALARIA_CONTROL_HAS_NETS'n = 2 then BedNet = 0;
else BedNet = .;

/* AGE GROUP (already categorized) */


AgeGroup = AGE_GROUP;

/* GENDER */
Gender = DEMOGRAFIC_GENDER;

/* LITERACY */
ReadAble = '[Link]'n;
WriteAble = '[Link]'n;
if ReadAble = 1 and WriteAble = 1 then Literate = 1;
else if ReadAble = 2 or WriteAble = 2 then Literate = 0;
else Literate = 0;

/* HOUSING VARIABLES */
MetalRoof = '[Link]'n;
DurableWall = ('[Link] bricks'n = 1 or '[Link] block'n = 1);

/* LABELS FOR CLARITY */


label MalariaPos = "Malaria RDT Positive"
IRSsprayed = "IRS Sprayed (2017)"
BedNet = "Household Owns Bed Net"
AgeGroup = "Age Group"
Gender = "Gender"
Literate = "Can Read and Write"
MetalRoof = "Zinc Roof"
DurableWall = "Durable Walls";

about:blank 1/23
5/15/25, 12:13 PM Program Summary - [Link]
keep MalariaPos IRSsprayed BedNet AgeGroup Gender Literate MetalRoof DurableWall;
run;

/* STEP 3: EXPORT CLEAN DATASET FOR ANALYSIS */


libname save '/home/u64151009/sasuser.v94/ZWhite/Files';

data save.clean_mopeia;
set mopeia;

/* STEP 4: DESCRIPTIVE ANALYSIS */

/* 4a. Distributions */
proc freq data=mopeia;
tables MalariaPos IRSsprayed BedNet AgeGroup Gender Literate MetalRoof DurableWall / missing;
title "Overall Distribution of Analysis Variables";
run;

/* 4b. Table 1: Stratified by IRSsprayed */


proc freq data=mopeia;
tables IRSsprayed*(MalariaPos BedNet AgeGroup Gender Literate MetalRoof DurableWall)
/ chisq norow nocol nopercent;
title "Table 1: Baseline Characteristics by IRSsprayed";
run;

/* STEP 5: BIVARIATE ASSOCIATIONS */


proc freq data=mopeia;
tables BedNet*MalariaPos / chisq relrisk;
title "Malaria and BedNet Ownership";
run;

proc freq data=mopeia;


tables IRSsprayed*MalariaPos / chisq relrisk;
title "Malaria and IRS Spraying";
run;

proc freq data=mopeia;


tables AgeGroup*MalariaPos Gender*MalariaPos Literate*MalariaPos
MetalRoof*MalariaPos DurableWall*MalariaPos / chisq;
title "Other Variables and Malaria Status";
run;

/* STEP 6: LOGISTIC REGRESSION */

/* 6a. Unadjusted: IRSsprayed only */


proc logistic data=mopeia descending plots(only)=roc;
class IRSsprayed (ref='0') / param=ref;
model MalariaPos = IRSsprayed / clodds=wald lackfit outroc=roc1;
title "Unadjusted Model: IRSsprayed → Malaria";
run;

/* 6b. Adjusted model */


proc logistic data=mopeia descending plots(only)=roc;
class IRSsprayed (ref='0')
BedNet (ref='0')
AgeGroup (ref='1')
Gender (ref='1')
Literate (ref='0')
MetalRoof (ref='0')
DurableWall (ref='0')
/ param=ref;
model MalariaPos = IRSsprayed BedNet AgeGroup Gender Literate MetalRoof DurableWall
/ clodds=wald lackfit outroc=roc2;
title "Adjusted Logistic Model: IRSsprayed & Covariates";
run;

/* STEP 7: INTERACTIONS & MODEL DIAGNOSTICS */

/* IRSsprayed × BedNet */
proc logistic data=mopeia descending;
class IRSsprayed (ref='0')
BedNet (ref='0')
AgeGroup (ref='1')
Gender (ref='1')
Literate (ref='0')
about:blank 2/23
5/15/25, 12:13 PM Program Summary - [Link]
MetalRoof (ref='0')
DurableWall (ref='0')
/ param=ref;
model MalariaPos = IRSsprayed|BedNet AgeGroup Gender Literate MetalRoof DurableWall
/ clodds=wald lackfit;
title "Interaction: IRSsprayed × BedNet";
run;

/* Goodness-of-fit and ROC */


ods graphics on;
ods select FitStatistics Association LackFitChiSq;
proc logistic data=mopeia descending plots(only)=roc(id=prob);
class IRSsprayed (ref='0')
BedNet (ref='0')
AgeGroup (ref='1')
Gender (ref='1')
Literate (ref='0')
MetalRoof (ref='0')
DurableWall (ref='0')
/ param=ref;

Log: [Link]

Notes (46)

1 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;


NOTE: ODS statements in the SAS Studio environment may disable some output features.
69
70 /* STEP 1: IMPORT DATA */
71 proc import
72 datafile="/home/u64151009/sasuser.v94/BioStat/Files/ngenirs-mopeia-cross-sectional-survey_2018.csv"
73 out=mopeia_raw
74 dbms=csv
75 replace;
76 getnames=yes;
77 guessingrows=max;
78 run;

NOTE: Unable to open parameter catalog: [Link] in update mode. Temporary parameter values will be saved to
[Link].
Name MAIN_CLINIC_DATA_SAMPLES_TEST_RESULT truncated to MAIN_CLINIC_DATA_SAMPLES_TEST_RE.
Name MAIN_MALARIA_CONTROL_MEASURES_LAST WEEK truncated to MAIN_MALARIA_CONTROL_MEASURES_LA.
Name PCD Reference File.2016 Spray Status truncated to PCD Reference File.2016 Spray St.
Name PCD Reference File.2017 Spray Status truncated to PCD Reference File.2017 Spray St.
Name PCD Reference [Link] Health Facility truncated to PCD Reference [Link] Healt.
Name PCD Reference [Link] To Closest Health Facility truncated to PCD Reference [Link] To Closest.
Problems were detected with provided names. See LOG.
79 /**********************************************************************
80 * PRODUCT: SAS
81 * VERSION: 9.4
82 * CREATOR: External File Interface
83 * DATE: 15MAY25
84 * DESC: Generated SAS Datastep Code
85 * TEMPLATE SOURCE: (None Specified.)
86 ***********************************************************************/
87 data WORK.MOPEIA_RAW ;
88 %let _EFIERR_ = 0; /* set the ERROR detection macro variable */
89 infile '/home/u64151009/sasuser.v94/BioStat/Files/ngenirs-mopeia-cross-sectional-survey_2018.csv' delimiter = ','
89 ! MISSOVER DSD lrecl=32767 firstobs=2 ;
90 informat DEMOGRAFIC_GENDER best32. ;
91 informat "[Link]"N best32. ;
92 informat MAIN_HAD_MALARIA best32. ;
93 informat "[Link]"N best32. ;
94 informat MAIN_WRITE best32. ;
95 informat "[Link]"N best32. ;
96 informat "[Link]"N best32. ;
97 informat "[Link] plastic"N best32. ;
98 informat "[Link] covered with sand"N best32. ;
99 informat "[Link] bark"N best32. ;
100 informat "[Link] and zinc"N best32. ;
101 informat "[Link] bricks"N best32. ;
102 informat "[Link]"N best32. ;
103 informat "[Link] block"N best32. ;
104 informat "[Link]"N best32. ;
105 informat MAIN_ADMIN_TIME best32. ;
106 informat MAIN_ROOF best32. ;
107 informat "[Link]"N best32. ;
108 informat "[Link]"N best32. ;
109 informat "[Link]"N best32. ;
110 informat "[Link]"N best32. ;
111 informat "[Link]"N best32. ;
112 informat "[Link]"N best32. ;
113 informat "[Link]"N best32. ;
114 informat MAIN_OCUPATION best32. ;

about:blank 3/23
5/15/25, 12:13 PM Program Summary - [Link]
115 informat "[Link]"N best32. ;
116 informat "[Link]"N best32. ;
117 informat "[Link]"N best32. ;
118 informat "[Link]"N best32. ;
119 informat "[Link]"N best32. ;
120 informat "[Link]"N best32. ;
121 informat "[Link]"N best32. ;
122 informat "[Link]"N best32. ;
123 informat PERM_ID $16. ;
124 informat CLUSTER_ID best32. ;
125 informat "[Link]"N best32. ;
126 informat MAIN_HF_TRAV_TIME best32. ;
127 informat MAIN_EDUCATION best32. ;
128 informat "[Link]"N best32. ;
129 informat "Ed.EP1"N best32. ;
130 informat "Ed.EP2"N best32. ;
131 informat "Ed.ESG1"N best32. ;
132 informat "Ed.ESG2"N best32. ;
133 informat "[Link]"N best32. ;
134 informat "[Link] basico"N best32. ;
135 informat "[Link] Medio"N best32. ;
136 informat "[Link]"N best32. ;
137 informat "[Link]"N best32. ;
138 informat "[Link]"N best32. ;
139 informat FAMILY_ID $12. ;
140 informat DEMOGRAFIC_DATE_VISIT anydtdtm40. ;
141 informat AGE_GROUP best32. ;
142 informat "Age.u5"N best32. ;
143 informat "Age.o5"N best32. ;
144 informat VILLAGE_ID best32. ;
145 informat "[Link]"N best32. ;
146 informat MAIN_MALARIA_CONTROL_HAS_NETS best32. ;
147 informat MAIN_MALARIA_CONTROL_NET_COUNT best32. ;
148 informat "[Link].48"N best32. ;
149 informat "[Link]"N best32. ;
150 informat "[Link]"N best32. ;
151 informat MAIN_CLINIC_DATA_SAMPLES_TEST_RE best32. ;
152 informat MAIN_MALARIA_CONTROL_NETS best32. ;
153 informat MAIN_MALARIA_CONTROL_MEASURES_LA best32. ;
154 informat "[Link]"N best32. ;
155 informat "[Link]"N best32. ;
156 informat "[Link]"N best32. ;
157 informat "[Link]"N best32. ;
158 informat "PCD Reference File.2016 Spray St"N $8. ;
159 informat "PCD Reference File.2017 Spray St"N $8. ;
160 informat "PCD Reference [Link]"N $25. ;
161 informat "[Link]"N best32. ;
162 informat "PCD Reference [Link]"N $6. ;
163 informat "[Link]"N best32. ;
164 informat "[Link]"N best32. ;
165 informat "[Link]"N best32. ;
166 informat "[Link]"N best32. ;
167 informat "PCD Reference [Link]"N best32. ;
168 informat "PCD Reference [Link]"N best32. ;
169 informat "PCD Reference [Link] Healt"N $30. ;
170 informat "PCD Reference [Link] To Closest"N best32. ;
171 format DEMOGRAFIC_GENDER best12. ;
172 format "[Link]"N best12. ;
173 format MAIN_HAD_MALARIA best12. ;
174 format "[Link]"N best12. ;
175 format MAIN_WRITE best12. ;
176 format "[Link]"N best12. ;
177 format "[Link]"N best12. ;
178 format "[Link] plastic"N best12. ;
179 format "[Link] covered with sand"N best12. ;
180 format "[Link] bark"N best12. ;
181 format "[Link] and zinc"N best12. ;
182 format "[Link] bricks"N best12. ;
183 format "[Link]"N best12. ;
184 format "[Link] block"N best12. ;
185 format "[Link]"N best12. ;
186 format MAIN_ADMIN_TIME best12. ;
187 format MAIN_ROOF best12. ;
188 format "[Link]"N best12. ;
189 format "[Link]"N best12. ;
190 format "[Link]"N best12. ;
191 format "[Link]"N best12. ;
192 format "[Link]"N best12. ;
193 format "[Link]"N best12. ;
194 format "[Link]"N best12. ;
195 format MAIN_OCUPATION best12. ;
196 format "[Link]"N best12. ;
197 format "[Link]"N best12. ;
198 format "[Link]"N best12. ;
199 format "[Link]"N best12. ;
200 format "[Link]"N best12. ;
201 format "[Link]"N best12. ;
202 format "[Link]"N best12. ;

about:blank 4/23
5/15/25, 12:13 PM Program Summary - [Link]
203 format "[Link]"N best12. ;
204 format PERM_ID $16. ;
205 format CLUSTER_ID best12. ;
206 format "[Link]"N best12. ;
207 format MAIN_HF_TRAV_TIME best12. ;
208 format MAIN_EDUCATION best12. ;
209 format "[Link]"N best12. ;
210 format "Ed.EP1"N best12. ;
211 format "Ed.EP2"N best12. ;
212 format "Ed.ESG1"N best12. ;
213 format "Ed.ESG2"N best12. ;
214 format "[Link]"N best12. ;
215 format "[Link] basico"N best12. ;
216 format "[Link] Medio"N best12. ;
217 format "[Link]"N best12. ;
218 format "[Link]"N best12. ;
219 format "[Link]"N best12. ;
220 format FAMILY_ID $12. ;
221 format DEMOGRAFIC_DATE_VISIT datetime. ;
222 format AGE_GROUP best12. ;
223 format "Age.u5"N best12. ;
224 format "Age.o5"N best12. ;
225 format VILLAGE_ID best12. ;
226 format "[Link]"N best12. ;
227 format MAIN_MALARIA_CONTROL_HAS_NETS best12. ;
228 format MAIN_MALARIA_CONTROL_NET_COUNT best12. ;
229 format "[Link].48"N best12. ;
230 format "[Link]"N best12. ;
231 format "[Link]"N best12. ;
232 format MAIN_CLINIC_DATA_SAMPLES_TEST_RE best12. ;
233 format MAIN_MALARIA_CONTROL_NETS best12. ;
234 format MAIN_MALARIA_CONTROL_MEASURES_LA best12. ;
235 format "[Link]"N best12. ;
236 format "[Link]"N best12. ;
237 format "[Link]"N best12. ;
238 format "[Link]"N best12. ;
239 format "PCD Reference File.2016 Spray St"N $8. ;
240 format "PCD Reference File.2017 Spray St"N $8. ;
241 format "PCD Reference [Link]"N $25. ;
242 format "[Link]"N best12. ;
243 format "PCD Reference [Link]"N $6. ;
244 format "[Link]"N best12. ;
245 format "[Link]"N best12. ;
246 format "[Link]"N best12. ;
247 format "[Link]"N best12. ;
248 format "PCD Reference [Link]"N best12. ;
249 format "PCD Reference [Link]"N best12. ;
250 format "PCD Reference [Link] Healt"N $30. ;
251 format "PCD Reference [Link] To Closest"N best12. ;
252 input
253 DEMOGRAFIC_GENDER
254 "[Link]"N
255 MAIN_HAD_MALARIA
256 "[Link]"N
257 MAIN_WRITE
258 "[Link]"N
259 "[Link]"N
260 "[Link] plastic"N
261 "[Link] covered with sand"N
262 "[Link] bark"N
263 "[Link] and zinc"N
264 "[Link] bricks"N
265 "[Link]"N
266 "[Link] block"N
267 "[Link]"N
268 MAIN_ADMIN_TIME
269 MAIN_ROOF
270 "[Link]"N
271 "[Link]"N
272 "[Link]"N
273 "[Link]"N
274 "[Link]"N
275 "[Link]"N
276 "[Link]"N
277 MAIN_OCUPATION
278 "[Link]"N
279 "[Link]"N
280 "[Link]"N
281 "[Link]"N
282 "[Link]"N
283 "[Link]"N
284 "[Link]"N
285 "[Link]"N
286 PERM_ID $
287 CLUSTER_ID
288 "[Link]"N
289 MAIN_HF_TRAV_TIME
290 MAIN_EDUCATION

about:blank 5/23
5/15/25, 12:13 PM Program Summary - [Link]
291 "[Link]"N
292 "Ed.EP1"N
293 "Ed.EP2"N
294 "Ed.ESG1"N
295 "Ed.ESG2"N
296 "[Link]"N
297 "[Link] basico"N
298 "[Link] Medio"N
299 "[Link]"N
300 "[Link]"N
301 "[Link]"N
302 FAMILY_ID $
303 DEMOGRAFIC_DATE_VISIT
304 AGE_GROUP
305 "Age.u5"N
306 "Age.o5"N
307 VILLAGE_ID
308 "[Link]"N
309 MAIN_MALARIA_CONTROL_HAS_NETS
310 MAIN_MALARIA_CONTROL_NET_COUNT
311 "[Link].48"N
312 "[Link]"N
313 "[Link]"N
314 MAIN_CLINIC_DATA_SAMPLES_TEST_RE
315 MAIN_MALARIA_CONTROL_NETS
316 MAIN_MALARIA_CONTROL_MEASURES_LA
317 "[Link]"N
318 "[Link]"N
319 "[Link]"N
320 "[Link]"N
321 "PCD Reference File.2016 Spray St"N $
322 "PCD Reference File.2017 Spray St"N $
323 "PCD Reference [Link]"N $
324 "[Link]"N
325 "PCD Reference [Link]"N $
326 "[Link]"N
327 "[Link]"N
328 "[Link]"N
329 "[Link]"N
330 "PCD Reference [Link]"N
331 "PCD Reference [Link]"N
332 "PCD Reference [Link] Healt"N $
333 "PCD Reference [Link] To Closest"N
334 ;
335 if _ERROR_ then call symputx('_EFIERR_',1); /* set ERROR detection macro variable */
336 run;

NOTE: The infile '/home/u64151009/sasuser.v94/BioStat/Files/ngenirs-mopeia-cross-sectional-survey_2018.csv' is:


Filename=/home/u64151009/sasuser.v94/BioStat/Files/ngenirs-mopeia-cross-sectional-survey_2018.csv,
Owner Name=u64151009,Group Name=oda,
Access Permission=-rw-r--r--,
Last Modified=14May2025:11:34:10,
File Size (bytes)=229211

NOTE: 805 records were read from the infile


'/home/u64151009/sasuser.v94/BioStat/Files/ngenirs-mopeia-cross-sectional-survey_2018.csv'.
The minimum record length was 267.
The maximum record length was 295.
NOTE: The data set WORK.MOPEIA_RAW has 805 observations and 81 variables.
NOTE: DATA statement used (Total process time):
real time 0.01 seconds
user cpu time 0.01 seconds
system cpu time 0.00 seconds
memory 10796.06k
OS Memory 35368.00k
Timestamp 05/15/2025 06:43:24 AM
Step Count 141 Switch Count 2
Page Faults 0
Page Reclaims 172
Page Swaps 0
Voluntary Context Switches 17
Involuntary Context Switches 0
Block Input Operations 0
Block Output Operations 1320

805 rows created in WORK.MOPEIA_RAW from /home/u64151009/sasuser.v94/BioStat/Files/ngenirs-mopeia-cross-sectional-survey_2018.csv.

NOTE: WORK.MOPEIA_RAW data set was successfully created.


NOTE: The data set WORK.MOPEIA_RAW has 805 observations and 81 variables.
NOTE: PROCEDURE IMPORT used (Total process time):
real time 1.37 seconds
user cpu time 1.34 seconds
system cpu time 0.00 seconds
memory 10796.06k
OS Memory 35752.00k

about:blank 6/23
5/15/25, 12:13 PM Program Summary - [Link]
Timestamp 05/15/2025 06:43:24 AM
Step Count 141 Switch Count 17
Page Faults 0
Page Reclaims 1428
Page Swaps 0
Voluntary Context Switches 136
Involuntary Context Switches 4
Block Input Operations 0
Block Output Operations 1400

337
338
339 /* STEP 2: CLEANING & VARIABLE CREATION */
340 data mopeia;
341 set mopeia_raw;
342
343 /* EXCLUSIONS: Remove missing RDT or spray status */
344 if missing('[Link]'n) or missing('PCD Reference File.2017 Spray St'n) then delete;
345
346 /* OUTCOME: Malaria (RDT positive) */
347 MalariaPos = '[Link]'n;
348
349 /* EXPOSURE: IRS in 2017 */
350 if 'PCD Reference File.2017 Spray St'n = 'Spray' then IRSsprayed = 1;
351 else if 'PCD Reference File.2017 Spray St'n = 'No_Spray' then IRSsprayed = 0;
352 else IRSsprayed = .;
353
354 /* BEDNET OWNERSHIP */
355 if 'MAIN_MALARIA_CONTROL_HAS_NETS'n = 1 then BedNet = 1;
356 else if 'MAIN_MALARIA_CONTROL_HAS_NETS'n = 2 then BedNet = 0;
357 else BedNet = .;
358
359 /* AGE GROUP (already categorized) */
360 AgeGroup = AGE_GROUP;
361
362 /* GENDER */
363 Gender = DEMOGRAFIC_GENDER;
364
365 /* LITERACY */
366 ReadAble = '[Link]'n;
367 WriteAble = '[Link]'n;
368 if ReadAble = 1 and WriteAble = 1 then Literate = 1;
369 else if ReadAble = 2 or WriteAble = 2 then Literate = 0;
370 else Literate = 0;
371
372 /* HOUSING VARIABLES */
373 MetalRoof = '[Link]'n;
374 DurableWall = ('[Link] bricks'n = 1 or '[Link] block'n = 1);
375
376 /* LABELS FOR CLARITY */
377 label MalariaPos = "Malaria RDT Positive"
378 IRSsprayed = "IRS Sprayed (2017)"
379 BedNet = "Household Owns Bed Net"
380 AgeGroup = "Age Group"
381 Gender = "Gender"
382 Literate = "Can Read and Write"
383 MetalRoof = "Zinc Roof"
384 DurableWall = "Durable Walls";
385
386 keep MalariaPos IRSsprayed BedNet AgeGroup Gender Literate MetalRoof DurableWall;
387 run;

NOTE: There were 805 observations read from the data set WORK.MOPEIA_RAW.
NOTE: The data set [Link] has 786 observations and 8 variables.
NOTE: DATA statement used (Total process time):
real time 0.00 seconds
user cpu time 0.00 seconds
system cpu time 0.00 seconds
memory 1875.93k
OS Memory 30512.00k
Timestamp 05/15/2025 06:43:24 AM
Step Count 142 Switch Count 2
Page Faults 0
Page Reclaims 203
Page Swaps 0
Voluntary Context Switches 10
Involuntary Context Switches 0
Block Input Operations 0
Block Output Operations 272

388
389
390 /* STEP 3: EXPORT CLEAN DATASET FOR ANALYSIS */
391 libname save '/home/u64151009/sasuser.v94/ZWhite/Files';
NOTE: Libref SAVE was successfully assigned as follows:
Engine: V9

about:blank 7/23
5/15/25, 12:13 PM Program Summary - [Link]
Physical Name: /home/u64151009/sasuser.v94/ZWhite/Files
392
393 data save.clean_mopeia;
394 set mopeia;
395
396
397 /* STEP 4: DESCRIPTIVE ANALYSIS */
398
399 /* 4a. Distributions */

NOTE: There were 786 observations read from the data set [Link].
NOTE: The data set SAVE.CLEAN_MOPEIA has 786 observations and 8 variables.
NOTE: DATA statement used (Total process time):
real time 0.01 seconds
user cpu time 0.01 seconds
system cpu time 0.00 seconds
memory 1059.06k
OS Memory 30128.00k
Timestamp 05/15/2025 06:43:24 AM
Step Count 143 Switch Count 1
Page Faults 0
Page Reclaims 107
Page Swaps 0
Voluntary Context Switches 34
Involuntary Context Switches 0
Block Input Operations 0
Block Output Operations 264

400 proc freq data=mopeia;


401 tables MalariaPos IRSsprayed BedNet AgeGroup Gender Literate MetalRoof DurableWall / missing;
402 title "Overall Distribution of Analysis Variables";
403 run;

NOTE: There were 786 observations read from the data set [Link].
NOTE: PROCEDURE FREQ used (Total process time):
real time 0.04 seconds
user cpu time 0.04 seconds
system cpu time 0.01 seconds
memory 1962.03k
OS Memory 30128.00k
Timestamp 05/15/2025 06:43:24 AM
Step Count 144 Switch Count 3
Page Faults 0
Page Reclaims 119
Page Swaps 0
Voluntary Context Switches 22
Involuntary Context Switches 2
Block Input Operations 0
Block Output Operations 288

404
405 /* 4b. Table 1: Stratified by IRSsprayed */
406 proc freq data=mopeia;
407 tables IRSsprayed*(MalariaPos BedNet AgeGroup Gender Literate MetalRoof DurableWall)
408 / chisq norow nocol nopercent;
409 title "Table 1: Baseline Characteristics by IRSsprayed";
410 run;

NOTE: There were 786 observations read from the data set [Link].
NOTE: PROCEDURE FREQ used (Total process time):
real time 0.13 seconds
user cpu time 0.13 seconds
system cpu time 0.00 seconds
memory 1694.56k
OS Memory 30388.00k
Timestamp 05/15/2025 06:43:24 AM
Step Count 145 Switch Count 5
Page Faults 0
Page Reclaims 185
Page Swaps 0
Voluntary Context Switches 39
Involuntary Context Switches 2
Block Input Operations 0
Block Output Operations 592

411
412
413 /* STEP 5: BIVARIATE ASSOCIATIONS */
414 proc freq data=mopeia;
415 tables BedNet*MalariaPos / chisq relrisk;
416 title "Malaria and BedNet Ownership";
417 run;

NOTE: There were 786 observations read from the data set [Link].
NOTE: PROCEDURE FREQ used (Total process time):

about:blank 8/23
5/15/25, 12:13 PM Program Summary - [Link]
real time 0.03 seconds
user cpu time 0.03 seconds
system cpu time 0.00 seconds
memory 1187.31k
OS Memory 30388.00k
Timestamp 05/15/2025 06:43:24 AM
Step Count 146 Switch Count 5
Page Faults 0
Page Reclaims 185
Page Swaps 0
Voluntary Context Switches 38
Involuntary Context Switches 2
Block Input Operations 0
Block Output Operations 544

418
419 proc freq data=mopeia;
420 tables IRSsprayed*MalariaPos / chisq relrisk;
421 title "Malaria and IRS Spraying";
422 run;

NOTE: There were 786 observations read from the data set [Link].
NOTE: PROCEDURE FREQ used (Total process time):
real time 0.03 seconds
user cpu time 0.03 seconds
system cpu time 0.01 seconds
memory 1184.65k
OS Memory 30388.00k
Timestamp 05/15/2025 06:43:24 AM
Step Count 147 Switch Count 5
Page Faults 0
Page Reclaims 185
Page Swaps 0
Voluntary Context Switches 40
Involuntary Context Switches 3
Block Input Operations 0
Block Output Operations 544

423
424 proc freq data=mopeia;
425 tables AgeGroup*MalariaPos Gender*MalariaPos Literate*MalariaPos
426 MetalRoof*MalariaPos DurableWall*MalariaPos / chisq;
427 title "Other Variables and Malaria Status";
428 run;

NOTE: There were 786 observations read from the data set [Link].
NOTE: PROCEDURE FREQ used (Total process time):
real time 0.11 seconds
user cpu time 0.11 seconds
system cpu time 0.00 seconds
memory 1340.93k
OS Memory 30388.00k
Timestamp 05/15/2025 06:43:24 AM
Step Count 148 Switch Count 5
Page Faults 0
Page Reclaims 185
Page Swaps 0
Voluntary Context Switches 36
Involuntary Context Switches 3
Block Input Operations 0
Block Output Operations 576

429
430
431 /* STEP 6: LOGISTIC REGRESSION */
432
433 /* 6a. Unadjusted: IRSsprayed only */
434 proc logistic data=mopeia descending plots(only)=roc;
435 class IRSsprayed (ref='0') / param=ref;
436 model MalariaPos = IRSsprayed / clodds=wald lackfit outroc=roc1;
437 title "Unadjusted Model: IRSsprayed → Malaria";
438 run;

NOTE: PROC LOGISTIC is modeling the probability that MalariaPos=1.


NOTE: Convergence criterion (GCONV=1E-8) satisfied.
NOTE: The Hosmer-Lemeshow test is not computed because the constructed number of groups (2) must be greater than the
degree-of-freedom correction (2).
NOTE: There were 786 observations read from the data set [Link].
NOTE: The data set WORK.ROC1 has 2 observations and 7 variables.
NOTE: PROCEDURE LOGISTIC used (Total process time):
real time 0.17 seconds
user cpu time 0.12 seconds
system cpu time 0.00 seconds
memory 11714.71k
OS Memory 38980.00k

about:blank 9/23
5/15/25, 12:13 PM Program Summary - [Link]
Timestamp 05/15/2025 06:43:24 AM
Step Count 149 Switch Count 3
Page Faults 0
Page Reclaims 2215
Page Swaps 0
Voluntary Context Switches 229
Involuntary Context Switches 3
Block Input Operations 0
Block Output Operations 1032

439
440 /* 6b. Adjusted model */
441 proc logistic data=mopeia descending plots(only)=roc;
442 class IRSsprayed (ref='0')
443 BedNet (ref='0')
444 AgeGroup (ref='1')
445 Gender (ref='1')
446 Literate (ref='0')
447 MetalRoof (ref='0')
448 DurableWall (ref='0')
449 / param=ref;
450 model MalariaPos = IRSsprayed BedNet AgeGroup Gender Literate MetalRoof DurableWall
451 / clodds=wald lackfit outroc=roc2;
452 title "Adjusted Logistic Model: IRSsprayed & Covariates";
453 run;

NOTE: PROC LOGISTIC is modeling the probability that MalariaPos=1.


NOTE: Convergence criterion (GCONV=1E-8) satisfied.
NOTE: There were 786 observations read from the data set [Link].
NOTE: The data set WORK.ROC2 has 105 observations and 7 variables.
NOTE: PROCEDURE LOGISTIC used (Total process time):
real time 0.17 seconds
user cpu time 0.12 seconds
system cpu time 0.00 seconds
memory 5609.62k
OS Memory 38724.00k
Timestamp 05/15/2025 06:43:24 AM
Step Count 150 Switch Count 3
Page Faults 0
Page Reclaims 738
Page Swaps 0
Voluntary Context Switches 531
Involuntary Context Switches 5
Block Input Operations 0
Block Output Operations 736

454
455
456 /* STEP 7: INTERACTIONS & MODEL DIAGNOSTICS */
457
458 /* IRSsprayed × BedNet */
459 proc logistic data=mopeia descending;
460 class IRSsprayed (ref='0')
461 BedNet (ref='0')
462 AgeGroup (ref='1')
463 Gender (ref='1')
464 Literate (ref='0')
465 MetalRoof (ref='0')
466 DurableWall (ref='0')
467 / param=ref;
468 model MalariaPos = IRSsprayed|BedNet AgeGroup Gender Literate MetalRoof DurableWall
469 / clodds=wald lackfit;
470 title "Interaction: IRSsprayed × BedNet";
471 run;

NOTE: PROC LOGISTIC is modeling the probability that MalariaPos=1.


NOTE: Convergence criterion (GCONV=1E-8) satisfied.
NOTE: Under full-rank parameterizations, Type 3 effect tests are replaced by joint tests. The joint test for an effect is a test
that all the parameters associated with that effect are zero. Such joint tests might not be equivalent to Type 3 effect tests
under GLM parameterization.
NOTE: The CLODDS= option does not compute odds ratios for effects involved in interactions or nestings.
Specify the ODDSRATIO statement to compute odds ratios for these effects.
NOTE: There were 786 observations read from the data set [Link].
NOTE: PROCEDURE LOGISTIC used (Total process time):
real time 0.18 seconds
user cpu time 0.12 seconds
system cpu time 0.01 seconds
memory 5966.21k
OS Memory 39488.00k
Timestamp 05/15/2025 06:43:25 AM
Step Count 151 Switch Count 1
Page Faults 0
Page Reclaims 732
Page Swaps 0
Voluntary Context Switches 223
Involuntary Context Switches 6

about:blank 10/23
5/15/25, 12:13 PM Program Summary - [Link]
Block Input Operations 0
Block Output Operations 472

472
473 /* Goodness-of-fit and ROC */
474 ods graphics on;
475 ods select FitStatistics Association LackFitChiSq;
476 proc logistic data=mopeia descending plots(only)=roc(id=prob);
477 class IRSsprayed (ref='0')
478 BedNet (ref='0')
479 AgeGroup (ref='1')
480 Gender (ref='1')
481 Literate (ref='0')
482 MetalRoof (ref='0')
483 DurableWall (ref='0')
484 / param=ref;
485 model MalariaPos = IRSsprayed BedNet AgeGroup Gender Literate MetalRoof DurableWall
486 / clodds=wald lackfit rsquare;
487 title "Model Diagnostics: Fit, ROC, HL Test";
488 run;

NOTE: PROC LOGISTIC is modeling the probability that MalariaPos=1.


NOTE: Convergence criterion (GCONV=1E-8) satisfied.
NOTE: There were 786 observations read from the data set [Link].
NOTE: PROCEDURE LOGISTIC used (Total process time):
real time 0.03 seconds
user cpu time 0.03 seconds
system cpu time 0.00 seconds
memory 2790.46k
OS Memory 37824.00k
Timestamp 05/15/2025 06:43:25 AM
Step Count 152 Switch Count 1
Page Faults 0
Page Reclaims 233
Page Swaps 0
Voluntary Context Switches 11
Involuntary Context Switches 2
Block Input Operations 0
Block Output Operations 296

489
490
491 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;
501

Results: [Link]

Overall Distribution of Analysis Variables


The FREQ Procedure

Malaria RDT Positive


Cumulative Cumulative
MalariaPos Frequency Percent Frequency Percent
0 487 61.96 487 61.96
1 299 38.04 786 100.00

IRS Sprayed (2017)


Cumulative Cumulative
IRSsprayed Frequency Percent Frequency Percent
0 387 49.24 387 49.24
1 399 50.76 786 100.00

Household Owns Bed Net


Cumulative Cumulative
BedNet Frequency Percent Frequency Percent
0 42 5.34 42 5.34
1 744 94.66 786 100.00

Age Group
Cumulative Cumulative
AgeGroup Frequency Percent Frequency Percent
1 384 48.85 384 48.85
2 59 7.51 443 56.36
3 343 43.64 786 100.00

Gender
Cumulative Cumulative
Gender Frequency Percent Frequency Percent
1 404 51.40 404 51.40

about:blank 11/23
5/15/25, 12:13 PM Program Summary - [Link]
Gender
Cumulative Cumulative
Gender Frequency Percent Frequency Percent
2 382 48.60 786 100.00

Can Read and Write


Cumulative Cumulative
Literate Frequency Percent Frequency Percent
0 468 59.54 468 59.54
1 318 40.46 786 100.00

Zinc Roof
Cumulative Cumulative
MetalRoof Frequency Percent Frequency Percent
0 634 80.66 634 80.66
1 152 19.34 786 100.00

Durable Walls
Cumulative Cumulative
DurableWall Frequency Percent Frequency Percent
0 417 53.05 417 53.05
1 369 46.95 786 100.00

Table 1: Baseline Characteristics by IRSsprayed


The FREQ Procedure

Frequency Table of IRSsprayed by MalariaPos


MalariaPos(Malaria RDT Positive)
IRSsprayed(IRS Sprayed (2017)) 0 1 Total
0 225 162 387
1 262 137 399
Total 487 299 786

Statistics for Table of IRSsprayed by MalariaPos

Statistic DF Value Prob


Chi-Square 1 4.7193 0.0298
Likelihood Ratio Chi-Square 1 4.7233 0.0298
Continuity Adj. Chi-Square 1 4.4054 0.0358
Mantel-Haenszel Chi-Square 1 4.7133 0.0299
Phi Coefficient -0.0775
Contingency Coefficient 0.0773
Cramer's V -0.0775

Fisher's Exact Test


Cell (1,1) Frequency (F) 225
Left-sided Pr <= F 0.0179
Right-sided Pr >= F 0.9877

Table Probability (P) 0.0056


Two-sided Pr <= P 0.0331

Sample Size = 786

Frequency Table of IRSsprayed by BedNet


BedNet(Household Owns Bed Net)
IRSsprayed(IRS Sprayed (2017)) 0 1 Total
0 23 364 387
1 19 380 399
Total 42 744 786

Statistics for Table of IRSsprayed by BedNet

Statistic DF Value Prob


Chi-Square 1 0.5420 0.4616
Likelihood Ratio Chi-Square 1 0.5424 0.4614
Continuity Adj. Chi-Square 1 0.3336 0.5636
Mantel-Haenszel Chi-Square 1 0.5413 0.4619
Phi Coefficient 0.0263
Contingency Coefficient 0.0262
Cramer's V 0.0263

Fisher's Exact Test


Cell (1,1) Frequency (F) 23
Left-sided Pr <= F 0.8145
Right-sided Pr >= F 0.2819

about:blank 12/23
5/15/25, 12:13 PM Program Summary - [Link]
Fisher's Exact Test
Table Probability (P) 0.0964
Two-sided Pr <= P 0.5269

Sample Size = 786

Frequency Table of IRSsprayed by AgeGroup


AgeGroup(Age Group)
IRSsprayed(IRS Sprayed (2017)) 1 2 3 Total
0 181 27 179 387
1 203 32 164 399
Total 384 59 343 786

Statistics for Table of IRSsprayed by AgeGroup

Statistic DF Value Prob


Chi-Square 2 2.1574 0.3400
Likelihood Ratio Chi-Square 2 2.1583 0.3399
Mantel-Haenszel Chi-Square 1 1.8234 0.1769
Phi Coefficient 0.0524
Contingency Coefficient 0.0523
Cramer's V 0.0524

Sample Size = 786

Frequency Table of IRSsprayed by Gender


Gender(Gender)
IRSsprayed(IRS Sprayed (2017)) 1 2 Total
0 190 197 387
1 214 185 399
Total 404 382 786

Statistics for Table of IRSsprayed by Gender

Statistic DF Value Prob


Chi-Square 1 1.6199 0.2031
Likelihood Ratio Chi-Square 1 1.6204 0.2030
Continuity Adj. Chi-Square 1 1.4433 0.2296
Mantel-Haenszel Chi-Square 1 1.6178 0.2034
Phi Coefficient -0.0454
Contingency Coefficient 0.0454
Cramer's V -0.0454

Fisher's Exact Test


Cell (1,1) Frequency (F) 190
Left-sided Pr <= F 0.1148
Right-sided Pr >= F 0.9106

Table Probability (P) 0.0254


Two-sided Pr <= P 0.2250

Sample Size = 786

Frequency Table of IRSsprayed by Literate


Literate(Can Read and Write)
IRSsprayed(IRS Sprayed (2017)) 0 1 Total
0 236 151 387
1 232 167 399
Total 468 318 786

Statistics for Table of IRSsprayed by Literate

Statistic DF Value Prob


Chi-Square 1 0.6562 0.4179
Likelihood Ratio Chi-Square 1 0.6563 0.4179
Continuity Adj. Chi-Square 1 0.5437 0.4609
Mantel-Haenszel Chi-Square 1 0.6553 0.4182
Phi Coefficient 0.0289
Contingency Coefficient 0.0289
Cramer's V 0.0289

Fisher's Exact Test


Cell (1,1) Frequency (F) 236
Left-sided Pr <= F 0.8113
Right-sided Pr >= F 0.2305

Table Probability (P) 0.0418

about:blank 13/23
5/15/25, 12:13 PM Program Summary - [Link]
Fisher's Exact Test
Two-sided Pr <= P 0.4249

Sample Size = 786

Frequency Table of IRSsprayed by MetalRoof


MetalRoof(Zinc Roof)
IRSsprayed(IRS Sprayed (2017)) 0 1 Total
0 309 78 387
1 325 74 399
Total 634 152 786

Statistics for Table of IRSsprayed by MetalRoof

Statistic DF Value Prob


Chi-Square 1 0.3259 0.5681
Likelihood Ratio Chi-Square 1 0.3259 0.5681
Continuity Adj. Chi-Square 1 0.2309 0.6308
Mantel-Haenszel Chi-Square 1 0.3255 0.5683
Phi Coefficient -0.0204
Contingency Coefficient 0.0204
Cramer's V -0.0204

Fisher's Exact Test


Cell (1,1) Frequency (F) 309
Left-sided Pr <= F 0.3154
Right-sided Pr >= F 0.7458

Table Probability (P) 0.0612


Two-sided Pr <= P 0.5886

Sample Size = 786

Frequency Table of IRSsprayed by DurableWall


DurableWall(Durable Walls)
IRSsprayed(IRS Sprayed (2017)) 0 1 Total
0 188 199 387
1 229 170 399
Total 417 369 786

Statistics for Table of IRSsprayed by DurableWall

Statistic DF Value Prob


Chi-Square 1 6.1285 0.0133
Likelihood Ratio Chi-Square 1 6.1360 0.0132
Continuity Adj. Chi-Square 1 5.7797 0.0162
Mantel-Haenszel Chi-Square 1 6.1207 0.0134
Phi Coefficient -0.0883
Contingency Coefficient 0.0880
Cramer's V -0.0883

Fisher's Exact Test


Cell (1,1) Frequency (F) 188
Left-sided Pr <= F 0.0081
Right-sided Pr >= F 0.9946

Table Probability (P) 0.0027


Two-sided Pr <= P 0.0151

Sample Size = 786

Malaria and BedNet Ownership


The FREQ Procedure

Frequency Table of BedNet by MalariaPos


Percent
Row Pct MalariaPos(Malaria RDT Positive)
Col Pct BedNet(Household Owns Bed Net) 0 1 Total
0 24 18 42
3.05 2.29 5.34
57.14 42.86
4.93 6.02
1 463 281 744
58.91 35.75 94.66
62.23 37.77
95.07 93.98
Total 487 299 786
61.96 38.04 100.00

Statistics for Table of BedNet by MalariaPos

about:blank 14/23
5/15/25, 12:13 PM Program Summary - [Link]
Statistic DF Value Prob
Chi-Square 1 0.4367 0.5087
Likelihood Ratio Chi-Square 1 0.4308 0.5116
Continuity Adj. Chi-Square 1 0.2475 0.6188
Mantel-Haenszel Chi-Square 1 0.4362 0.5090
Phi Coefficient -0.0236
Contingency Coefficient 0.0236
Cramer's V -0.0236

Fisher's Exact Test


Cell (1,1) Frequency (F) 24
Left-sided Pr <= F 0.3065
Right-sided Pr >= F 0.7960

Table Probability (P) 0.1026


Two-sided Pr <= P 0.5171

Odds Ratio and Relative Risks


Statistic Value 95% Confidence Limits
Odds Ratio 0.8092 0.4315 1.5176
Relative Risk (Column 1) 0.9182 0.7025 1.2002
Relative Risk (Column 2) 1.1347 0.7907 1.6284

Sample Size = 786

Malaria and IRS Spraying


The FREQ Procedure

Frequency Table of IRSsprayed by MalariaPos


Percent
Row Pct MalariaPos(Malaria RDT Positive)
Col Pct IRSsprayed(IRS Sprayed (2017)) 0 1 Total
0 225 162 387
28.63 20.61 49.24
58.14 41.86
46.20 54.18
1 262 137 399
33.33 17.43 50.76
65.66 34.34
53.80 45.82
Total 487 299 786
61.96 38.04 100.00

Statistics for Table of IRSsprayed by MalariaPos

Statistic DF Value Prob


Chi-Square 1 4.7193 0.0298
Likelihood Ratio Chi-Square 1 4.7233 0.0298
Continuity Adj. Chi-Square 1 4.4054 0.0358
Mantel-Haenszel Chi-Square 1 4.7133 0.0299
Phi Coefficient -0.0775
Contingency Coefficient 0.0773
Cramer's V -0.0775

Fisher's Exact Test


Cell (1,1) Frequency (F) 225
Left-sided Pr <= F 0.0179
Right-sided Pr >= F 0.9877

Table Probability (P) 0.0056


Two-sided Pr <= P 0.0331

Odds Ratio and Relative Risks


Statistic Value 95% Confidence Limits
Odds Ratio 0.7263 0.5440 0.9696
Relative Risk (Column 1) 0.8854 0.7929 0.9887
Relative Risk (Column 2) 1.2191 1.0189 1.4588

Sample Size = 786

Other Variables and Malaria Status


The FREQ Procedure

Frequency Table of AgeGroup by MalariaPos


Percent
Row Pct MalariaPos(Malaria RDT Positive)
Col Pct AgeGroup(Age Group) 0 1 Total
1 176 208 384
22.39 26.46 48.85
45.83 54.17
36.14 69.57
2 38 21 59
4.83 2.67 7.51
64.41 35.59

about:blank 15/23
5/15/25, 12:13 PM Program Summary - [Link]
Table of AgeGroup by MalariaPos
MalariaPos(Malaria RDT Positive)
AgeGroup(Age Group) 0 1 Total
7.80 7.02
3 273 70 343
34.73 8.91 43.64
79.59 20.41
56.06 23.41
Total 487 299 786
61.96 38.04 100.00

Statistics for Table of AgeGroup by MalariaPos

Statistic DF Value Prob


Chi-Square 2 87.7617 <.0001
Likelihood Ratio Chi-Square 2 90.6092 <.0001
Mantel-Haenszel Chi-Square 1 87.5837 <.0001
Phi Coefficient 0.3341
Contingency Coefficient 0.3169
Cramer's V 0.3341

Sample Size = 786

Frequency Table of Gender by MalariaPos


Percent
Row Pct MalariaPos(Malaria RDT Positive)
Col Pct Gender(Gender) 0 1 Total
1 244 160 404
31.04 20.36 51.40
60.40 39.60
50.10 53.51
2 243 139 382
30.92 17.68 48.60
63.61 36.39
49.90 46.49
Total 487 299 786
61.96 38.04 100.00

Statistics for Table of Gender by MalariaPos

Statistic DF Value Prob


Chi-Square 1 0.8619 0.3532
Likelihood Ratio Chi-Square 1 0.8623 0.3531
Continuity Adj. Chi-Square 1 0.7308 0.3926
Mantel-Haenszel Chi-Square 1 0.8608 0.3535
Phi Coefficient -0.0331
Contingency Coefficient 0.0331
Cramer's V -0.0331

Fisher's Exact Test


Cell (1,1) Frequency (F) 244
Left-sided Pr <= F 0.1963
Right-sided Pr >= F 0.8418

Table Probability (P) 0.0381


Two-sided Pr <= P 0.3780

Sample Size = 786

Frequency Table of Literate by MalariaPos


Percent
Row Pct MalariaPos(Malaria RDT Positive)
Col Pct Literate(Can Read and Write) 0 1 Total
0 288 180 468
36.64 22.90 59.54
61.54 38.46
59.14 60.20
1 199 119 318
25.32 15.14 40.46
62.58 37.42
40.86 39.80
Total 487 299 786
61.96 38.04 100.00

Statistics for Table of Literate by MalariaPos

Statistic DF Value Prob


Chi-Square 1 0.0869 0.7681
Likelihood Ratio Chi-Square 1 0.0870 0.7681
Continuity Adj. Chi-Square 1 0.0484 0.8259
Mantel-Haenszel Chi-Square 1 0.0868 0.7683
Phi Coefficient -0.0105
Contingency Coefficient 0.0105
Cramer's V -0.0105

Fisher's Exact Test


Cell (1,1) Frequency (F) 288
Left-sided Pr <= F 0.4134

about:blank 16/23
5/15/25, 12:13 PM Program Summary - [Link]
Fisher's Exact Test
Right-sided Pr >= F 0.6438

Table Probability (P) 0.0572


Two-sided Pr <= P 0.8224

Sample Size = 786

Frequency Table of MetalRoof by MalariaPos


Percent
Row Pct MalariaPos(Malaria RDT Positive)
Col Pct MetalRoof(Zinc Roof) 0 1 Total
0 379 255 634
48.22 32.44 80.66
59.78 40.22
77.82 85.28
1 108 44 152
13.74 5.60 19.34
71.05 28.95
22.18 14.72
Total 487 299 786
61.96 38.04 100.00

Statistics for Table of MetalRoof by MalariaPos

Statistic DF Value Prob


Chi-Square 1 6.6110 0.0101
Likelihood Ratio Chi-Square 1 6.8097 0.0091
Continuity Adj. Chi-Square 1 6.1414 0.0132
Mantel-Haenszel Chi-Square 1 6.6026 0.0102
Phi Coefficient -0.0917
Contingency Coefficient 0.0913
Cramer's V -0.0917

Fisher's Exact Test


Cell (1,1) Frequency (F) 379
Left-sided Pr <= F 0.0061
Right-sided Pr >= F 0.9965

Table Probability (P) 0.0026


Two-sided Pr <= P 0.0119

Sample Size = 786

Frequency Table of DurableWall by MalariaPos


Percent
Row Pct MalariaPos(Malaria RDT Positive)
Col Pct DurableWall(Durable Walls) 0 1 Total
0 258 159 417
32.82 20.23 53.05
61.87 38.13
52.98 53.18
1 229 140 369
29.13 17.81 46.95
62.06 37.94
47.02 46.82
Total 487 299 786
61.96 38.04 100.00

Statistics for Table of DurableWall by MalariaPos

Statistic DF Value Prob


Chi-Square 1 0.0030 0.9565
Likelihood Ratio Chi-Square 1 0.0030 0.9565
Continuity Adj. Chi-Square 1 0.0000 1.0000
Mantel-Haenszel Chi-Square 1 0.0030 0.9566
Phi Coefficient -0.0019
Contingency Coefficient 0.0019
Cramer's V -0.0019

Fisher's Exact Test


Cell (1,1) Frequency (F) 258
Left-sided Pr <= F 0.5078
Right-sided Pr >= F 0.5508

Table Probability (P) 0.0586


Two-sided Pr <= P 1.0000

Sample Size = 786

Unadjusted Model: IRSsprayed → Malaria


The LOGISTIC Procedure

Model Information
Data Set [Link]
Response Variable MalariaPos Malaria RDT Positive

about:blank 17/23
5/15/25, 12:13 PM Program Summary - [Link]
Model Information
Number of Response Levels 2
Model binary logit
Optimization Technique Fisher's scoring

Number of Observations Read 786


Number of Observations Used 786

Response Profile
Ordered Total
Value MalariaPos Frequency
1 1 299
2 0 487

Probability modeled is MalariaPos=1.

Class Level Information


Class Value Design Variables
IRSsprayed 0 0
1 1

Model Convergence Status


Convergence criterion (GCONV=1E-8) satisfied.

Model Fit Statistics


Criterion Intercept Only Intercept and Covariates
AIC 1046.222 1043.498
SC 1050.889 1052.832
-2 Log L 1044.222 1039.498

Testing Global Null Hypothesis: BETA=0


Test Chi-Square DF Pr > ChiSq
Likelihood Ratio 4.7233 1 0.0298
Score 4.7193 1 0.0298
Wald 4.7075 1 0.0300

Type 3 Analysis of Effects


Wald
Effect DF Chi-Square Pr > ChiSq
IRSsprayed 1 4.7075 0.0300

Analysis of Maximum Likelihood Estimates


Standard Wald
Parameter DF Estimate Error Chi-Square Pr > ChiSq
Intercept 1 -0.3285 0.1030 10.1640 0.0014
IRSsprayed 1 1 -0.3199 0.1474 4.7075 0.0300

Association of Predicted Probabilities and Observed Responses


Percent Concordant 29.1 Somers' D 0.080
Percent Discordant 21.2 Gamma 0.159
Percent Tied 49.7 Tau-a 0.038
Pairs 145613 c 0.540

Odds Ratio Estimates and Wald Confidence Intervals


Effect Unit Estimate 95% Confidence Limits
IRSsprayed 1 vs 0 1.0000 0.726 0.544 0.970

about:blank 18/23
5/15/25, 12:13 PM Program Summary - [Link]

Adjusted Logistic Model: IRSsprayed & Covariates


The LOGISTIC Procedure

Model Information
Data Set [Link]
Response Variable MalariaPos Malaria RDT Positive
Number of Response Levels 2
Model binary logit
Optimization Technique Fisher's scoring

Number of Observations Read 786


Number of Observations Used 786

Response Profile
Ordered Total
Value MalariaPos Frequency
1 1 299
2 0 487

Probability modeled is MalariaPos=1.

Class Level Information


Class Value Design Variables
IRSsprayed 0 0
1 1
BedNet 0 0
1 1
AgeGroup 1 0 0
2 1 0
3 0 1
Gender 1 0
2 1
Literate 0 0
1 1
MetalRoof 0 0
1 1
DurableWall 0 0
1 1

Model Convergence Status


Convergence criterion (GCONV=1E-8) satisfied.

Model Fit Statistics


Criterion Intercept Only Intercept and Covariates
AIC 1046.222 953.584
SC 1050.889 995.587

about:blank 19/23
5/15/25, 12:13 PM Program Summary - [Link]
Model Fit Statistics
Criterion Intercept Only Intercept and Covariates
-2 Log L 1044.222 935.584

Testing Global Null Hypothesis: BETA=0


Test Chi-Square DF Pr > ChiSq
Likelihood Ratio 108.6377 8 <.0001
Score 103.5159 8 <.0001
Wald 94.6308 8 <.0001

Type 3 Analysis of Effects


Wald
Effect DF Chi-Square Pr > ChiSq
IRSsprayed 1 7.7782 0.0053
BedNet 1 0.1771 0.6739
AgeGroup 2 85.5661 <.0001
Gender 1 0.0201 0.8873
Literate 1 0.5300 0.4666
MetalRoof 1 7.9114 0.0049
DurableWall 1 0.1119 0.7379

Analysis of Maximum Likelihood Estimates


Standard Wald
Parameter DF Estimate Error Chi-Square Pr > ChiSq
Intercept 1 0.7105 0.3841 3.4217 0.0643
IRSsprayed 1 1 -0.4459 0.1599 7.7782 0.0053
BedNet 1 1 -0.1458 0.3465 0.1771 0.6739
AgeGroup 2 1 -0.7627 0.2949 6.6865 0.0097
AgeGroup 3 1 -1.6031 0.1737 85.2263 <.0001
Gender 2 1 -0.0249 0.1758 0.0201 0.8873
Literate 1 1 -0.1318 0.1810 0.5300 0.4666
MetalRoof 1 1 -0.6290 0.2236 7.9114 0.0049
DurableWall 1 1 0.0576 0.1720 0.1119 0.7379

Association of Predicted Probabilities and Observed Responses


Percent Concordant 70.8 Somers' D 0.436
Percent Discordant 27.2 Gamma 0.444
Percent Tied 2.0 Tau-a 0.206
Pairs 145613 c 0.718

Odds Ratio Estimates and Wald Confidence Intervals


Effect Unit Estimate 95% Confidence Limits
IRSsprayed 1 vs 0 1.0000 0.640 0.468 0.876
BedNet 1 vs 0 1.0000 0.864 0.438 1.705
AgeGroup 2 vs 1 1.0000 0.466 0.262 0.831
AgeGroup 3 vs 1 1.0000 0.201 0.143 0.283
Gender 2 vs 1 1.0000 0.975 0.691 1.377
Literate 1 vs 0 1.0000 0.877 0.615 1.250
MetalRoof 1 vs 0 1.0000 0.533 0.344 0.826
DurableWall 1 vs 0 1.0000 1.059 0.756 1.484

about:blank 20/23
5/15/25, 12:13 PM Program Summary - [Link]

Partition for the Hosmer and Lemeshow Test


MalariaPos = 1 MalariaPos = 0
Group Total Observed Expected Observed Expected
1 84 20 11.90 64 72.10
2 87 14 15.61 73 71.39
3 82 14 18.30 68 63.70
4 83 21 21.62 62 61.38
5 82 25 27.76 57 54.24
6 91 37 43.16 54 47.84
7 76 37 39.27 39 36.73
8 89 55 49.92 34 39.08
9 112 76 71.46 36 40.54

Hosmer and Lemeshow Goodness-of-Fit Test


Chi-Square DF Pr > ChiSq
12.2842 7 0.0916

Interaction: IRSsprayed × BedNet


The LOGISTIC Procedure

Model Information
Data Set [Link]
Response Variable MalariaPos Malaria RDT Positive
Number of Response Levels 2
Model binary logit
Optimization Technique Fisher's scoring

Number of Observations Read 786


Number of Observations Used 786

Response Profile
Ordered Total
Value MalariaPos Frequency
1 1 299
2 0 487

Probability modeled is MalariaPos=1.

Class Level Information


Class Value Design Variables
IRSsprayed 0 0
1 1
BedNet 0 0
1 1
AgeGroup 1 0 0
2 1 0

about:blank 21/23
5/15/25, 12:13 PM Program Summary - [Link]
Class Level Information
Class Value Design Variables
3 0 1
Gender 1 0
2 1
Literate 0 0
1 1
MetalRoof 0 0
1 1
DurableWall 0 0
1 1

Model Convergence Status


Convergence criterion (GCONV=1E-8) satisfied.

Model Fit Statistics


Criterion Intercept Only Intercept and Covariates
AIC 1046.222 951.526
SC 1050.889 998.195
-2 Log L 1044.222 931.526

Testing Global Null Hypothesis: BETA=0


Test Chi-Square DF Pr > ChiSq
Likelihood Ratio 112.6960 9 <.0001
Score 107.1334 9 <.0001
Wald 97.2028 9 <.0001

Joint Tests
Wald
Effect DF Chi-Square Pr > ChiSq
IRSsprayed 1 1.6502 0.1989
BedNet 1 1.1326 0.2872
IRSsprayed*BedNet 1 3.9238 0.0476
AgeGroup 2 84.4534 <.0001
Gender 1 0.0418 0.8379
Literate 1 0.7165 0.3973
MetalRoof 1 7.5786 0.0059
DurableWall 1 0.0660 0.7972

Note: Under full-rank parameterizations, Type 3 effect tests are replaced by joint tests. The joint test for an effect is a test that all the parameters associated with that effect are zero. Such joint tests might not be equivalent to
Type 3 effect tests under GLM parameterization.

Analysis of Maximum Likelihood Estimates


Standard Wald
Parameter DF Estimate Error Chi-Square Pr > ChiSq
Intercept 1 0.0834 0.5116 0.0266 0.8705
IRSsprayed 1 1 0.8951 0.6968 1.6502 0.1989
BedNet 1 1 0.5354 0.5031 1.1326 0.2872
IRSsprayed*BedNet 1 1 1 -1.4180 0.7158 3.9238 0.0476
AgeGroup 2 1 -0.7634 0.2963 6.6363 0.0100
AgeGroup 3 1 -1.5960 0.1740 84.1154 <.0001
Gender 2 1 -0.0361 0.1765 0.0418 0.8379
Literate 1 1 -0.1541 0.1821 0.7165 0.3973
MetalRoof 1 1 -0.6179 0.2244 7.5786 0.0059
DurableWall 1 1 0.0444 0.1727 0.0660 0.7972

Association of Predicted Probabilities and Observed Responses


Percent Concordant 71.0 Somers' D 0.440
Percent Discordant 27.0 Gamma 0.449
Percent Tied 2.0 Tau-a 0.208
Pairs 145613 c 0.720

Odds Ratio Estimates and Wald Confidence Intervals


Effect Unit Estimate 95% Confidence Limits
AgeGroup 2 vs 1 1.0000 0.466 0.261 0.833
AgeGroup 3 vs 1 1.0000 0.203 0.144 0.285
Gender 2 vs 1 1.0000 0.965 0.682 1.363
Literate 1 vs 0 1.0000 0.857 0.600 1.225
MetalRoof 1 vs 0 1.0000 0.539 0.347 0.837
DurableWall 1 vs 0 1.0000 1.045 0.745 1.466

about:blank 22/23
5/15/25, 12:13 PM Program Summary - [Link]

Partition for the Hosmer and Lemeshow Test


MalariaPos = 1 MalariaPos = 0
Group Total Observed Expected Observed Expected
1 79 17 10.83 62 68.17
2 79 14 13.71 65 65.29
3 78 13 15.99 65 62.01
4 83 19 21.85 64 61.15
5 79 24 24.64 55 54.36
6 66 26 28.89 40 37.11
7 81 37 39.98 44 41.02
8 84 41 43.64 43 40.36
9 89 59 54.34 30 34.66
10 68 49 45.12 19 22.88

Hosmer and Lemeshow Goodness-of-Fit Test


Chi-Square DF Pr > ChiSq
8.6184 8 0.3755

Model Diagnostics: Fit, ROC, HL Test


The LOGISTIC Procedure

Model Fit Statistics


Criterion Intercept Only Intercept and Covariates
AIC 1046.222 953.584
SC 1050.889 995.587
-2 Log L 1044.222 935.584

Association of Predicted Probabilities and Observed Responses


Percent Concordant 70.8 Somers' D 0.436
Percent Discordant 27.2 Gamma 0.444
Percent Tied 2.0 Tau-a 0.206
Pairs 145613 c 0.718

Hosmer and Lemeshow Goodness-of-Fit Test


Chi-Square DF Pr > ChiSq
12.2842 7 0.0916

about:blank 23/23

You might also like