Research Article  Open Access
Xiaobao Yang, Mei Huan, Bingfeng Si, Liang Gao, Hongwei Guo, "Crossing at a Red Light: Behavior of Cyclists at Urban Intersections", Discrete Dynamics in Nature and Society, vol. 2012, Article ID 490810, 12 pages, 2012. https://doi.org/10.1155/2012/490810
Crossing at a Red Light: Behavior of Cyclists at Urban Intersections
Abstract
To investigate the relationship between cyclist violation and waiting duration, the redlight running behavior of nonmotorized vehicles is examined at signalized intersections. Violation waiting duration is collected by video cameras and it is assigned as censored and uncensored data to distinguish between normal crossing and redlight running. A proportional hazardbased duration model is introduced, and variables revealing personal characteristics and traffic conditions are used to describe the effects of internal and external factors. Empirical results show that the redlight running behavior of cyclist is time dependent. Cyclist’s violating behavior represents positive duration dependence, that the longer the waiting time elapsed, the more likely cyclists would end the wait soon. About 32% of cyclists are at high risk of violation and low waiting time to cross the intersections. About 15% of all the cyclists are generally nonrisk takers who can obey the traffic rules after waiting for 95 seconds. The human factors and external environment play an important role in cyclists’ violation behavior. Minimizing the effects of unfavorable condition in traffic planning and designing may be an effective measure to enhance traffic safety.
1. Introduction
Urban traffic problem is widely recognized as one of the main maladies of life in large cities. Many scholars have paid much attention to the traffic problem [1, 2]. A Mix of nonmotorized and motorized vehicles is an important traffic type in China. Some surveys show that the nonmotorized vehicle is one of the most widely used traffic tools in Chinese daily travel activity [3]. At present, cycling still has heavy proportion among all travel modes in China, and in Tianjin, as large as over 60% [4]. Even in the developed countries, bicycle travel is recognized as a lowenergy consumption and as being healthy to the users and it does not damage the health of others. Meanwhile, the electric bike (ebike) has emerged as a popular mode of transportation in many large cities during recent years. Ebike use has rapidly expanded in China, in the process changing the mode split of many cities. Currently, China produces over 20 million ebikes yearly, up from a few thousands a decade ago [5]. Ebikes in China are defined as electric twowheelers with relatively low speeds and weights compared to a motorcycle. Both bicyclestyle ebikes (with functioning pedals) and scooterstyle ebikes (with many features of gasoline scooters) are classified as bicycles and are given access to bicycle infrastructure.
However, the growing popularity of cycling traffic also entails safety concerns as observed in accident statistics. Accident analysis reveals that over 60% of fatal crashes involving cyclists result from violation of traffic rules [6]. One typical type of rule violation behavior is violation behavior at red period. Because of the poor law enforcement and peoples’ low safety awareness, violation behavior at red period is rather prevalent and represents a substantial safety problem in Chinese urban intersections. Especially, electric bicycles with relatively high speed are likely to increase the risk of traffic incident.
Previous research on drivers and pedestrians also points to several variables of interest regarding violation behaviors. Keegan and O’Mahony gave reports about pedestrians’ streetcrossing behavior influenced by travel distance and waiting time [7]. Other researchers paid much attention to the influences of personal features on the streetcrossing behavior [8–10]. Some useful reviews of the existing research on pedestrian streetcrossing behavior in urban areas can be found in Ishaque and Noland [11] and Papadimitriou et al. [12].
Unfortunately, only a few studies have investigated the violation behavior of cyclists. Johnson et al. identified three distinct types of violated cyclists that are exposed to different levels of risk: racers, impatients, and runners [13]. Johnson et al. used video recordings to analyze urban commuter cyclists’ violation behaviors in Melbourne, VIC, Australia [14]. The field observation approach was also used by researchers studying the influence of cycle paths on accident numbers [15].
In this paper, a hazardbased duration approach is adopted to describe the cyclist violation behavior at signalized intersections. The hazardbased duration models have been used extensively in biometrics and reliability engineering for decades [16]. Duration models can be used to determine causality in duration data and they are also useful tools in the field of transportation [17–21]. These models represent a type of analytical methods to describe the duration of a certain state and how various factors have affected the duration. More importantly, duration models can deal with not only uncensored data but censored data. For example, the exact waiting duration reflecting cyclist endurance cannot be observed if cyclists could wait until the permission of traffic rules. However, most statistical models are unable to analyze these uncensored data. Accordingly, cyclists’ waiting times are modeled by a proportional hazardbased duration model. The covariates relevant to traffic conditions and personal features are investigated to capture the influenced factors of cyclist behavior. The results give the time when cyclists are easy to violate traffic rules and the significantly influential factors on waiting behavior. The findings will imply some effective countermeasures for improving the road safety of urban intersections.
2. Method
2.1. Duration Model
The variable of interest in duration model is the survival time that elapsed from the beginning of an event until its end. The waiting time of a cyclist at red light can be regarded as the waiting duration that starts when a cyclist arrives at the intersection at the red period and ends when the cyclist begins to cross the intersection.
Let be a nonnegative variable representing the waiting duration of a cyclist at a signalized intersection. Let denote the probability density function of and the cumulative distribution be Let denote the probability that the waiting duration does not end prior to is called survival function or survivor probability. The survival function is defined to be the probability that the waiting time of a cyclist at a red light is longer than some specific time, .
In the survival analysis, can be characterized by a hazard function, . The hazard function under study is the instantaneous rate at which the waiting duration will end in an infinitesimally small time period, , after time , given that the duration time has lasted to time
Note that the waiting time of a cyclist is influenced by various factors. The influential factors can be defined as a vector of explanatory variables, . To accommodate the effects of these influential factors is a main objective of this paper. Thus the proportional hazard form is introduced, which specifies the effects of explanatory variables to be multiplicative on a hazard function where is called the baseline hazard function and can be interpreted as the hazard function when all covariates are ignored. is a known function to represent the effects of explanatory variables, is a vector of estimable coefficients for . In this paper, a typical specification with , which was proposed by Cox [22], is used. The specification is convenient since it guarantees the positivity of the hazard function without placing constraints on the signs of the elements of . The Cox proportional hazard model is
Combining (2.3) and (2.5), the survival function can be written as where represents the baseline cumulative hazard function. Thus, the covariates can be incorporated into the survival function.
2.2. Model Estimation
The main interest of this paper is to identify from the covariates a subset of variables that affects the hazard more significantly, and consequently, the waiting duration time at a signalized intersection. We are concerned with the regression coefficients. If is zero, the corresponding covariate is not related to the waiting time. If is not zero, it represents the magnitude of the effect of on hazard when the other covariates are considered simultaneously.
To estimate the coefficients, , a partial likelihood method is adopted. Suppose that of the duration times from cyclists is observed and distinct. Let be the ordered distinct duration times with corresponding covariates . Let be the risk set at time . consists of all cyclists whose duration times are at least . For the particular duration time , conditionally on the risk set , the probability is Each distinct duration time contributes a factor and hence the partial likelihood function is and the logpartial likelihood is
The overall goodness of fit of the model estimation is determined by the likelihood ratio (LR) statistics, which is specified as where is the logpartial likelihood for null model with all the regression coefficients are set as zero and is the logpartial likelihood at convergence with regression coefficients. The Cox proportional hazard model has been widely cited in the literature. For the estimation of and other detailed discussion of this model see, Lee and Wang [16] and Bhat [17].
2.3. Covariate Selection
The covariate selection takes into account the previous research [13, 14] and arguments regarding the effects of the exogenous variables on cyclist crossing behavior. The practical effects on waiting behavior and the feasibility of data acquisition are considered in the covariate selection. Two broad sets of variables are considered as covariates: personal characteristics and traffic conditions. Personal characteristics involve age and gender. The selected covariates of traffic conditions can determine the effects on the waiting time and traffic volume. The following covariates, as shown in Table 1, are adopted to construct the duration model.

3. Survey and Data
3.1. Site Survey Design
To record cyclists’ waiting durations, the whole redlight period of a signal cycle was observed as a data collection unit. Only the cyclists who arrived in the redlight period were defined as a valid sample. The waiting duration was from the time a cyclist arrived at the crossing location to the time he/she began to cross. It can be classified into two kinds: uncensored data and censored data. The uncensored data was defined as the waiting duration which ended within the redlight period (violating crossing). Otherwise, the waiting duration was called as the censored data as long as it ended within the greenlight period (normal crossing). For the censored data, it is unknown about the exact waiting duration which can reflect the endurance of waiting time for cyclists.
The site survey was conducted at three selected signalized intersections near Jiaotong University in Beijing, China. Data collections were done by placing video cameras at each location. The survey periods included peak hour (7:30 a.m.–9:30 a.m.) and offpeak hour (10:00 a.m.–4:00 p.m.). The survey area covered the zebra crossing and a part of traffic lanes so that the cyclist crossing behavior and the corresponding traffic conditions can be monitored clearly. Some additional explanations are needed for the site survey. (a)The signals were old traditional person heads so that the influence of type of signals could be neglected [7]. The selected sites had similar characteristics involved geometric, and traffic conditions, traffic control. (b)The survey was conducted in good weather and the absence of pointsmen. Cyclists were unobtrusively observed.
3.2. Descriptive Statistics
Of the 459 valid observations, 295 (64.27%) cyclists violated the traffic regulations. The average waiting time of all samples was 25.16 seconds, with a standard deviation of 27.13 seconds. The average waiting time of the violating crossing was 15.71 seconds while the average waiting time of the normal crossing is 43.14 seconds. The maximum waiting duration was 116 seconds while the minimum was 0 second. The latter means people cross the street without any wait. This descriptive statistic cannot reflect the exact waiting behavior due to the neglect of the censored data. The estimation of the waiting duration with censored data will be discussed later.
4. Empirical Results
The results are discussed in two sections. The overall results are presented in the first section including model fit statistics and survival probability estimation. The second subsection presents the effects of covariates.
4.1. Overall Results
(1) Model fit statistics: the LR statistic of the estimated model clearly indicates the overall goodness of fit (the LR statistic is 3201.0, which is greater than the chisquared statistic with 8 degrees of freedom at any reasonable level of significance). The significant level corresponding to each covariate is given by value in Table 2. From the results, most of the included covariates are statistically significant at the 0.10 level of significance. It means that these covariates are significantly related to violation behavior. Only gender has relative low significant level. It is partly because that the female rate (24.2%) in the sample is relative low. The significance level of each covariate suggests that the importance of covariate should be interpreted carefully.

(2) Survival probability: Figure 1 gives the survival probability calculated by the duration model, which represents the probability of complying with the traffic rules while waiting at the signalized intersections. The survival probability for estimated model presents a general decline trend with elapsed waiting duration. The survival probability can be divided into three parts according to the gradient. Firstly, a sharp decline for the short duration indicates that there are a number of cyclists would violate to cross without any delay. Especially, about 32 percent of cyclists can be defined as risk takers since they show high violation inclination and very low waiting endurance (<3 seconds). Then, the probability decreases smoothly from 3 seconds to 95 seconds. This steady reduction reflects that the number of cyclist violations is increasing continuously. The declining trend of the survival probability indicates that the redlight running behavior of most cyclists is time dependent. It means that cyclists are easy to end waiting duration and violate the traffic rules with the elapsed duration. Note that about half of the observed cyclists cannot endure 29 seconds or longer. Finally, there are 15 percent of cyclists who wait and wait longer, and they are generally nonrisk takers.
(3) Figure 1 also gives the comparison between the estimated survival probability and the observed survival probability. Here, the estimated results are calculated by the Cox proportional hazard model; while the observed results are calculated by the nonparametric approach in which the covariate effects are not considered. The detailed discussion of the nonparametric approach can refer to the work of Lee and Wang [16]. The results show that there are some differences between them though the general shape is the same between the two results. Specifically, compared to the estimated results, the observed survival probability is smaller until about 24.0?s, larger thereafter. This difference is expected to be the covariate effects, at least partly. The observed results indicate the waiting time under the specific condition for individual sample, while the estimated results indicate the waiting time under the average condition for all the samples. The estimated survival probability reflects the characteristics of the waiting time which has an average value for every variable. Any change of the variables could influence the estimated results. The effects of variables are discussed in the next subsection.
4.2. Analysis of Covariate Effects
According to (2.5), the effects of the explanatory variables can be interpreted by the signs of the coefficients in a rather straightforward fashion. If the coefficient is negative, it implies that an increase in the corresponding variable decreases the hazard rate, or equivalently, increases the waiting duration. With regard to the magnitude of the variable effects, when a variable changes by one unit, the hazard would change by .
To assess the effects of the explanatory variables on the duration time, a function of hazard ratio (HR) can be obtained by dividing both sides of (2.5) by , yield where the x’s are covariates for the th cyclist; is a vector of the coefficients which has been estimated by using the Cox proportional hazard model. The left side of (4.1) is a function of hazard ratio (HR) and the right side is a linear function of the covariates and their respective coefficients. The HR can represent the multiple relations between the hazard under the covariate effects and the hazard when all variables are ignored.
If the covariates are standardized about the mean and the model used is where and is the average of the jth covariate for all cyclists, the left side of (4.2) is the logarithm of the relative hazard ratio (RHR). RHR represents the hazard ratio for a cyclist with a given set of values to that for a cyclist which has an average value for every covariate. If RHR is more than one, it means that the covariate effects can increase the hazard and so the variables are favorable. That is to say, the waiting time in such a favorable condition is less than the average level of the survey sample. On the contrary, the unfavorable variable corresponds to a low hazard. Therefore, a cyclist in the unfavorable condition would have longer waiting time than that in the favorable condition.
In order to make a quantitative analysis on the effects of covariates, the relative hazard for each variable is calculated by considering favorable or unfavorable values of that variable, assuming that other variables are at their average value. The favorable or unfavorable values of that variable are given on the basis of the hazard with the value of the variable. The value of the variable with the low hazard is regarded as the favorable condition. Take the age as an example, old people are defined as the favorable condition since old people have lower violation risk than young people. The assumed conditions and corresponding RHRs and HRs are shown in Table 3. The RHRs for three continuous covariates are shown in Figure 2.

(a)
(b)
(c)
The effect of age (AG) indicates that older cyclists have longer waiting time. This is partly because older cyclists have stronger risk consciousness of traffic violations. In addition, older cyclists’ trip purposes are seldom related to work or school so they are not in a hurry. Note that this conclusion is a statistical result; traffic violations involved older cyclists are also common sometimes.
The effect of gender (GEN) indicates that male cyclists have shorter waiting time and higher tendency to disobey the traffic rules. They are 1.38 times more likely than females to have shorter waiting times. Hamed reported that male pedestrians are 2.61 times more likely than females to have shorter waiting times [23], and other qualitatively similar results were obtained by Tiwari et al. [24].
The effect of covariate nonmotor vehicle type (NT) indicates that cyclists of electric bike have shorter waiting time and higher tendency to disobey the traffic rules. They are 1.60 times more likely than humanpowered bicyclists to have shorter waiting times.
The waiting time of cyclists would increase with the bigger number of other cyclists that are waiting for a green, (WN) when arrives (see Figure 2(a)). Otherwise, the waiting time decreases with the bigger number of other cyclists that are crossing against the red light (CN) when arrives (see Figure 2(b)). This is caused by two reasons. First, many people may consider that the more people cross together, the safer they would be. They take it for granted that drivers must yield to a group of people more often than one person. Second, the conformity psychology would work well in dense cyclist environments.
The effect of covariate TC (twice crossing) shows that the cyclists of twice crossing have higher hazard and shorter waiting time. Cyclists who are apt to twice crossing behavior have little or no patience to wait at a red light. They are 2.35 times more likely than onetime crossing cyclists to have shorter waiting times.
The effect of covariate MV (motor vehicle volume) indicates that heavy traffic can increase waiting time or decrease the risk of cyclist violations (see Figure 2(c)).This is because that the larger motor vehicle volume is, the smaller the average time gap between successive cars is.
The characteristics of travel time also have impacts on cyclists’ redlight running behavior. It indicates that cyclists are at high risk level of traffic violation in peak hour. The cyclists who travel in peak hour are 1.44 times more than those who travel in offpeak hour to end waiting duration and cross illegally. In peak hour period, both cyclists and drivers are in a hurry to the destination related to work or school, so the heavy mixed traffic with impatient cyclists and drivers would cause traffic accidents easily.
5. Conclusions
This paper applies a proportional hazardbased duration model to study the cyclist crossing behavior at signalized intersections by using data acquired in Beijing, China. The crossing behavior is examined by modeling the duration between the arrival at survey area and the start to cross the intersection. If cyclists violate the traffic rules to cross the intersection, their waiting times are recorded as uncensored data, while the waiting durations of normal crossing are recorded as censored data.
The paper provides several important insights into the determinants of the regularity and frequency of cyclist crossing behavior, especially the relation between violation behavior and waiting duration. First, the results indicate that the crossing behavior of cyclists is time dependent, as well as the risk of traffic violation. Cyclists’ crossing behavior presents positive duration dependence, which also implies a “snowballing” effect. It means the longer the time has elapsed since the start of the waiting duration, the more likely cyclists will end the wait soon. Such positive duration dependence also indicates that longer waiting time would increase the risk of cyclist violation. Second, some crucial time points deserve our concern: 3 seconds and 95 seconds. The 3 seconds indicate cyclists who are at high risk of violating crossing the street and low waiting time, and they account for 32% of the sample in the study. The duration of 95 seconds reflects the cyclists’ endurance. About 15 percent of all the cyclists can obey the traffic rules after waiting for 95 seconds. These people are generally nonrisk takers. Third, the human factors and the external environment play an important role in redlight running behavior. Various factors in the unfavorable condition could increase the risk of traffic violation, as well as traffic accidents. The effects of covariates can help to modify cyclists’ crossing behavior. Specifically, rational traffic planning and designing should fully consider cyclist behavioral characteristics. More importantly, minimizing the effects of unfavorable condition involved human factors may be an effective measure to obtain conscious cooperation and behavioral changes of cyclists. Finally, it is noted that, for different cities, the model should be estimated by using the specified field data. Additionally, the explanatory variables can be chosen flexibly according to the research aim and the traffic reality.
In terms of the future work, more parameters under different situations should be taken into account. Next, some engineering solutions should be proposed to improve the safe crossing behavior of cyclists in urban traffic environment. In addition, from the viewpoint of cyclist prevention, the interaction between cyclists and motor vehicles could be analyzed based on such crossing behavior. Findings from this paper may partly supplement previous research which helps us in inspiration. It is also hoped that these findings may give better understanding of cyclist behavioral characteristics at signalized intersections and help to plan and design proper facilities for nonmotorized vehicles.
Acknowledgments
This work was supported by the National Basic Research Program of China (Grant no. 2012CB725400), the National Natural Science Foundation of China (Grant no. 70901005, 71131001), and the Fundamental Research Funds for the Central Universities (Grant no. 2012YJS067).
References
 W. Wang, W. Zhang, H. Guo, H. Bubb, and K. Ikeuchi, “A safetybased approaching behavioural model with various driving characteristics,” Transportation Research C, vol. 19, no. 6, pp. 1202–1214, 2011. View at: Publisher Site  Google Scholar
 C. K. Wong and Y. Y. Lee, “Convergence study of minimizing the nonconvex total delay using the lanebased optimization method for signalcontrolled junctions,” Discrete Dynamics in Nature and Society, vol. 2012, Article ID 858731, 14 pages, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 K. J. Krizek, “Estimating the economic benefits of bicycling and bicycle facilities: an interpretive review and proposed methods,” in Proceedings of the TRB 2004 Annual Meeting, 2004, CDROM 042070. View at: Google Scholar
 D. Wang, T. Feng, and C. Liang, “Research on bicycle conversion factors,” Transportation Research A, vol. 42, no. 8, pp. 1129–1139, 2008. View at: Publisher Site  Google Scholar
 J. X. Weinert, C. Ma, X. Yang, and C. R. Cherry, “Electric twowheelers in China: effect on travel behavior, mode shift, and user safety perceptions in a mediumsized city,” Transportation Research Record, no. 2038, pp. 62–68, 2007. View at: Publisher Site  Google Scholar
 CRTASR, China Road Traffic Accidents Statistics Report, Traffic Administration Bureau of China State Security Ministry, Beijing, China, 2004.
 O. Keegan and M. O'Mahony, “Modifying pedestrian behaviour,” Transportation Research A, vol. 37, no. 10, pp. 889–901, 2003. View at: Publisher Site  Google Scholar
 T. Rosenbloom, D. Nemrodov, and H. Barkan, “For heaven's sake follow the rules: pedestrians' behavior in an ultraorthodox and a nonorthodox city,” Transportation Research F, vol. 7, no. 6, pp. 395–404, 2004. View at: Publisher Site  Google Scholar
 K. Pfeffer, H. P. Fagbemi, and S. Stennet, “Adult Pedestrian behavior when accompanying children on the route to school,” Traffic Injury Prevention, vol. 11, no. 2, pp. 188–193, 2010. View at: Publisher Site  Google Scholar
 R. Zhou and W. J. Horrey, “Predicting adolescent pedestrians' behavioral intentions to follow the masses in risky crossing situations,” Transportation Research F, vol. 13, no. 3, pp. 153–163, 2010. View at: Publisher Site  Google Scholar
 M. M. Ishaque and R. B. Noland, “Behavioural issues in pedestrian speed choice and street crossing behaviour: a review,” Transport Reviews, vol. 28, no. 1, pp. 61–85, 2008. View at: Publisher Site  Google Scholar
 E. Papadimitriou, G. Yannis, and J. Golias, “A critical assessment of pedestrian behaviour models,” Transportation Research F, vol. 12, no. 3, pp. 242–255, 2009. View at: Publisher Site  Google Scholar
 M. Johnson, J. Charlton, and J. Oxley, “Cyclists and red lights: a study of behavior of commuter cyclists in Melbourne,” in Proceedings of the Australasian Road Safety Research, Policing and Education Conference, 2008. View at: Publisher Site  Google Scholar
 M. Johnson, S. Newstead, J. Charlton, and J. Oxley, “Riding through red lights: the rate, characteristics and risk factors of noncompliant urban commuter cyclists,” Accident Analysis and Prevention, vol. 43, no. 1, pp. 323–328, 2011. View at: Publisher Site  Google Scholar
 R. O. Phillips, T. Bjørnskau, R. Hagman, and F. Sagberg, “Reduction in carbicycle conflict at a roadcycle path intersection: evidence of road user adaptation?” Transportation Research F, vol. 14, no. 2, pp. 87–95, 2011. View at: Publisher Site  Google Scholar
 E. T. Lee and J. W. Wang, Statistical Methods for Survival Data Analysis, John Wiley & Sons, New York, NY, USA, 2003. View at: Publisher Site
 C. R. Bhat, “Duration modeling,” in Handbook of Transport Modelling, D. A. Hensher and K. J. Button, Eds., pp. 91–111, Elsevier Science, Oxford, UK, 2000. View at: Google Scholar
 C. R. Bhat, T. Frusti, H. Zhao, S. Schönfelder, and K. W. Axhausen, “Intershopping duration: an analysis using multiweek data,” Transportation Research B, vol. 38, no. 1, pp. 39–60, 2004. View at: Publisher Site  Google Scholar
 X. B. Yang, M. Huan, H. W. Guo, and L. Gao, “Car travel time estimation near a bus stop with nonmotorized vehicles,” International Journal of Computational Intelligence Systems, vol. 6, no. 4, pp. 1350–1357, 2011. View at: Google Scholar
 W. Wang, H. Guo, Z. Gao, and H. Bubb, “Individual differences of pedestrian behaviour in midblock crosswalk and intersection,” International Journal of Crashworthiness, vol. 16, no. 1, pp. 1–9, 2011. View at: Publisher Site  Google Scholar
 H. W. Guo, W. H. Wang, W. W. Guo, X. B. Jiang, and H. Bubb, “Reliability analysis of pedestrian safety crossing in urban traffic environment,” Safety Science, vol. 50, no. 4, pp. 968–973, 2012. View at: Publisher Site  Google Scholar
 D. R. Cox, “Regression models and lifetables,” Journal of the Royal Statistical Society B, vol. 34, no. 2, pp. 187–220, 1972. View at: Google Scholar  Zentralblatt MATH
 M. M. Hamed, “Analysis of pedestrians' behavior at pedestrian crossings,” Safety Science, vol. 38, no. 1, pp. 63–82, 2001. View at: Publisher Site  Google Scholar
 G. Tiwari, S. Bangdiwala, A. Saraswat, and S. Gaurav, “Survival analysis: pedestrian risk exposure at signalized intersections,” Transportation Research F, vol. 10, no. 2, pp. 77–89, 2007. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2012 Xiaobao Yang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.