manuscript submission site
eISSN : 2233-6842
Sorry.
You are not permitted to access the full text of articles.
If you have any questions about permissions,
please contact the Society.
죄송합니다.
회원님은 논문 이용 권한이 없습니다.
권한 관련 문의는 학회로 부탁 드립니다.
[ Original articles ] | |
Journal of Exercise Nutrition & Biochemistry - Vol. 21, No. 3, pp.43-49 | |
Abbreviation: J Exerc Nutr Biochem | |
ISSN: 2233-6842 (Online) | |
Print publication date 30 Sep 2017 Online publication date 30 Sep 2017 | |
Received 01 May 2017 Revised 24 Jul 2017 Accepted 28 Aug 2017 | |
DOI: https://doi.org/10.20463/jenb.2017.0016 | |
Development of an anaerobic threshold (HRLT, HRVT) estimation equation using the heart rate threshold (HRT) during the treadmill incremental exercise test | |
Joo-ho Ham^{1} ; Hun-Young Park^{1}^{, 3} ; Youn-ho Kim^{2} ; Sang-kon Bae^{2} ; Byung-hoon Ko^{2} ; Sang-seok Nam^{1}
| |
1.Department of Sports Medicine, Kyung Hee University, Yongin-si, Republic of Korea | |
2.Samsung Advanced Institute of Technology, Suwon-si, Republic of Korea | |
3.Performance Activity and Performance Institute, Konkuk University, Seoul, Republic of Korea | |
Correspondence to : ^{*}Sang-seok Nam Department of Sports Medicine, Kyung Hee University, 1732, Deogyeong-daero, Giheung-gu, Yongin-si, Gyeonggi-do, 17104, Republic of Korea. Tel : +82-31-201-3759 Fax: +82-31-201-3760 Email : captain@khu.ac.kr | |
©2017 The Korean Society for Exercise Nutrition ©2017 Joo-ho Ham et al.; License Journal of Exercise Nutrition and Biochemistry. 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 orginal work is properly cited. | |
Funding Information ▼ |
The purpose of this study was to develop a regression model to estimate the heart rate at the lactate threshold (HR_{LT}) and the heart rate at the ventilatory threshold (HR_{VT}) using the heart rate threshold (HRT), and to test the validity of the regression model.
We performed a graded exercise test with a treadmill in 220 normal individuals (men: 112, women: 108) aged 20–59 years. HRT, HR_{LT}, and HR_{VT} were measured in all subjects. A regression model was developed to estimate HR_{LT} and HR_{VT} using HRT with 70% of the data (men: 79, women: 76) through randomization (7:3), with the Bernoulli trial. The validity of the regression model developed with the remaining 30% of the data (men: 33, women: 32) was also examined.
Based on the regression coefficient, we found that the independent variable HRT was a significant variable in all regression models. The adjusted R^{2} of the developed regression models averaged about 70%, and the standard error of estimation of the validity test results was 11 bpm, which is similar to that of the developed model.
These results suggest that HRT is a useful parameter for predicting HR_{LT} and HR_{VT}.
Keywords: Heart rate threshold, Anaerobic threshold, Lactate threshold, Ventilation threshold, Treadmill, Estimation equation |
The term “customized exercise” is recently gaining popularity. The application of customized exercise means prescribing the frequency, intensity, type, and time of exercise as appropriate for each individual. Many sports scientists argue that the anaerobic threshold (AT) is the most useful measure of personalized exercise intensity^{1, 2, 3}.
AT is the time point at which the rate of energy supply from anaerobic metabolism increases during exercise, and it consists of the lactate threshold (LT) and the ventilatory threshold (VT)^{4, 5}. Therefore, to use AT as a parameter of customized exercise intensity, the LT or VT should be determined. However, LT determination requires blood collection and VT measurement requires expensive equipment and manpower, making them difficult to measure outside of the laboratory. Therefore, some sports scientists have investigated ways to easily predict AT.
Conconi et al.^{6} proposed a method of estimating AT using the heart rate threshold (HRT). HRT is the deflection point of the heart rate increase during incremental exercise, and many scholars have suggested that HRT can be as useful as LT and VT in estimating the AT^{2, 7, 8}. Nam et al.^{2} estimated the AT using HRT during the incremental bicycle ergometer exercise. They reported an adjusted R^{2} of about 72% and emphasized that the estimation method of AT using HRT is practical. However, few studies have suggested a linear relationship between AT and HRT, as the study by Nam et al.^{2}.
Therefore, the purpose of this study was to develop a regression equation that estimates AT using HRT for the treadmill maximal exercise load test.
This study included 220 healthy male and female subjects aged 20–59 years. We confirmed the health of the subjects through the physical activity readiness questionnaire (PAR-Q)^{9} and blood pressure measurement. All subjects were informed about the purpose and process of the study, and voluntarily agreed to participate in the study. All procedures followed were in accordance with the ethical standards of the responsible committee on human experimentation and with the Helsinki Declaration.
Male (n = 112) and female (n = 108) subjects were divided in a ratio of 7 : 3 using the Bernoulli trial. Approximately 70% of the divided data (men: 79, women: 76) were used in the development of the AT estimation formula with HRT, and about 30% of the data (men: 33, women: 32) were used for the validity test. The characteristics of the subjects are presented in Table 1.
Sex | n | Age (years) | Height (cm) | Weight (kg) | BMI (kg/m_{2}) | |
---|---|---|---|---|---|---|
Regression model | Male | 79 | 37.5 ± 12.4 | 173.5 ± 6.5 | 72.3 ± 10.1 | 24.0 ± 2.6 |
Female | 76 | 37.0 ± 13.0 | 161.2 ± 5.5 | 59.0 ± 7.1 | 22.7 ± 2.6 | |
Validity data | Male | 33 | 35.8 ± 10.8 | 173.5 ± 7.4 | 72.5 ± 10.3 | 24.0 ± 2.5 |
Female | 32 | 37.3 ± 11.6 | 161.5 ± 5.5 | 60.1 ± 8.2 | 23.1 ± 3.1 |
We instructed all subjects to abstain from drinking and to get adequate sleep from 3 days before the measurement. Subjects who arrived in the laboratory rested for 30 min before starting the test.
A treadmill was used in the maximal exercise test, and we determined the HRT, LT, and VT after the test. We used the heart rate at the determined LT or VT time point as the dependent variable, and HRT as the independent variable.
We used about 70% of the data of all subjects for the development of the regression equation and about 30% of the data of all subjects for the validity test of the developed regression equation.
For all tests, the laboratory temperature was 21 ± 2°C and the humidity was 50 ± 5%.
For the test, subjects wore the mask of a respiratory gas analyzer (Sensor Medics Vmax229, Milan, Italy) on the face and a heart rate monitor (Polar 610i, Kempele, Finland) on the chest area. We measured the stability parameters of the subjects for 2 min in a sitting position and performed an incremental maximal exercise test using a treadmill (TaeHa TH-6000, Seoul, Korea) until exhaustion.
The protocol of the incremental maximal exercise test was a modification of the Bruce protocol, with the aim to achieve two purposes. First, the protocol should have a small incremental loading between stages so that the responses of the body during exercise can be closely observed. Second, the protocol should be available for use in outdoor fields and indoor tracks. We performed a pilot test >10 times to modify the protocol. These pilot tests included 1–10 subjects aged 20–50 years. The modified Bruce protocol was started at 3.6 km/h, the speed was increased by 1.2 km/h every 2 min, and the slope was 0% at all speeds (Table 2).
Time (min) | Speed (km/h) | Grade (%) |
---|---|---|
0 – 2 | 3.6 | 0 |
2 – 4 | 4.8 | |
4 – 6 | 6.0 | |
6 – 8 | 7.2 | |
8 – 10 | 8.4 | |
10 – 12 | 9.6 | |
12 – 14 | 10.8 | |
14 – 16 | 12.0 | |
16 – 18 | 13.2 | |
18 – 20 | 14.4 | |
20 – 22 | 15.6 | |
22 – 24 | 16.8 | |
24 – 26 | 18.0 | |
26 – 28 | 19.2 |
We finished the exercise test when three or more of the following occurred: (i) when the subject requested to stop the test; (ii) when the current heart rate exceeded 90% of the predicted maximal heart rate; (iii) when the respiratory exchange rate was >1.15; (iv) when the oxygen intake did not change even when the exercise intensity increased^{9}.
Oxygen uptake (VO_{2}), carbon dioxide production (VCO_{2}), respiratory exchange rate, and ventilation (VE) were measured in real time during exercise. The variables were measured every 1 min and the mean value was obtained. The blood lactate concentration and heart rate were measured as described below.
Analysis of blood lactate concentration during exercise was performed every 1 min from the start to the end of the exercise test. For the analysis, blood samples were collected from the subjects’ index finger using a heparinized capillary tube through the “finger-tip” method. The blood samples were analyzed using an automatic blood lactate analyzer (YSI-1500; Yellowspring Instruments, Yellowspring, OH, USA). We determined the LT using the result of blood lactate analysis.
Measurement of heart rate during exercise was performed every 5 s from the start to the end of the exercise test, using a wireless heart rate monitor. The variables were measured every 1 min and the mean value was obtained. We used these data to determine the maximal heart rate (HRmax) and HRT.
The data of blood lactate concentration and respiratory gas during the incremental maximal exercise test were used to determine LT and VT.
LT was determined to be the point at which the blood lactate concentration increased sharply from the resting point^{3}.
VT was determined according to certain criteria: (i) the point at which VE/ VCO_{2} does not increase but VE/VO_{2} increases; (ii) the point at which VE and VCO_{2} increase abruptly; (iii) the point at which VE increases relative to VO_{2}; and (iv) the point at which VCO_{2} increases relative to VO_{2}^{3}.
After LT and VT were determined, we defined the heart rate at the time point of LT and VT as HR_{LT} and HR_{VT}, respectively.
HRT during the incremental maximal exercise test was determined using the maximal distance (Dmax) method^{2, 10}.
The procedure for determining the HRT with the Dmax method follows a certain protocol: (i) Draw a hypothetical straight line connecting both ends of the heart rate trajectory. (ii) Measure the distance from the hypothetical straight line to each heart rate point. (iii) The inflection point is determined as the point having the largest measured distance. (iv) The inflection point is named the HRT. It has been reported that the determination rate of the inflection point with the Dmax method is very high^{10}.
However, in this study, HRT was defined as the shortest distance from the virtual straight line. This is because of two reasons. First, according to some studies on HRT^{10}, the blood lactate concentration was 3–4 mmoL/l^{11, 12}, which corresponds to the onset of blood lactate accumulation. Therefore, Nam et al.^{2} pointed out that early detection of HRT and continuation of exercise were difficult. Second, Hofmann et al.^{11} reported that the HR at lactate turn point 2 (2–4 mmoL/L) was largely different among subjects, whereas the HR at the lactate turn point 1 (0–2 mmoL/L) was not.
The HRT was defined in this study as the first inflection point after setting the point at which the heart rate continuously increases (Figure 1).
The descriptive statistics for each variable were calculated. Different statistical methods were applied according to the purpose of analysis. Statistical analysis was conducted using SPSS 24.0 for Windows and the Excel 2016 program. The level of significance was set at <5%.
In this study, the detected HRT showed a distinct linearity with the dependent variable (HR_{LT} or HR_{VT}). Therefore, we included only HRT as an independent variable considering the simplicity of the regression model.
We used the Bernoulli trial to divide the total subjects in a ratio of about 7:3 (men: 79, women: 76), and developed a regression model for each male and female subject using about 70% of the data. In the development of the regression model, the independent variable was HRT and the dependent variable was HR_{LT} or HR_{VT}. In addition, we rigorously conformed to the basic assumptions of the regression model (linearity, independency, continuity, normality, homoscedasticity, autocorrelation, and outlier).
The validity of the regression model was tested using approximately 30% of the data (men: 33, women: 32) of the total subjects who had already been divided through the Bernoulli trial and were not included in the development of the regression model. We calculated the predicted value of HRLT or HRVT by substituting the data for the validity test into the regression equation. Finally, we calculated the residuals, the mean error (ME, %), and the standard error of estimation (SEE, bpm), using the difference between the predicted value and the actual measured value of HRLT or HRVT, according to the following formulas^{2}:
ME: mean error (%)
HRreal: actual measured value of HR_{LT} or HR_{VT}
HRpred.: predicted value of H_{RLT} or HR_{VT}
Formula 1. The calculation formula for the ME rate for the validity test.
SEE: standard error of estimation (bpm)
HRreal: actual measured value of HR_{LT} or HR_{VT}
HRpred.: predicted value of HR_{LT} or HR_{VT}
Formula 2. The calculation formula of the SEE for the validity test.
The outlier data in the regression model were identified when the absolute value of the standardized residual was 3 or more. In this study, one outlier datum was deleted from the HR_{LT} prediction model of female subjects.
We analyzed the maximal exercise capacity of the subjects. The maximal oxygen uptake of male subjects was about 50 mL·kg^{-1}·min^{-1}, the maximal oxygen uptake of female subjects was about 40 mL·kg^{-1}·min^{-1}, the maximal heart rate of male subjects was 188 bpm, and the maximal heart rate of female subjects was 185 bpm (Table 3).
Items | Male model | Female model | ||
---|---|---|---|---|
HR_{LT} | HR_{VT} | HR_{LT} | HR_{VT} | |
VO_{2}max (mL·kg^{-1}·min^{-1}) | 49.73 ± 7.04 | 49.73 ± 7.04 | 39.75 ± 6.67 | 39.68 ± 6.68 |
VO_{2} (mL·kg^{-1}·min^{-1}) | 29.41 ± 8.70 | 31.08 ± 8.47 | 23.52 ± 6.59 | 24.02 ± 7.06 |
%VO_{2}max (%) | 58.72 ± 13.13 | 61.98 ± 11.94 | 58.82 ± 10.72 | 60.05 ± 11.33 |
HRmax (bpm) | 188.28 ± 15.04 | 188.28 ± 15.04 | 185.11 ± 11.88 | 184.96 ± 11.89 |
HR (bpm) | 134.77 ± 19.36 | 138.73 ± 19.97 | 135.56 ± 21.33 | 135.84 ± 21.13 |
%HR (bpm) | 71.71 ± 9.28 | 73.73 ± 8.92 | 73.06 ± 9.18 | 73.30 ± 9.37 |
We also analyzed the relative ratios of AT (LT or VT) to the maximal athletic performance. The oxygen consumption at AT was 59–62% of the maximal oxygen uptake, and the heart rate at AT was 72–74% of the maximal heart rate (Table 3).
We analyzed the relative ratio of HRT to the maximal heart rate. The HRT of male subjects was 68% of the maximal heart rate and the HRT of female subjects was 67% of the maximal heart rate (Table 4).
Regression model | HRT (bpm) | %HRmax (%) | |
---|---|---|---|
Male model | HR_{LT} | 127.2 ± 18.7 | 67.6 ± 9.0 |
HR_{VT} | 127.2 ± 18.7 | 67.6 ± 9.0 | |
Female model | HR_{LT} | 124.2 ± 19.2 | 67.0 ± 8.7 |
HR_{VT} | 124.4 ± 19.2 | 67.2 ± 8.7 |
We tested the significance of the model using the F-test for each regression model developed. As a result, we confirmed that all regression models were statistically significant (Table 5).
Regression | F-value | Sig. | t-value | Sig. | |
---|---|---|---|---|---|
Male model | HR_{LT} | 175.431 | 0.000 | 13.245 | 0.000 |
HR_{VT} | 201.299 | 0.000 | 14.188 | 0.000 | |
Female model | HR_{LT} | 183.351 | 0.000 | 13.541 | 0.000 |
HR_{VT} | 154.513 | 0.000 | 12.430 | 0.000 |
We used the t-test for each regression model to test the significance of the regression coefficients of the independent variable (HRT). As a result, we found that the regression coefficient of the independent variable was statistically significant in all regression models (Table 5).
For each regression model developed, we calculated the adjusted coefficient of determination (adjusted R^{2}, adj R^{2}) and the SEE. The mean explanatory power of the regression models was about 69.9%, and the mean SEE was about 11.23 bpm (Table 6).
Regression model | F-value | Sig. | t-value | Sig. | |
---|---|---|---|---|---|
Male model | HR_{LT} | 0.691 | 10.76 | HRLT = 0.861 × HRT + 25.309 | 0.000 |
HR_{VT} | 0.720 | 10.57 | HRVT = 0.906 × HRT + 23.544 | 0.000 | |
Female model | HR_{LT} | 0.709 | 11.52 | HRLT = 0.936 × HRT + 19.311 | 0.000 |
HR_{VT} | 0.675 | 12.05 | HRVT = 0.902 × HRT + 23.583 | 0.000 | |
Mean | 0.699 | 11.23 | - |
The final regression equations are presented in (Table 6).
We assessed the validity of the regression model developed using the unused data (30% of the total) in the regression model development. In all regression models, ME was within ±3% for both male and female subjects, and SEE was similar to the regression model developed (about 11 bpm) (Table 7).
Regression model | HRT (bpm) | %HRmax (%) | |
---|---|---|---|
Male model | HR_{LT} | 1.86 | 11.34 |
HR_{VT} | 0.46 | 10.47 | |
Female model | HR_{LT} | -2.61 | 10.88 |
HR_{VT} | -1.64 | 10.72 | |
Mean | -0.48 | 10.85 |
In linear regression analysis, it is recommended that outliers are removed because they increase the error of the estimates. Generally, outlier determinations in regression analysis use the absolutely value of standardized residual. In this study, only one outlier was observed in the HR_{VT} prediction model of women. This result implies that the linearity between the independent variable and the dependent variable is clear. Moreover, it is consistent with the previous studys’ claim that there is a relationship between HRT and AT^{6, 8}
Nam et al.^{2} performed bicycle ergometer tests on subjects similar to those of the present study, and reported that the VO_{2}max and HRmax of male subjects were 43 mL·kg^{-1}·min^{-1} and 182 bpm, respectively, and those of female subjects were 34 mL·kg^{-1}·min^{-1} and 177 bpm, respectively. However, in this study in which a treadmill was used, the VO_{2}max and HRmax of male subjects were 50 mL·kg^{-1}·min^{-1} and 188 bpm, respectively, and those of female subjects were 40 mL·kg^{-1}·min^{-1} and 185 bpm, respectively. The maximal exercise capacity measured in this study was higher than that of Nam et al.'s study^{2}, and we considered that this was due to differences in the types of exercise (treadmill vs. bike). In other words, treadmill running involves more wholebody motion than bicycle exercise; thus, the maximal exercise capacity was higher in treadmill running than in bicycle exercise. However, many studies claimed^{13, 14} that maximal exercise capacity is not a good indicator of the fitness level. Therefore, it is necessary to develop a reasonable athletic index.
In this study, VO_{2 }at AT was slightly higher than reported by Nam et al.^{2}, and heart rate at AT was similar to that of Nam et al.’s study^{2}. The results of these comparisons show that VO at AT is high when VO_{2}max is high, and heart rate at AT remains the same when HRmax is high. This suggests that AT is a reasonable indicator of exercise intensity^{15} but results may vary depending on the type of exercise. Therefore, AT should be used with caution. If these problems are fundamental because of the type of exercise, we could not explain these problems solely with the results of this study. Therefore, further study is needed on the timing of AT expression according to the type of exercise (or type of test). Kang^{16} measured the AT using an arm ergometer and a bicycle ergometer, respectively. They found a significant difference between the test types in both LT and VT. In addition, Kang^{16} explained that this difference was due to the combination of various factors such as skeletal muscle mobilization, distance between the muscle and heart, isometric kinetic factor, muscle fiber component, motor unit mobilization type, total capillary cross-sectional area, concentration of aerobic enzyme, and mechanical pressure of blood vessels. These explanations suggest that there may be differences in the timing of AT expression according to the type of exercise.
On the other hand, we also considered the possibility that this difference occurred in the measurement process of the AT. However, as the measurement method of LT and VT^{3} is traditionally clear, the possibility of error in the measurement method is not high.
Compared with a previous study^{2}, the relative level of HRT was similar but the absolute level was somewhat higher in this study. Further research is needed to determine if the results are due to differences in exercise types, as described above.
Significance of the independent variable (HRT)
The regression models of HR_{LT} and HR_{VT} developed in this study were all statistically significant. This is because we strictly adhered to the statistical baseline assumptions of the regression analysis. Moreover, we found that HRT was a statistically significant independent variable in all regression models. These results provide evidence that HRT is a significant predictor of AT (HR_{LT} or HR_{VT}).
The average adjusted R^{2} of the regression models developed in this study was about 70%. Considering that the adjusted R^{2} has a lower value than R^{2}, it is important that the adjusted R^{2} of this study is about 70%. The reason we chose adjusted R^{2} in this study is to avoid the problem of increasing R^{2} with a larger sample size. The adjusted R^{2} of 70% emphasizes that HRT is useful as a predictor of HR_{LT} and HR_{VT}.
The results of this study demonstrate that HRT can be useful in estimating AT during treadmill exercise in male and female subjects aged 20–50 years.
In the present study, the validity of the regression model was analyzed using the data of subjects who were not used in the development of the regression model. The mean error rate was within ±3%. This result is similar to that of Nam et al.’s study^{2} with similar subjects. The SEE of validity test was about 11 bpm, which is similar to the error of the regression model developed in this study (about 11 bpm). Therefore, we confirmed that the regression model developed in this study has statistical significance, good performance, and validity.
The results of this study confirm that HRT is a statistically significant variable in all regression models. Moreover, the adjusted R^{2} of the developed regression models is about 70% on average. In addition, the validity test showed that the SEE was about 11 bpm, which is similar to that of the regression model developed in this study. In summary, HRT is a useful parameter for predicting HR_{LT} and HR_{VT}.
This study was supported by the Samsung Advanced Institute of Technology.
1. | Krista, S, & Brian, RM, Anaerobic threshold: The concept and methods of measurement, Can J Appl Physiol, (2003), 28, p299-323. |
2. | Nam, SS, Ham, JH, Park, HY, Bae, SK, Ko, BH, Kim, YH, Sunoo, S, Development of anaerobic threshold (HRLT, HRVT) estimation equation using heart rate threshold (HRT) during maximal bike exercise test, J Korean Soci Sports Sci, (2015), 24, p1165-77. |
3. | Wasserman, K, Whipp, BJ, Koyal, SN, & Beaver, W, Anaerobic threshold and respiratory gas exchange during exercise, J Appl Physiol, (1973), 35, p236-43. |
4. | Wasserman, K, Determinants and detection of anaerobic threshold and consequences of exercise above it, Circulation, (1987), 76, p29-39. |
5. | Wasserman, K, The anaerobic threshold measurement in exercise testing, Clin Chest Med, (1984), 5, p77-88. |
6. | Conconi, F, Ferrari, M, Ziglio, PG, Droghetti, P, & Codeca, L, Determination of the anaerobic threshold by noninvasive field test in runners, J Appl Physiol, (1982), 52, p869-73. |
7. | Kim, KJ, Determination of anaerobic threshold by the kinetics of heart rate during graded rowing exercise, Korean J Physical Edu, (1996), 35, p349-57. |
8. | Ko, SS, Lee, DS, Chung, JS, Relationship between Heart Rate Threshold and Ventilatory Threshold for Determination of Anaerobic Threshold in College Students, Korean J Growth Develop, (2013), 21, p307-313. |
9. | American College of Sports Medicine, ACSM’s guidelines for exercise testing and prescription, ACSM, (2014), p24. |
10. | Kara, M, Gokbel, H, Bediz, C, Ergene, N, Ucok, K, & Yysal, H, Determination of the heart rate deflection point by the Dmax method, J. Sports Med. Phys. Fitness, (1996), 36, p31-4. |
11. | Hofmann, P, Pokan, R, Preidler, K, Leitner, H, Szolar, D, Eber, B, & Schwaberger, G, Relationship between heart rate threshold, lactate turn point and myocardial function, Int J Sports Med, (1994), 15, p232-37. |
12. | Ko, SS, Validity of the Heart Rate Threshold by Dmax Method as the Lactate Threshold, JKSSPE, (2007), 12, p149-56. |
13. | Bruce, RA, Blackmon, JR, Jones, JW, & Strait, G, Exercising testing in adult normal subjects and cardiac patients, Pediatrics, (1963), 32, p742-56. |
14. | Costill, DL, Thomason, H, & Roberts, E, Fractional utilization of the aerobic capacity during distance running, Med Sci Sports Exerc, (1973), 5, p248-52. |
15. | Davis, JA, Anaerobic threshold: review of the concept and directions for future research, Med Sci Sports Exerc, (1985), 17, p6-21. |
16. | Kang, KS, Analysis of anaerobic threshold and plasma catecholamine during arm ergometer and bicycle ergometer exercise, Graduate School of Education in Hanman University, (2001). |