Summer activity patterns among teenage girls: harmonic shape invariant modeling to estimate circadian cycles
© Ogbagaber et al.; licensee BioMed Central Ltd. 2012
Received: 13 February 2012
Accepted: 15 April 2012
Published: 6 May 2012
Physical activity as measured by activity counts over short time intervals across a 24 h period are often used to assess circadian variation. We are interested in characterizing circadian patterns in activity among adolescents and examining how these patterns vary by obesity status. New statistical approaches are needed to examine how factors affect different features of the circadian pattern and to make appropriate covariate adjustments when the outcomes are longitudinal count data.
We develop a statistical model for longitudinal or repeated activity count data that is used to examine differences in the overall activity level, amplitude (defined as the difference between the lowest and highest activity level over a 24 hour period), and phase shift. Using seven days of continuous activity monitoring, we characterize the circadian patterns and compare them between obese and non-obese adolescent girls.
We find a statistically significant phase delay in adolescent girls who were obese compared with their non-obese counterparts. After the appropriate adjustment for measured potential confounders, we did not find differences in mean activity level between the two groups.
New statistical methodology was developed to identify a phase delay in obese compared with non-obese adolescents. This new approach for analyzing longitudinal circadian rhythm count data provides a useful statistical technique to add to the repertoire for those analyzing circadian rhythm data.
The purpose of this study was to characterize circadian rhythms patterns of activity among adolescent girls ages 15 to 16. Numerous health benefits have been associated with adolescent lifestyles that include moderate to vigorous physical activity (PA) . Furthermore, adolescent PA tracks into adulthood and relates to adult obesity [2–4]. However, most adolescent girls do not meet recommendations for daily PA  and levels of PA decrease in girls during adolescence [6–8]. Changes in sleep health and circadian regulation have also been associated with decreased physical activity as well as health problems (e.g., obesity and diabetes), and decrements in quality of life and neurocognitive function [9–14].
In general the development of methodology to characterize changes in physical activity and to estimate circadian rhythms and the timing of sleep has the potential to provide insight into the causes of disease and could serve as a tool to evaluate treatment outcomes. In this particular case, characterizing circadian patterns in adolescent girls may provide guidance in targeting particular PA lifestyles or times of day most conducive to intervention. Although there is evidence for objective differences in the daily quantity of vigorous physical activity engaged in by normal weight versus overweight children and adolescents [15, 16], less is known about how PA in overweight and normal adolescents may differ on other dimensions, such as time of day for peak activity, differences in periods of peak and minimal activity, and how PA fits within their daily rhythms.
The NEXT Generation Health Study provides a unique resource for characterization of adolescent PA and also provides an opportunity to examine the role of obesity on the circadian rhythm of PA. Specifically, are circadian PA patterns different across optimal weight and overweight groups?
In order to address these questions, we developed a new statistical modeling approach that allows us to characterize the effect of obesity on important features of the circadian rhythm. The proposed methodology allowed for analyses with and without adjustments for demographic variables. We proposed a shape invariant model for the effects of covariates on the circadian rhythm patterns in longitudinal activity (count) measurements. This modeling approach assumes that all subjects have the same underlying circadian pattern with differences across individuals and subgroups being reflected by changes in the mean, amplitude and phase shift of the underlying circadian pattern in activity that can also be viewed as changes in sleep period timing.
The study of circadian rhythms is common in the biological and social sciences  and is becoming increasingly important for health as circadian clock genes  have been identified in both neural and non-neural tissue. Coordination of these clock genes through the body may be critical for metabolic function, immunity and tissue repair as well as neurocognitive function. While the health effects of circadian regulation have only been studied in recent years, there is a longer history of interest in estimating the effect of important groupings or covariates on features of the circadian rhythm. These features can be characterized by an overall mean, amplitude (defined here as the distance from the lowest to the highest point in the rhythm), and phase shift (i.e., shifting of the whole pattern).
Others have proposed approaches for statistical analysis in ciradian rhythm longitudinal data [19, 20]. In this paper, we develop an approach that allows investigators to compare circadian patterns on longitudinal count data in terms of amplitude, phase-shift, and overall mean, while adjusting for confounding factors.
Wang, Kee and Brown  and Albert and Hunsberger  proposed a regression-based approach for analyzing longitudinal continuous circadian rhythm data where individual variations were incorporated through random effects added to the mean, amplitude, and phase shift. The adaptation of this approach to the analysis of longitudinal circadian rhythm count data (as compared to continuous data) requires the development of new statistical methodology, which we develop in this paper.
Materials and methods
We adapted Albert and Hunsberger’s approach to the Poisson regression framework. The goal was to develop a simple approach to estimating circadian rhythms that can be easily implemented by practitioners and can be used for analyzing longitudinal activity data.
The analysis is based on data obtained from the NEXT Generation Health Study. Data was collected by the Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD) of the Nation Institutes of Health (NIH) in the summer of 2010. Analyses were performed using SAS version 9.2 software  and STATA version 11 . Model fitting was done using PROC NLMIXED in SAS.
A nationally-representative cohort of U.S students in grade 10 was recruited using a multistage stratified design. Primary sampling units consisted of school districts or groups of school districts stratified across the nine U.S. Census divisions. Within this sampling framework 137 schools were selected and formally recruited; 80 (58.4%) agreed to participate. Tenth-grade classes were randomly selected within each recruited school and 3,796 students were recruited to participate; youth assent and parental consent were obtained from 2,619 (69.0%) students. Of those who consented, 2,519 (96.24%) completed the Wave 1 survey which was administered by trained research assistants. In a US nationally representative sample, the subsample of African American youth might not be sufficient to support robust statistical comparisons across race/ethnicity; therefore, additional African-American students were recruited into the sample. The prevalence of Hispanic youth in a national sample of this age group did not make it necessary to oversample Hispanic youth.
The study protocol was reviewed and approved by the Institutional Review Board of the Eunice Kennedy Shriver National Institute of Child Health and Human Development.
Demographic information was provided by the adolescent and her parents. Race/ethnicity of the adolescent was classified into four categories: Hispanic, African-American, Asian, and White. In the analysis there were only three categories: Hispanic, African-American, and Other. Since there were only 2 Asians in the first classification, we combined them with the “White” category to create “Other”. The adolescent provided an estimate of family socioeconomic status using the Family Affluence Scale (FAS), a validated measure of socioeconomic status . The adolescent answered four questions regarding material conditions of their household: “Does your family own a car, van or truck” (none = 0, 1 = 1, 2 or more = 2); “ Do you have your own bedroom for yourself” (no = 0, yes = 1); “ During the past 12 months, how many times did you travel away on holiday with your family”(not at all = 0, once = 1, twice = 2, more than twice = 3); and “ How many computers does your family own” (none = 0, one = 1, two = 2, more than two = 3). The FAS was used as a continuous scale in analyses with possible scores ranging from 0 to 9.
Adolescents’ primary caretakers reported their education level using a seven point scale: less than a high school diploma; high school diploma; GED; some college or technical school; associate’s degree; bachelor’s degree; or graduate degree. And, if applicable, reported education level of another parent or guardian providing support for the adolescent whether living with the adolescent or living separately.
Height was assessed using a portable stadiometer places on a level, hard surface. Students removed shoes and measures were taken at least twice. If the first two measures were not within ± 1.0 cm of each other, measurement was repeated. Failure to obtain two measures ± 1.0 cm of each other resulted in measures being repeated by a supervisor. Weight was measured using a Healthometer Model 498 KL Digital Scale on a hard, level surface. Students removed heavy objects from their pockets and extra outerwear (e.g., sweater, sweat shirt, or jacket). If the first two measurements were not within ± 0.2 kg, a third measurement was taken. Mean height and weight values of the two measures meeting the criteria were used to calculate Body Mass Index (BMI = wt(kg)/ht(m)2). Weight status (underweight, normal weight, overweight, and obese) were determined from BMI-for-age percentiles for each gender using the CDC 2000 growth chart ). Underweight was defined as a BMI below the 5th percentile; normal weight was ≥ 5th but < 85th percentile; overweight was when BMI was ≥ 85th but < 95th percentile; and obese was ≥ 95th percentile. For these analyses, overweight and obese were combined and are labeled ‘obese’ in subsequent descriptions and normal weight (‘non-obese’) was the comparison group.
Based on weight status, 281 overweight/obese and 286 normal weight adolescents were recruited from 40 of the participating schools for additional assessment procedures. Underweight adolescents were excluded. For one of the additional assessments, adolescents in both groups were asked to wear an Actiwatch2 on their non-dominant wrist for 24 hours/day for seven consecutive days. This device has been validated as a measure of activity in youth. The Actiwatch2 recorded motion in 30-sec epochs. Although the primary purpose of the Actiwatch2 is to assess sleep patterns, because it can easily be worn throughout the day, it provides an acceptable measure of activity during the day 23 .
A subsample of 96 female adolescents was selected for analysis in this paper (19 obese and 77 non-obese participants). The criteria for selection included: 1) having worn the Actigraph2 for a full 7 days; 2) the entire period of assessment was during the summer when school is not in session in order to avoid the effect of differences in school schedules; and 3) female (to avoid gender differences in activity level). Activity counts collected in the 30-sec epochs were summed over 15 min periods to provide 96 observations/day quantifying activity. This resulted in 672 observations within each student for 7 days.
We propose a shape-invariant Poisson model for circadian rhythm count data, which is an adaptation of Albert and Hunsberger  to the longitudinal activity example. Albert and Hunsberger developed a random effects shape invariant model that incorporates differences between individual circadian rhythms of cortisol, a highly predictable biological rhythm, by allowing important features, mean, amplitude and shift parameters, to vary based on fixed effect covariates and individual random effects. A direct extension to the Poisson outcome (or more generally to a generalized linear model outcome) would require the development of new non-standard software, and would be difficult to implement for the practitioner. We propose a simple two-stage approach that is simple for the practitioner to implement.
In our analysis, we chose K = 2 (two harmonic terms) since this provided a flexible pattern that nicely characterized our data (i.e. mean curves from fitted models are close to the empirical means). Analyses with K = 3 provided similar results (data not shown). We used an inverse logit transformation for , which restricts the phase shift parameter to be within the interval [0, 1], corresponding to characterizing the shift in terms of the percentage of a 24 hour period. For simplicity, the time variable of a 24-hour period (t) is rescaled to [0, 1]. For a typical day, one segment corresponds to 15 minutes from 12:00 AM to 12:15 AM (i.e., t = 1/96). The second segment corresponds to the time from 12:15 AM to 12:30 AM (i.e., t = 2/96), etc. There are 4 such segments in 1 hour and 96 segments in 24 hours. The distribution of mean, amplitude, and phase-shift parameters can be compared across populations using simple statistical methods such as two-sample tests and regression techniques.
Although the modeling framework assumes Poisson counts given individual parameters of intercept, amplitude, and phase-shift, the proposed estimation procedures can easily be altered by using a negative binomial model. However, the effect of covariates on changes on harmonic parameters using Poisson regression should be robust to the assumption of no overdispersion.
We propose a two-stage approach for estimation of the shape invariant model for longitudinal activity count data. In the first stage, subject-specific model parameters are estimated by iterating between fitting individual Poisson regression models and using a nonlinear optimizer which involve the 2 steps described below. In the second stage, we regress the individual model parameters on important covariates.
The first stage can be implemented as follows
We iterate between Step 1 and 2 until we have convergence.
In the second stage, we regress individual estimates obtained from stage 1 on covariates such as BMI. For the dichotomous single BMI covariate this simplifies to conducting a t-test on the three components characterizing the circadian pattern. A regression model can be applied when adjusting for confounding factors. Estimation of , , and D using the proposed two-step procedure may induce dependence in estimates across individuals, a violation of a key assumption in standard regression or the t-test. We investigate if any induced dependence between , , and D might affect the type I error rate of the statistical test.
A simulation was conducted to examine whether the two-side t-test in stage 2 results in valid hypothesis tests (i.e., an alpha-level procedure). More specifically, we were interested if the t-test was rejecting the null hypothesis of no significant BMI differences 5% of the time when BMI has no effect on the circadian pattern. We simulated 2000 circadian rhythm data sets. The simulation verified that the nominal is contained within a 95% confidence interval for the rejection rate for each parameter. That is, : (0.0395, 0.058), : (0.046, 0.066), : (0.045, 0.065). Since, in each case, the interval contains 0.05, the two-stage approach appears to be a valid test.
N (%) or Mean (SD)
2: African American
Race (redefined) *
2: African American
BMI Group *
0: Low BMI Group
1: High BMI Group
25.45 (5.45); N= 96
Family Affluence **†
5.05 (1.61); N= 96
Parental Education 1 **††
3.70 (1.79); N= 92
Parental Education 2 **†††
3.36 (1.73); N= 61
All individuals contributed 96 segments of 15 minute activity count data for 7 days of the week. The two stage analysis was performed to assess the association between obesity (obese and non-obese) and the circadian rhythm of activity levels.
A visual inspection of Figure 2 suggests that there was a shift in activity among the obese students. Specifically, obese adolescents tend to go to sleep and wake up later than non-obese adolescents. There also appears to be an overall lower mean activity level in obese as compared to non-obese adolescents. The proposed shape invariant statistical model for longitudinal activity count data can be used to formally test these empirical observations.
Differences in average circadian parameters between non-obese and obese participants
Difference Estimate (SE)
We compared the difference in circadian activity patterns for obese and non-obese participants after adjusting for race, family affluence, and parental education. Differences in the phase-shift between groups remained statistically significant (P = 0.016), while differences in the mean and amplitude were not ( P = 0.19 and P = 0.75, respectively). In addition to the previous analysis where we dichotomizing BMI into two categories of obesity status, we performed an analysis treating BMI as continuous(17.78 kg/m 2 to 47.98 kg/m 2 ). The second stage of this analysis used a regression model as compared with a t-test. After adjusting for race, family affluence, and parental education, the phase shift increased with BMI ( P = 0.0072) but BMI was not related to either the mean amplitude ( P = 0.39 and P = 0.62, respectively).
We developed a new statistical approach to examine the effect of various factors on the circadian patterns in longitudinal count data. The modeling approach was developed with the aim of determining whether circadian patterns in activity are different between obese and non-obese teenagers. Our approach allowed us to focus the comparison on the mean, amplitude, and phase-shift of the circadian patterns. We found that, after adjusting for potentially confounding factors such as parental education and income as well as race, there was a statistically significant phase delay in the circadian timing of sleep and activity for obese versus normal weight adolescent girls. These results raise important scientific questions regarding the contribution of circadian phase abnormalities, changes in sleep time and timing and patterns of activity on overweight and obesity.
Although we found a statistically significant relationship between the phase-shift in activity and obesity, we cannot determine whether overweight girls develop a tendency to begin the day later or the shift in circadian patterns results in subsequent gains in weight. Change in sleep patterns (shifting later in day) have been associated with sleep deprivation  and sleep patterns have been related to weight gain in adolescents  and TV watching . Much less is known about potential effects of shifts circadian patterns for physical activity.
Future work should examine whether these shifts in circadian patterns have implications for other health outcomes. For example, previous work has related adolescent morningness, or a tendency to prefer morning activities, with positive mental, social and physical health and adolescent eveningness with problem behaviors [31, 32]. In the current study, the shift in the circadian pattern is also associated with less physical activity and relationships have been found between levels of adolescent physical activity and indicators of social, behavioral and physical health .
The new statistical methodology allowed us to focus attention on important features of the circadian pattern and find this interesting association. The larger overall activity and phase shift among obese as compared with non-obese adolescents raises important questions about the association between circadian phase, activity and sleep patterns and obesity. An overall increase in activity in normal weight adolescent girls could be attributed to a healthier pattern of activity and increased capacity to engage in a range of activities, thereby contributing to normal weight. There are several explanations for the phase shift between obese and non-obese adolescents. A mis-match between circadian phase and the timing of sleep may be present in this sample of obese adolescents that could result in abnormalities in the hormones that regulate appetite, energy metabolism and adipose deposition. It is also possible that decreased energy expenditure during the day and night and the phase shift could be an indicator of abnormalities or differences in the metabolic rate. It is not clear whether this phase shift is a precursor or a manifestation of obesity or a shared marker of a metabolic type or abnormality. Also, it is possible that the delayed activity pattern among obese girls may be associated with factors that were not controlled for in this analysis. For example, in the NEXT study we did not collect information about whether participants had summer jobs or scheduled activities or whether they could sleep on an ad hoc basis. As this appears to confirm that this approach has some promise, it will be important to evaluate other contributors to circadian phase and lifestyle differences in subsequent studies.
This work was supported by the Intramural Research Program of the Eunice Kennedy Shriver National Institute of Child Health and Human Development and the National Heart Lung and Blood Institute.
- Pate RR, Pratt M, Blair SN, Haskell WL, Macera CA, Bouchard C, Buchner D, Ettinger W, Heath GW, King AC, Kriska A, Leon AS, Marcus BH, Morris J, Paffenbarger RS, Patrick K, Pollock ML, Rippe JM, Sallis J, Wilmore JH: Physical activity and public health: a recommendation from the Centers for Disease Control and Prevention and the American College of Sports Medicine. JAMA 1995, 273:402–407.PubMedView ArticleGoogle Scholar
- Hallal PC, Victora CG, Azevedo RM, Wells JCK: Adolescent physical activity and health. A systematic review. Sports Med 2006, 36:1019–1030.Google Scholar
- Strong WB, Malina RM, Bumke CJ, et al.: Evidence based physical activity for school-age children. J Pediatr 2005, 146:732–737.PubMedView ArticleGoogle Scholar
- Pietiläinen KH, Kaprio J, Borg P, Plasqui G, Yki-Järvinen H, Kujala UM, Rose RJ, Westerterp KR, Rissanen A: Physical inactivity and obesity: a vicious circle. Obesity 2008, 16:409–414.PubMedView ArticleGoogle Scholar
- Borraccino A, Lemma P, Iannotti R, Zambon A, Dalmasso P, Lazzeri G, Giacchi M, Cavallo F: Socio-economic effects on meeting PA guidelines: comparisons among 32 countries. Med Sci Sports Exerc 2009, 41:749–756.PubMedView ArticleGoogle Scholar
- Brodersen NH, Steptoe A, Boniface DR, Wardle J: Trends in physical activity and sedentary behaviour in adolescence: ethnic and socioeconomic differences. Br J Sports Med 2007,41(3):140–144.PubMedView ArticleGoogle Scholar
- Nader PR, Bradley RH, Houts RM, McRitchie SL, O’Brien M: Moderate-to-vigorous PA from ages 9 to 15 years. JAMA 2008,300(3):295–305.PubMedView ArticleGoogle Scholar
- Nelson MC, Neumark-Stzainer D, Hannan PJ, Sirard JR, Story M: Longitudinal and secular trends in PA and sedentary behavior during adolescence. Pediatrics 2006,118(6):e1627-e1634.PubMedView ArticleGoogle Scholar
- Carskadon MA, Acebo C, Jenni OG: Regulation of Adolescent Sleep: Implications for Behavior. Ann N Y Aca Sci 2004, 1021:276–291.View ArticleGoogle Scholar
- Beebe DW: Cognitive, behavioral, and functional consequences of inadequate sleep in children and adolescents. Pediatr Clin North Am 2011 Jun,58(3):649–665. doi:10.1016/j.pcl.2011.03.002.PubMedView ArticleGoogle Scholar
- McKnight-Eily LR, Eaton DK, Richard Lowry, Croft JB, Presley-Cantrell L, Perry GS: Relationships between hours of sleep and health-risk behaviors in US adolescent students. Prev. Med 2011, doi:10.1016/j.ypmed.2011.06.020.Google Scholar
- Owens JA, Belon K, Moss P: Impact of delaying school start time on adolescent sleep, mood, and behavior. Arch Pediatr Adolesc Med 2010 Jul,164(7):608–614.PubMedView ArticleGoogle Scholar
- Taheri S: The link between short sleep duration and obesity: we should recommend more sleep to prevent obesity. Arch Dis Child 2006,91(11):881–884.PubMedView ArticleGoogle Scholar
- Taveras EM, Rifas-Shiman SL, Oken E, Gunderson EP, Gillman MW: Short sleep duration in infancy and risk of childhood overweight. Arch Pediatr Adolesc Med 2008, 162:305–311.PubMedView ArticleGoogle Scholar
- Dorsey KB, Herrin J, Krumholz HM: Patterns of moderate and vigorous physical activity in obese and overweight compared with non-overweight children. Int J Pediatr Obes 2011, 6:e547-e555.PubMedView ArticleGoogle Scholar
- Patrick K, Norman GJ, Calfas KJ, Sallis JF, Zabinski MF, Rupp J, Cella J: Diet, physical activity, and sedentary behaviors as risk factors for overweight in adolescence. Arch Pediatr Adolesc Med 2004, 158:385–390.PubMedView ArticleGoogle Scholar
- Wehr TA, Goodwin FK: eds: Circadian Rhythms in Psychiatry. Boxwood Press, Pacific Grove, CA; 1983.Google Scholar
- Green CB, Takahashih JS, Bass J: Circadian rhythms and metabolism. Cell 2008, 134:728–742.PubMedView ArticleGoogle Scholar
- Brown EN, Czeisler CA: The statistical analysis of circadian phase and amplitute in constant-routine core-temperature data. J Biological Rhythms 1992, 7:177–202.View ArticleGoogle Scholar
- Wang J, Xian H, Licis A, Deych E, Ding J, McLeland J, Toedebusch C, Li T, Duntley S, Shannon W: Measuring the Impact of Apnea and Obesity on Circadian Activity Patterns Using Functional Linear Modeling of Actigraphy Data. J Circadian Rhythms 2011, 9:11.PubMedView ArticleGoogle Scholar
- Wang Y, Ke C, Brown MB: Shape-invariant modeling of circadian rhythms with random effects and smoothing spline ANOVA decompositions. Biometrics 2003, 59:804–812.PubMedView ArticleGoogle Scholar
- Albert PS, Hunsberger S: On analyzing circadian rhythms data using nonlinear mixed models with harmonic terms. Biometrics 2005, 61:1115–1120.PubMedView ArticleGoogle Scholar
- SAS User’s Guide: Statistics. SAS Institute, Cary, NC; 2008.Google Scholar
- StataCorp: Stata Statistical Software: Release 11. StataCorp LP, College Station, TX; 2009.Google Scholar
- Currie C, Molcho M, Boyce W, Holstein BR, Torsheim TR, Richter M: Researching health inequalities in adolescents: The development of the Health Behaviour in School-Aged Children (HBSC) Family Affluence Scale. Social Science & Medicine 2008,66(6):1429–1436.View ArticleGoogle Scholar
- Kuczmarski RJ, Ogden CL, Guo SS, Grummer-Strawn LM, Flegal KM, Mei Z, Wei R, Curtin LR, Roche AF, Johnson CL: 2000 CDC Growth Charts for the United States: methods and development. Vital Health Stat 11 2002, 246:1–190.PubMedGoogle Scholar
- Long AC, Palerno TM, Manees AM: Brief report: using actigraphy to compare physical activity levels in adolescents with chronic pain and healthy adolescents. J Pediatr Psychology 2008, 33:660–665.View ArticleGoogle Scholar
- Brown SA, McGue M, Maggs J, Schulenberg J, Hingson R, Swartzwelder S, Martin C, Chung T, Tapert SF, Sher K, Winters KC, Lowman C, Murphy S: A developmental perspective on alcohol and youths 16 to 20 years of age. Pediatrics 2008, 121:S290.PubMedView ArticleGoogle Scholar
- Chaput JP, Klingenberg L, Sjodin A: Do all sedentary activities lead to weight gain: sleep does not. Curr Opin Clin Nutr Metab Care 2010, 13:601–607.PubMedView ArticleGoogle Scholar
- Ortega FB, Chillón P, Ruiz JR, Delgado M, Albers U, Alvarez-Granda JL, Marcos A, Moreno LA, Castillo MJ: Sleep patterns in Spanish adolescents: associations with TV watching and leisure-time physical activity. Eur J Appl Physiol 2010, 110:563–573.PubMedView ArticleGoogle Scholar
- Lange L, Randler C: Morningness-eveningness and behavior problems in adolescents. Sleep Biol Rhythms 2011, 9:12–18.View ArticleGoogle Scholar
- Randler C: Association between morningness-eveningness and mental and physical health in adolescents. Psychol Health Med 2011, 16:29–38.PubMedView ArticleGoogle Scholar
- Iannotti RJ, Kogan MD, Janssen I, Boyce WF: Patterns of adolescent physical activity, screen-based media use and positive and negative health indicators in the U.S. and Canada. J Adolescent Health 2009, 44:493–499.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.