Detec%ng and Analyzing Urban Regions with High Impact of Weather Change on Transport

Size: px
Start display at page:

Download "Detec%ng and Analyzing Urban Regions with High Impact of Weather Change on Transport"

Transcription

1 Detec%ng and Analyzing Urban Regions with High Impact of Weather Change on Transport Ye Ding, Yanhua Li, Ke Deng, Haoyu Tan, Mingxuan Yuan, Lionel M. Ni Presenta;on by Karan Somaiah Napanda, Suchithra Balakrishnan, Zhaoning Su

2 URBAN compu%ng connects urban sensing, data management, data analy%c and service providing *Urban compu;ng with taxis, MSRA *A Cloud-Based Knowledge Discovery System for Monitoring Fine-Grained Air Quality, MSRA

3 URBAN compu%ng connects urban sensing, data management, data analy%c and service providing Focus Traffic conges;on Energy consump;on Pollu;on Based on data traffic flow human mobility geographical data infer the real-;me and fine grained air quality informa;on throughout a city iden;fy the hot spots of moving vehicles in an urban area propose a framework, called DRoF, to discover regions of different func;ons in a city try to sense the refueling behavior and citywide petrol consump;on in real-;me

4 The impact of inclement weather to traffic May slow down the traffic May cause conges;ons due to low visibility and high demand of May influence the transport performance,

5 the impact of inclement weather to traffic *Changes in Climate and Weather Relevant on US Transport, The impact of climate change and weather on transport: An overview of empirical fin

6 Mo%va%on *IBM Smarter Traveler Traffic Predic;on So^ware *Google Opera;ng System how can we identify those regions being highly influenced by weather change on transport?

7 Challenges lacking of effec%ve traffic monitoring system in city-wide scale To enable weather-traffic index throughout a city and factor analysis extract traffic informa;on from numerous taxis driving on roads due to its availability, wide-coverage and low-cost. A taxi tracking system con;nuously record the informa;on including loca;on, speed, occupancy status, and orienta;on of the taxis

8 Challenges lacking of effec%ve traffic monitoring system in city-wide scale Voronoi diagram

9 Challenges lacking of effec%ve traffic monitoring system in city-wide scale Equal-sized rectangles Ø some cells are highly dense and in others are highly sparse Voronoi diagram Ø seeds are the road intersec;ons Ø every cell include at least one road intersec;on and a number of roads connected to this intersec;on Ø if several road intersec;ons are very close to each other, for example within 50 meters, they are grouped together as a complex intersec;on *The Voronoi diagrams par;;ons in Shanghai. The under layer represents the road networks.

10 Ques%ons The average of taxi can reflect the traffic informa%on? In the area of residence community, everyone have their own car so they don t need the taxi By calcula%ng the average driving speeds of all taxis in each Voronoi cell can reflect the average speed? However maybe some cells have a lower average speed of taxis is just because people are likely get on or get off from the taxis here

11 Challenges how to disclose the key factors behind the weather-traffic index density of roads number of road intersec;ons number of POIs(points of interest) traffic volume average age of the household density of buildings and more in the surrounding regions

12 The Goal of this paper To develope weather-traffic index (WTI) system The first is to set up a weather-traffic index throughout a city, which indicates the impact of weather to traffic from light to heavy. The second is to reveal the key factors behind the weather-traffic index throughout the city and their rela;ve weights. Previous works mainly focus on the analysis of road segments; on the contrary, this paper is the first study on local traffic-weather sensi;vity throughout a city and the inves;ga;on to reveal the key factors behind the sensi;vity

13 OVERVIEW 1. Data Prepara;on 2. Weather-traffic Index Establishment 3. Factor Analysis

14 1. DATA PREPARATION i. Regional Par;;oning Voronoi Diagram ii. Source Data Road Network Traffic Data Regional features Weather Report Data

15 DATA PREPARATION - Regional Par%%oning VORONOI DIAGRAM Par;;oning of a plane into cells/ regions based on the distance between seeds Shape and size of each cells is different from each other SEEDS- road intersec;ons ROAD- INTERSECTION- ORIENTED Par;;oning Proper;es Even distribu;on of road networks Potray the rela;on of weather and traffic

16 DATA PREPARATION Source Data i. Road Networks ii. Traffic Data iii. Regional Features iv. Weather Report Data i. Road Networks - G(V, E) E- set of road segments V- set of road intersec;ons E- type, length, speed limit, two end points V- loca;on (la;tude and longitude), type

17 DATA PREPARATION Source Data Time Mean Speed/ Average Speed- Traffic Parameter Arithme;c mean of individual spot speeds that are recorded over a selected ;me period Extracted from taxi trajectories in each cell ii. Traffic Data Traffic Parameters of interest are extracted Average Speed Quan;ty measures Quality assessment measures Movement measures Composi;on/ Classifica;on measures Split into 7 classes <10 km/h km/h km/h km/h km/h km/h >110 km/h The average speed of one road segment is subject to the traffic parameter of that road segment only, not comparable with other road segments

18 DATA PREPARATION Source Data iv. Regional Features Four Categories a) Points of Interest b) Structure c) Density d) Community iii. Weather Report Data State of atmosphere Degree to which it is hot or cold, wet or dry, calm or stormy, clear or cloudy

19 2. WEATHER TRAFFIC INDEX ESTABLISHMENT Input Weather data and Traffic data Indicates the impact of weather to traffic in different cells Given a cell g, its value in the weather traffic index is the correla;on between traffic and weather, denoted as ρ(g). ρ(g) takes a value from the discrete range [1, 2, 3, 4, 5]

20 WEATHER TRAFFIC INDEX ESTABLISHMENT Correla%on Detec%on Correla;on between traffic speed F t and weather F w. Classifier is trained with F t and F w Input- Weather as a feature vector Output- one of the seven speed classes Inference accuracy Correla;on between F t and F w in that cell Inference accuracy Correla;on between F t and F w in that cell Weakness of this method, Does not consider other reasons that affect traffic

21 WEATHER TRAFFIC INDEX ESTABLISHMENT Correla%on Detec%on Granger Causality test whether one ;me series is useful in forecas;ng another A ;me series X is said to Granger-cause Y if it can be shown that those X values provide sta;s;cally significant informa;on about future values of Y. A variable X that evolves over ;me Granger-cause another evolving variable Y if predic;ons of the value of Y based on its own past values and on the past values of X are beuer than predic;ons of Y based only on its own past values.

22 WEATHER TRAFFIC INDEX ESTABLISHMENT Correla%on Detec%on Traffic Predic;on models are trained separately for different ;me Weather-traffic index value ρ(g) is assigned to each cell to indicate the extent to which traffic is impacted by weather Cells are organised in ascending order of traffic predic;on accuracy improvement and then divided into k- equal sized subsets k- quan;les show the correla;on between traffic and weather from weak to strong

23 WEATHER TRAFFIC INDEX ESTABLISHMENT Traffic Predic%on Classifica;on Problem Methods used Support Vector Machine (SVM) Logis;c Regression (LOGIT) Perceptron 10-fold cross valida;on SVM is applied and LOGIT and Perceptron is used to verify the output of SVM

24 3. FACTOR ANALYSIS AssumpRon weather traffic index of all cells have been certainly assigned This method iden;fies the key factors and their weights contribu;ng to the weather- traffic indices of cells. Discloses what regional features make the traffic in some cells vulnerable to inclement weather Steps: 1. Key Factor Verifica;on by Index Inference (KFVII) 2. Weight Es;ma;on of Regional Features

25 FACTOR ANALYSIS Key Factor Verifica%on by Index Inference (KFVII) Intui;on Weather- traffic index of one region can be inferred from the indices of its closely located (or adjacent) cells Given a set of regional features F r, if the inference accuracy is sa;sfactory with F r as input, it indicates that such set of regional features are the key features This model is not symmetric Naïve Bayes classifier

26 FACTOR ANALYSIS Key Factor Verifica%on by Index Inference (KFVII) Marginal DistribuRon Probability distribu;on of the regions contained in a similarity subset Probability of one region being the index of i given one of its adjacent regions with index j, if two regions have a certain similarity Cosine Similarity is used to describe the similarity m uv between two regions g u and g v Similarity ranges of k- groups b0 is minimum similarity Bk is maximum similarity

27 FACTOR ANALYSIS Key Factor Verifica%on by Index Inference (KFVII) Pairs of cells in the group of [b i-1, b i ] are in marginal distribu;on matrix B i Rows of B i - weather traffic indices of g u Columns of B i - weather traffic indices of g v p ij - probability that the weather-traffic index of g v is ρ j when the weather- traffic index of g u is ρ i

28 FACTOR ANALYSIS Key Factor Verifica%on by Index Inference (KFVII) Index- Index Inference From marginal distribu;on, weather traffic index of a par;cular cell is inferred from adjacent cell using naïve Bayes classifier. Given a cell g u, the marginal distribu;on allows naïve Bayes classifier to infer which value the weather- traffic index of g u is most likely to be, based on the weathertraffic indices of its adjacent cells ρ(g 1 ), ρ(g 2 ),

29 FACTOR ANALYSIS Weight Es%ma%on of Regional Features Consider a regional feature with nontrivial impact to weather-traffic index, 1. Remove this regional features 2. Using KFVII method test if remaining features have high overall accuracy 3. If not, then removed regional feature is very important 4. δ(f i r) weight of the regional feature F i r by compu;ng the difference of the inference accuracy with and without F i r

30 EMPIRICAL STUDY Data Sets Road Network Data of Shanghai Taxi Trajectory Data Weather Report Data (of the same period of ;me as taxi trajectory data) Regional Informa;on Data

31 More about the Data Road Network Data Obtained from the government 7 levels of roads: Na;onal Expressway City Expressway Regular Highway Large Avenue Primary Way Secondary Way Regular Road Minor Roads Major Roads Weather Report Data Obtained from Weather Underground (wunderground.com) Includes 14 weather features For data alignment, recorded during the same ;me frame as taxi trajectory.

32 More about the Data Trajectory Data Spa;al-Temporal points Provides addi;onal informa;on including average speed. Regional Informa;on Data Real Estate Data & POI (Points of Interest) data. Real estate (soufun.com) provides informa;on about loca;on, price and residen;al communi;es. Dianping.com provides local merchant informa;on and other aurac;ons around the locality

33 Weather Traffic Index What is the WTI Average speed is inferred with/without weather. Difference indicates sensi;vity of the cell at that ;me slot. Robustness of the model: Accuracy of traffic predic;on is test with 3 models SVM gives 0.5 Perceptron and Logis;c regression give low accuracy, hence they will be used to test the robustness of WTI (with/without weather). A Posi;ve value indicates high impact or weather and a nega;ve value indicates low impact.

34 Robustness Evalua;on The above diagram shows most of the cells have an average accuracy of 0.5 for SVM. The above diagram gives changes in the predic;on accuracy which is the least for SVM. Thus, SVM can be used as the base model and Logis;c Regression and Perceptron are used to test. The changes in predic;on accuracy are close to mean for SVM and Logis;c Regression and it s a regular distribu;on where as for Perceptron is greater.

35 Valida%on The five labelled regions indicate high impact of weather on traffic. Most of them are tourist aurac;ons. Not all tourist aurac;ons have high impact due to inclement weather. eg: Xin;andi is not influenced by weather

36 Valida%on The rain aspect has been separated to show that it has significant impact on traffic. Example area 1: It shows that weather and rain both have significant impact on traffic. Rain label gives more informa;on about rain affec;ng traffic than weather affec;ng traffic. Example area 2: Average speed is more on rainy days.

37 Effec%veness Effec;veness verifies the regional features and es;mates their weights. Reciprocal method is used to evaluate It is considerably a fairer method than maximum likelihood since it widens the gap between different likelihoods. For this, Naive Bayes, Ar;ficial Neural Networks and Random Guess is used. Bayes performs the best, followed by ANN and then finally Random Guess.

38 Factor Analysis Four categories of regional features are tested: POI, structure, density and community. Figure shows that community regional features are more important than structure regional features. House Age & Number of Neighboring Cells are important factors.

39 Conclusion This paper fills the gap in the study on the impact of weather on traffic. The effec;veness of the system is verified. Regional factor analysis has significant impact: Regional House Age is an important factor The factors retrieved from this study can be applied to other ci;es to generate same knowledge. Further developments are expected since research in this field is always emerging.

40 Ques%ons?

Classifica(on and predic(on omics style. Dr Nicola Armstrong Mathema(cs and Sta(s(cs Murdoch University

Classifica(on and predic(on omics style. Dr Nicola Armstrong Mathema(cs and Sta(s(cs Murdoch University Classifica(on and predic(on omics style Dr Nicola Armstrong Mathema(cs and Sta(s(cs Murdoch University Classifica(on Learning Set Data with known classes Prediction Classification rule Data with unknown

More information

CS 6140: Machine Learning Spring What We Learned Last Week 2/26/16

CS 6140: Machine Learning Spring What We Learned Last Week 2/26/16 Logis@cs CS 6140: Machine Learning Spring 2016 Instructor: Lu Wang College of Computer and Informa@on Science Northeastern University Webpage: www.ccs.neu.edu/home/luwang Email: luwang@ccs.neu.edu Sign

More information

CS 6140: Machine Learning Spring What We Learned Last Week. Survey 2/26/16. VS. Model

CS 6140: Machine Learning Spring What We Learned Last Week. Survey 2/26/16. VS. Model Logis@cs CS 6140: Machine Learning Spring 2016 Instructor: Lu Wang College of Computer and Informa@on Science Northeastern University Webpage: www.ccs.neu.edu/home/luwang Email: luwang@ccs.neu.edu Assignment

More information

CS 6140: Machine Learning Spring 2016

CS 6140: Machine Learning Spring 2016 CS 6140: Machine Learning Spring 2016 Instructor: Lu Wang College of Computer and Informa?on Science Northeastern University Webpage: www.ccs.neu.edu/home/luwang Email: luwang@ccs.neu.edu Logis?cs Assignment

More information

Bias/variance tradeoff, Model assessment and selec+on

Bias/variance tradeoff, Model assessment and selec+on Applied induc+ve learning Bias/variance tradeoff, Model assessment and selec+on Pierre Geurts Department of Electrical Engineering and Computer Science University of Liège October 29, 2012 1 Supervised

More information

STAD68: Machine Learning

STAD68: Machine Learning STAD68: Machine Learning Russ Salakhutdinov Department of Statistics! rsalakhu@utstat.toronto.edu! h0p://www.cs.toronto.edu/~rsalakhu/ Lecture 1 Evalua;on 3 Assignments worth 40%. Midterm worth 20%. Final

More information

Predic've Analy'cs for Energy Systems State Es'ma'on

Predic've Analy'cs for Energy Systems State Es'ma'on 1 Predic've Analy'cs for Energy Systems State Es'ma'on Presenter: Team: Yingchen (YC) Zhang Na/onal Renewable Energy Laboratory (NREL) Andrey Bernstein, Rui Yang, Yu Xie, Anthony Florita, Changfu Li, ScoD

More information

Logis&c Regression. Robot Image Credit: Viktoriya Sukhanova 123RF.com

Logis&c Regression. Robot Image Credit: Viktoriya Sukhanova 123RF.com Logis&c Regression These slides were assembled by Eric Eaton, with grateful acknowledgement of the many others who made their course materials freely available online. Feel free to reuse or adapt these

More information

COMP 562: Introduction to Machine Learning

COMP 562: Introduction to Machine Learning COMP 562: Introduction to Machine Learning Lecture 20 : Support Vector Machines, Kernels Mahmoud Mostapha 1 Department of Computer Science University of North Carolina at Chapel Hill mahmoudm@cs.unc.edu

More information

DART Tutorial Sec'on 1: Filtering For a One Variable System

DART Tutorial Sec'on 1: Filtering For a One Variable System DART Tutorial Sec'on 1: Filtering For a One Variable System UCAR The Na'onal Center for Atmospheric Research is sponsored by the Na'onal Science Founda'on. Any opinions, findings and conclusions or recommenda'ons

More information

Regression.

Regression. Regression www.biostat.wisc.edu/~dpage/cs760/ Goals for the lecture you should understand the following concepts linear regression RMSE, MAE, and R-square logistic regression convex functions and sets

More information

VISUAL EXPLORATION OF SPATIAL-TEMPORAL TRAFFIC CONGESTION PATTERNS USING FLOATING CAR DATA. Candra Kartika 2015

VISUAL EXPLORATION OF SPATIAL-TEMPORAL TRAFFIC CONGESTION PATTERNS USING FLOATING CAR DATA. Candra Kartika 2015 VISUAL EXPLORATION OF SPATIAL-TEMPORAL TRAFFIC CONGESTION PATTERNS USING FLOATING CAR DATA Candra Kartika 2015 OVERVIEW Motivation Background and State of The Art Test data Visualization methods Result

More information

Linear Regression and Correla/on. Correla/on and Regression Analysis. Three Ques/ons 9/14/14. Chapter 13. Dr. Richard Jerz

Linear Regression and Correla/on. Correla/on and Regression Analysis. Three Ques/ons 9/14/14. Chapter 13. Dr. Richard Jerz Linear Regression and Correla/on Chapter 13 Dr. Richard Jerz 1 Correla/on and Regression Analysis Correla/on Analysis is the study of the rela/onship between variables. It is also defined as group of techniques

More information

Linear Regression and Correla/on

Linear Regression and Correla/on Linear Regression and Correla/on Chapter 13 Dr. Richard Jerz 1 Correla/on and Regression Analysis Correla/on Analysis is the study of the rela/onship between variables. It is also defined as group of techniques

More information

An Introduc+on to Sta+s+cs and Machine Learning for Quan+ta+ve Biology. Anirvan Sengupta Dept. of Physics and Astronomy Rutgers University

An Introduc+on to Sta+s+cs and Machine Learning for Quan+ta+ve Biology. Anirvan Sengupta Dept. of Physics and Astronomy Rutgers University An Introduc+on to Sta+s+cs and Machine Learning for Quan+ta+ve Biology Anirvan Sengupta Dept. of Physics and Astronomy Rutgers University Why Do We Care? Necessity in today s labs Principled approach:

More information

UVA CS / Introduc8on to Machine Learning and Data Mining

UVA CS / Introduc8on to Machine Learning and Data Mining UVA CS 4501-001 / 6501 007 Introduc8on to Machine Learning and Data Mining Lecture 13: Probability and Sta3s3cs Review (cont.) + Naïve Bayes Classifier Yanjun Qi / Jane, PhD University of Virginia Department

More information

Founda'ons of Large- Scale Mul'media Informa'on Management and Retrieval. Lecture #3 Machine Learning. Edward Chang

Founda'ons of Large- Scale Mul'media Informa'on Management and Retrieval. Lecture #3 Machine Learning. Edward Chang Founda'ons of Large- Scale Mul'media Informa'on Management and Retrieval Lecture #3 Machine Learning Edward Y. Chang Edward Chang Founda'ons of LSMM 1 Edward Chang Foundations of LSMM 2 Machine Learning

More information

Cri$ques Ø 5 cri&ques in total Ø Each with 6 points

Cri$ques Ø 5 cri&ques in total Ø Each with 6 points Cri$ques Ø 5 cri&ques in total Ø Each with 6 points 1 Distributed Applica$on Alloca$on in Shared Sensor Networks Chengjie Wu, You Xu, Yixin Chen, Chenyang Lu Shared Sensor Network Example in San Francisco

More information

STA 4273H: Sta-s-cal Machine Learning

STA 4273H: Sta-s-cal Machine Learning STA 4273H: Sta-s-cal Machine Learning Russ Salakhutdinov Department of Computer Science! Department of Statistical Sciences! rsalakhu@cs.toronto.edu! h0p://www.cs.utoronto.ca/~rsalakhu/ Lecture 1 Evalua:on

More information

Least Mean Squares Regression. Machine Learning Fall 2017

Least Mean Squares Regression. Machine Learning Fall 2017 Least Mean Squares Regression Machine Learning Fall 2017 1 Lecture Overview Linear classifiers What func?ons do linear classifiers express? Least Squares Method for Regression 2 Where are we? Linear classifiers

More information

PSAAP Project Stanford

PSAAP Project Stanford PSAAP Project QMU @ Stanford Component Analysis and rela:on to Full System Simula:ons 1 What do we want to predict? Objec:ve: predic:on of the unstart limit expressed as probability of unstart (or alterna:vely

More information

DART Tutorial Part IV: Other Updates for an Observed Variable

DART Tutorial Part IV: Other Updates for an Observed Variable DART Tutorial Part IV: Other Updates for an Observed Variable UCAR The Na'onal Center for Atmospheric Research is sponsored by the Na'onal Science Founda'on. Any opinions, findings and conclusions or recommenda'ons

More information

Least Squares Parameter Es.ma.on

Least Squares Parameter Es.ma.on Least Squares Parameter Es.ma.on Alun L. Lloyd Department of Mathema.cs Biomathema.cs Graduate Program North Carolina State University Aims of this Lecture 1. Model fifng using least squares 2. Quan.fica.on

More information

Structural Equa+on Models: The General Case. STA431: Spring 2013

Structural Equa+on Models: The General Case. STA431: Spring 2013 Structural Equa+on Models: The General Case STA431: Spring 2013 An Extension of Mul+ple Regression More than one regression- like equa+on Includes latent variables Variables can be explanatory in one equa+on

More information

Sta$s$cal sequence recogni$on

Sta$s$cal sequence recogni$on Sta$s$cal sequence recogni$on Determinis$c sequence recogni$on Last $me, temporal integra$on of local distances via DP Integrates local matches over $me Normalizes $me varia$ons For cts speech, segments

More information

Computer Vision. Pa0ern Recogni4on Concepts. Luis F. Teixeira MAP- i 2014/15

Computer Vision. Pa0ern Recogni4on Concepts. Luis F. Teixeira MAP- i 2014/15 Computer Vision Pa0ern Recogni4on Concepts Luis F. Teixeira MAP- i 2014/15 Outline General pa0ern recogni4on concepts Classifica4on Classifiers Decision Trees Instance- Based Learning Bayesian Learning

More information

Class Notes. Examining Repeated Measures Data on Individuals

Class Notes. Examining Repeated Measures Data on Individuals Ronald Heck Week 12: Class Notes 1 Class Notes Examining Repeated Measures Data on Individuals Generalized linear mixed models (GLMM) also provide a means of incorporang longitudinal designs with categorical

More information

Predicate abstrac,on and interpola,on. Many pictures and examples are borrowed from The So'ware Model Checker BLAST presenta,on.

Predicate abstrac,on and interpola,on. Many pictures and examples are borrowed from The So'ware Model Checker BLAST presenta,on. Predicate abstrac,on and interpola,on Many pictures and examples are borrowed from The So'ware Model Checker BLAST presenta,on. Outline. Predicate abstrac,on the idea in pictures 2. Counter- example guided

More information

Machine Learning and Data Mining. Bayes Classifiers. Prof. Alexander Ihler

Machine Learning and Data Mining. Bayes Classifiers. Prof. Alexander Ihler + Machine Learning and Data Mining Bayes Classifiers Prof. Alexander Ihler A basic classifier Training data D={x (i),y (i) }, Classifier f(x ; D) Discrete feature vector x f(x ; D) is a con@ngency table

More information

Effect of Environmental Factors on Free-Flow Speed

Effect of Environmental Factors on Free-Flow Speed Effect of Environmental Factors on Free-Flow Speed MICHAEL KYTE ZAHER KHATIB University of Idaho, USA PATRICK SHANNON Boise State University, USA FRED KITCHENER Meyer Mohaddes Associates, USA ABSTRACT

More information

Machine learning for Dynamic Social Network Analysis

Machine learning for Dynamic Social Network Analysis Machine learning for Dynamic Social Network Analysis Manuel Gomez Rodriguez Max Planck Ins7tute for So;ware Systems UC3M, MAY 2017 Interconnected World SOCIAL NETWORKS TRANSPORTATION NETWORKS WORLD WIDE

More information

CSCI 360 Introduc/on to Ar/ficial Intelligence Week 2: Problem Solving and Op/miza/on

CSCI 360 Introduc/on to Ar/ficial Intelligence Week 2: Problem Solving and Op/miza/on CSCI 360 Introduc/on to Ar/ficial Intelligence Week 2: Problem Solving and Op/miza/on Professor Wei-Min Shen Week 13.1 and 13.2 1 Status Check Extra credits? Announcement Evalua/on process will start soon

More information

Graphical Models. Lecture 3: Local Condi6onal Probability Distribu6ons. Andrew McCallum

Graphical Models. Lecture 3: Local Condi6onal Probability Distribu6ons. Andrew McCallum Graphical Models Lecture 3: Local Condi6onal Probability Distribu6ons Andrew McCallum mccallum@cs.umass.edu Thanks to Noah Smith and Carlos Guestrin for some slide materials. 1 Condi6onal Probability Distribu6ons

More information

Exploring the Patterns of Human Mobility Using Heterogeneous Traffic Trajectory Data

Exploring the Patterns of Human Mobility Using Heterogeneous Traffic Trajectory Data Exploring the Patterns of Human Mobility Using Heterogeneous Traffic Trajectory Data Jinzhong Wang April 13, 2016 The UBD Group Mobile and Social Computing Laboratory School of Software, Dalian University

More information

Mining Classification Knowledge

Mining Classification Knowledge Mining Classification Knowledge Remarks on NonSymbolic Methods JERZY STEFANOWSKI Institute of Computing Sciences, Poznań University of Technology COST Doctoral School, Troina 2008 Outline 1. Bayesian classification

More information

CSE546: SVMs, Dual Formula5on, and Kernels Winter 2012

CSE546: SVMs, Dual Formula5on, and Kernels Winter 2012 CSE546: SVMs, Dual Formula5on, and Kernels Winter 2012 Luke ZeClemoyer Slides adapted from Carlos Guestrin Linear classifiers Which line is becer? w. = j w (j) x (j) Data Example i Pick the one with the

More information

Short introduc,on to the

Short introduc,on to the OXFORD NEUROIMAGING PRIMERS Short introduc,on to the An General Introduction Linear Model to Neuroimaging for Neuroimaging Analysis Mark Jenkinson Mark Jenkinson Janine Michael Bijsterbosch Chappell Michael

More information

Correla'on. Keegan Korthauer Department of Sta's'cs UW Madison

Correla'on. Keegan Korthauer Department of Sta's'cs UW Madison Correla'on Keegan Korthauer Department of Sta's'cs UW Madison 1 Rela'onship Between Two Con'nuous Variables When we have measured two con$nuous random variables for each item in a sample, we can study

More information

Gene Regulatory Networks II Computa.onal Genomics Seyoung Kim

Gene Regulatory Networks II Computa.onal Genomics Seyoung Kim Gene Regulatory Networks II 02-710 Computa.onal Genomics Seyoung Kim Goal: Discover Structure and Func;on of Complex systems in the Cell Identify the different regulators and their target genes that are

More information

BBM406 - Introduc0on to ML. Spring Ensemble Methods. Aykut Erdem Dept. of Computer Engineering HaceDepe University

BBM406 - Introduc0on to ML. Spring Ensemble Methods. Aykut Erdem Dept. of Computer Engineering HaceDepe University BBM406 - Introduc0on to ML Spring 2014 Ensemble Methods Aykut Erdem Dept. of Computer Engineering HaceDepe University 2 Slides adopted from David Sontag, Mehryar Mohri, Ziv- Bar Joseph, Arvind Rao, Greg

More information

Ensemble of Climate Models

Ensemble of Climate Models Ensemble of Climate Models Claudia Tebaldi Climate Central and Department of Sta7s7cs, UBC Reto Knu>, Reinhard Furrer, Richard Smith, Bruno Sanso Outline Mul7 model ensembles (MMEs) a descrip7on at face

More information

2.1 Traffic Stream Characteristics. Time Space Diagram and Measurement Procedures Variables of Interest

2.1 Traffic Stream Characteristics. Time Space Diagram and Measurement Procedures Variables of Interest 2.1 Traffic Stream Characteristics Time Space Diagram and Measurement Procedures Variables of Interest Traffic Stream Models 2.1 Traffic Stream Characteristics Time Space Diagram Speed =100km/h = 27.78

More information

Priors in Dependency network learning

Priors in Dependency network learning Priors in Dependency network learning Sushmita Roy sroy@biostat.wisc.edu Computa:onal Network Biology Biosta2s2cs & Medical Informa2cs 826 Computer Sciences 838 hbps://compnetbiocourse.discovery.wisc.edu

More information

MAVEN OBSERVATIONS OF SIDING SPRING AND MARS. Roger Yelle and MAVEN Team

MAVEN OBSERVATIONS OF SIDING SPRING AND MARS. Roger Yelle and MAVEN Team MAVEN OBSERVATIONS OF SIDING SPRING AND MARS Roger Yelle and MAVEN Team APP/CSS Mee=ng, 11 August 2014 Some Per=nent Facts The closest approach occurs on 19 Oct 2014, at 18:28 UT at a distance of 136,000

More information

Molinas. June 15, 2018

Molinas. June 15, 2018 ITT8 SAMBa Presentation June 15, 2018 ling Data The data we have include: Approx 30,000 questionnaire responses each with 234 questions during 1998-2017 A data set of 60 questions asked to 500,000 households

More information

CSE 473: Ar+ficial Intelligence. Hidden Markov Models. Bayes Nets. Two random variable at each +me step Hidden state, X i Observa+on, E i

CSE 473: Ar+ficial Intelligence. Hidden Markov Models. Bayes Nets. Two random variable at each +me step Hidden state, X i Observa+on, E i CSE 473: Ar+ficial Intelligence Bayes Nets Daniel Weld [Most slides were created by Dan Klein and Pieter Abbeel for CS188 Intro to AI at UC Berkeley. All CS188 materials are available at hnp://ai.berkeley.edu.]

More information

Sta$s$cal Significance Tes$ng In Theory and In Prac$ce

Sta$s$cal Significance Tes$ng In Theory and In Prac$ce Sta$s$cal Significance Tes$ng In Theory and In Prac$ce Ben Cartere8e University of Delaware h8p://ir.cis.udel.edu/ictir13tutorial Hypotheses and Experiments Hypothesis: Using an SVM for classifica$on will

More information

Linear Regression with mul2ple variables. Mul2ple features. Machine Learning

Linear Regression with mul2ple variables. Mul2ple features. Machine Learning Linear Regression with mul2ple variables Mul2ple features Machine Learning Mul4ple features (variables). Size (feet 2 ) Price ($1000) 2104 460 1416 232 1534 315 852 178 Mul4ple features (variables). Size

More information

Traffic Demand Forecast

Traffic Demand Forecast Chapter 5 Traffic Demand Forecast One of the important objectives of traffic demand forecast in a transportation master plan study is to examine the concepts and policies in proposed plans by numerically

More information

Outline. What is Machine Learning? Why Machine Learning? 9/29/08. Machine Learning Approaches to Biological Research: Bioimage Informa>cs and Beyond

Outline. What is Machine Learning? Why Machine Learning? 9/29/08. Machine Learning Approaches to Biological Research: Bioimage Informa>cs and Beyond Outline Machine Learning Approaches to Biological Research: Bioimage Informa>cs and Beyond Robert F. Murphy External Senior Fellow, Freiburg Ins>tute for Advanced Studies Ray and Stephanie Lane Professor

More information

Transport in the extratropical UTLS as revealed by chemical tracers during the START08 campaign

Transport in the extratropical UTLS as revealed by chemical tracers during the START08 campaign Transport in the extratropical UTLS as revealed by chemical tracers during the START08 campaign Jasna V. Pi*man, Laura L. Pan Na#onal Center for Atmospheric Research INTRODUCTION Climate is sensi4ve to

More information

Estimating Large Scale Population Movement ML Dublin Meetup

Estimating Large Scale Population Movement ML Dublin Meetup Deutsche Bank COO Chief Data Office Estimating Large Scale Population Movement ML Dublin Meetup John Doyle PhD Assistant Vice President CDO Research & Development Science & Innovation john.doyle@db.com

More information

Experimental Designs for Planning Efficient Accelerated Life Tests

Experimental Designs for Planning Efficient Accelerated Life Tests Experimental Designs for Planning Efficient Accelerated Life Tests Kangwon Seo and Rong Pan School of Compu@ng, Informa@cs, and Decision Systems Engineering Arizona State University ASTR 2015, Sep 9-11,

More information

Boos$ng Can we make dumb learners smart?

Boos$ng Can we make dumb learners smart? Boos$ng Can we make dumb learners smart? Aarti Singh Machine Learning 10-601 Nov 29, 2011 Slides Courtesy: Carlos Guestrin, Freund & Schapire 1 Why boost weak learners? Goal: Automa'cally categorize type

More information

Generative Model (Naïve Bayes, LDA)

Generative Model (Naïve Bayes, LDA) Generative Model (Naïve Bayes, LDA) IST557 Data Mining: Techniques and Applications Jessie Li, Penn State University Materials from Prof. Jia Li, sta3s3cal learning book (Has3e et al.), and machine learning

More information

Be ready to take a few more Notes (Chapter 10 Section 3) Chapter 10 Quiz next Tuesday. If you were absent yesterday - graph review page - Exit ticket

Be ready to take a few more Notes (Chapter 10 Section 3) Chapter 10 Quiz next Tuesday. If you were absent yesterday - graph review page - Exit ticket February 6 Be ready to take a few more Notes (Chapter 10 Section 3) Chapter 10 Quiz next Tuesday If you were absent yesterday - graph review page - Exit ticket Chapter 10 Section 3 Acceleration measures

More information

TRAVEL TIME RELIABILITY ON EXPRESSWAY NETWORK UNDER UNCERTAIN ENVIRONMENT OF SNOWFALL AND TRAFFIC REGULATION

TRAVEL TIME RELIABILITY ON EXPRESSWAY NETWORK UNDER UNCERTAIN ENVIRONMENT OF SNOWFALL AND TRAFFIC REGULATION TRAVEL TIME RELIABILITY ON EXPRESSWAY NETWORK UNDER UNCERTAIN ENVIRONMENT OF SNOWFALL AND TRAFFIC REGULATION Hiroshi Wakabayashi Faculty of Urban Science, Meijo University, 4-3-3, Nijigaoka, Kani-City,

More information

Computer Vision. Pa0ern Recogni4on Concepts Part I. Luis F. Teixeira MAP- i 2012/13

Computer Vision. Pa0ern Recogni4on Concepts Part I. Luis F. Teixeira MAP- i 2012/13 Computer Vision Pa0ern Recogni4on Concepts Part I Luis F. Teixeira MAP- i 2012/13 What is it? Pa0ern Recogni4on Many defini4ons in the literature The assignment of a physical object or event to one of

More information

Spatio-temporal models

Spatio-temporal models Spatio-temporal models Involve a least a three dimensional representation of one or more key attribute variation in planar (X-Y) space and through time. (a 4 th dimension could also be use, like Z for

More information

Paul Waddell Professor, City and Regional Planning University of California, Berkeley Director, Urban Analy;cs Lab President, UrbanSim Inc.

Paul Waddell Professor, City and Regional Planning University of California, Berkeley Director, Urban Analy;cs Lab President, UrbanSim Inc. 2017 Philomathia Forum Urban Data Science and Simula9on For Metropolitan Sustainability Paul Waddell Professor, City and Regional Planning University of California, Berkeley Director, Urban Analy;cs Lab

More information

KNMI Climate Explorer A tool for climate analysis [ and seasonal prediction]

KNMI Climate Explorer A tool for climate analysis [ and seasonal prediction] KNMI Climate Explorer A tool for climate analysis [ and seasonal prediction] MedCOF Training Workshop Madrid, 26-30 October 2015 Jonathan Eden Royal Netherlands Meteorological Institute (KNMI) An introduc+on

More information

DART Ini)al Condi)ons for a Refined Grid CAM- SE Forecast of Hurricane Katrina. Kevin Raeder (IMAGe) Colin Zarzycki (ASP)

DART Ini)al Condi)ons for a Refined Grid CAM- SE Forecast of Hurricane Katrina. Kevin Raeder (IMAGe) Colin Zarzycki (ASP) DART Ini)al Condi)ons for a Refined Grid CAM- SE Forecast of Hurricane Katrina Kevin Raeder (IMAGe) Colin Zarzycki (ASP) 1 Mo)va)on Thousands of processors on current supercomputers. - > new CAM dynamical

More information

UVA CS 4501: Machine Learning. Lecture 6: Linear Regression Model with Dr. Yanjun Qi. University of Virginia

UVA CS 4501: Machine Learning. Lecture 6: Linear Regression Model with Dr. Yanjun Qi. University of Virginia UVA CS 4501: Machine Learning Lecture 6: Linear Regression Model with Regulariza@ons Dr. Yanjun Qi University of Virginia Department of Computer Science Where are we? è Five major sec@ons of this course

More information

The Changing Landscape of Land Administration

The Changing Landscape of Land Administration The Changing Landscape of Land Administration B r e n t J o n e s P E, PLS E s r i World s Largest Media Company No Journalists No Content Producers No Photographers World s Largest Hospitality Company

More information

Statistical Learning Learning From Examples

Statistical Learning Learning From Examples Statistical Learning Learning From Examples We want to estimate the working temperature range of an iphone. We could study the physics and chemistry that affect the performance of the phone too hard We

More information

Natural Language Processing with Deep Learning. CS224N/Ling284

Natural Language Processing with Deep Learning. CS224N/Ling284 Natural Language Processing with Deep Learning CS224N/Ling284 Lecture 4: Word Window Classification and Neural Networks Christopher Manning and Richard Socher Classifica6on setup and nota6on Generally

More information

Mining Classification Knowledge

Mining Classification Knowledge Mining Classification Knowledge Remarks on NonSymbolic Methods JERZY STEFANOWSKI Institute of Computing Sciences, Poznań University of Technology SE lecture revision 2013 Outline 1. Bayesian classification

More information

Machine Learning for NLP

Machine Learning for NLP Machine Learning for NLP Linear Models Joakim Nivre Uppsala University Department of Linguistics and Philology Slides adapted from Ryan McDonald, Google Research Machine Learning for NLP 1(26) Outline

More information

Short- term solar forecas/ng with sta/s/cal models and combina/on of models

Short- term solar forecas/ng with sta/s/cal models and combina/on of models Short- term solar forecas/ng with sta/s/cal models and combina/on of models Philippe Lauret, Mathieu David PIMENT Laboratory, University of La Réunion L. Mazorra Aguiar University of Las Palmas de Gran

More information

Encapsulating Urban Traffic Rhythms into Road Networks

Encapsulating Urban Traffic Rhythms into Road Networks Encapsulating Urban Traffic Rhythms into Road Networks Junjie Wang +, Dong Wei +, Kun He, Hang Gong, Pu Wang * School of Traffic and Transportation Engineering, Central South University, Changsha, Hunan,

More information

Slides modified from: PATTERN RECOGNITION AND MACHINE LEARNING CHRISTOPHER M. BISHOP

Slides modified from: PATTERN RECOGNITION AND MACHINE LEARNING CHRISTOPHER M. BISHOP Slides modified from: PATTERN RECOGNITION AND MACHINE LEARNING CHRISTOPHER M. BISHOP Predic?ve Distribu?on (1) Predict t for new values of x by integra?ng over w: where The Evidence Approxima?on (1) The

More information

Diagnosing New York City s Noises with Ubiquitous Data

Diagnosing New York City s Noises with Ubiquitous Data Diagnosing New York City s Noises with Ubiquitous Data Dr. Yu Zheng yuzheng@microsoft.com Lead Researcher, Microsoft Research Chair Professor at Shanghai Jiao Tong University Background Many cities suffer

More information

Op#mized Tropical Cyclone Retrievals from ASCAT, OceanSAT- 2 and QuikSCAT

Op#mized Tropical Cyclone Retrievals from ASCAT, OceanSAT- 2 and QuikSCAT Op#mized Tropical Cyclone Retrievals from ASCAT, OceanSAT- 2 and QuikSCAT Bryan W. S#les, Alex Fore, W. Lee Poulsen, Svetla- Hristova Veleva Jet Propulsion Laboratory, California Ins#tute of Technology

More information

REGRESSION AND CORRELATION ANALYSIS

REGRESSION AND CORRELATION ANALYSIS Problem 1 Problem 2 A group of 625 students has a mean age of 15.8 years with a standard devia>on of 0.6 years. The ages are normally distributed. How many students are younger than 16.2 years? REGRESSION

More information

Achieving a Readable Style Part 2: Sentence Structure

Achieving a Readable Style Part 2: Sentence Structure Achieving a Readable Style Part 2: Sentence Structure A wri&ng workshop presented by BACTER (Bringing Advanced Computa&onal Techniques to Energy Research) University of Wisconsin Madison October 2009 Last

More information

Fi#h Mee(ng of the Inter-agency and Expert Group on Sustainable Development Goals Indicators O"awa, Canada March 2017

Fi#h Mee(ng of the Inter-agency and Expert Group on Sustainable Development Goals Indicators Oawa, Canada March 2017 Fi#h Mee(ng of the Inter-agency and Expert Group on Sustainable Development Goals Indicators O"awa, Canada 28 31 March 2017 Report of the IAEG-SDGs Working Group on Geospa(al Informa(on co-chairs: Ms Marie

More information

Development and Evaluation of an Automated System for Empirical Forecasting of Flares

Development and Evaluation of an Automated System for Empirical Forecasting of Flares Development and Evaluation of an Automated System for Empirical Forecasting of Flares Sung-Hong Park1 Collaborators: J. Baek2, S. Bong2, Y. Yuan3, J. Jing3 1 National Observatory of Athens (NOA) 2 Korea

More information

1. Introduc9on 2. Bivariate Data 3. Linear Analysis of Data

1. Introduc9on 2. Bivariate Data 3. Linear Analysis of Data Lecture 3: Bivariate Data & Linear Regression 1. Introduc9on 2. Bivariate Data 3. Linear Analysis of Data a) Freehand Linear Fit b) Least Squares Fit c) Interpola9on/Extrapola9on 4. Correla9on 1. Introduc9on

More information

CIV3703 Transport Engineering. Module 2 Transport Modelling

CIV3703 Transport Engineering. Module 2 Transport Modelling CIV3703 Transport Engineering Module Transport Modelling Objectives Upon successful completion of this module you should be able to: carry out trip generation calculations using linear regression and category

More information

Announcements. Topics: Work On: - sec0ons 1.2 and 1.3 * Read these sec0ons and study solved examples in your textbook!

Announcements. Topics: Work On: - sec0ons 1.2 and 1.3 * Read these sec0ons and study solved examples in your textbook! Announcements Topics: - sec0ons 1.2 and 1.3 * Read these sec0ons and study solved examples in your textbook! Work On: - Prac0ce problems from the textbook and assignments from the coursepack as assigned

More information

Managing Growth: Integrating Land Use & Transportation Planning

Managing Growth: Integrating Land Use & Transportation Planning Managing Growth: Integrating Land Use & Transportation Planning Metro Vancouver Sustainability Community Breakfast Andrew Curran Manager, Strategy June 12, 2013 2 Integrating Land Use & Transportation

More information

Figure 8.2a Variation of suburban character, transit access and pedestrian accessibility by TAZ label in the study area

Figure 8.2a Variation of suburban character, transit access and pedestrian accessibility by TAZ label in the study area Figure 8.2a Variation of suburban character, transit access and pedestrian accessibility by TAZ label in the study area Figure 8.2b Variation of suburban character, commercial residential balance and mix

More information

Texas A&M University

Texas A&M University Texas A&M University CVEN 658 Civil Engineering Applications of GIS Hotspot Analysis of Highway Accident Spatial Pattern Based on Network Spatial Weights Instructor: Dr. Francisco Olivera Author: Zachry

More information

Rela%ons and Their Proper%es. Slides by A. Bloomfield

Rela%ons and Their Proper%es. Slides by A. Bloomfield Rela%ons and Their Proper%es Slides by A. Bloomfield What is a rela%on Let A and B be sets. A binary rela%on R is a subset of A B Example Let A be the students in a the CS major A = {Alice, Bob, Claire,

More information

CS 188: Artificial Intelligence Spring Announcements

CS 188: Artificial Intelligence Spring Announcements CS 188: Artificial Intelligence Spring 2011 Lecture 12: Probability 3/2/2011 Pieter Abbeel UC Berkeley Many slides adapted from Dan Klein. 1 Announcements P3 due on Monday (3/7) at 4:59pm W3 going out

More information

Caesar s Taxi Prediction Services

Caesar s Taxi Prediction Services 1 Caesar s Taxi Prediction Services Predicting NYC Taxi Fares, Trip Distance, and Activity Paul Jolly, Boxiao Pan, Varun Nambiar Abstract In this paper, we propose three models each predicting either taxi

More information

Urban Computing Using Big Data to Solve Urban Challenges

Urban Computing Using Big Data to Solve Urban Challenges Urban Computing Using Big Data to Solve Urban Challenges Dr. Yu Zheng Lead Researcher, Microsoft Research Asia Chair Professor at Shanghai Jiaotong University http://research.microsoft.com/en-us/projects/urbancomputing/default.aspx

More information

Par$$oned Elias- Fano indexes. Giuseppe O)aviano Rossano Venturini

Par$$oned Elias- Fano indexes. Giuseppe O)aviano Rossano Venturini Par$$oned Elias- Fano indexes Giuseppe O)aviano Rossano Venturini Inverted indexes Core data structure of Informa$on Retrieval Documents are sequences of terms 1: [it is what it is not] 2: [what is a]

More information

A Framework for Protec/ng Worker Loca/on Privacy in Spa/al Crowdsourcing

A Framework for Protec/ng Worker Loca/on Privacy in Spa/al Crowdsourcing A Framework for Protec/ng Worker Loca/on Privacy in Spa/al Crowdsourcing Nov 12 2014 Hien To, Gabriel Ghinita, Cyrus Shahabi VLDB 2014 1 Mo/va/on Ubiquity of mobile users 6.5 billion mobile subscrip/ons,

More information

I. M. Schoeman North West University, South Africa. Abstract

I. M. Schoeman North West University, South Africa. Abstract Urban Transport XX 607 Land use and transportation integration within the greater area of the North West University (Potchefstroom Campus), South Africa: problems, prospects and solutions I. M. Schoeman

More information

West meets East: Monitoring and modeling urbanization in China Land Cover-Land Use Change Program Science Team Meeting April 3, 2012

West meets East: Monitoring and modeling urbanization in China Land Cover-Land Use Change Program Science Team Meeting April 3, 2012 West meets East: Monitoring and modeling urbanization in China Land Cover-Land Use Change Program Science Team Meeting April 3, 2012 Annemarie Schneider Center for Sustainability and the Global Environment,

More information

City monitoring with travel demand momentum vector fields: theoretical and empirical findings

City monitoring with travel demand momentum vector fields: theoretical and empirical findings City monitoring with travel demand momentum vector fields: theoretical and empirical findings Xintao Liu 1, Joseph Y.J. Chow 2 1 Department of Civil Engineering, Ryerson University, Canada 2 Tandon School

More information

Whether you are driving or walking, if you come to a flooded road, Turn Around Don't Drown

Whether you are driving or walking, if you come to a flooded road, Turn Around Don't Drown Whether you are driving or walking, if you come to a flooded road, Turn Around Don't Drown You will not know the depth of the water nor will you know the condition of the road under the water. Many people

More information

Pilot Evaluation of a UAS Detect-and-Avoid System s Effectiveness in Remaining Well Clear!

Pilot Evaluation of a UAS Detect-and-Avoid System s Effectiveness in Remaining Well Clear! National Aeronautics and Space Administration Pilot Evaluation of a UAS Detect-and-Avoid System s Effectiveness in Remaining Well Clear! Presenter: Confesor Santiago Co-author: Eric Mueller Affiliation:

More information

Seman&cs with Dense Vectors. Dorota Glowacka

Seman&cs with Dense Vectors. Dorota Glowacka Semancs with Dense Vectors Dorota Glowacka dorota.glowacka@ed.ac.uk Previous lectures: - how to represent a word as a sparse vector with dimensions corresponding to the words in the vocabulary - the values

More information

Ensemble Data Assimila.on and Uncertainty Quan.fica.on

Ensemble Data Assimila.on and Uncertainty Quan.fica.on Ensemble Data Assimila.on and Uncertainty Quan.fica.on Jeffrey Anderson, Alicia Karspeck, Tim Hoar, Nancy Collins, Kevin Raeder, Steve Yeager Na.onal Center for Atmospheric Research Ocean Sciences Mee.ng

More information

About places and/or important events Landmarks Maps How the land is, hills or flat or mountain range Connected to maps World Different countries

About places and/or important events Landmarks Maps How the land is, hills or flat or mountain range Connected to maps World Different countries What do you think you know about geography? About places and/or important events Landmarks Maps How the land is, hills or flat or mountain range Connected to maps World Different countries What do you

More information

Unsupervised Learning of Discrimina4ve Rela4ve Visual A9ributes

Unsupervised Learning of Discrimina4ve Rela4ve Visual A9ributes Unsupervised Learning of Discrimina4ve Rela4ve Visual A9ributes Shugao Ma, Stan Sclaroff, Nazli Ikizler- Cinbis * Department of Computer Science, Boston University * Department of Computer Engineering,

More information

Prediction of Monthly Rainfall of Nainital Region using Artificial Neural Network (ANN) and Support Vector Machine (SVM)

Prediction of Monthly Rainfall of Nainital Region using Artificial Neural Network (ANN) and Support Vector Machine (SVM) Vol- Issue-3 25 Prediction of ly of Nainital Region using Artificial Neural Network (ANN) and Support Vector Machine (SVM) Deepa Bisht*, Mahesh C Joshi*, Ashish Mehta** *Department of Mathematics **Department

More information