3.1) Selection of schools
A comprehensive list of operational primary (n = 784) and elementary (n = 1719) schools across the seven provinces was acquired from the National Department of Education, serving as a sampling frame for the selection of thirty schools along with two backups. The School Malaria Survey (SMS) aimed to assess malaria prevalence among children aged 5–16 years; thus, pairs of elementary (grades 1 and 2) and primary schools (grades 3–8) were chosen to encompass this age range. The corresponding elementary school within the primary school's catchment area was selected. To account for the altitude effect on malaria prevalence, a stratified sampling approach was employed. Catchment areas of the schools were mapped using ArcGIS 10.3, and these delineated areas were grouped by malaria incidence categories. Hence, three incidence categories were defined: low (< 1), moderate (1–10), and high (> 10 per 1000). A weighted random sampling, utilizing a random number table in STATA, was then performed within each incidence category to select the schools.
3.2) Selection of students, and households for RCD
Within each pair of schools, 200 students were randomly sampled from a list of all students. In case a sampling frame contained less than 200 individuals, all available students were selected. Then, sampled students were tested for malaria as described below. Those who tested positive using mRDTs were chosen for a subsequent RCD household survey. Additional households were randomly selected from the group of students who tested negative, bringing the total to ten households per school. Every present household member was then tested for malaria.
3.3) Data collection forms
In the SMS, three data collection forms were utilized: a school summary form, a malaria test register, and a student questionnaire. The first two forms were filled out once per school by the team leader and research nurse, while the third form, an electronic structured questionnaire was administered to the selected students. Additionally, during the RCD survey, three electronic questionnaires were employed, aimed at the household head and other members. These questionnaires, adapted from the Malaria Indicator Survey set [16], included a household questionnaire (one per household), a treatment-seeking questionnaire (one per household member with a recent febrile illness), and a prevalence form (one per household member tested for malaria), all programmed in Open Data Kit (ODK) at PNGIMR and administered via tablet computers.
3.4) Survey implementation procedures
The SMS and RCD were conducted from July to November 2019 by three PNGIMR field teams operating concurrently at different locations. Each team comprised two nursing officers, a scientific officer, and one or two research assistants, all of whom underwent comprehensive training on the project background, survey protocol, and sampling techniques.
Upon arriving at a selected school, an informational meeting with parents and the school board was arranged to explain the study and its methods. RCD interviews were held with household heads, while all household members provided blood samples and those with a recent fever episode were interviewed on treatment seeking. The survey spanned 6–7 days per school (including RCD households). Handheld GPS devices (Garmin) were used to log the geo-coordinates and elevations of the surveyed schools and households.
3.5) Testing of malaria and anaemia
Trained nurses collected blood samples by finger-prick from study participants (i.e., selected students and household members) aged six months or older. One thick and one thin smear were prepared on the same glass slide for malaria diagnosis by light microscopy. In addition, an mRDT (CareStart Malaria HRP2/pLDH (Pf/pan) Combo Test, Access Bio) was performed, and a microcuvette sample was prepared to measure Hb levels using a handheld HemoCue Hb 201 + analyser (HemoCue, Angelholm, Sweden).
During the surveys, positive mRDT results were classified as P. falciparum (control line with HRP2 test-line only), non-P. falciparum (control line with pLDH test-line only), or P. falciparum or mixed infection (control line with two test lines for HRP2 and pLDH) [17]. If an invalid RDT result was encountered (i.e., no control line, irrespective of the two test lines), the test was repeated.
The axillary temperature of participants was measured with an electronic thermometer. In addition, children aged (2–9) years were examined for splenomegaly and graded according to the Hackett grading system [18].
The nursing officers provided treatment for mRDT positive individuals following the national treatment protocol. They also offered treatment for other minor ailments encountered or gave referral advice.
Light microscopy was done at the PNGIMR in Madang following established procedures for research studies [5, 19]. Two qualified microscopists had examined blood slides independently and were blinded to results of each other. Slides were considered positive for malaria if judged positive by at least two microscopists. A third senior microscopist examined slides with discrepant results.
3.6) Passive case detection in the nearby health facilities
Data on malaria cases reported in the health facilities (HFs) between 2017 and 2019 was sourced from the electronic Health Information System (eNHIS). Since the obtained dataset does not provide information on malaria cases at the aid post level, our analysis was limited to the nearest health centres and urban clinics. Nearest HFs to the surveyed schools and households were identified using shapefiles of the surveyed sites and health facilities in the Highlands, and a travel friction raster, as described in [3]. The gdistance package in R was employed to calculate cost distance of travel time (minutes) from surveyed sites to HFs.
Only cases confirmed by mRDTs or microscopy at the HFs were included in the analyses. Malaria incidence within the catchment areas of surveyed schools and households was computed for the general population, as well as for the children under five years and children at school age (5–14 years). The population projections for the three years, based on 2011 Census and growth rates reported by the National Statistics Office of PNG were utilized in incidence estimations. Given that the annual incidence rates exhibit over-dispersion, negative binomial models were employed to estimate the average incidence rates, along with 95% confidence interval.
3.7) Data management and analysis
Data were collected electronically using the ODK Collect application installed on tablet computers. Completed and verified collected forms were uploaded directly to the project server at the Swiss TPH using the local mobile phone network (Digicel). ODK Briefcase v1.4.9 was used to download the data and export it for analysis in STATA/IC 14.2 (StataCorp LLC, College Station, TX, USA).
Weighted statistics with 95% confidence limits were calculated using the survey design command set in Stata (svy). Schools and households were established as the primary sampling units, while the seven provinces were treated as strata. We calculated separate sampling weights for schools and household surveys based on the inverse of an observation's probability of selection. The overall selection probability was determined from the eight combination categories of altitude and incidence within the primary schools' sampling frame. For each school, the selection probability of a RCD household was computed among tested students. Since all members of the sampled family were eligible, individual-level weights were set equal to the household weights to which an individual belonged.
Malaria prevalence was calculated separately for tested students and household members as a proportion of positive results detected by light microscopy. WHO adjustments of Hb measurements, considering sex and altitude [20] were applied prior to determining the anaemia status of each tested individual and calculation of prevalence among the study groups. Prevalence of malaria and anaemia were age-standardised using the overall age distribution of the PNG population reported from the 2011 Census by the National Statistics Office (NSO) of PNG.
Contingency tables and odds ratios with 95% confidence intervals were generated separately for the students and household groups to explore the relative risk of malaria infection by travel history, ownership, and use of bed nets. Where appropriate, Microsoft Excel graphs of percentage and cumulative frequencies were used to compare the participants' groups. Further, to account for the relationship of altitude and malaria transmission, four altitudinal categories were used in the analyses in line with previous publications [3]: <1200m, 1200-1600m, 1600–2000, and above 2000m.
3.7.4) LLIN ownership and use
Bed net ownership and use were calculated following standard procedures [21]. The population percentage with access to an LLIN was evaluated at the household level by creating an intermediate variable named the potential users of LLINs [22]. Hence, the number of LLIN sleeping spaces, assuming two per LLIN, was calculated as potential users, which was then divided by the number of people sleeping in the household the previous night. However, the access proportion was adjusted if the potential use exceeded the number of people in the house the last night, i.e., to converge the access proportion to one. Net use among the people with access to LLINs was calculated by dividing the number of people using an LLIN by the total population with access (multiplying the weighted proportion with access by the total population).