Weighted rank regression for clustered data analysis

We consider ranked-based regression models for clustered data analysis. A weighted Wilcoxon rank method is proposed to take account of within-cluster correlations and varying cluster sizes. The asymptotic normality of the resulting estimators is established. A method to estimate covariance of the estimators is also given, which can bypass estimation of the density function. Simulation studies are carried out to compare different estimators for a number of scenarios on the correlation structure, presence/absence of outliers and different correlation values. The proposed methods appear to perform well, in particular, the one incorporating the correlation in the weighting achieves the highest efficiency and robustness against misspecification of correlation structure and outliers. A real example is provided for illustration.

Rank-based regression for analysis of repeated measures

We consider rank-based regression models for repeated measures. To account for possible withinsubject correlations, we decompose the total ranks into between- and within-subject ranks and obtain two different estimators based on between- and within-subject ranks. A simple perturbation method is then introduced to generate bootstrap replicates of the estimating functions and the parameter estimates. This provides a convenient way for combining the corresponding two types of estimating function for more efficient estimation.

Analysing commercial catch and effort data from a penaeid trawl fishery: A comparison of linear models, mixed models, and generalised estimating equations approaches

Statistical methods are often used to analyse commercial catch and effort data to provide standardised fishing effort and/or a relative index of fish abundance for input into stock assessment models. Achieving reliable results has proved difficult in Australia's Northern Prawn Fishery (NPF), due to a combination of such factors as the biological characteristics of the animals, some aspects of the fleet dynamics, and the changes in fishing technology. For this set of data, we compared four modelling approaches (linear models, mixed models, generalised estimating equations, and generalised linear models) with respect to the outcomes of the standardised fishing effort or the relative index of abundance. We also varied the number and form of vessel covariates in the models. Within a subset of data from this fishery, modelling correlation structures did not alter the conclusions from simpler statistical models. The random-effects models also yielded similar results. This is because the estimators are all consistent even if the correlation structure is mis-specified, and the data set is very large. However, the standard errors from different models differed, suggesting that different methods have different statistical efficiency. We suggest that there is value in modelling the variance function and the correlation structure, to make valid and efficient statistical inferences and gain insight into the data. We found that fishing power was separable from the indices of prawn abundance only when we offset the impact of vessel characteristics at assumed values from external sources. This may be due to the large degree of confounding within the data, and the extreme temporal changes in certain aspects of individual vessels, the fleet and the fleet dynamics.

Calculating the Fickian diffusivity for a lattice-based random walk with agents and obstacles of different shapes and sizes

Random walk models are often used to interpret experimental observations of the motion of biological cells and molecules. A key aim in applying a random walk model to mimic an in vitro experiment is to estimate the Fickian diffusivity (or Fickian diffusion coefficient),D. However, many in vivo experiments are complicated by the fact that the motion of cells and molecules is hindered by the presence of obstacles. Crowded transport processes have been modeled using repeated stochastic simulations in which a motile agent undergoes a random walk on a lattice that is populated by immobile obstacles. Early studies considered the most straightforward case in which the motile agent and the obstacles are the same size. More recent studies considered stochastic random walk simulations describing the motion of an agent through an environment populated by obstacles of different shapes and sizes. Here, we build on previous simulation studies by analyzing a general class of lattice-based random walk models with agents and obstacles of various shapes and sizes. Our analysis provides exact calculations of the Fickian diffusivity, allowing us to draw conclusions about the role of the size, shape and density of the obstacles, as well as examining the role of the size and shape of the motile agent. Since our analysis is exact, we calculateDdirectly without the need for random walk simulations. In summary, we find that the shape, size and density of obstacles has a major influence on the exact Fickian diffusivity. Furthermore, our results indicate that the difference in diffusivity for symmetric and asymmetric obstacles is significant.

Does psychological stress mediate social deprivation in tooth loss?

It is unclear which theoretical dimension of psychological stress affects health status. We hypothesized that both distress and coping mediate the relationship between socio-economic position and tooth loss. Cross-sectional data from 2915 middle-aged adults evaluated retention of < 20 teeth, behaviors, psychological stress, and sociodemographic characteristics. Principal components analysis of the Perceived Stress Scale (PSS) extracted 'distress' (a = 0.85) and 'coping' (a =0.83) factors, consistent with theory. Hierarchical entry of explanatory variables into age- and sex-adjusted logistic regression models estimated odds ratios (OR) and 95% confidence intervals [95% CI] for retention of < 20 teeth. Analysis of the separate contributions of distress and coping revealed a significant main effect of coping (OR = 0.7 [95% CI = 0.7-0.8]), but no effect for distress (OR = 1.0 [95% CI = 0.9-1.1]) or for the interaction of coping and distress. Behavior and psychological stress only modestly attenuated socio-economic inequality in retention of < 20 teeth, providing evidence to support a mediating role of coping.

Exploratory spatial analysis of social and environmental factors associated with the incidence of Ross River Virus in Brisbane, Australia

We used geographic information systems and a spatial analysis approach to explore the pattern of Ross River virus (RRV) incidence in Brisbane, Australia. Climate, vegetation and socioeconomic data in 2001 were obtained from the Australian Bureau of Meteorology, the Brisbane City Council and the Australian Bureau of Statistics, respectively. Information on the RRV cases was obtained from the Queensland Department of Health. Spatial and multiple negative binomial regression models were used to identify the socioeconomic and environmental determinants of RRV transmission. The results show that RRV activity was primarily concentrated in the northeastern, northwestern, and southeastern regions in Brisbane. Multiple negative binomial regression models showed that the spatial pattern of RRV disease in Brisbane seemed to be determined by a combination of local ecologic, socioeconomic, and environmental factors.

Long-term trends in fine particle number concentrations in the urban atmosphere of Brisbane : the relevance of traffic emissions and new particle formation

The measurement of submicrometre (< 1.0 m) and ultrafine particles (diameter < 0.1 m) number concentration have attracted attention since the last decade because the potential health impacts associated with exposure to these particles can be more significant than those due to exposure to larger particles. At present, ultrafine particles are not regularly monitored and they are yet to be incorporated into air quality monitoring programs. As a result, very few studies have analysed their long-term and spatial variations in ultrafine particle concentration, and none have been in Australia. To address this gap in scientific knowledge, the aim of this research was to investigate the long-term trends and seasonal variations in particle number concentrations in Brisbane, Australia. Data collected over a five-year period were analysed using weighted regression models. Monthly mean concentrations in the morning (6:00-10:00) and the afternoon (16:00-19:00) were plotted against time in months, using the monthly variance as the weights. During the five-year period, submicrometre and ultrafine particle concentrations increased in the morning by 105.7% and 81.5% respectively whereas in the afternoon there was no significant trend. The morning concentrations were associated with fresh traffic emissions and the afternoon concentrations with the background. The statistical tests applied to the seasonal models, on the other hand, indicated that there was no seasonal component. The spatial variation in size distribution in a large urban area was investigated using particle number size distribution data collected at nine different locations during different campaigns. The size distributions were represented by the modal structures and cumulative size distributions. Particle number peaked at around 30 nm, except at an isolated site dominated by diesel trucks, where the particle number peaked at around 60 nm. It was found that ultrafine particles contributed to 82%-90% of the total particle number. At the sites dominated by petrol vehicles, nanoparticles (< 50 nm) contributed 60%-70% of the total particle number, and at the site dominated by diesel trucks they contributed 50%. Although the sampling campaigns took place during different seasons and were of varying duration these variations did not have an effect on the particle size distributions. The results suggested that the distributions were rather affected by differences in traffic composition and distance to the road. To investigate the occurrence of nucleation events, that is, secondary particle formation from gaseous precursors, particle size distribution data collected over a 13 month period during 5 different campaigns were analysed. The study area was a complex urban environment influenced by anthropogenic and natural sources. The study introduced a new application of time series differencing for the identification of nucleation events. To evaluate the conditions favourable to nucleation, the meteorological conditions and gaseous concentrations prior to and during nucleation events were recorded. Gaseous concentrations did not exhibit a clear pattern of change in concentration. It was also found that nucleation was associated with sea breeze and long-range transport. The implications of this finding are that whilst vehicles are the most important source of ultrafine particles, sea breeze and aged gaseous emissions play a more important role in secondary particle formation in the study area.

Teacher-child relationship quality for young children with parent reported language concerns

Previous research has demonstrated the importance of the qualities of the teacher-child relationship on children’s development. Close teacher-child relationships are especially important for children at risk. Positive relationships have been shown to have beneficial effects on children’s social and academic development (Birch & Ladd, 1997; Pianta & Stuhlman, 2004). Children with language difficulties are likely to face increased risks with regard to long term social and academic outcomes. The purpose of the current research was to gain greater understanding of the qualities of teacher-child relationships for young children with parent reported language concerns. The research analyses completed for this thesis involved the use of data from the public-access database of Growing Up in Australia: The Longitudinal Study of Australian Children (LSAC). LSAC is a longitudinal study involving a nationally representative sample of 10,000 Australian children. Data are being collected biennially from 2004 (Wave 1 data collection) until 2010 (Wave 4 data collection). LSAC has a cross-sequential research design involving two cohorts, an infant cohort (0-1 year at age of recruitment) and a kindergarten cohort (4-5 years at age of recruitment). Two studies are reported in this thesis using data for the LSAC Kindergarten Cohort which had 4983 child participants at recruitment. Study 1 used Wave 1 data to identify the differences between teacher-child relationship qualities for children with parent reported language concerns and their peers. Children identified by parents for whom concerns were held about their receptive and expressive language, as measured by items from the Parents’ Evaluation of Developmental Status (PEDS) (Glascoe, 2000) were the target (at risk) group in the study (n = 210). A matched case control group of peers (n = 210), matched on the child characteristics of sex, age, cultural and linguistic differences (CALD), and socio-economic positioning (SEP), were the comparison group for this analysis. Teacher-child relationship quality was measured by teacher reports on the Closeness and Conflict scales from the short version of the Student-Teacher Relationship Scale (STRS) (Pianta, 2001). There were statistically significant differences in the levels of closeness and conflict between the two groups. The target group had relationships with their teachers that had lower levels of closeness and higher levels of conflict than the control group. Study 2 reports analyses that examined the stability of the qualities of the teacher-child relationships at Wave 1 (4-5 years) and the qualities of the teacher-child relationships at Wave 2 (6-7 years). This time frame crosses the period of the children’s transition to school. The study examined whether early patterns in the qualities of the teacher-child relationship for children with parent reported language concerns at Wave 1 predicted the qualities of the teacher-child relationship outcomes in the early years of formal school. The sample for this study consisted of the group of children identified with PEDS language concerns at Wave 1 who also had teacher report data at Wave 2 (n = 145). Teacher-child relationship quality at Wave 1 and Wave 2 was again measured by the STRS scales of Closeness and Conflict. Results from multiple regression models indicated that teacher-child relationship quality at Wave 1 significantly contributed to the prediction of the quality of the teacher-child relationship at Wave 2, beyond other predictor variables included in the regression models. Specifically, Wave 1 STRS Closeness scores were the most significant predictor for STRS Closeness scores at Wave 2, while Wave 1 STRS Conflict scores were the only significant predictor for Wave 2 STRS Conflict outcomes. These results indicate that the qualities of the teacher-child relationship experienced prior to school by children with parent reported language concerns remained stable across transitions into formal schooling at which time the child had a different teacher. The results of these studies provide valuable insight into the nature of teacher-child relationship quality for young children with parent reported language concerns. These children experienced teacher-child relationships of a lower quality when compared with peers and, additionally, the qualities of these relationships prior to formal schooling were predictive of the qualities of the relationships in the early years of formal schooling. This raises concerns, given the increased risks of poorer social and academic outcomes already faced by children with language difficulties, that these early teacher-child relationships have an impact on future teacher-child relationships. Results of these studies are discussed with these considerations in mind and also discussed in terms of the implications for educational theory, policy and practice.

Space-time variation of malaria incidence in Yunnan province, China

Abstract Background Understanding spatio-temporal variation in malaria incidence provides a basis for effective disease control planning and monitoring. Methods Monthly surveillance data between 1991 and 2006 for Plasmodium vivax and Plasmodium falciparum malaria across 128 counties were assembled for Yunnan, a province of China with one of the highest burdens of malaria. County-level Bayesian Poisson regression models of incidence were constructed, with effects for rainfall, maximum temperature and temporal trend. The model also allowed for spatial variation in county-level incidence and temporal trend, and dependence between incidence in June–September and the preceding January–February. Results Models revealed strong associations between malaria incidence and both rainfall and maximum temperature. There was a significant association between incidence in June–September and the preceding January–February. Raw standardised morbidity ratios showed a high incidence in some counties bordering Myanmar, Laos and Vietnam, and counties in the Red River valley. Clusters of counties in south-western and northern Yunnan were identified that had high incidence not explained by climate. The overall trend in incidence decreased, but there was significant variation between counties. Conclusion Dependence between incidence in summer and the preceding January–February suggests a role of intrinsic host-pathogen dynamics. Incidence during the summer peak might be predictable based on incidence in January–February, facilitating malaria control planning, scaled months in advance to the magnitude of the summer malaria burden. Heterogeneities in county-level temporal trends suggest that reductions in the burden of malaria have been unevenly distributed throughout the province.

Using Trauma Injury Severity Score (TRISS) variables to predict length of hospital stay following trauma in New Zealand

Aim – To develop and assess the predictive capabilities of a statistical model that relates routinely collected Trauma Injury Severity Score (TRISS) variables to length of hospital stay (LOS) in survivors of traumatic injury. Method – Retrospective cohort study of adults who sustained a serious traumatic injury, and who survived until discharge from Auckland City, Middlemore, Waikato, or North Shore Hospitals between 2002 and 2006. Cubic-root transformed LOS was analysed using two-level mixed-effects regression models. Results – 1498 eligible patients were identified, 1446 (97%) injured from a blunt mechanism and 52 (3%) from a penetrating mechanism. For blunt mechanism trauma, 1096 (76%) were male, average age was 37 years (range: 15-94 years), and LOS and TRISS score information was available for 1362 patients. Spearman’s correlation and the median absolute prediction error between LOS and the original TRISS model was ρ=0.31 and 10.8 days, respectively, and between LOS and the final multivariable two-level mixed-effects regression model was ρ=0.38 and 6.0 days, respectively. Insufficient data were available for the analysis of penetrating mechanism models. Conclusions – Neither the original TRISS model nor the refined model has sufficient ability to accurately or reliably predict LOS. Additional predictor variables for LOS and other indicators for morbidity need to be considered.

Contemporary New Zealand coefficients for the Trauma Injury Severity Score: TRISS(NZ)

Aims – To develop local contemporary coefficients for the Trauma Injury Severity Score in New Zealand, TRISS(NZ), and to evaluate their performance at predicting survival against the original TRISS coefficients. Methods – Retrospective cohort study of adults who sustained a serious traumatic injury, and who survived until presentation at Auckland City, Middlemore, Waikato, or North Shore Hospitals between 2002 and 2006. Coefficients were estimated using ordinary and multilevel mixed-effects logistic regression models. Results – 1735 eligible patients were identified, 1672 (96%) injured from a blunt mechanism and 63 (4%) from a penetrating mechanism. For blunt mechanism trauma, 1250 (75%) were male and average age was 38 years (range: 15-94 years). TRISS information was available for 1565 patients of whom 204 (13%) died. Area under the Receiver Operating Characteristic (ROC) curves was 0.901 (95%CI: 0.879-0.923) for the TRISS(NZ) model and 0.890 (95% CI: 0.866-0.913) for TRISS (P<0.001). Insufficient data were available to determine coefficients for penetrating mechanism TRISS(NZ) models. Conclusions – Both TRISS models accurately predicted survival for blunt mechanism trauma. However, TRISS(NZ) coefficients were statistically superior to TRISS coefficients. A strong case exists for replacing TRISS coefficients in the New Zealand benchmarking software with these updated TRISS(NZ) estimates.

Temperature enhanced effects of ozone on cardiovascular mortality in 95 large US communities, 1987-2000 - assessment using the NMMAPS data

A few studies examined interactive effects between air pollution and temperature on health outcomes. This study is to examine if temperature modified effects of ozone and cardiovascular mortality in 95 large US cities. A nonparametric and a parametric regression models were separately used to explore interactive effects of temperature and ozone on cardiovascular mortality during May and October, 1987-2000. A Bayesian meta-analysis was used to pool estimates. Both models illustrate that temperature enhanced the ozone effects on mortality in the northern region, but obviously in the southern region. A 10-ppb increment in ozone was associated with 0.41 % (95% posterior interval (PI): -0.19 %, 0.93 %), 0.27 % (95% PI: -0.44 %, 0.87 %) and 1.68 % (95% PI: 0.07 %, 3.26 %) increases in daily cardiovascular mortality corresponding to low, moderate and high levels of temperature, respectively. We concluded that temperature modified effects of ozone, particularly in the northern region.

An investigation into how work-related road safety can be enhanced

Despite the facts that vehicle incidents continue to be the most common mechanism for Australian compensated fatalities and that employers have statutory obligations to provide safe workplaces, very few organisations are proactively and comprehensively managing their work-related road risks. Unfortunately, limited guidance is provided in the existing literature to assist practitioners in managing work-related road risks. The current research addresses this gap in the literature. To explore how work-related road safety can be enhanced, three studies were conducted. Study one explored the effectiveness of a range of risk management initiatives and whether comprehensive risk management practices were associated with safety outcomes. Study two explored barriers to, and facilitators for, accepting risk management initiatives. Study three explored the influence of organisational factors on road safety outcomes to identify optimal work environments for managing road risks. To maximise the research sample and increase generalisability, the studies were designed to allow data collection to be conducted simultaneously drawing upon the same sample obtained from four Australian organisations. Data was collected via four methods. A structured document review of published articles was conducted to identify what outcomes have been observed in previously investigated work-related road safety initiatives. The documents reviewed collectively assessed the effectiveness of 19 work-related road safety initiatives. Audits of organisational practices and process operating within the four researched organisations were conducted to identify whether organisations with comprehensive work-related road risk management practices and processes have better safety outcomes than organisations with limited risk management practices and processes. Interviews were conducted with a sample of 24 participants, comprising 16 employees and eight managers. The interviews were conducted to identify what barriers and facilitators within organisations are involved in implementing work-related road safety initiatives and whether differences in fleet safety climate, stage of change and safety ownership relate to work-related road safety outcomes. Finally, questionnaires were administered to a sample of 679 participants. The questionnaires were conducted to identify which initiatives are perceived by employees to be effective in managing work-related road risks and whether differences in fleet safety climate, stage of change and safety ownership relate to work-related road safety outcomes. Seven research questions were addressed in the current research project. The key findings with respect to each of the research questions are presented below. Research question one: What outcomes have been observed in previously investigated work-related road safety initiatives? The structured document review indicated that initiatives found to be positively associated with occupational road safety both during and after the intervention period included: a pay rise; driver training; group discussions; enlisting employees as community road safety change agents; safety reminders; and group and individual rewards. Research question two: Which initiatives are perceived by employees to be effective in managing work-related road risks? Questionnaire findings revealed that employees believed occupational road risks could best be managed through making vehicle safety features standard, providing practical driver skills training and through investigating serious vehicle incidents. In comparison, employees believed initiatives including signing a promise card commitment to drive safely, advertising the organisation’s phone number on vehicles and consideration of driving competency in staff selection process would have limited effectiveness in managing occupational road safety. Research question three: Do organisations with comprehensive work-related road risk management practices and processes have better safety outcomes than organisations with limited risk management practices and processes? The audit identified a difference among the organisations in their management of work-related road risks. Comprehensive risk management practices were associated with employees engaging in overall safer driving behaviours, committing less driving errors, and experiencing less fatigue and distraction issues when driving. Given that only four organisations participated in this research, these findings should only be considered as preliminary. Further research should be conducted to explore the relationship between comprehensiveness of risk management practices and road safety outcomes with a larger sample of organisations. Research question four: What barriers and facilitators within organisations are involved in implementing work-related road safety initiatives? The interviews identified that employees perceived six organisational characteristics as potential barriers to implementing work-related road safety initiatives. These included: prioritisation of production over safety; complacency towards work-related road risks; insufficient resources; diversity; limited employee input in safety decisions; and a perception that road safety initiatives were an unnecessary burden. In comparison, employees perceived three organisational characteristics as potential facilitators to implementing work-related road safety initiatives. These included: management commitment; the presence of existing systems that could support the implementation of initiatives; and supportive relationships. Research question five: Do differences in fleet safety climate relate to work-related road safety outcomes? The interviews and questionnaires identified that organisational climates with high management commitment, support for managing work demands, appropriate safety rules and safety communication were associated with employees who engaged in safer driving behaviours. Regression analyses indicated that as participants’ perceptions of safety climate increased, the corresponding likelihood of them engaging in safer driving behaviours increased. Fleet safety climate was perceived to influence road safety outcomes through several avenues. Some of these included: the allocation of sufficient resources to manage occupational road risks; fostering a supportive environment of mutual responsibility; resolving safety issues openly and fairly; clearly communicating to employees that safety is the top priority; and developing appropriate work-related road safety policies and procedures. Research question six: Do differences in stage of change relate to work-related road safety outcomes? The interviews and questionnaires identified that participants’ perceptions of initiative effectiveness were found to vary with respect to their individual stage of readiness, with stage-matched initiatives being perceived most effective. In regards to safety outcomes, regression analyses identified that as participants’ progress through the stages of change, the corresponding likelihood of them being involved in vehicle crashes decreases. Research question seven: Do differences in safety ownership relate to work-related road safety outcomes? The interviews and questionnaires revealed that management of road risks is often given less attention than other areas of health and safety management in organisations. In regards to safety outcomes, regression analyses identified that perceived authority and perceived shared ownership both emerged as significant independent predictors of self-reported driving behaviours pertaining to fatigue and distractions. The regression models indicated that as participants’ perceptions of the authority of the person managing road risks increases, and perceptions of shared ownership of safety tasks increases, the corresponding likelihood of them engaging in driving while fatigued or multitasking while driving decreases. Based on the findings from the current research, the author makes several recommendations to assist practitioners in developing proactive and comprehensive approaches to managing occupational road risks. The author also suggests several avenues for future research in the area of work-related road safety.

Different cognitive profiles for single compared with recurrent fallers without Dementia

Relationships between self-reported retrospective falls and cognitive measures (executive function, reaction time, processing speed, working memory, visual attention) were examined in a population based sample of older adults (n = 658). Two of the choice reaction time tests involved inhibiting responses to either targets of a specific color or location with hand and foot responses. Potentially confounding demographic variables, medical conditions and postural sway were controlled for in logistic regression models, excluding participants with possible cognitive impairment. A factor analysis of cognitive measures extracted factors measuring reaction time, accuracy and inhibition, and visual search. Single fallers did not differ from non-fallers in terms of health, sway or cognitive function, except that they performed worse on accuracy and inhibition. In contrast, recurrent fallers performed worse than non-fallers on all measures. Results suggest that occasional falls in late life may be associated with subtle age-related changes in the pre-frontal cortex leading to failures of executive control, whereas recurrent falling may result from more advanced brain ageing that is associated with generalized cognitive decline.

Confidence to cook vegetables and the buying habits of Australian households

Cooking skills are emphasized in nutrition promotion but their distribution among population subgroups and relationship to dietary behavior is researched by few population-based studies. This study examined the relationships between confidence to cook, sociodemographic characteristics, and household vegetable purchasing. This cross-sectional study of 426 randomly selected households in Brisbane, Australia, used a validated questionnaire to assess household vegetable purchasing habits and the confidence to cook of the person who most often prepares food for these households. The mutually adjusted odds ratios (ORs) of lacking confidence to cook were assessed across a range of demographic subgroups using multiple logistic regression models. Similarly, mutually adjusted mean vegetable purchasing scores were calculated using multiple linear regression for different population groups and for respondents with varying confidence levels. Lacking confidence to cook using a variety of techniques was more common among respondents with less education (OR 3.30; 95% confidence interval [CI] 1.01 to 10.75) and was less common among respondents who lived with minors (OR 0.22; 95% CI 0.09 to 0.53) and other adults (OR 0.43; 95% CI 0.24 to 0.78). Lack of confidence to prepare vegetables was associated with being male (OR 2.25; 95% CI 1.24 to 4.08), low education (OR 6.60; 95% CI 2.08 to 20.91), lower household income (OR 2.98; 95% CI 1.02 to 8.72) and living with other adults (OR 0.53; 95% CI 0.29 to 0.98). Households bought a greater variety of vegetables on a regular basis when the main chef was confident to prepare them (difference: 18.60; 95% CI 14.66 to 22.54), older (difference: 8.69; 95% CI 4.92 to 12.47), lived with at least one other adult (difference: 5.47; 95% CI 2.82 to 8.12) or at least one minor (difference: 2.86; 95% CI 0.17 to 5.55). Cooking skills may contribute to socioeconomic dietary differences, and may be a useful strategy for promoting fruit and vegetable consumption, particularly among socioeconomically disadvantaged groups.