Interpersonal and Intrapersonal Variabilities in Daily Activity-Travel Patterns: A Networked Spatiotemporal Analysis

Interpersonal and intrapersonal variabilities are two important perspectives to understand daily travel behaviors, while only a small number of studies incorporate them for understanding human dynamics. This paper employed a network analysis approach to detecting daily activity-travel patterns of 680 Beijing’s residents within a week and then used a multilevel multinomial logit model to analyze the intrapersonal variability in patterns and the socioeconomic linkages behind them. Results suggest that most activity-travel patterns have significant day-to-day intrapersonal and interpersonal variabilities. This suggests that the application of a typical day of activity-travel behaviors to measure and represent a week’s or even longer-term behaviors may be biased, due to the existence of day-to-day intrapersonal variability. This study also provides a hint for the selection of days of a week to conduct a diary survey for activity pattern mining or travel demand modeling.


Introduction
With a shift from road capacity expansion to travel demand management in transportation planning, it is necessary to better understand the variability of travel behavior [1]. Behavioral variability includes interpersonal variability and intrapersonal variability. Most travel studies focus on explaining the variability of travel behavior of different individuals (or groups) in terms of socioeconomic attributes, built environment, and perceived accessibility [2][3][4][5][6]. In other words, the complexity of daily activity-travel patterns is subject to several spatiotemporal patterns determined by a set of individual characteristics of the traveler. The other one is intrapersonal variability, which describes the travel behavioral variability of individuals (or groups) from trip to trip, day to day, and week to week [7][8][9][10].
Most existing research on behavioral variability focused on interpersonal variability, while intrapersonal variability are under-researched. Individual daily activity-travel patterns were often assumed to be stable in the short term, and prediction of multi-day behaviors was usually achieved by choosing the travel behavioral characteristics of "typical day" to represent multi-day features [10]. In the reality, people's daily travel patterns may significantly vary across days and weeks. Ignorance of day-to-day variability can easily lead to errors of travel demand model estimation and model interpretation [11]. Researchers found that sociodemographic characteristics which have been used to capture interpersonal variability hold large difference in explaining the level of intrapersonal variability [2,10], which shows the importance of intrapersonal variability analysis.
There are two major types of research on intrapersonal variability. One is based on the characteristics of travel behavior based on trip frequency, trip chain, travel time, time allocation, route choice, mode choice, and activity space [12][13][14][15]. However, these simplified measures do not regard the behavioral trajectory as a whole, and the integrity of travel behavior was neglected. The other is to measure the intrapersonal variability based on behavioral integrity, for example, based on activity chain patterns [16,17] and the measure of spatiotemporal prism [1,18]. Some attempts emphasize on the regularity of individual travel behavior, representing the repeated activity-travel episodes over time [19][20][21][22]. However, most holistic measurement of intrapersonal variability cannot well consider the characteristics of its sequence, dynamics, spatial-temporal interaction, and multi-dimensional attributes in pattern mining. This is mainly due to the lack of multi-day mobility datasets and appropriate trajectory clustering methods. Traditionally, the collection of mobility data relies on the Origin-Destination (OD) data collection technology through travel surveys, including the in-person/in-home interviews since the 1970s, the computer-aided telephone interviewing (CATI) technique in the 1990s, and current Web, mail, GPS, and automatic fare collection (AFC) technologies [23,24]. Particularly, development of internet communication technology (ICT) and wearable devices (e.g., cell phone and GPS devices) make it possible for large-scale, long time-series, and fine-grained spatial-temporal trajectories, researchers have been able to observe behavioral variability at larger temporal scale, cover multi-dimensional attributes. To deal with these datasets, Zhang and Thill [25] proposed a spatiotemporal behavioral trajectory pattern mining method from the field of network science. This method not only measures the sequence of activities effectively, but also ensures its temporal concurrence. However, they fail to consider the day-to-day variability of activity-travel patterns and do not check the factors resulting in such interpersonal and intrapersonal variabilities.
To fill these gaps, this study thus focuses on two research questions: (1) How to extract variant activity-travel patterns from daily trajectory data and examine their intrapersonal variabilities? (2) Do the occurrence likelihoods of different activity-travel patterns significantly vary across the day of a week (i.e., intrapersonal variability) and between individual (i.e., interpersonal variability)? To examine these questions, we used the GPS-based activity-travel dairy data of suburban residents in 2012 Beijing. We firstly applied Zhang and Thill's network analysis approach [25] to capturing day-to-day (intrapersonal) variability of activity-travel patterns in a week, thus detecting the typology of representative activity-travel patterns. We then developed a multilevel multinomial logit model to estimate the intrapersonal variability of different activity-travel patterns in a week and the relationship between socioeconomic attributes between the clusters of patterns.

Activity-Travel Pattern Mining
Activity-travel pattern mining can inform transportation policies with the information of target segments for tailor-made policy measures. To gain insights into the interpersonal and intrapersonal variability in activity-travel patterns, most of the research followed a two-step approach [26,27]. First is to find the interdependent activity-travel pattern with various pattern recognition methods. Second is to elucidate the relationship between these interdependent patterns and the corresponding factors with correlation analysis. These help transportation planners design better policies and envision more accurate forecasts of individual's travel demands.
Current patterns mining methods for spatiotemporal behavior trajectories mainly include hot-spot detection, clustering based on trajectory similarity, and sequential alignment method [28][29][30][31]. First, hotspot detection refers to the analysis of spatial and temporal density distribution of point or line elements in trajectory data [32][33][34], which decomposes the space-time paths into several discrete events and aggregates them according to their spatial and temporal distribution and density characteristics. However, it often needs to split a whole behavior trajectory into several fragments, so only fragmented behavior pattern features can be obtained.
In addition to the hot-spot detection, there are some other qualitative and quantitative technologies to extract space-time behavioral patterns from trajectories data. As summarized by Thériault et al. [35], there are two complementary GIS methods used to analyze the dynamics of entities in space, including deductive methods that search for qualitative reasoning behind the spatial dynamics and inductive methods that examine the properties and trends of the distribution of the entities in space using spatial statistical toolkits, such as the point clustering analysis, minimum convex polygon detection, centrographic analysis, and spatial autocorrelation analysis. Thériault et al. [35] further developed an approach combining both deductive and inductive methods to understand the qualitative reasoning behind the point changes in spatiotemporal evolution of entities, as well as to measure the overall spatial patterns and their evolution of a set of entities. Similar qualitative or quantitative measures have been used to categorize varying behavioral patterns in space and time as well as to track the changing pattern of activity space [12,36,37].
Second, there are many types of trajectory-based similarity measurement and dissimilarity measurement (e.g., distance). One of the most commonly used methods is to extract key variables from spatiotemporal behavior trajectories and measure vector distances, followed by cluster analysis [16,20]. Common distance measures include Euclidean distance, Mahalanobis distance, Hamming distance, minimum outer rectangular distance, etc. [33,38,39]. However, these distance-based approaches are often difficult to consider the similarity of trajectory's temporal dimension. They are not suitable for trajectory data with different lengths, acquisition frequencies and time scales, and often neglect the sequence, or the order of activities [22].
The third genre of pattern mining relies on the sequential alignment method (SAM), which measures the distance between behavior trajectories directly, including dynamic time planning (DTW), longest common subsequence (LCS), and edit distance-based alignment method. A key advantage of SAM is its capability to consider the multiple attributes and the compositional and sequential characteristics of human activity patterns [39]. Early sequential alignment measurements focused mainly on activity type alignment of trajectory sequence, without considering other multi-dimensional attributes of activity, such as location, timing, duration, travel mode, arrival time at the destination, accompany person, and existence of secondary activities, while recent analysis began to consider these attributes, forming a multi-dimensional sequential alignment method [2,29,[39][40][41]. Although the SAM can effectively measure the sequence similarity of activity events, it ignores the temporal concurrence of activity occurrence, and is difficult to measure the spatial interaction between activities. Moreover, the computational cost is high, especially for large dataset, and the computational efficiency is low [25].

Interpersonal and Intrapersonal Variabilities
Many interpersonal variability studies focused on the associations between sociodemographic attributes or built environment and specific individual's mobility features, such as route choice, travel mode, destination choice, health-engagement activity. For example, Abdel-Aty and Kitamura [42] explored interpersonal variability in relation to route choice, the results showed that travel time, traffic safety, and roadway characteristics are of significance on route choice. Larsen et al. [43] found that gender, trip distance, land use mix level, and presence of street trees are positively associated to children's travel mode. Yang et al. [44] concluded individual occupation as an important factor which influences the work destination choice. Moudon et al. [45] found strong correlation between sociodemographic attributes, built environment and individual's health-enhancement walking activity. Jahre et al. [46] studied the type of bike use and its relation to sociodemographic attributes, finding that ethnicity correlated to the bike use. In addition to urban activity-travel behavior, studies also focused on the rural travel activity. Fan et al. [47] studied the influence of built environment and sociodemographic characteristics on active commuting in America rural area, finding that sociodemographic factors explained more variance in active commuting than physical environmental factors. Some researchers also investigated interpersonal variability from the perspective of household characteristics, and perceived accessibility variables [5,48,49]. These studies provide interesting insights into the interpersonal variability, but the intrapersonal variability, or day-to-day variability, was ignored and left uncovered.
There are two main causes of the intrapersonal variability in travel behavior: (1) From the perspective of time geography, individual behavior is driven by needs and desires, which are determined by a set of constraints. The needs, desires and the constraints on individuals are different every day [1,50]. (2) The variation of daytime behavior is the feedback to the built environment. Changes in urban environment, transportation system, population or socio-economic structure at both macro and micro levels induce the temporal variation of activity-travel patterns. The dynamics of behavior may come from the randomness associated with travelers' values, perceptions, attitudes, needs, preferences, and decision-making processes [1].
At present, the research on intrapersonal variability mainly focuses on individual urban travel behavior in trip frequency, trip chain, departure time, route choice, time allocation, activity space, and other behavioral indicators by using longitudinal travel survey data. Pas [8] used the dataset of an activity diary survey to measure intrapersonal variability by examining the daily trip frequency (the difference between individual daily trip frequency and average trip frequency over a period). Hatcher and Mahmassani [51] used data from a sample of commuters in Austin, Texas, to examine whether the departure time or route for a given day was different from the previous day, as well as the deviation from the median departure time and the most commonly used route. Taking these as indices to measure intrapersonal variability of individuals, they studied the dayto-day variability of trip frequency, departure time and morning commuting route choice. Chikaraishi et al. [52] decomposed the total variation of departure time choice into five components: spatial variation, temporal variation, inter-household variation, inter-individual variation, and intra-individual variation. Moreover, they used multi-level analysis approach to analyze the proportion of each type of variation. Susilo and Axhausen [53] explored the stability and variability of individual daily activity-travel-location patterns based on Herfindahl-Hirschman index.
As the emergence of multiple-day daily data, a limited number of recent studies has attempted to investigate both behavioral variability in a comprehensive and integrated framework. For example, Dharmowijoyo et al. [5] studied day-to-day interpersonal and intrapersonal variability of individuals' activity spaces with four-consecutive-day Household Activity-Travel Survey Database of Jakarta Metropolitan Area. Their results showed that intrapersonal variability is more important than interpersonal in terms of total squared distance of individual out-of-home activity locations to the centroid of activity locations. Zhang et al. [10] investigated multi-day activity-travel patterns sampling based on single-day data with the Mobidrive six-week travel diary dataset and reported that interpersonal variability observed in cross-sectional single-day data of a group of people can be used to generate the day-to-day intrapersonal variability. Egu and Bonnel [6] explored day-to-day variability of transit usage on a multi-month scale with smart card data. Findings suggested that it is possible to correlate this intrapersonal variability to interpersonal variability using cluster analysis. All the above studies show substantial relationships between interpersonal and intrapersonal variability, and more importantly, highlight the importance of combining interpersonal and intrapersonal variations to gain indepth understandings on behavioral variability studies.

Method
The research methods of this article include two parts. We first applied the network analysis approach proposed by Zhang and Thill [25] to convert individual-based trajectory data into individual-by-individual networked data and used the community detection algorithm (i.e., the Louvain method) to cluster individuals into groups by their cohesive activity-travel patterns. Different from Zhang and Thill's analysis, this study continued to compare the intrapersonal (day-to-day) variability of daily activity-travel patterns in a week. In order to statistically test the intrapersonal and interpersonal variations in daily activity-travel patterns, we used a multilevel multinomial logit model to examine whether the activity-travel patterns significantly vary across days of a week and their associations with individual's socioeconomic attributes.

A Network Analysis Approach
Based on the social affiliation theory of Georg Simmel, Zhang and Thill [25] proposed a new method of activity-travel patterns mining. They constructed the interpersonal relationship network by calculating the strength (or degree of similarity) between individuals and used the community detection algorithm in the field of complex network analysis to conduct pattern mining and visualization of individual's spatiotemporal activity-travel trajectories.
This method consists of two steps. Firstly, we need to transform the spatiotemporal activity-travel trajectory (or time-geography paths) into a two-mode network based on individual and activity-travel event. Secondly, based on the spatiotemporal relationships among the behavioral events formed under different narrative modes (e.g., activity, trip, tour, sequence, and a composite of previous four narrative modes), a two-mode network is transformed into an individual-based one-mode network by constructing the adjacency matrix between individuals. Figure 1 shows the process of transforming spatiotemporal activity-travel trajectory into individual-by-individual network based on four different narrative modes and their spatiotemporal connections.
Specifically, Equation (1) represents the calculation of individual-by-individual adjacency matrix based on the shared duration of the conjoint events between individuals i and j, as defined in Equation (2), and the spatial distance between the two events as defined in Equation (3). By Equations (1)-(3), the space-time trajectory is transformed into the adjacency matrix of the network between individuals, the element value of the adjacency matrix ranges from 0 to 1. The higher the value, the higher the similarity between two individuals.
, , , where , is the k-th (k = 1…, K) common conjoint event of individual i and j, Given the narrative mode n and K common conjoint event shared by each pair individuals i and j, we can compute the adjacency matrix = ∑ ( , ) . The above equations are subjective to the calculation processes when n is a narrative mode of activity, trip or tour events. When n is a sequence, the calculation equations have subtle changes: the time relationship is no longer the overlapping time of the k-th common conjoint event, but the minimum duration of the k-th common conjoint event in the similar sequence. When n is a composite event, the aggregate adjacency matrix of four narrative perspectives is calculated according to the specific weighted index of each narrative perspective. In summary, the adjacency matrix by individual is determined by two parameters: the type of narrative mode n and the spatial interaction parameter . In this study, we focused on the activity event, and the parameter is set at 0, because this study mainly looks at interpersonal and intrapersonal variabilities between different activity-travel patterns. If we consider the spatial interaction, it can divide individuals of the same activity-travel pattern into several groups in different locations, but it is not the focus of this study.
Finally, based on the adjacency matrix derived above, we used the Louvain method [54] to partition the individual-by-individual network into communities, in order to maximize the quality score, or the modularity Q, as defined below: where represents the weight of linkage strength between individual nodes i and j in the network, = ∑ is the sum of the edge weights attached to vertex i, is the community to which the individual i is assigned, the δ function δ(u, v) is 1 if u = v and 0 otherwise and = ∑ . Intuitively, this modularity-maximization method aims to search for the best partition or community property, in which the interconnections or similarities inside the communities are much larger than those between communities. While there are many modularity-maximization methods, this study selected the Louvain method because (1) its algorithm is fast and highly efficient, particularly for a large data set; and (2) it is non-parametric and unsupervised, with the number of partitions determined by the algorithm [55][56][57][58].
After the analysis of community detection, we can derive an outcome of partition (i.e., a set of communities with similar activity-travel patterns). By using such a clustering outcome, we can simply tabulate the day-to-day variability and check whether the detected communities (i.e., groups of activity-travel patterns) vary across days within a week (e.g., between the weekday and weekend). Furthermore, we developed a multilevel model, as introduced in the next section, to estimate whether the intrapersonal and interpersonal variabilities in activity-travel patterns are significant in statistics.

Multilevel Multinomial Logit Model
In this section we developed a two-level multinomial logit model to estimate the intrapersonal variability (i.e., whether the detected category of activity-travel patterns varies across days in a week) and interpersonal variability (i.e., whether the detected category varies with individual's socioeconomic characteristics). We used a multinomial-logit form of model because the dependent variable is the class of detected communities or patterns, which is a multiple categorical variable [59,60]. Furthermore, we developed a two-level regression model, because the independent variables are measured at two levels, with day variables nested within each individual. Table 1 defines and describes both dependent and independent variables used in the model. The specification of a two-level multinomial logit model is written as follows: Level 2: = + (Individual Attriibutes) + Here, in Level-1′s equation, is the odd ratio of an individual j has a Type-m activity-travel pattern on Day i (i=Monday, …, Sunday). m=1…, M-1, given there are M-type of communities detected by the previous community detection analysis and the Type-M community is set as the reference class. (Day) is a dummy day variable, and (Day) equals one when the day is i, otherwise zero.
is the constant; is the coefficient of a Day variable; is the error term. At Level-2, we further assume the constant terms vary with different individuals and are associated with individual socioeconomic attributes. The model is finally estimated by using the maximum likelihood estimation method.  1 In fact, 9 types of communities are found in the case study in the next section. However, considering that community size affects representativeness, 7 types of communities with community size over 100 are selected as the dependent variable. 2 Hukou is a system of household registration used in mainland China, which officially identifies a person as a resident of a district.

Data and Variables
The research data comes from a one-week activity-travel diary survey in the Shangdi-Qinghe area of Beijing from September to December 2012. This survey recorded individual residents' activity and travel trajectory data within seven consecutive days by wearable GPS devises, in combination with an interactive website survey which collected residents' socioeconomic attributes and activity diary for each day of the week. The questionnaire of the activity diary includes the investigation of starting and ending time of each activity, activity details, travel mode, activity and travel companion(s), activity-travel flexibility, etc.
Finally, the survey collected the data of 680 residents from 23 neighborhoods and 19 representative companies, including 456 neighborhood samples and 224 company samples. Although the residents were asked to report their diaries of all the seven days in a week, some of them only recalled several days less than a week. After cleaning the missing data, we finally have 3793 person days for the analysis. Table 1 shows the description and statistics of intrapersonal and interpersonal characteristics. Figure 2 is a time-geography representation of the trajectory data. The horizontal coordinates represent two-dimensional space (such as longitude and latitude), while the vertical coordinates represent the temporal axis. Each trajectory displays a person's activity travel within seven days, and colored lines represent different types of activity events. This study focuses on one-week individual spatial-temporal activity-travel pattern.

Activity-Travel Patterns and Their Variabilities Across Days in a Week
The network analysis approach divided the daily trajectories of 680 residents into nine communities or clusters, with a modularity value of 0.44. As defined in Equation (4), this value represents the quality of the partition detected by the community detection algorithm (i.e., the Louvain method here). When the value is between 0.3 and 0.7, the partition is deemed as well-detected [61]. Because two communities only have a small number of person-day trajectories, we focused on the rest of 7 communities with more than 100 trajectories in the communities. Figure 3 visualizes varying activity patterns by day of the seven communities, including the plots of daily activity sequence (Figure 3a), sequence index plot (b), state distribution graph (c), polar area chart (d), and trip pair distribution (e). Colors in these plots indicate 10 types of daily activities, including at-home activities, housework, meals, shopping, working, school, leisure (including exercise, entertainment, and travelling), private affair (including family or friend visiting, running errands, and taking care of the elderly and children), medical, and other activities.
As shown in Figure 3, different communities have varying activity-travel patterns. Among them, the activity-travel pattern as shown in Community 1 is a dominant type of pattern, consisting of 2282 person-days of trajectories and accounting for 60.16% of the total observations. Residents of the Community 1 have a typical day dominant by working activities (hereinafter referred to the working pattern). This pattern is characterized by working activities in the daytime between the breakfast and lunch and between the lunch and the dinner, along with leisure and recreational activities after 6 p.m.
By contrast, the primary activity-travel pattern in Community 2 is the private affairoriented pattern (i.e., the private-affair pattern), containing leisure and private affairs activities (both peaked after 6 p.m.) except for working and having meals. The dominant activity of the pattern is private affair. This Community has 390 person days, comprising 10.28% of the total observations. Community 3 is primarily the leisure-oriented pattern (i.e., the leisure pattern), which is characterized by multiple types of activities and dominated by leisure, recreational (peaked at 9 p.m.) and shopping (peaked at 3 p.m.) activities. This community contains 269 person days, accounting for 7.09% of the total observations. In addition, Community 4 is dominated by working and housework activities (i.e., the working-housework pattern). Communities 5-7 have a relatively lower percentage of person-day trajectories, with a significant pattern dominated by other activities, school-related activities, and medical activities, respectively.
We further calculated the proportional distribution of the 7 activity-travel patterns detected across a week, as shown in Table 2. The distribution significantly varies across days of a week. For example, the working pattern as detected in Community 1 is not distributed evenly from Monday to Sunday, but more likely to occur on the weekday while less likely on the weekend. By contrast, the private-affair pattern, the leisure pattern and the school-related pattern are more likely to occur on the weekend. As to the leisure pattern, the proportional gap between the weekday and the weekend is significant. In addition, there is a lower probability to see the private-affair pattern on Thursday, the leisure pattern on Tuesday, and the school pattern on Monday. Similarly, the occurrence likelihood of the other-activities pattern is the lowest on Saturday but the highest on Sunday. The occurrence likelihood of the medical-activity pattern appears periodic, with relatively higher proportions on Monday, Wednesday, Friday, and Sunday. These findings demonstrate that the occurrence likelihood of specific activity-travel patterns may vary between weekday and weekend, as well as between the day of weekday/weekend. In addition, we calculated the proportion of different types of activity-travel patterns during the weekday (five consecutive days), the weekend (two consecutive days), and a full week (seven consecutive days), to see to what extent individuals present more than one pattern in different time periods. In the weekday, the proportion of individuals with two types of activity-travel patterns is the highest (45.78%), followed by one type (33.73%) and three different patterns (18.37%), respectively. Less than 3% individuals have more than three types of activity-travel patterns in the weekday. In the weekend, more individuals (53.24%) have two types of patterns. For an entire week, two types (40.82%) of activity-travel patterns are the most found among the samples, followed by the three types (29.08%).  Table 3 shows the estimated effects of intrapersonal (i.e., the observed day of a week) and interpersonal (individual socioeconomics) variations on the occurrence likelihood of the seven activity-travel patterns, based on the multilevel multinomial logit model. Based on the significance test, we can find that not all of activity-travel patterns have significant intrapersonal or day-to-day variability, such as the pattern of working-housework, while some significant day-to-day variations exist in the patterns of private-affair, leisure, school, other-activities and medical-activities. Particularly, compared to Monday, the patterns of private-affair, leisure, and school-related are more likely to occur on the weekend, and the patterns of other-activities and medical-activity are more likely to occur on Sunday. On the one hand, these statistical findings generally correspond to the descriptive analysis in last section. The assumption of using a typical day to measure the activitytravel pattern may result in significant biased findings. These findings are consistent with some existing studies (e.g., [2,10]). On the other hand, the activity-travel patterns have no significant variations in the weekday, but significant variations in the weekend, after controlling for socioeconomic features. This finding provides a hint for the selection of days of a week to conduct a diary survey of activity pattern or travel demand. In general, if researchers want to estimate the activity-travel pattern or the travel demand in the research area, the selection of at least three days within a week as the typical days is necessary: one from the weekday and two of the weekend. To investigate the interpersonal variability of activity-travel patterns, we need to look at the individual socioeconomic correlates of the occurrence likelihood of a specific pattern. Findings suggest that after controlling for the intrapersonal variations, the occurrence likelihoods of most activity-travel patterns vary between individuals; the interpersonal variability prevalently exists. For example, the married individuals are more likely to have private-affair pattern, even compared to the occurrence of working pattern. Furthermore, the groups of the male, aged between 40 and 50, medium education level and employed status have lower probabilities to present a private-affair pattern, compared to the other groups. Similarly, the occurrence of the leisure pattern significantly varies with gender, employment status, income level, and marital status. Interestingly, the workinghousework pattern is preferable for people with age between 30 and 50, while the medicalactivity pattern is preferable for people with age above 50, hukou in Beijing and higheducation level. These findings also demonstrate that the interpersonal variability varies between the activity-travel pattern.

Discussion and Conclusions
This paper investigates the intrapersonal and interpersonal variabilities in activitytravel pattern mining. We developed a method incorporating a network analysis approach (e.g., the community detection algorithm) with a multilevel multinomial logit model to measure and estimate intrapersonal and interpersonal variabilities. Particularly, we firstly used the network analysis approach to detect activity-travel patterns within a week of 680 residents in Shangdi-Qinghe area, Beijing. We then adopted a multilevel multinomial logit model to analyze the day-to-day variability of different patterns and their associations with days of the week and individual's socioeconomic attributes.
In this study we transformed the person-day trajectory data into network data and then used the Louvain algorithm to partition the trajectories into clusters or communities. The algorithm detects seven large communities of activity-travel patterns, which are characterized as the working pattern, private-affair pattern, leisure pattern, working-housework pattern, other-activities pattern, school pattern, and medical-activity pattern. Modeling findings suggest that most activity-travel patterns have significant intrapersonal or day-to-day variability while some have not, such as the working-housework pattern. The occurrence probabilities of most activity-travel patterns vary between weekday and weekend but have no significant difference between Monday to Friday. The patterns of privateaffair, leisure and school-related are more likely to occur in the weekend, while the patterns of other-activities and medical-activity are preferable on Saturday. These findings suggest that application of a typical day of activity-travel behaviors to represent a week's or even longer-term behaviors may be biased, due to the existence of day-to-day intrapersonal variability. Further, finding provides a hint for the selection of days of a week to conduct a diary survey of activity pattern or travel demand.
In addition, the multilevel multinomial logit modeling results reveal significant interpersonal variations in activity-travel patterns. Findings suggest that women have a higher level of interpersonal variabilities than men, whose activity-travel patterns tend to be more diversified, and the female's time-allocation differences of different patterns are larger, which may be related to their social constraints.
This study also demonstrates the advantage of an incorporation of network analysis approach and multilevel multinomial logit model to spatiotemporal behavior pattern mining. This approach provides a feasible framework for integrating interpersonal and intrapersonal variabilities, which yielding new insights into activity travel pattern and linkage to the individual's socioeconomic attributes. The analytical framework thus helps researchers and policymakers better understand different activity-travel patterns and the underlying variabilities.
There are some limitations that deserve further investigations. First, although this study relying on a week's activity-travel data can reveal both intrapersonal and interpersonal variabilities in activity-travel patterns, it fails to capture the robust variabilities which often occur in the longer term. Future work should try to collect a longer-term dataset, for example, by using one-month or longer-period mobile phone signaling data. By comparing behavioral data in the short run and long run, we can further examine whether the individual residents have seasonal variations in activity-travel arrangement. Second, because this study only focused on a typical suburban area in Beijing, modeling findings might be biased. Future studies should investigate and compare other areas in Beijing and even other cities around the world.

Data Availability Statement:
The data is available from the authors upon reasonable request.