Estimation of Large Truck Volume Using Single Loop Detector Data

Size: px
Start display at page:

Download "Estimation of Large Truck Volume Using Single Loop Detector Data"

Transcription

1 Estmaton of Large Truck Volume Usng Sngle Loop Detector Data Yunlong Zhang*, Ph.D., P.E. Assstant Professor Department of Cvl Engneerng Texas A&M Unversty 3136 TAMU College Staton TX, Phone: (979) Emal: Zhru Ye Department of Cvl Engneerng Texas A&M Unversty 3136 TAMU College Staton, TX Phone (979) Emal: Yuanchang Xe Department of Cvl Engneerng Texas A&M Unversty 3116 TAMU College Staton, TX Phone (979) Emal: July, 2007 *Correspondng Author Word Count: (10 Tables and Fgures)*250 = 6676

2 Estmaton of Large Truck Volume Usng Sngle Loop Detector Data Yunlong Zhang, Zhru Ye, and Yuanchang Xe ABSTRACT Ths paper presents a methodology to estmate the volume of large trucks by usng sngle loop outputs. An Unscented Kalman Flterng (UKF) algorthm s used to generate speed estmaton based on sngle loop data, and the mean effectve vehcle length of each tme nterval s thus obtaned wth the loop outputs and speed data provded. We nvestgated the bmodal vehcle length dstrbuton and further separated t nto short and long vehcle dstrbutons. Based on the estmated mean effectve vehcle lengths and vehcle length dstrbutons, an algorthm s proposed to estmate large truck volume. Data collected from Texas Transportaton Insttute s vehcle detecton test beds are used for ths study. The results show that the proposed method produces good performances n large truck volume estmaton under varous traffc condtons. Key Words: Trucks; Volume Estmaton; Sngle Loop Detectors 2

3 INTRODUCTION Truck traffc has been ncreasng steadly on Amerca s hghways due to the growng economy and ncreased nternatonal trade. Trucks have dfferent characterstcs from passenger cars and thus are a key consderaton n many stages of a hghway project. Trucks are crtcal to transportaton plannng, pavement desgn, geometrc desgn, and traffc operatons. Wth the rapd growth of truck volumes on hghways, safety wth respect to trucks has also receved more and more attenton. Vehcle classfcaton nformaton can be obtaned from vehcle detecton equpment wth a specfc classfcaton mechansm. Vsual, axle, and presence sensors are three common types of detectors used for truck dentfcaton. The presence sensors, such as loop detectors, are the most wdely used detectors n the U.S. due to ther cost-effectveness. Dual-loop (double loop) detectors can measure vehcle length nformaton and act as a crude vehcle classfer (1). For sngle loop detectors, however, the vehcle length and classfcaton nformaton s not avalable n the output. Snce sngle loops are a very common type of detector n the feld that provde onlne traffc data, an effectve estmaton method for vehcle length and truck traffc based on sngle loop output would be very useful for real-tme traffc management and control Research n the estmaton of large truck volume from sngle loop outputs has attracted more and more attenton n recent years. A study n the state of Washngton (2) proposed a pattern dscrmnaton method to estmate large truck volume usng sngle loop data. One assumpton used n that study s that vehcle speeds n each fve-mnute tme perod are 3

4 consstent and speed varance can be gnored. The developed method s hence sutable for stable traffc condtons wthout large speed varatons. Another study n Calforna ( 3 ) used lane-to-lane speed correlaton to estmate large truck volume usng data from sngle loop detectors. They assumed that the percentage of large trucks n the leftmost lane can be gnored, and vehcle speeds over dfferent lanes tend to be synchronzed. However, due to wde ranges of traffc and geometrc condtons, varatons of speeds and truck lane dstrbuton can be very sgnfcant. The assumptons made n the pror studes may not be vald for some types or locatons of hghways or for some traffc condtons. Ths drawback lmts the applcaton scope and may affect the effectveness of the approach. Most recently, an Artfcal Neural Network (ANN) was proposed for vehcle classfcaton usng sngle loop outputs (4). Four length-based vehcle classfcaton categores were used n ths study. The results showed that ths method performs well for two categores wth the shortest and the largest vehcle lengths. The performance of an ANN s drectly related to the tranng dataset that s used, and overpredcton based on lmted tranng data s often assocated wth ANNs. Ths paper provdes a new methodology to estmate the volume of large trucks from sngle loop detector outputs. In ths study, a Large Truck (LT) or a long vehcle s defned as a truck wth a length greater than or equal to m (40 feet) as n the studes of Wang and Nhan (2) and Kwon et al. (3). An accurate sngle loop speed estmaton algorthm usng an Unscented Kalman Flter (UKF) s employed to generate speed data, and the Mean Effectve Vehcle Length (MEVL), whch s equal to the sum of average vehcle length and sngle loop detector length, can then be calculated at each tme nterval. The paper then nvestgates the overall 4

5 vehcle length dstrbuton as well as short and long vehcle length dstrbutons. An algorthm s then proposed to estmate LT volumes based on the MEVL and vehcle length dstrbutons. Fnally, estmaton results are presented and evaluated wth data collected from two vehcle detecton test beds. DATA FOR THE STUDY The data we used are collected from two freeway locatons, State Hghway 6 (SH6) near the FM 60 (Unversty Drve) Brdge n College Staton, Texas, and Interstate Hghway 35 (IH-35) near the 47 th street ext n Austn, Texas. Texas Transportaton Insttute (TTI) has two test beds at those locatons. SH6 has two lanes n each drecton whle IH-35 has four at the test bed locaton. Data are collected usng Peek ADR 6000 detectors, each of whch ncludes two nductve loops embedded n the pavement. Ths relatvely new ADR 6000 detector adopts the nductve loop technology and generates accurate outputs of vehcle presence tme, speed, count, vehcle length, and vehcle classfcaton (5). Count and occupancy data, whch are typcal outputs of sngle loop detectors, can be extracted from vehcle presence tmes (the duraton a vehcle passng a sngle loop) wth a certan tme nterval for speed and vehcle length estmaton purposes. Vehcle classfcaton, vehcle length, and speed data are also avalable for every tme nterval from the detector outputs and are used to compare wth estmaton results to assess the performances of the estmaton algorthms. 5

6 Traffc volume on SH6 s moderate, wth daly volume around 13,000 vehcles for the rght lane, and 10,000 vehcles for the left lane. Most vehcles were travelng at hgh speeds at ths locaton. Traffc volume s hgh on IH-35, wth an average daly volume over 27,000 vehcles for the rghtmost lane, and about 15,000 vehcles for the leftmost lane. Congeston repeatedly occurred durng mornng and afternoon peak hours. Durng those congested perods, vehcles traveled at very low speeds and many experenced stop and go stuatons. The rghtmost lanes of both locatons are selected for ths study because they experenced the worst traffc condtons and contaned the hghest LT percentages. Data were collected for one weekday from each test bed (January 27, 2004 on northbound SH6 and October 27, 2004 on southbound IH-35) and compled nto 30-second tme ntervals, whch s a common type of tme nterval used by sngle loop detectors. Large Truck Classfcatons As has been mentoned before, whether a vehcle s a LT or not s determned by ts length n ths paper. Thus, n the LT volume estmaton usng sngle loop outputs, vehcles are smply classfed nto two types: short vehcles (less than m) and large trucks (greater than or equal to m). Ths way of classfcaton s dfferent from the Vehcle Classfcaton Scheme F (6), whch classfes vehcles nto 15 types wth 13 types explctly defned based on the number of axles, axle spacng, and other features of vehcles. In the Scheme F that are used by most state DOTs, classes 5 to 13 are classfed as trucks. Trucks of class 5 (sngle-unt trucks) have two axles and those of classes 6 to 13 have three axles or more. 6

7 Wth data collected from Peek systems, we are able to dscover the dfference between LTs and trucks. We conducted a comparson of LTs (based on length) and trucks wth three axles and up, usng two-day data collected from the test beds. Results of the comparson are shown n Table 1. It s found that the percentages of LTs based on vehcle length are only slghtly less than those of trucks wth three or more axles based on the Scheme F. Therefore, based on our data, the LT classfcaton based on vehcle length n ths study s reasonable and the LT results can be used to represent all trucks excludng sngle unt trucks from the Scheme F classfcaton method. METHODOLOGY In order to determne LT volume from sngle loop detector data, the most crtcal parameter to be determned n ths study s the MEVL n a tme nterval. Snce both vehcle lengths and speeds are unknown, one way to obtan the MEVL s to employ some algorthms to frst estmate the average speed at each tme nterval. Speed Estmaton Usng the UKF For tme nterval wth duratont, assume N vehcles pass over a loop detector, vehcle j has speed s j, and the effectve vehcle length (EVL), defned as the vehcle length plus sngle loop length, s L j. The occupancy s then defned by 7

8 O 1 T N = j=1 L s j j + ε (1) where ε s a zero-mean nose of occupancy. Daley (7) has statstcally transformed ths equaton to an aggregated level. Assume perfect measurement ( ε = 0 ) for a loop detector, Equaton 1 s transformed to N O T s = L 1 / s [ 2 2 σ s ] + 1 (2) where s s the mean speed, σ 2 s s the speed varance and L s the MEVL at th tme nterval. A number of research efforts have been conducted to estmate speeds usng sngle loop outputs (8; 9; 10). However, most prevous studes of speed estmaton cannot produce very accurate estmates under varous traffc condtons due to varous smplfcatons and lnear assumptons. As descrbed n Equaton 2, the problem of speed estmaton s a nonlnear system. An Extended Kalman Flterng (EKF) technque was proposed to solve ths problem through lnearzaton (7). However, ths flterng technque has ts weaknesses n dealng wth nonlnear systems (11, 12). Partcularly, ts nablty to make good predctons of the system state, observatons, and assocated covarance matrces when the system and/or observaton models are nonlnear affects predcton accuracy.. Recently, Ye et al. (13) proposed a UKF method for speed estmaton. Ths method s able to generate hghly accurate speed estmates under both normal and congested traffc condtons. The speed estmaton s not the man focus of ths paper and t can be also 8

9 ndependent from the LT estmaton methodology proposed n ths paper. However, to make ths paper self-suffcent, the UKF s brefly dscussed here. The fundamental part of the UKF s the Unscented Transformaton (UT) (12). It uses a determnstc algorthm to sample a set of weghted sgma ponts, whch are appled to parameterze the means and covarances of probablty dstrbutons. For a d-dmensonal random varable x wth mean x and covarance P xx, ( 2d + 1) weghted sgma ponts are drawn to approxmate the random varable. The ( 2d + 1) sgma ponts { χ,..., χ } are determned by 0 2 d χ = x, W = κ/(d+ κ) 0 0 χ = x + ( ( d + κ) P ), W = 1/[2( d + κ)], = 1,..., d xx where χ = x ( ( d + κ) P ), W = 1/[2( d + κ)], = 1,..., d + d xx +d κ R provdes an extra degree of freedom to fne-tune the hgher order moments of the approxmaton, P xx s the covarance of random varable x, ( ( d + κ) P xx ) s the th row or column of the matrx root of ( d + κ P, and W s the weght assocated wth the th sgma ) xx pont. After sgma ponts are selected, the UKF wll undergo two steps at each tme nterval: the tme update step and the measurement update step. The tme update step predcts the state and covarance, and the measurement step corrects the state and covarance wth consderaton of the most recent observaton. The operaton process of the UKF s shown n Fgure 1. The calculated sgma ponts are used n both tme update and measurement update steps to approxmated mean, covarance and cross-varance of both state varable s and measurement varable N /O. For more detals about the UKF algorthm, refer to (13). 9

10 Estmaton of the MEVL As the mean speed s s estmated usng the UKF at each tme nterval, the MEVL for the tme nterval ( L )s obtaned by rearrangng Equaton 2: O T s L = N 1 / s [ 2 2 σ s ] + 1 (3) O and N are occupancy and count for the th tme nterval and are drectly from the sngle loop output. Once the MEVL s determned, t wll be used for further LT volume estmaton. Estmaton of LT Volume Vehcle Length Dstrbuton In the study by Wang and Nhan (2), a large sample of vehcle lengths was collected from dual-loop detectors and t was found that there were two peaks n the dstrbuton of vehcle length. After separatng vehcle lengths less than m (40ft) from the lengths greater than m, two vehcle length dstrbutons were obtaned and each ft a normal dstrbuton very well. In ths paper, a sample of vehcle lengths from IH-35 was chosen for the study of vehcle length dstrbuton. The Kernel densty of vehcle lengths s shown n the part a) of Fgure 2. The densty plot dsplays a bmodal dstrbuton of vehcle lengths. Small vehcle lengths are hghly concentrated around 5 m and long vehcle lengths are dstrbuted around 20 m. In order to nvestgate the dstrbutons of small and long vehcles, the sample was separated nto two groups, usng m as the separaton pont. To examne f both samples ft normal dstrbutons, 10

11 normal probablty plots were utlzed. Parts b) and c) of Fgure 2 show the Q-Q plots of short and long vehcle lengths. A Q-Q plot can be used to determne graphcally f a dataset comes from a normal dstrbuton. Obvously, the dstrbuton of short vehcle lengths s hghly rght skewed and the dstrbuton of long vehcle lengths s left skewed. The dstrbutons of both short and long vehcle lengths are not normally dstrbuted. The methodology of LT volume estmaton n ths paper does not hnge upon the normal dstrbuton assumptons. No matter what knd of dstrbuton vehcle lengths ft for, we assume that the mean values of short and long vehcle lengths are u 1 and u 2. It s also assumed that the short vehcle length dstrbuton X and the long vehcle length dstrbuton Y are ndependent. Then the average vehcle length at each tme nterval wll be subjected to a dstrbuton wth the expectaton u = p1 * u1 + p2 * u2, where p 1 and p 2 denote the proportons of short and long vehcles. At th tme nterval, assume that N vehcle are detected, and suppose that n 1 vehcles are short vehcles and n 2 vehcles are LTs. Then, X 1, X 2,..., X n form a random 1 sample from the dstrbuton X, and Y 1, Y2,..., Yn form the other random sample from the 2 dstrbuton Y. To compute the average vehcle length ( τ ), we smply take the expectaton X + X + X + Y + Y + + Y τ = E N n1 u1 + n2 u2 = E N n1 n2 = u1 + u2 N N 1 2 n1 1 2 n2 where n 1 +n 2 =N, and the average vehcle length ( τ ) equals to the MEVL mnus sngle loop (4) length. 11

12 Two samples of hstorcal vehcle length data from IH-35 and SH6 are gathered to generate vehcle length dstrbutons of both short vehcles and LTs. Mean length ( u 1 ) of short vehcles s around 4.88 m at SH6 and 4.72 m at IH-35; mean LT length ( u 2 ) s about m at SH6 and m at IH-35. IH 35 s an urban freeway whle SH6 s a rural area. The truck percentages are sgnfcantly dfferent as shown n Table 1. However, u 1 and u 2 are about the same. Further nvestgaton showed that at both locatons, u 1 and u 2 for dfferent tme perods pretty much stay the same even though truck volumes and percentage may vary sgnfcantly between tme perods. Therefore, u 1 and u 2 are treated as known parameters n our method and can be determned based on samples from hstorcal vehcle length data from a smlar hghway n the area. Estmaton of LT Volume Gven the MEVL value estmated n the prevous secton and the number of vehcles at a certan tme nterval, we can calculate the number of LTs durng ths tme perod by rearrangng Equaton 4: n 2 = N( τ u1) ( u u ) 2 1 (5) where u 1 and u 2 are mean short vehcle and LT lengths correspondng to a specfc hghway or locaton. Note that n 2 s rounded to the nearest nteger. 12

13 RESULTS AND DISCUSSIONS One set of weekday data from each TTI s test bed, as descrbed n the second secton, s processed and analyzed. The UKF algorthm s frstly used to generate speed nformaton. Then the average vehcle length at each tme nterval s calculated through Equaton 3. Comparsons between estmated and observed speed data, of 30 second tme ntervals n a day are shown n parts a) and b) of Fgures 3 and 4. The results demonstrate that the UKF does well n speed estmaton under both normal and congested condtons. The Root Mean Square Error (RMSE) s used for the evaluaton of speed estmaton. The RMSE s square root of Mean Square Error (MSE), whch can capture both the varance of errors and the bas of estmates. The RMSE values are shown n parts b) of Fgures 3 and 4. The estmated MEVL and vehcle length estmaton errors for all 30 seconds ntervals are shown n parts c) and d) of Fgures 3 and 4. It should be noted that the estmaton errors of the MEVL shown n the part d) s of Fgures 3 and 4 have some larger values and varatons durng nghttme and durng daytme when traffc congeston exsts, ndcatng potental estmaton accuracy ssues for both tme perods. The observed average vehcle lengths between mdnght and 6 a.m. are generally longer than m from the part c) s of Fgures 3 and 4, a result of very hgh percentage of LTs n early mornng hours when few passenger cars are on the road. Wth speed and MEVL estmated, the LT volume at each tme nterval s obtaned. For better llustraton, LT volume data are aggregated nto 1-hour ntervals. Fgure 5 shows both observed and estmated LT volumes. The estmaton curves capture the observaton curves very 13

14 well durng most hours of the day. A summary of hourly-based estmaton results of LT volume for the entre day s descrbed n Table 2. In ths table, two Measurements of Effectveness (MOEs) are used: the Mean Absolute Error (MAE) and Mean Absolute Percentage Error (MAPE). They are defned as MAE = MAPE = = 1 24 nˆ = 1 ( ) est. nˆ n ( ) est ( ) obs.. n ( ) n obs. ( ) obs. (6) where ( ) ˆ est. n s the estmated LT volume durng th hour and n s the observed LT volume. ( ) obs. The overall estmated LT percentages are very close (wthn 0.5%) to observed percentages wth slght underestmatons. By comparng the estmated LT percentages wth truck percentages that are computed based on vehcle classfcatons and shown n Table 1, the estmaton errors are less than 1% for both data sets. To further explore the factors affectng estmaton results, the effects of LT percentage and traffc volume level are nvestgated usng the SH6 data. Frstly, a study of estmaton errors under dfferent traffc volume levels s conducted. Part a) of Fgure 6 llustrates hourly traffc volumes throughout the experment whle part c) of ths fgure shows hourly Absolute Percentage Errors (APEs). Based on parts a) and c) of ths fgure, t s observed that APE values durng nghttme when traffc volume s relatvely low are generally not much smaller than those durng daytme when traffc volume s relatvely hgh, ndcatng that overall traffc volume does not seem to affect LT estmaton accuracy. Secondly, hourly-based LT percentages 14

15 are also calculated from the observed dataset and shown n the part b) of Fgure 6. The LT percentages are very hgh between mdnght and 6 a.m., and the maxmum percentage s close to 40%. Nevertheless, the relatve larger estmaton errors n the part c) are not generated durng those hgh LT percentage hours. Thus, the LT percentage tself does not show sgnfcant effects on the LT volume estmaton ether accordng to our experment. From Fgure 6(c), there are perods of tme durng the day when there are more fluctuatons n APE. For example, at the 22:00 hour, the APE reaches 20%. However, a closer lookng at Fgure 5 reveals that the estmaton error for the hour s only three LTs, 18 estmated versus 15 observed, and obvously acceptable. In addton, the fluctuatons llustrated n Fgure 6(c) are centered around the MAPE lne, ndcatng no system error trend wth respect to any partcular factor. Ths seems to ndcate that the estmaton results are not sgnfcantly affected by ether LT percentage or overall traffc volume. To further nvestgate the performance of the algorthm under varous traffc condtons, partcularly under congested condtons, data of two consecutve days (November 9 and 10, 2004) from IH-35 ste were analyzed. The ste has regular AM and PM recurrent congeston perods. Fgure 7 presents LT estmaton results along wth speed estmaton results and Fgure 8 llustrate the LT estmaton error dstrbuton. For ths analyss, LTs are estmated at 15-mnutes nterval, an nterval commonly used n real-tme traffc flow predcton. It can be seen that agan MAPE lne centers around zero, ndcatng no systematc estmaton errors. A closer look nto the error dstrbuton of the LT estmaton at 15-mnute nterval n Fgure 8 show that there are only 6 ntervals that have error over 15 LTs per nterval, or 1 LT per mnute on average. About 90% of 15

16 the ntervals have errors smaller than 10 Lts per nterval. About 80% of the ntervals have errors of less than 8 LTs, or about 1 LT every 2 mnutes. That s certanly good accuracy for practcal applcatons. From Fgure 7, wth the excepton of early mornng perods durng whch overall traffc s very lght but LT percentage s very hgh, LT estmaton accuracy dd not seem to show any trend wth respect to overall traffc volume, LT volume, or LT percentage. Durng the early mornng perods (when LT% s at the hghest), LT estmaton accuracy s the best. The worst estmaton accuracy occurred durng the perods of congeston. A further nvestgaton reveal that for the AM peak on the frst day and both the AM peak and PM peak on the second day the speed estmaton results are not as good as the rest of the perods. Concdently, those are the same three perods that have sgnfcant LT estmaton errors. Ths ndcates that traffc congeston s a factor nfluencng the LT volume estmaton results. Traffc congeston affects speed estmaton and consequently contrbutes to the accuracy of LT estmaton. CONCLUSIONS Ths paper frstly uses an accurate speed estmaton method to estmate speed and mean effectve vehcle length based on sngle loop detector outputs. A vehcle length sample s selected to nvestgate length dstrbutons, ncludng short and long vehcle length dstrbutons. Ths paper does not assume that ether short or long vehcle length dstrbuton ft a normal dstrbuton and proposes an algorthm for LT volume estmaton based on vehcle length estmaton results. Data collected from two Texas Transportaton Insttute s vehcle detecton test beds are used for ths 16

17 study. The estmated large truck volumes are compared wth feld data. The results show that the proposed method has good performances under varous traffc condtons. Fnally, the factors that nfluence the estmaton accuraces are explored. It s found that whle LT percentage and traffc volume generally do not sgnfcantly affect estmaton accuracy, traffc congeston s a factor affectng the results of large truck volume estmaton for the congested perods manly due to the fact that speed estmaton under congested flow condtons are typcally less accurate. Real tme large truck volume can be easly estmated usng ths proposed method and provde useful nformaton for traffc management and control n an Intellgent Transportaton System (ITS). The results can also be valuable to detaled and more relable analyses such as level of servce (LOS) analyss of a freeway segment or safety study of a locaton at a gven demand level. It should be ponted out that the UFK speed estmaton method s used n the paper for the purpose of provdng relable speed nputs for the LT estmaton methodology. Any accurate speed estmaton method usng sngle loop data can work wth the LT estmaton methodology proposed n ths paper. Ths study uses ADR sensor that provdes more relable and accurate measures of flow and occupancy than a normal sngle loop detector. ADR data are also used as ground truth data for valdaton. Actual sngle loop detectors have more measurement error such as mssng data or erroneous readngs. Because of ths, further nvestgaton and valdaton wth actual sngle loop and ground truth data are needed n the future. 17

18 ACKNOWLEDGEMENT The authors would lke to acknowledge TTI for ther support. Specal thanks go to Dan R. Mddleton for provdng detector data. REFERENCES 1 Cofman, B., and Ergueta, E. Improved Vehcle Redentfcaton and Travel Tme Measurement on Congested Freeways. In Proceedngs of the 81st Annual Meetng Transportaton Research Board, Washngton D.C., Wang, Y., and Nhan, N. L. Dynamc Estmaton of Freeway Large Truck Volume Based on Sngle-loop Measurements. Journal of Intellgent Transportaton Systems, Vol (8), No. 3,,2004, pp Kwon, J., Varaya, P., and Skabardons, A. Estmaton of Truck Traffc Volume from Sngle Loop Detector Usng Lane-to-Lane Speed Correlaton. In Transportaton Research Record: Journal of the Transportaton Research Record, No. 1232, Natonal Research Councl, Washngton, D.C., 2002, pp Zhang, G., Wang, Y., and We, H. Artfcal Neural Network Method for Vehcle Classfcaton Usng Sngle-Loop Outputs.. In Proceedngs of the 85th Annual Meetng of the Transportaton Research Board, Washngton D.C., Mddleton, D.R., and R.T. Parker. Summary Report on Evaluaton of Vehcle Detecton Systems. PSR S. Texas Transportaton Insttute, College Staton, TX. January Wyman, J.H., Gary, A. B., and Robert, I. S. Feld Evaluaton of FHWA Vehcle Classfcaton Categores. Mane Department of Transportaton. Fnal Report for Contract DTFH ME-03 for USDOT, Daley, D.J. A Statstcal Algorthm for Estmatng Speed from Sngle Loop Volume and Occupancy Measurements. Transportaton Research Part B, 33B(5), 1999, pp Athol, P. Interdependence of Certan Operatonal Characterstcs wthn a Movng Traffc Stream. Hghway Research Record 72, HRB, Natonal Research Councl, Washngton D.C., 1965, pp

19 9 Hall, F. An nterpretaton of Speed-flow-concentraton Relatonshps Usng Catastrophe Theory. Transportaton Research Part A, Vol.21A, No.3, 1987, pp Ln, W., Dahlgren, J., and Huo, H. Enhancement of Vehcle Speed Estmaton wth Sngle Loop Detectors. Transportaton Research Record 1870, Natonal Research Councl, Washngton, D.C., 2004, pp Juler, S., and Uhlmann, J. A General Method for Approxmatng Nonlnear Transformatons of Probablty Dstrbutons. Oxford Unversty: Oxford, Juler, S., and Uhlmann, J. A New Extenson of the Kalman Flter to Nonlnear Systems. Internatonal Symposum of Aerospace/Defense Sensng, Smulaton, and Controls, Orlando, FL, 1997, pp Ye, Z., Zhang, Y., and Mddleton, D. An Unscented Kalman Flter Method for Speed Estmaton Usng Sngle Loop Detector Data. In Transportaton Research Record: Journal of the Transportaton Research Record, No. 1968, Natonal Research Councl, Washngton, D.C., 2006, pp

20 LIST OF TABLES TABLE 1 Comparson of Truck Classfcaton Results Based on Scheme F and LTs TABLE 2 Results of Hourly-Based Estmaton of LT Volume LIST OF FIGURES FIGURE 1 Operaton process of the UKF FIGURE 2 Vehcle length dstrbutons FIGURE 3 Speed and vehcle length nformaton of SH6 (Jan. 27, 2004) FIGURE 4 Speed and vehcle length nformaton of IH-35 (Oct. 27, 2004) FIGURE 5 Observed and estmated LT volumes FIGURE 6 Hourly-based traffc volumes, LT percentages and estmaton errors on SH FIGURE 7 Speed and LT Volume Estmaton Results for IH-35 (Nov. 9-10, 2004) FIGURE 8 IH-35 LT Estmaton Error Dstrbuton (Nov. 9-10, 2004)

21 TABLE 1 Comparson of Truck Classfcaton Results Based on Scheme F and LTs Daly Traffc Volume (Veh.) Trucks (Based on Scheme F) Number of Percentage Trucks (%) (Veh.) LT (Based on Length) Number of Percentage LT (%) (Veh.) IH SH

22 TABLE 2 Results of Hourly-Based Estmaton of LT Volume Traffc Number of LT Total LT MAE of MAPE of Volume Real Estmated Error LT LT (veh.) (veh.) (veh.) (veh.) (veh.) (%) SH (6.26%) 795 (6.10%) IH (11.50%) 3064 (11.07%)

23 Intal estmates of state and covarance Calculate sgma ponts Tme Update 1) Project the state ahead 2) Project the error covarance ahead Measurement Update 1) Compute the Kalman gan 2) Update estmate wth measurement 3) Update the error FIGURE 1 Operaton process of the UKF. 23

24 FIGURE 2 Vehcle length dstrbutons. a) Kernel densty of vehcle lengths. b) Q-Q plot of short vehcle lengths. c) Q-Q plot of long vehcle lengths. 24

25 FIGURE 3 Speed and vehcle length nformaton of SH6 (Jan. 27, 2004). a) Observed speed. b) Estmated speed usng the UKF. c) Observed vehcle lengths. d) Estmaton error of vehcle length. 25

26 FIGURE 4 Speed and vehcle length nformaton of IH-35 (Oct. 27, 2004). a) Observed speed. b) Estmated speed usng the UKF. c) Observed vehcle lengths. d) Estmaton error of vehcle length. 26

27 FIGURE 5 Observed and estmated LT volumes. a) Hourly-based LT volumes on SH6. b) Hourly-based LT volumes on IH

28 FIGURE 6 Hourly-based traffc volumes, LT percentages and estmaton errors on SH6. a) Traffc volume. b) LT percentages. c) APE. 28

29 FIGURE 7 Speed and LT Volume Estmaton Results for IH-35 (Nov. 9-10, 2004). 29

30 LT Estmaton Error Dstrbuton Frequency Estmaton Error (vehs/15 mnutes) Error (LTs/15-mn) <= % <= % <= % >15 6 ntervals Max Error/Interval 22 FIGURE 8 IH-35 LT Estmaton Error Dstrbuton (Nov. 9-10, 2004). 30

Negative Binomial Regression

Negative Binomial Regression STATGRAPHICS Rev. 9/16/2013 Negatve Bnomal Regresson Summary... 1 Data Input... 3 Statstcal Model... 3 Analyss Summary... 4 Analyss Optons... 7 Plot of Ftted Model... 8 Observed Versus Predcted... 10 Predctons...

More information

Parameter Estimation for Dynamic System using Unscented Kalman filter

Parameter Estimation for Dynamic System using Unscented Kalman filter Parameter Estmaton for Dynamc System usng Unscented Kalman flter Jhoon Seung 1,a, Amr Atya F. 2,b, Alexander G.Parlos 3,c, and Klto Chong 1,4,d* 1 Dvson of Electroncs Engneerng, Chonbuk Natonal Unversty,

More information

Comparison of Regression Lines

Comparison of Regression Lines STATGRAPHICS Rev. 9/13/2013 Comparson of Regresson Lnes Summary... 1 Data Input... 3 Analyss Summary... 4 Plot of Ftted Model... 6 Condtonal Sums of Squares... 6 Analyss Optons... 7 Forecasts... 8 Confdence

More information

Chapter 13: Multiple Regression

Chapter 13: Multiple Regression Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to

More information

A Particle Filter Algorithm based on Mixing of Prior probability density and UKF as Generate Importance Function

A Particle Filter Algorithm based on Mixing of Prior probability density and UKF as Generate Importance Function Advanced Scence and Technology Letters, pp.83-87 http://dx.do.org/10.14257/astl.2014.53.20 A Partcle Flter Algorthm based on Mxng of Pror probablty densty and UKF as Generate Importance Functon Lu Lu 1,1,

More information

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers Psychology 282 Lecture #24 Outlne Regresson Dagnostcs: Outlers In an earler lecture we studed the statstcal assumptons underlyng the regresson model, ncludng the followng ponts: Formal statement of assumptons.

More information

Statistical Evaluation of WATFLOOD

Statistical Evaluation of WATFLOOD tatstcal Evaluaton of WATFLD By: Angela MacLean, Dept. of Cvl & Envronmental Engneerng, Unversty of Waterloo, n. ctober, 005 The statstcs program assocated wth WATFLD uses spl.csv fle that s produced wth

More information

Turbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH

Turbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH Turbulence classfcaton of load data by the frequency and severty of wnd gusts Introducton Oscar Moñux, DEWI GmbH Kevn Blebler, DEWI GmbH Durng the wnd turbne developng process, one of the most mportant

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

More information

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth

More information

Parking Demand Forecasting in Airport Ground Transportation System: Case Study in Hongqiao Airport

Parking Demand Forecasting in Airport Ground Transportation System: Case Study in Hongqiao Airport Internatonal Symposum on Computers & Informatcs (ISCI 25) Parkng Demand Forecastng n Arport Ground Transportaton System: Case Study n Hongqao Arport Ln Chang, a, L Wefeng, b*, Huanh Yan 2, c, Yang Ge,

More information

See Book Chapter 11 2 nd Edition (Chapter 10 1 st Edition)

See Book Chapter 11 2 nd Edition (Chapter 10 1 st Edition) Count Data Models See Book Chapter 11 2 nd Edton (Chapter 10 1 st Edton) Count data consst of non-negatve nteger values Examples: number of drver route changes per week, the number of trp departure changes

More information

/ n ) are compared. The logic is: if the two

/ n ) are compared. The logic is: if the two STAT C141, Sprng 2005 Lecture 13 Two sample tests One sample tests: examples of goodness of ft tests, where we are testng whether our data supports predctons. Two sample tests: called as tests of ndependence

More information

A LINEAR PROGRAM TO COMPARE MULTIPLE GROSS CREDIT LOSS FORECASTS. Dr. Derald E. Wentzien, Wesley College, (302) ,

A LINEAR PROGRAM TO COMPARE MULTIPLE GROSS CREDIT LOSS FORECASTS. Dr. Derald E. Wentzien, Wesley College, (302) , A LINEAR PROGRAM TO COMPARE MULTIPLE GROSS CREDIT LOSS FORECASTS Dr. Derald E. Wentzen, Wesley College, (302) 736-2574, wentzde@wesley.edu ABSTRACT A lnear programmng model s developed and used to compare

More information

Chapter 11: Simple Linear Regression and Correlation

Chapter 11: Simple Linear Regression and Correlation Chapter 11: Smple Lnear Regresson and Correlaton 11-1 Emprcal Models 11-2 Smple Lnear Regresson 11-3 Propertes of the Least Squares Estmators 11-4 Hypothess Test n Smple Lnear Regresson 11-4.1 Use of t-tests

More information

2016 Wiley. Study Session 2: Ethical and Professional Standards Application

2016 Wiley. Study Session 2: Ethical and Professional Standards Application 6 Wley Study Sesson : Ethcal and Professonal Standards Applcaton LESSON : CORRECTION ANALYSIS Readng 9: Correlaton and Regresson LOS 9a: Calculate and nterpret a sample covarance and a sample correlaton

More information

Composite Hypotheses testing

Composite Hypotheses testing Composte ypotheses testng In many hypothess testng problems there are many possble dstrbutons that can occur under each of the hypotheses. The output of the source s a set of parameters (ponts n a parameter

More information

Lecture 6: Introduction to Linear Regression

Lecture 6: Introduction to Linear Regression Lecture 6: Introducton to Lnear Regresson An Manchakul amancha@jhsph.edu 24 Aprl 27 Lnear regresson: man dea Lnear regresson can be used to study an outcome as a lnear functon of a predctor Example: 6

More information

Statistics for Managers Using Microsoft Excel/SPSS Chapter 13 The Simple Linear Regression Model and Correlation

Statistics for Managers Using Microsoft Excel/SPSS Chapter 13 The Simple Linear Regression Model and Correlation Statstcs for Managers Usng Mcrosoft Excel/SPSS Chapter 13 The Smple Lnear Regresson Model and Correlaton 1999 Prentce-Hall, Inc. Chap. 13-1 Chapter Topcs Types of Regresson Models Determnng the Smple Lnear

More information

Department of Quantitative Methods & Information Systems. Time Series and Their Components QMIS 320. Chapter 6

Department of Quantitative Methods & Information Systems. Time Series and Their Components QMIS 320. Chapter 6 Department of Quanttatve Methods & Informaton Systems Tme Seres and Ther Components QMIS 30 Chapter 6 Fall 00 Dr. Mohammad Zanal These sldes were modfed from ther orgnal source for educatonal purpose only.

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

Uncertainty as the Overlap of Alternate Conditional Distributions

Uncertainty as the Overlap of Alternate Conditional Distributions Uncertanty as the Overlap of Alternate Condtonal Dstrbutons Olena Babak and Clayton V. Deutsch Centre for Computatonal Geostatstcs Department of Cvl & Envronmental Engneerng Unversty of Alberta An mportant

More information

This column is a continuation of our previous column

This column is a continuation of our previous column Comparson of Goodness of Ft Statstcs for Lnear Regresson, Part II The authors contnue ther dscusson of the correlaton coeffcent n developng a calbraton for quanttatve analyss. Jerome Workman Jr. and Howard

More information

A Robust Method for Calculating the Correlation Coefficient

A Robust Method for Calculating the Correlation Coefficient A Robust Method for Calculatng the Correlaton Coeffcent E.B. Nven and C. V. Deutsch Relatonshps between prmary and secondary data are frequently quantfed usng the correlaton coeffcent; however, the tradtonal

More information

Statistics for Economics & Business

Statistics for Economics & Business Statstcs for Economcs & Busness Smple Lnear Regresson Learnng Objectves In ths chapter, you learn: How to use regresson analyss to predct the value of a dependent varable based on an ndependent varable

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

RBF Neural Network Model Training by Unscented Kalman Filter and Its Application in Mechanical Fault Diagnosis

RBF Neural Network Model Training by Unscented Kalman Filter and Its Application in Mechanical Fault Diagnosis Appled Mechancs and Materals Submtted: 24-6-2 ISSN: 662-7482, Vols. 62-65, pp 2383-2386 Accepted: 24-6- do:.428/www.scentfc.net/amm.62-65.2383 Onlne: 24-8- 24 rans ech Publcatons, Swtzerland RBF Neural

More information

Appendix B: Resampling Algorithms

Appendix B: Resampling Algorithms 407 Appendx B: Resamplng Algorthms A common problem of all partcle flters s the degeneracy of weghts, whch conssts of the unbounded ncrease of the varance of the mportance weghts ω [ ] of the partcles

More information

Uncertainty and auto-correlation in. Measurement

Uncertainty and auto-correlation in. Measurement Uncertanty and auto-correlaton n arxv:1707.03276v2 [physcs.data-an] 30 Dec 2017 Measurement Markus Schebl Federal Offce of Metrology and Surveyng (BEV), 1160 Venna, Austra E-mal: markus.schebl@bev.gv.at

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

x = , so that calculated

x = , so that calculated Stat 4, secton Sngle Factor ANOVA notes by Tm Plachowsk n chapter 8 we conducted hypothess tests n whch we compared a sngle sample s mean or proporton to some hypotheszed value Chapter 9 expanded ths to

More information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

Topic- 11 The Analysis of Variance

Topic- 11 The Analysis of Variance Topc- 11 The Analyss of Varance Expermental Desgn The samplng plan or expermental desgn determnes the way that a sample s selected. In an observatonal study, the expermenter observes data that already

More information

Lecture 12: Classification

Lecture 12: Classification Lecture : Classfcaton g Dscrmnant functons g The optmal Bayes classfer g Quadratc classfers g Eucldean and Mahalanobs metrcs g K Nearest Neghbor Classfers Intellgent Sensor Systems Rcardo Guterrez-Osuna

More information

Department of Statistics University of Toronto STA305H1S / 1004 HS Design and Analysis of Experiments Term Test - Winter Solution

Department of Statistics University of Toronto STA305H1S / 1004 HS Design and Analysis of Experiments Term Test - Winter Solution Department of Statstcs Unversty of Toronto STA35HS / HS Desgn and Analyss of Experments Term Test - Wnter - Soluton February, Last Name: Frst Name: Student Number: Instructons: Tme: hours. Ads: a non-programmable

More information

Chapter 14 Simple Linear Regression

Chapter 14 Simple Linear Regression Chapter 4 Smple Lnear Regresson Chapter 4 - Smple Lnear Regresson Manageral decsons often are based on the relatonshp between two or more varables. Regresson analss can be used to develop an equaton showng

More information

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method Appled Mathematcal Scences, Vol. 7, 0, no. 47, 07-0 HIARI Ltd, www.m-hkar.com Comparson of the Populaton Varance Estmators of -Parameter Exponental Dstrbuton Based on Multple Crtera Decson Makng Method

More information

Statistics for Business and Economics

Statistics for Business and Economics Statstcs for Busness and Economcs Chapter 11 Smple Regresson Copyrght 010 Pearson Educaton, Inc. Publshng as Prentce Hall Ch. 11-1 11.1 Overvew of Lnear Models n An equaton can be ft to show the best lnear

More information

Tracking with Kalman Filter

Tracking with Kalman Filter Trackng wth Kalman Flter Scott T. Acton Vrgna Image and Vdeo Analyss (VIVA), Charles L. Brown Department of Electrcal and Computer Engneerng Department of Bomedcal Engneerng Unversty of Vrgna, Charlottesvlle,

More information

Homework Assignment 3 Due in class, Thursday October 15

Homework Assignment 3 Due in class, Thursday October 15 Homework Assgnment 3 Due n class, Thursday October 15 SDS 383C Statstcal Modelng I 1 Rdge regresson and Lasso 1. Get the Prostrate cancer data from http://statweb.stanford.edu/~tbs/elemstatlearn/ datasets/prostate.data.

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

Regularized Discriminant Analysis for Face Recognition

Regularized Discriminant Analysis for Face Recognition 1 Regularzed Dscrmnant Analyss for Face Recognton Itz Pma, Mayer Aladem Department of Electrcal and Computer Engneerng, Ben-Guron Unversty of the Negev P.O.Box 653, Beer-Sheva, 845, Israel. Abstract Ths

More information

CHAPTER 4 SPEECH ENHANCEMENT USING MULTI-BAND WIENER FILTER. In real environmental conditions the speech signal may be

CHAPTER 4 SPEECH ENHANCEMENT USING MULTI-BAND WIENER FILTER. In real environmental conditions the speech signal may be 55 CHAPTER 4 SPEECH ENHANCEMENT USING MULTI-BAND WIENER FILTER 4.1 Introducton In real envronmental condtons the speech sgnal may be supermposed by the envronmental nterference. In general, the spectrum

More information

Supporting Information

Supporting Information Supportng Informaton The neural network f n Eq. 1 s gven by: f x l = ReLU W atom x l + b atom, 2 where ReLU s the element-wse rectfed lnear unt, 21.e., ReLUx = max0, x, W atom R d d s the weght matrx to

More information

1. Inference on Regression Parameters a. Finding Mean, s.d and covariance amongst estimates. 2. Confidence Intervals and Working Hotelling Bands

1. Inference on Regression Parameters a. Finding Mean, s.d and covariance amongst estimates. 2. Confidence Intervals and Working Hotelling Bands Content. Inference on Regresson Parameters a. Fndng Mean, s.d and covarance amongst estmates.. Confdence Intervals and Workng Hotellng Bands 3. Cochran s Theorem 4. General Lnear Testng 5. Measures of

More information

Temperature. Chapter Heat Engine

Temperature. Chapter Heat Engine Chapter 3 Temperature In prevous chapters of these notes we ntroduced the Prncple of Maxmum ntropy as a technque for estmatng probablty dstrbutons consstent wth constrants. In Chapter 9 we dscussed the

More information

Topic 23 - Randomized Complete Block Designs (RCBD)

Topic 23 - Randomized Complete Block Designs (RCBD) Topc 3 ANOVA (III) 3-1 Topc 3 - Randomzed Complete Block Desgns (RCBD) Defn: A Randomzed Complete Block Desgn s a varant of the completely randomzed desgn (CRD) that we recently learned. In ths desgn,

More information

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient 58:080 Expermental Engneerng 1 OBJECTIVE Lab 2e Thermal System Response and Effectve Heat Transfer Coeffcent Warnng: though the experment has educatonal objectves (to learn about bolng heat transfer, etc.),

More information

Copyright 2017 by Taylor Enterprises, Inc., All Rights Reserved. Adjusted Control Limits for P Charts. Dr. Wayne A. Taylor

Copyright 2017 by Taylor Enterprises, Inc., All Rights Reserved. Adjusted Control Limits for P Charts. Dr. Wayne A. Taylor Taylor Enterprses, Inc. Control Lmts for P Charts Copyrght 2017 by Taylor Enterprses, Inc., All Rghts Reserved. Control Lmts for P Charts Dr. Wayne A. Taylor Abstract: P charts are used for count data

More information

RELIABILITY ASSESSMENT

RELIABILITY ASSESSMENT CHAPTER Rsk Analyss n Engneerng and Economcs RELIABILITY ASSESSMENT A. J. Clark School of Engneerng Department of Cvl and Envronmental Engneerng 4a CHAPMAN HALL/CRC Rsk Analyss for Engneerng Department

More information

STAT 511 FINAL EXAM NAME Spring 2001

STAT 511 FINAL EXAM NAME Spring 2001 STAT 5 FINAL EXAM NAME Sprng Instructons: Ths s a closed book exam. No notes or books are allowed. ou may use a calculator but you are not allowed to store notes or formulas n the calculator. Please wrte

More information

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 31 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 6. Rdge regresson The OLSE s the best lnear unbased

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

Assessment of Site Amplification Effect from Input Energy Spectra of Strong Ground Motion

Assessment of Site Amplification Effect from Input Energy Spectra of Strong Ground Motion Assessment of Ste Amplfcaton Effect from Input Energy Spectra of Strong Ground Moton M.S. Gong & L.L Xe Key Laboratory of Earthquake Engneerng and Engneerng Vbraton,Insttute of Engneerng Mechancs, CEA,

More information

Power law and dimension of the maximum value for belief distribution with the max Deng entropy

Power law and dimension of the maximum value for belief distribution with the max Deng entropy Power law and dmenson of the maxmum value for belef dstrbuton wth the max Deng entropy Bngy Kang a, a College of Informaton Engneerng, Northwest A&F Unversty, Yanglng, Shaanx, 712100, Chna. Abstract Deng

More information

Boostrapaggregating (Bagging)

Boostrapaggregating (Bagging) Boostrapaggregatng (Baggng) An ensemble meta-algorthm desgned to mprove the stablty and accuracy of machne learnng algorthms Can be used n both regresson and classfcaton Reduces varance and helps to avod

More information

PHYS 450 Spring semester Lecture 02: Dealing with Experimental Uncertainties. Ron Reifenberger Birck Nanotechnology Center Purdue University

PHYS 450 Spring semester Lecture 02: Dealing with Experimental Uncertainties. Ron Reifenberger Birck Nanotechnology Center Purdue University PHYS 45 Sprng semester 7 Lecture : Dealng wth Expermental Uncertantes Ron Refenberger Brck anotechnology Center Purdue Unversty Lecture Introductory Comments Expermental errors (really expermental uncertantes)

More information

Andreas C. Drichoutis Agriculural University of Athens. Abstract

Andreas C. Drichoutis Agriculural University of Athens. Abstract Heteroskedastcty, the sngle crossng property and ordered response models Andreas C. Drchouts Agrculural Unversty of Athens Panagots Lazards Agrculural Unversty of Athens Rodolfo M. Nayga, Jr. Texas AMUnversty

More information

Chapter 3 Describing Data Using Numerical Measures

Chapter 3 Describing Data Using Numerical Measures Chapter 3 Student Lecture Notes 3-1 Chapter 3 Descrbng Data Usng Numercal Measures Fall 2006 Fundamentals of Busness Statstcs 1 Chapter Goals To establsh the usefulness of summary measures of data. The

More information

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U) Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of

More information

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur Analyss of Varance and Desgn of Experment-I MODULE VII LECTURE - 3 ANALYSIS OF COVARIANCE Dr Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur Any scentfc experment s performed

More information

Traffic Signal Timing: Basic Principles. Development of a Traffic Signal Phasing and Timing Plan. Two Phase and Three Phase Signal Operation

Traffic Signal Timing: Basic Principles. Development of a Traffic Signal Phasing and Timing Plan. Two Phase and Three Phase Signal Operation Traffc Sgnal Tmng: Basc Prncples 2 types of sgnals Pre-tmed Traffc actuated Objectves of sgnal tmng Reduce average delay of all vehcles Reduce probablty of accdents by mnmzng possble conflct ponts Objectves

More information

Module 9. Lecture 6. Duality in Assignment Problems

Module 9. Lecture 6. Duality in Assignment Problems Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept

More information

Using the estimated penetrances to determine the range of the underlying genetic model in casecontrol

Using the estimated penetrances to determine the range of the underlying genetic model in casecontrol Georgetown Unversty From the SelectedWorks of Mark J Meyer 8 Usng the estmated penetrances to determne the range of the underlyng genetc model n casecontrol desgn Mark J Meyer Neal Jeffres Gang Zheng Avalable

More information

BOOTSTRAP METHOD FOR TESTING OF EQUALITY OF SEVERAL MEANS. M. Krishna Reddy, B. Naveen Kumar and Y. Ramu

BOOTSTRAP METHOD FOR TESTING OF EQUALITY OF SEVERAL MEANS. M. Krishna Reddy, B. Naveen Kumar and Y. Ramu BOOTSTRAP METHOD FOR TESTING OF EQUALITY OF SEVERAL MEANS M. Krshna Reddy, B. Naveen Kumar and Y. Ramu Department of Statstcs, Osmana Unversty, Hyderabad -500 007, Inda. nanbyrozu@gmal.com, ramu0@gmal.com

More information

Continuous vs. Discrete Goods

Continuous vs. Discrete Goods CE 651 Transportaton Economcs Charsma Choudhury Lecture 3-4 Analyss of Demand Contnuous vs. Dscrete Goods Contnuous Goods Dscrete Goods x auto 1 Indfference u curves 3 u u 1 x 1 0 1 bus Outlne Data Modelng

More information

Correlation and Regression. Correlation 9.1. Correlation. Chapter 9

Correlation and Regression. Correlation 9.1. Correlation. Chapter 9 Chapter 9 Correlaton and Regresson 9. Correlaton Correlaton A correlaton s a relatonshp between two varables. The data can be represented b the ordered pars (, ) where s the ndependent (or eplanator) varable,

More information

Global Sensitivity. Tuesday 20 th February, 2018

Global Sensitivity. Tuesday 20 th February, 2018 Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance Calculation for Conductors of Arbitrary Shape CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors

More information

DETERMINATION OF UNCERTAINTY ASSOCIATED WITH QUANTIZATION ERRORS USING THE BAYESIAN APPROACH

DETERMINATION OF UNCERTAINTY ASSOCIATED WITH QUANTIZATION ERRORS USING THE BAYESIAN APPROACH Proceedngs, XVII IMEKO World Congress, June 7, 3, Dubrovn, Croata Proceedngs, XVII IMEKO World Congress, June 7, 3, Dubrovn, Croata TC XVII IMEKO World Congress Metrology n the 3rd Mllennum June 7, 3,

More information

Basically, if you have a dummy dependent variable you will be estimating a probability.

Basically, if you have a dummy dependent variable you will be estimating a probability. ECON 497: Lecture Notes 13 Page 1 of 1 Metropoltan State Unversty ECON 497: Research and Forecastng Lecture Notes 13 Dummy Dependent Varable Technques Studenmund Chapter 13 Bascally, f you have a dummy

More information

Durban Watson for Testing the Lack-of-Fit of Polynomial Regression Models without Replications

Durban Watson for Testing the Lack-of-Fit of Polynomial Regression Models without Replications Durban Watson for Testng the Lack-of-Ft of Polynomal Regresson Models wthout Replcatons Ruba A. Alyaf, Maha A. Omar, Abdullah A. Al-Shha ralyaf@ksu.edu.sa, maomar@ksu.edu.sa, aalshha@ksu.edu.sa Department

More information

Predictive Analytics : QM901.1x Prof U Dinesh Kumar, IIMB. All Rights Reserved, Indian Institute of Management Bangalore

Predictive Analytics : QM901.1x Prof U Dinesh Kumar, IIMB. All Rights Reserved, Indian Institute of Management Bangalore Sesson Outlne Introducton to classfcaton problems and dscrete choce models. Introducton to Logstcs Regresson. Logstc functon and Logt functon. Maxmum Lkelhood Estmator (MLE) for estmaton of LR parameters.

More information

4DVAR, according to the name, is a four-dimensional variational method.

4DVAR, according to the name, is a four-dimensional variational method. 4D-Varatonal Data Assmlaton (4D-Var) 4DVAR, accordng to the name, s a four-dmensonal varatonal method. 4D-Var s actually a drect generalzaton of 3D-Var to handle observatons that are dstrbuted n tme. The

More information

Low default modelling: a comparison of techniques based on a real Brazilian corporate portfolio

Low default modelling: a comparison of techniques based on a real Brazilian corporate portfolio Low default modellng: a comparson of technques based on a real Brazlan corporate portfolo MSc Gulherme Fernandes and MSc Carlos Rocha Credt Scorng and Credt Control Conference XII August 2011 Analytcs

More information

Pop-Click Noise Detection Using Inter-Frame Correlation for Improved Portable Auditory Sensing

Pop-Click Noise Detection Using Inter-Frame Correlation for Improved Portable Auditory Sensing Advanced Scence and Technology Letters, pp.164-168 http://dx.do.org/10.14257/astl.2013 Pop-Clc Nose Detecton Usng Inter-Frame Correlaton for Improved Portable Audtory Sensng Dong Yun Lee, Kwang Myung Jeon,

More information

Computing MLE Bias Empirically

Computing MLE Bias Empirically Computng MLE Bas Emprcally Kar Wa Lm Australan atonal Unversty January 3, 27 Abstract Ths note studes the bas arses from the MLE estmate of the rate parameter and the mean parameter of an exponental dstrbuton.

More information

Notes on Frequency Estimation in Data Streams

Notes on Frequency Estimation in Data Streams Notes on Frequency Estmaton n Data Streams In (one of) the data streamng model(s), the data s a sequence of arrvals a 1, a 2,..., a m of the form a j = (, v) where s the dentty of the tem and belongs to

More information

January Examinations 2015

January Examinations 2015 24/5 Canddates Only January Examnatons 25 DO NOT OPEN THE QUESTION PAPER UNTIL INSTRUCTED TO DO SO BY THE CHIEF INVIGILATOR STUDENT CANDIDATE NO.. Department Module Code Module Ttle Exam Duraton (n words)

More information

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran

More information

Chapter 8 Indicator Variables

Chapter 8 Indicator Variables Chapter 8 Indcator Varables In general, e explanatory varables n any regresson analyss are assumed to be quanttatve n nature. For example, e varables lke temperature, dstance, age etc. are quanttatve n

More information

Resource Allocation with a Budget Constraint for Computing Independent Tasks in the Cloud

Resource Allocation with a Budget Constraint for Computing Independent Tasks in the Cloud Resource Allocaton wth a Budget Constrant for Computng Independent Tasks n the Cloud Wemng Sh and Bo Hong School of Electrcal and Computer Engneerng Georga Insttute of Technology, USA 2nd IEEE Internatonal

More information

Week 9 Chapter 10 Section 1-5

Week 9 Chapter 10 Section 1-5 Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,

More information

CONTRAST ENHANCEMENT FOR MIMIMUM MEAN BRIGHTNESS ERROR FROM HISTOGRAM PARTITIONING INTRODUCTION

CONTRAST ENHANCEMENT FOR MIMIMUM MEAN BRIGHTNESS ERROR FROM HISTOGRAM PARTITIONING INTRODUCTION CONTRAST ENHANCEMENT FOR MIMIMUM MEAN BRIGHTNESS ERROR FROM HISTOGRAM PARTITIONING N. Phanthuna 1,2, F. Cheevasuvt 2 and S. Chtwong 2 1 Department of Electrcal Engneerng, Faculty of Engneerng Rajamangala

More information

x i1 =1 for all i (the constant ).

x i1 =1 for all i (the constant ). Chapter 5 The Multple Regresson Model Consder an economc model where the dependent varable s a functon of K explanatory varables. The economc model has the form: y = f ( x,x,..., ) xk Approxmate ths by

More information

The Study of Teaching-learning-based Optimization Algorithm

The Study of Teaching-learning-based Optimization Algorithm Advanced Scence and Technology Letters Vol. (AST 06), pp.05- http://dx.do.org/0.57/astl.06. The Study of Teachng-learnng-based Optmzaton Algorthm u Sun, Yan fu, Lele Kong, Haolang Q,, Helongang Insttute

More information

j) = 1 (note sigma notation) ii. Continuous random variable (e.g. Normal distribution) 1. density function: f ( x) 0 and f ( x) dx = 1

j) = 1 (note sigma notation) ii. Continuous random variable (e.g. Normal distribution) 1. density function: f ( x) 0 and f ( x) dx = 1 Random varables Measure of central tendences and varablty (means and varances) Jont densty functons and ndependence Measures of assocaton (covarance and correlaton) Interestng result Condtonal dstrbutons

More information

Generalized Linear Methods

Generalized Linear Methods Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set

More information

A Bayes Algorithm for the Multitask Pattern Recognition Problem Direct Approach

A Bayes Algorithm for the Multitask Pattern Recognition Problem Direct Approach A Bayes Algorthm for the Multtask Pattern Recognton Problem Drect Approach Edward Puchala Wroclaw Unversty of Technology, Char of Systems and Computer etworks, Wybrzeze Wyspanskego 7, 50-370 Wroclaw, Poland

More information

ANSWERS CHAPTER 9. TIO 9.2: If the values are the same, the difference is 0, therefore the null hypothesis cannot be rejected.

ANSWERS CHAPTER 9. TIO 9.2: If the values are the same, the difference is 0, therefore the null hypothesis cannot be rejected. ANSWERS CHAPTER 9 THINK IT OVER thnk t over TIO 9.: χ 2 k = ( f e ) = 0 e Breakng the equaton down: the test statstc for the ch-squared dstrbuton s equal to the sum over all categores of the expected frequency

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

Lecture 9: Linear regression: centering, hypothesis testing, multiple covariates, and confounding

Lecture 9: Linear regression: centering, hypothesis testing, multiple covariates, and confounding Recall: man dea of lnear regresson Lecture 9: Lnear regresson: centerng, hypothess testng, multple covarates, and confoundng Sandy Eckel seckel@jhsph.edu 6 May 8 Lnear regresson can be used to study an

More information

Feb 14: Spatial analysis of data fields

Feb 14: Spatial analysis of data fields Feb 4: Spatal analyss of data felds Mappng rregularly sampled data onto a regular grd Many analyss technques for geophyscal data requre the data be located at regular ntervals n space and/or tme. hs s

More information

8 Derivation of Network Rate Equations from Single- Cell Conductance Equations

8 Derivation of Network Rate Equations from Single- Cell Conductance Equations Physcs 178/278 - Davd Klenfeld - Wnter 2015 8 Dervaton of Network Rate Equatons from Sngle- Cell Conductance Equatons We consder a network of many neurons, each of whch obeys a set of conductancebased,

More information

Lecture 9: Linear regression: centering, hypothesis testing, multiple covariates, and confounding

Lecture 9: Linear regression: centering, hypothesis testing, multiple covariates, and confounding Lecture 9: Lnear regresson: centerng, hypothess testng, multple covarates, and confoundng Sandy Eckel seckel@jhsph.edu 6 May 008 Recall: man dea of lnear regresson Lnear regresson can be used to study

More information

DERIVATION OF THE PROBABILITY PLOT CORRELATION COEFFICIENT TEST STATISTICS FOR THE GENERALIZED LOGISTIC DISTRIBUTION

DERIVATION OF THE PROBABILITY PLOT CORRELATION COEFFICIENT TEST STATISTICS FOR THE GENERALIZED LOGISTIC DISTRIBUTION Internatonal Worshop ADVANCES IN STATISTICAL HYDROLOGY May 3-5, Taormna, Italy DERIVATION OF THE PROBABILITY PLOT CORRELATION COEFFICIENT TEST STATISTICS FOR THE GENERALIZED LOGISTIC DISTRIBUTION by Sooyoung

More information

Estimating the Fundamental Matrix by Transforming Image Points in Projective Space 1

Estimating the Fundamental Matrix by Transforming Image Points in Projective Space 1 Estmatng the Fundamental Matrx by Transformng Image Ponts n Projectve Space 1 Zhengyou Zhang and Charles Loop Mcrosoft Research, One Mcrosoft Way, Redmond, WA 98052, USA E-mal: fzhang,cloopg@mcrosoft.com

More information

LINEAR REGRESSION ANALYSIS. MODULE VIII Lecture Indicator Variables

LINEAR REGRESSION ANALYSIS. MODULE VIII Lecture Indicator Variables LINEAR REGRESSION ANALYSIS MODULE VIII Lecture - 7 Indcator Varables Dr. Shalabh Department of Maematcs and Statstcs Indan Insttute of Technology Kanpur Indcator varables versus quanttatve explanatory

More information

Feature Selection & Dynamic Tracking F&P Textbook New: Ch 11, Old: Ch 17 Guido Gerig CS 6320, Spring 2013

Feature Selection & Dynamic Tracking F&P Textbook New: Ch 11, Old: Ch 17 Guido Gerig CS 6320, Spring 2013 Feature Selecton & Dynamc Trackng F&P Textbook New: Ch 11, Old: Ch 17 Gudo Gerg CS 6320, Sprng 2013 Credts: Materal Greg Welch & Gary Bshop, UNC Chapel Hll, some sldes modfed from J.M. Frahm/ M. Pollefeys,

More information