pacman::p_load(tmap, sf, sp,
performance,
ggpubr,
tidyverse, reshape2,
DT, knitr, kableExtra)
1 Overview
Spatial interaction modelling (or gravity models) is one of the most widely used analytical tools in studying interactions between social and economic agents observed in geographical space.
Spatial Interaction Models (SIMs) are mathematical models for estimating flows between spatial entities developed by Alan Wilson in the late 1960s and early 1970, with considerable uptake and refinement for transport modelling since then Boyce and Williams (2015). Spatial interaction models, in general terms, describe and explain flow or movement between places, based on (1) their spatial separation; (2) their complementarity; and (3) other intervening opportunities or spatial structural elements serve to augment or diminish the expected flow.
There are four main types of traditional SIMs:
- Unconstrained
- Production-constrained
- Attraction-constrained
- Doubly-constrained
Ordinary least square (OLS), log-normal, Poisson and negative binomial (NB) regression methods have been used extensively to calibrate OD flow models by processing flow data as different types of dependent variables. In this exercise, we will explore how to calibrate SIM by using the four regression methods in R environment.
2 The Packages
| Package | Usage |
|---|---|
| sf | For importing, integrating, processing and transforming geospatial data |
| tmap | Choropleth mapping; thematic maps |
| sp | For handling spatial objects |
| performance | For model performance measurement |
| ggpubr | For visualisation |
| tidyverse, reshape2 | For importing, integrating, wrangling and visualising data |
| DT, knitr, kableExtra | For building tables |
3 The Data
This exercise is a continuation of Hands-on Exercise 9a: Processing and Visualising Flow Data and the following data will be used:
| Name | Details |
|---|---|
| od_data.rds | Weekday morning peak passenger flows at planning subzone level. |
| mpsz.rds | URA Master Plan 2019 Planning Subzome boundary in simple feature tibble data frame format. |
mpsz <- read_rds("data/rds/mpsz.rds")
mpszSimple feature collection with 328 features and 6 fields
Geometry type: MULTIPOLYGON
Dimension: XY
Bounding box: xmin: 2667.538 ymin: 21448.47 xmax: 55941.94 ymax: 50256.33
Projected CRS: SVY21 / Singapore TM
First 10 features:
SUBZONE_N SUBZONE_C PLN_AREA_N PLN_AREA_C REGION_N
1 MARINA EAST MESZ01 MARINA EAST ME CENTRAL REGION
2 INSTITUTION HILL RVSZ05 RIVER VALLEY RV CENTRAL REGION
3 ROBERTSON QUAY SRSZ01 SINGAPORE RIVER SR CENTRAL REGION
4 JURONG ISLAND AND BUKOM WISZ01 WESTERN ISLANDS WI WEST REGION
5 FORT CANNING MUSZ02 MUSEUM MU CENTRAL REGION
6 MARINA EAST (MP) MPSZ05 MARINE PARADE MP CENTRAL REGION
7 SENTOSA SISZ01 SOUTHERN ISLANDS SI CENTRAL REGION
8 CITY TERMINALS BMSZ17 BUKIT MERAH BM CENTRAL REGION
9 ANSON DTSZ10 DOWNTOWN CORE DT CENTRAL REGION
10 STRAITS VIEW SVSZ01 STRAITS VIEW SV CENTRAL REGION
REGION_C geometry
1 CR MULTIPOLYGON (((33222.98 29...
2 CR MULTIPOLYGON (((28481.45 30...
3 CR MULTIPOLYGON (((28087.34 30...
4 WR MULTIPOLYGON (((14557.7 304...
5 CR MULTIPOLYGON (((29542.53 31...
6 CR MULTIPOLYGON (((35279.55 30...
7 CR MULTIPOLYGON (((26879.04 26...
8 CR MULTIPOLYGON (((27891.15 28...
9 CR MULTIPOLYGON (((29201.07 28...
10 CR MULTIPOLYGON (((31269.21 28...od <- read_rds("data/rds/od_data.rds")
od# A tibble: 20,849 × 3
# Groups: ORIGIN_SZ [309]
ORIGIN_SZ DESTIN_SZ MORNING_PEAK
<chr> <chr> <dbl>
1 AMSZ01 AMSZ01 2435
2 AMSZ01 AMSZ02 10418
3 AMSZ01 AMSZ03 15578
4 AMSZ01 AMSZ04 2862
5 AMSZ01 AMSZ05 7938
6 AMSZ01 AMSZ06 2357
7 AMSZ01 AMSZ07 1513
8 AMSZ01 AMSZ08 2551
9 AMSZ01 AMSZ09 2358
10 AMSZ01 AMSZ10 291
# ℹ 20,839 more rows4 Computing Distance Matrix
In spatial interaction, a distance matrix is a table that shows the distance between pairs of locations. In this section, we will conduct additional data preparation necessary to compute distance matrix to run spatial interaction modelling.
4.1 Converting from sf data.table to SpatialPolygonsDataFrame
There are at least two ways to compute the required distance matrix. One is based on sf and the other is based on sp.
Past experience shown that computing distance matrix by using sf function took relatively longer time that sp method especially the data set is large. In view of this, sp method is used in the code chunks below.
First as.Spatial() will be used to convert mpsz from sf tibble data frame to SpatialPolygonsDataFrame of sp object as shown in the code chunk below.
mpsz_sp <- as(mpsz, "Spatial")
mpsz_spclass : SpatialPolygonsDataFrame
features : 328
extent : 2667.538, 55941.94, 21448.47, 50256.33 (xmin, xmax, ymin, ymax)
crs : +proj=tmerc +lat_0=1.36666666666667 +lon_0=103.833333333333 +k=1 +x_0=28001.642 +y_0=38744.572 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs
variables : 6
names : SUBZONE_N, SUBZONE_C, PLN_AREA_N, PLN_AREA_C, REGION_N, REGION_C
min values : ADMIRALTY, AMSZ01, ANG MO KIO, AM, CENTRAL REGION, CR
max values : YUNNAN, YSSZ09, YISHUN, YS, WEST REGION, WR 4.2 Computing Distance Matrix
Next, spDists() of sp package will be used to compute the Euclidean distance between the centroids of the planning subzones.
dist <- spDists(mpsz_sp,
longlat = FALSE)head(dist, n=c(10, 10)) [,1] [,2] [,3] [,4] [,5] [,6] [,7]
[1,] 0.000 3926.0025 3939.108 20252.96 2989.9839 1431.330 6391.342
[2,] 3926.003 0.0000 305.737 16513.86 951.8314 5254.066 4975.002
[3,] 3939.108 305.7370 0.000 16412.06 1045.9088 5299.849 4669.295
[4,] 20252.964 16513.8648 16412.062 0.00 17450.3044 21665.795 15469.566
[5,] 2989.984 951.8314 1045.909 17450.30 0.0000 4303.232 5226.873
[6,] 1431.330 5254.0664 5299.849 21665.79 4303.2323 0.000 7707.091
[7,] 6391.342 4975.0021 4669.295 15469.57 5226.8731 7707.091 0.000
[8,] 4716.445 3176.1592 2873.497 16142.96 3341.2116 6103.071 1893.049
[9,] 3593.009 2345.1741 2074.691 16922.15 2264.2014 5007.197 3068.627
[10,] 2434.850 3455.5791 3277.921 18597.24 2890.8696 3699.242 4009.437
[,8] [,9] [,10]
[1,] 4716.445 3593.009 2434.850
[2,] 3176.159 2345.174 3455.579
[3,] 2873.497 2074.691 3277.921
[4,] 16142.960 16922.153 18597.236
[5,] 3341.212 2264.201 2890.870
[6,] 6103.071 5007.197 3699.242
[7,] 1893.049 3068.627 4009.437
[8,] 0.000 1200.264 2532.383
[9,] 1200.264 0.000 1709.443
[10,] 2532.383 1709.443 0.000Notice that the output dist is a matrix object class of R. Also notice that the column heanders and row headers are not labeled with the planning subzone codes. We need to carry out additional data preparation steps to tidy this data up.
4.3 Labelling Column and Row Headers of Distance Matrix
First, we will create a list sorted according to the the distance matrix by planning sub-zone code.
sz_names <- mpsz$SUBZONE_CNext we will attach SUBZONE_C to row and column for distance matrix matching ahead
colnames(dist) <- paste0(sz_names)
rownames(dist) <- paste0(sz_names)4.4 Pivoting Distance Value by SUBZONE_C
Next, we will pivot the distance matrix into a long table by using the row and column subzone codes as show in the code chunk below.
distPair <- melt(dist) %>%
rename(dist = value)
head(distPair, 10) Var1 Var2 dist
1 MESZ01 MESZ01 0.000
2 RVSZ05 MESZ01 3926.003
3 SRSZ01 MESZ01 3939.108
4 WISZ01 MESZ01 20252.964
5 MUSZ02 MESZ01 2989.984
6 MPSZ05 MESZ01 1431.330
7 SISZ01 MESZ01 6391.342
8 BMSZ17 MESZ01 4716.445
9 DTSZ10 MESZ01 3593.009
10 SVSZ01 MESZ01 2434.850Notice that the within zone distance is 0.
4.5 Updating Intra-zonal Distances
In this section, we are going to append a constant value to replace the intra-zonal distance of 0. First, we will select and find out the minimum value of the distance by using summary().
distPair %>%
filter(dist > 0) %>%
summary() Var1 Var2 dist
MESZ01 : 327 MESZ01 : 327 Min. : 173.8
RVSZ05 : 327 RVSZ05 : 327 1st Qu.: 7019.9
SRSZ01 : 327 SRSZ01 : 327 Median :11731.2
WISZ01 : 327 WISZ01 : 327 Mean :12029.5
MUSZ02 : 327 MUSZ02 : 327 3rd Qu.:16203.5
MPSZ05 : 327 MPSZ05 : 327 Max. :44892.8
(Other):105294 (Other):105294 Next, a constant distance value of 50m is added into intra-zones distance.
distPair$dist <- ifelse(distPair$dist == 0,
50, distPair$dist)
distPair %>%
summary() Var1 Var2 dist
MESZ01 : 328 MESZ01 : 328 Min. : 50
RVSZ05 : 328 RVSZ05 : 328 1st Qu.: 6969
SRSZ01 : 328 SRSZ01 : 328 Median :11706
WISZ01 : 328 WISZ01 : 328 Mean :11993
MUSZ02 : 328 MUSZ02 : 328 3rd Qu.:16189
MPSZ05 : 328 MPSZ05 : 328 Max. :44893
(Other):105616 (Other):105616 The code chunk below is used to rename the origin and destination fields.
distPair <- distPair %>%
rename(orig = Var1,
dest = Var2)Lastly, the code chunk below is used to save the dataframe into rds file format
write_rds(distPair, "data/rds/distPair.rds") 5 Preparing the Flow Data
Firstly, we will compute the total passenger trip between and within planning subzones by using the code chunk below. The output is called flow_data.
flow_data <- od %>%
group_by(ORIGIN_SZ, DESTIN_SZ) %>%
summarize(TRIPS = sum(MORNING_PEAK)) Let’s look at some of the data rows.
head(flow_data, 10)# A tibble: 10 × 3
# Groups: ORIGIN_SZ [1]
ORIGIN_SZ DESTIN_SZ TRIPS
<chr> <chr> <dbl>
1 AMSZ01 AMSZ01 2435
2 AMSZ01 AMSZ02 10418
3 AMSZ01 AMSZ03 15578
4 AMSZ01 AMSZ04 2862
5 AMSZ01 AMSZ05 7938
6 AMSZ01 AMSZ06 2357
7 AMSZ01 AMSZ07 1513
8 AMSZ01 AMSZ08 2551
9 AMSZ01 AMSZ09 2358
10 AMSZ01 AMSZ10 2915.1 Separating intra-flow from passenger volume df
Code chunk below is used to add two new fields in flow_data dataframe namely FlowNoIntra and offset.
flow_data$FlowNoIntra <- ifelse(
flow_data$ORIGIN_SZ == flow_data$DESTIN_SZ,
0,
flow_data$TRIPS)
flow_data$offset <- ifelse(
flow_data$ORIGIN_SZ == flow_data$DESTIN_SZ,
0.000001,
1)5.2 Combining passenger volume data with distance value
Before we can join flow_data and distPair, we need to convert data value type of ORIGIN_SZ and DESTIN_SZ fields of flow_data dataframe into factor data type.
flow_data$ORIGIN_SZ <- as.factor(flow_data$ORIGIN_SZ)
flow_data$DESTIN_SZ <- as.factor(flow_data$DESTIN_SZ)Now, left_join() of dplyr will be used to flow_data dataframe and distPair dataframe. The output is called flow_data1.
flow_data1 <- flow_data %>%
left_join (distPair,
by = c("ORIGIN_SZ" = "orig",
"DESTIN_SZ" = "dest"))6 Preparing Origin and Destination Attributes
6.1 Importing population data
Firstly, we will import the population data.
pop <- read_csv("data/aspatial/pop.csv")Next, we will do a left_join to pop data frame with mpsz. The output will be a sf object where each polygon in mpsz will be assigned a population value.
pop <- pop %>%
left_join(mpsz,
by = c("PA" = "PLN_AREA_N",
"SZ" = "SUBZONE_N")) %>%
select(1:6) %>%
rename(SZ_NAME = SZ,
SZ = SUBZONE_C)Next, we will need to do another left_join() with flow_data1 that we have prepared earlier to prepare both origin and destination attributes.
flow_data1 <- flow_data1 %>%
left_join(pop,
by = c(ORIGIN_SZ = "SZ")) %>% #<< ORIGIN_SZ
rename(ORIGIN_AGE7_12 = AGE7_12,
ORIGIN_AGE13_24 = AGE13_24,
ORIGIN_AGE25_64 = AGE25_64) %>%
select(-c(PA, SZ_NAME))flow_data1 <- flow_data1 %>%
left_join(pop,
by = c(DESTIN_SZ = "SZ")) %>% #<< DESTIN_SZ
rename(DESTIN_AGE7_12 = AGE7_12,
DESTIN_AGE13_24 = AGE13_24,
DESTIN_AGE25_64 = AGE25_64) %>%
select(-c(PA, SZ_NAME))We will called the output data file SIM_data. it is in rds data file format.
write_rds(flow_data1, "data/rds/flow_data_6-9.rds")
SIM_data <- read_rds("data/rds/flow_data_6-9.rds")7 Calibrating Spatial Interaction Models
In this section, we explore how to calibrate Spatial Interaction Models by using Poisson Regression method.
7.1 Visualising the Dependent Variables
Firstly, let us plot the distribution of the dependent variable (i.e. TRIPS) by using histogram method by using the code chunk below.
Show the code
ggplot(data = SIM_data,
aes(x = TRIPS)) +
geom_histogram() +
theme(
plot.title = element_text(face= 'bold', size = 10),
panel.grid.major = element_line(colour = "#ede5de", linetype = 1, linewidth = 0.5),
panel.grid.minor = element_line(colour = "#ede5de", linetype = 1, linewidth= 0.5),
plot.background = element_rect(fill="#E4D5C9",colour="#E4D5C9"),
panel.border = element_blank(),
panel.background = element_blank(),
axis.ticks = element_blank(),
axis.title.y = element_text(hjust=1, angle=0, size = 8),
axis.title.x = element_text(size = 8),
strip.text = element_text(face= 'bold'),
strip.background = element_rect(color="#E4D5C9", fill="#E4D5C9")
)
Notice that the distribution is highly skewed and not resemble bell shape or also known as normal distribution.
Next, let us visualise the relation between the dependent variable and one of the key independent variable in Spatial Interaction Model, namely distance.
Show the code
ggplot(data = SIM_data,
aes(x = dist,
y = TRIPS)) +
geom_point() +
geom_smooth(method = lm) +
theme(
plot.title = element_text(face= 'bold', size = 10),
panel.grid.major = element_line(colour = "#ede5de", linetype = 1, linewidth = 0.5),
panel.grid.minor = element_line(colour = "#ede5de", linetype = 1, linewidth= 0.5),
plot.background = element_rect(fill="#E4D5C9",colour="#E4D5C9"),
panel.border = element_blank(),
panel.background = element_blank(),
axis.ticks = element_blank(),
axis.title.y = element_text(hjust=1, angle=0, size = 8),
axis.title.x = element_text(size = 8),
strip.text = element_text(face= 'bold'),
strip.background = element_rect(color="#E4D5C9", fill="#E4D5C9")
)
Notice that their relationship hardly resemble linear relationship.
On the other hand, if we plot the scatter plot by using the log transformed version of both variables, we can see that their relationship is more resemble linear relationship.
Show the code
ggplot(data = SIM_data,
aes(x = log(dist),
y = log(TRIPS))) +
geom_point() +
geom_smooth(method = lm) +
theme(
plot.title = element_text(face= 'bold', size = 10),
panel.grid.major = element_line(colour = "#ede5de", linetype = 1, linewidth = 0.5),
panel.grid.minor = element_line(colour = "#ede5de", linetype = 1, linewidth= 0.5),
plot.background = element_rect(fill="#E4D5C9",colour="#E4D5C9"),
panel.border = element_blank(),
panel.background = element_blank(),
axis.ticks = element_blank(),
axis.title.y = element_text(hjust=1, angle=0, size = 8),
axis.title.x = element_text(size = 8),
strip.text = element_text(face= 'bold'),
strip.background = element_rect(color="#E4D5C9", fill="#E4D5C9")
)
7.2 Log Transformation
Since Poisson Regression is based of log and log 0 is undefined, it is important for us to ensure that no 0 values in the explanatory variables.
In the code chunk below, summary() of Base R is used to compute the summary statistics of all variables in SIM_data data frame.
summary(SIM_data) ORIGIN_SZ DESTIN_SZ TRIPS FlowNoIntra
Length:20849 Length:20849 Min. : 1 Min. : 0.0
Class :character Class :character 1st Qu.: 16 1st Qu.: 15.0
Mode :character Mode :character Median : 85 Median : 79.0
Mean : 1149 Mean : 954.7
3rd Qu.: 438 3rd Qu.: 403.0
Max. :339548 Max. :216167.0
offset dist ORIGIN_AGE7_12 ORIGIN_AGE13_24
Min. :0.000001 Min. : 50 Min. : 0.0 Min. : 0
1st Qu.:1.000000 1st Qu.: 3361 1st Qu.: 50.0 1st Qu.: 100
Median :1.000000 Median : 6155 Median : 510.0 Median : 1130
Mean :0.986043 Mean : 6979 Mean : 887.6 Mean : 1952
3rd Qu.:1.000000 3rd Qu.: 9884 3rd Qu.:1360.0 3rd Qu.: 2950
Max. :1.000000 Max. :26136 Max. :6340.0 Max. :16380
ORIGIN_AGE25_64 DESTIN_AGE7_12 DESTIN_AGE13_24 DESTIN_AGE25_64
Min. : 0 Min. : 0.0 Min. : 0 Min. : 0
1st Qu.: 730 1st Qu.: 10.0 1st Qu.: 60 1st Qu.: 630
Median : 5730 Median : 510.0 Median : 1100 Median : 5710
Mean : 9088 Mean : 852.4 Mean : 1890 Mean : 8805
3rd Qu.:14180 3rd Qu.:1200.0 3rd Qu.: 2920 3rd Qu.:13650
Max. :74610 Max. :6340.0 Max. :16380 Max. :74610 The print report above reveals that variables ORIGIN_AGE7_12, ORIGIN_AGE13_24, ORIGIN_AGE25_64, DESTIN_AGE7_12, DESTIN_AGE13_24, DESTIN_AGE25_64 consist of 0 values.
In view of this, code chunk below will be used to replace zero values to 0.99.
SIM_data$DESTIN_AGE7_12 <- ifelse(
SIM_data$DESTIN_AGE7_12 == 0,
0.99,
SIM_data$DESTIN_AGE7_12)
SIM_data$DESTIN_AGE13_24 <- ifelse(
SIM_data$DESTIN_AGE13_24 == 0,
0.99,
SIM_data$DESTIN_AGE13_24)
SIM_data$DESTIN_AGE25_64 <- ifelse(
SIM_data$DESTIN_AGE25_64 == 0,
0.99,
SIM_data$DESTIN_AGE25_64)
SIM_data$ORIGIN_AGE7_12 <- ifelse(
SIM_data$ORIGIN_AGE7_12 == 0,
0.99,
SIM_data$ORIGIN_AGE7_12)
SIM_data$ORIGIN_AGE13_24 <- ifelse(
SIM_data$ORIGIN_AGE13_24 == 0,
0.99,
SIM_data$ORIGIN_AGE13_24)
SIM_data$ORIGIN_AGE25_64 <- ifelse(
SIM_data$ORIGIN_AGE25_64 == 0,
0.99,
SIM_data$ORIGIN_AGE25_64)Let’s summarise the new data:
summary(SIM_data) ORIGIN_SZ DESTIN_SZ TRIPS FlowNoIntra
Length:20849 Length:20849 Min. : 1 Min. : 0.0
Class :character Class :character 1st Qu.: 16 1st Qu.: 15.0
Mode :character Mode :character Median : 85 Median : 79.0
Mean : 1149 Mean : 954.7
3rd Qu.: 438 3rd Qu.: 403.0
Max. :339548 Max. :216167.0
offset dist ORIGIN_AGE7_12 ORIGIN_AGE13_24
Min. :0.000001 Min. : 50 Min. : 0.99 Min. : 0.99
1st Qu.:1.000000 1st Qu.: 3361 1st Qu.: 50.00 1st Qu.: 100.00
Median :1.000000 Median : 6155 Median : 510.00 Median : 1130.00
Mean :0.986043 Mean : 6979 Mean : 887.84 Mean : 1952.55
3rd Qu.:1.000000 3rd Qu.: 9884 3rd Qu.:1360.00 3rd Qu.: 2950.00
Max. :1.000000 Max. :26136 Max. :6340.00 Max. :16380.00
ORIGIN_AGE25_64 DESTIN_AGE7_12 DESTIN_AGE13_24 DESTIN_AGE25_64
Min. : 0.99 Min. : 0.99 Min. : 0.99 Min. : 0.99
1st Qu.: 730.00 1st Qu.: 10.00 1st Qu.: 60.00 1st Qu.: 630.00
Median : 5730.00 Median : 510.00 Median : 1100.00 Median : 5710.00
Mean : 9087.96 Mean : 852.59 Mean : 1890.03 Mean : 8805.49
3rd Qu.:14180.00 3rd Qu.:1200.00 3rd Qu.: 2920.00 3rd Qu.:13650.00
Max. :74610.00 Max. :6340.00 Max. :16380.00 Max. :74610.00 Notice that all the 0 values have been replaced by 0.99.
7.3 Unconstrained Spatial Interaction Model
In this section, we will calibrate an unconstrained spatial interaction model by using glm() function. The explanatory variables are origin population by different age cohort, destination population by different age cohort (i.e. ORIGIN_AGE25_64) and distance between origin and destination in km (i.e. dist).
uncSIM <- glm(formula = TRIPS ~
log(ORIGIN_AGE25_64) +
log(DESTIN_AGE25_64) +
log(dist),
family = poisson(link = "log"),
data = SIM_data,
na.action = na.exclude)
uncSIM
Call: glm(formula = TRIPS ~ log(ORIGIN_AGE25_64) + log(DESTIN_AGE25_64) +
log(dist), family = poisson(link = "log"), data = SIM_data,
na.action = na.exclude)
Coefficients:
(Intercept) log(ORIGIN_AGE25_64) log(DESTIN_AGE25_64)
10.51560 0.26534 0.02819
log(dist)
-0.72945
Degrees of Freedom: 20848 Total (i.e. Null); 20845 Residual
Null Deviance: 102700000
Residual Deviance: 60630000 AIC: 607700007.3.1 R-squared function
In order to measure how much variation of the trips can be accounted by the model we will write a function to calculate R-Squared value as shown below.
CalcRSquared <- function(observed,estimated){
r <- cor(observed,estimated)
R2 <- r^2
R2
}Next, we will compute the R-squared of the unconstrained SIM by using the code chunk below.
CalcRSquared(uncSIM$data$TRIPS, uncSIM$fitted.values)[1] 0.1951635r2_mcfadden(uncSIM)# R2 for Generalized Linear Regression
R2: 0.409
adj. R2: 0.4097.4 Origin Constrained Spatial Interaction Model
In this section, we will calibrate an origin constrained SIM. For origin constrained SIM, only explanatory variables representing the attractiveness at the destinations will be used. This is because such models emphasize the limitations or capacities of the origins rather than the demand or attractiveness of the destinations. The capacity or limitation at the origin sites determines the potential for generating interactions or flows.
orcSIM <- glm(formula = TRIPS ~
ORIGIN_SZ +
log(DESTIN_AGE25_64) +
log(dist),
family = poisson(link = "log"),
data = SIM_data,
na.action = na.exclude)
summary(orcSIM)
Call:
glm(formula = TRIPS ~ ORIGIN_SZ + log(DESTIN_AGE25_64) + log(dist),
family = poisson(link = "log"), data = SIM_data, na.action = na.exclude)
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 1.249e+01 3.153e-03 3963.125 < 2e-16 ***
ORIGIN_SZAMSZ02 1.059e+00 3.700e-03 286.348 < 2e-16 ***
ORIGIN_SZAMSZ03 5.866e-01 3.791e-03 154.729 < 2e-16 ***
ORIGIN_SZAMSZ04 -4.110e-02 4.266e-03 -9.634 < 2e-16 ***
ORIGIN_SZAMSZ05 -1.968e-01 4.852e-03 -40.565 < 2e-16 ***
ORIGIN_SZAMSZ06 3.212e-01 4.402e-03 72.972 < 2e-16 ***
ORIGIN_SZAMSZ07 -1.203e+00 7.440e-03 -161.703 < 2e-16 ***
ORIGIN_SZAMSZ08 -8.599e-01 6.930e-03 -124.080 < 2e-16 ***
ORIGIN_SZAMSZ09 1.910e-01 4.569e-03 41.811 < 2e-16 ***
ORIGIN_SZAMSZ10 4.601e-01 4.014e-03 114.626 < 2e-16 ***
ORIGIN_SZAMSZ11 -1.468e+00 9.049e-03 -162.172 < 2e-16 ***
ORIGIN_SZAMSZ12 -1.559e+00 9.123e-03 -170.923 < 2e-16 ***
ORIGIN_SZBDSZ01 1.002e+00 3.660e-03 273.818 < 2e-16 ***
ORIGIN_SZBDSZ02 4.894e-01 4.227e-03 115.786 < 2e-16 ***
ORIGIN_SZBDSZ03 9.607e-01 3.745e-03 256.566 < 2e-16 ***
ORIGIN_SZBDSZ04 1.711e+00 3.273e-03 522.847 < 2e-16 ***
ORIGIN_SZBDSZ05 6.961e-01 3.738e-03 186.199 < 2e-16 ***
ORIGIN_SZBDSZ06 9.371e-01 3.768e-03 248.676 < 2e-16 ***
ORIGIN_SZBDSZ07 -1.008e+00 6.985e-03 -144.324 < 2e-16 ***
ORIGIN_SZBDSZ08 -8.857e-01 6.728e-03 -131.649 < 2e-16 ***
ORIGIN_SZBKSZ01 -3.528e-01 5.370e-03 -65.706 < 2e-16 ***
ORIGIN_SZBKSZ02 4.323e-01 4.257e-03 101.557 < 2e-16 ***
ORIGIN_SZBKSZ03 7.078e-01 3.957e-03 178.873 < 2e-16 ***
ORIGIN_SZBKSZ04 -7.670e-02 4.915e-03 -15.607 < 2e-16 ***
ORIGIN_SZBKSZ05 -9.955e-02 4.705e-03 -21.159 < 2e-16 ***
ORIGIN_SZBKSZ06 2.934e-02 4.872e-03 6.021 1.73e-09 ***
ORIGIN_SZBKSZ07 7.640e-01 3.684e-03 207.356 < 2e-16 ***
ORIGIN_SZBKSZ08 6.920e-02 4.368e-03 15.844 < 2e-16 ***
ORIGIN_SZBKSZ09 4.699e-02 4.572e-03 10.277 < 2e-16 ***
ORIGIN_SZBLSZ01 -1.352e+00 1.055e-02 -128.111 < 2e-16 ***
ORIGIN_SZBLSZ02 -2.144e+00 1.587e-02 -135.047 < 2e-16 ***
ORIGIN_SZBLSZ03 -3.269e+00 3.422e-02 -95.539 < 2e-16 ***
ORIGIN_SZBLSZ04 -1.816e+00 1.618e-02 -112.215 < 2e-16 ***
ORIGIN_SZBMSZ01 1.304e-01 4.006e-03 32.562 < 2e-16 ***
ORIGIN_SZBMSZ02 -1.204e+00 5.746e-03 -209.516 < 2e-16 ***
ORIGIN_SZBMSZ03 -3.589e-01 4.572e-03 -78.504 < 2e-16 ***
ORIGIN_SZBMSZ04 1.720e-02 4.033e-03 4.263 2.01e-05 ***
ORIGIN_SZBMSZ05 -1.327e+00 6.195e-03 -214.135 < 2e-16 ***
ORIGIN_SZBMSZ06 -1.826e+00 9.733e-03 -187.649 < 2e-16 ***
ORIGIN_SZBMSZ07 -4.338e-01 4.557e-03 -95.191 < 2e-16 ***
ORIGIN_SZBMSZ08 -4.782e-01 4.532e-03 -105.519 < 2e-16 ***
ORIGIN_SZBMSZ09 -1.151e+00 5.838e-03 -197.209 < 2e-16 ***
ORIGIN_SZBMSZ10 -1.229e+00 6.220e-03 -197.635 < 2e-16 ***
ORIGIN_SZBMSZ11 -8.715e-01 5.405e-03 -161.247 < 2e-16 ***
ORIGIN_SZBMSZ12 -1.183e+00 7.269e-03 -162.746 < 2e-16 ***
ORIGIN_SZBMSZ13 -7.287e-02 4.475e-03 -16.285 < 2e-16 ***
ORIGIN_SZBMSZ14 -6.203e-01 5.321e-03 -116.562 < 2e-16 ***
ORIGIN_SZBMSZ15 -3.498e-01 4.818e-03 -72.594 < 2e-16 ***
ORIGIN_SZBMSZ16 -1.315e+00 6.208e-03 -211.798 < 2e-16 ***
ORIGIN_SZBMSZ17 -1.780e+00 9.537e-03 -186.685 < 2e-16 ***
ORIGIN_SZBPSZ01 1.962e-01 4.512e-03 43.469 < 2e-16 ***
ORIGIN_SZBPSZ02 1.751e-01 5.029e-03 34.818 < 2e-16 ***
ORIGIN_SZBPSZ03 4.120e-01 4.665e-03 88.322 < 2e-16 ***
ORIGIN_SZBPSZ04 4.980e-01 4.102e-03 121.392 < 2e-16 ***
ORIGIN_SZBPSZ05 5.585e-01 3.803e-03 146.871 < 2e-16 ***
ORIGIN_SZBPSZ06 -1.016e+00 6.738e-03 -150.756 < 2e-16 ***
ORIGIN_SZBPSZ07 -8.616e-01 6.688e-03 -128.817 < 2e-16 ***
ORIGIN_SZBSSZ01 -1.342e-02 4.384e-03 -3.061 0.0022 **
ORIGIN_SZBSSZ02 3.874e-01 3.946e-03 98.168 < 2e-16 ***
ORIGIN_SZBSSZ03 2.957e-01 3.898e-03 75.877 < 2e-16 ***
ORIGIN_SZBTSZ01 2.073e-02 4.326e-03 4.793 1.64e-06 ***
ORIGIN_SZBTSZ02 -9.811e-01 6.128e-03 -160.096 < 2e-16 ***
ORIGIN_SZBTSZ03 -9.794e-02 4.573e-03 -21.419 < 2e-16 ***
ORIGIN_SZBTSZ04 -8.325e-01 7.975e-03 -104.395 < 2e-16 ***
ORIGIN_SZBTSZ05 -1.568e+00 8.698e-03 -180.279 < 2e-16 ***
ORIGIN_SZBTSZ06 -7.266e-01 5.746e-03 -126.454 < 2e-16 ***
ORIGIN_SZBTSZ07 -1.974e+00 9.180e-03 -215.080 < 2e-16 ***
ORIGIN_SZBTSZ08 -1.178e+00 7.125e-03 -165.399 < 2e-16 ***
ORIGIN_SZCCSZ01 -1.751e+00 1.131e-02 -154.914 < 2e-16 ***
ORIGIN_SZCHSZ01 -1.183e+00 9.458e-03 -125.098 < 2e-16 ***
ORIGIN_SZCHSZ02 -6.343e-01 6.796e-03 -93.334 < 2e-16 ***
ORIGIN_SZCHSZ03 5.968e-01 4.679e-03 127.546 < 2e-16 ***
ORIGIN_SZCKSZ01 3.745e-01 4.096e-03 91.428 < 2e-16 ***
ORIGIN_SZCKSZ02 8.855e-01 4.092e-03 216.399 < 2e-16 ***
ORIGIN_SZCKSZ03 9.006e-01 3.736e-03 241.036 < 2e-16 ***
ORIGIN_SZCKSZ04 1.365e+00 3.791e-03 359.919 < 2e-16 ***
ORIGIN_SZCKSZ05 1.077e+00 4.395e-03 245.092 < 2e-16 ***
ORIGIN_SZCKSZ06 1.292e+00 4.175e-03 309.531 < 2e-16 ***
ORIGIN_SZCLSZ01 -4.418e-01 5.804e-03 -76.113 < 2e-16 ***
ORIGIN_SZCLSZ02 -1.653e+00 1.067e-02 -154.865 < 2e-16 ***
ORIGIN_SZCLSZ03 -4.454e-01 5.413e-03 -82.283 < 2e-16 ***
ORIGIN_SZCLSZ04 8.566e-01 3.622e-03 236.520 < 2e-16 ***
ORIGIN_SZCLSZ05 -1.701e+00 1.058e-02 -160.702 < 2e-16 ***
ORIGIN_SZCLSZ06 8.556e-01 3.509e-03 243.842 < 2e-16 ***
ORIGIN_SZCLSZ07 -1.857e-01 4.475e-03 -41.498 < 2e-16 ***
ORIGIN_SZCLSZ08 3.329e-01 4.867e-03 68.389 < 2e-16 ***
ORIGIN_SZCLSZ09 -1.770e+00 1.390e-02 -127.288 < 2e-16 ***
ORIGIN_SZDTSZ01 -1.730e+00 6.972e-03 -248.070 < 2e-16 ***
ORIGIN_SZDTSZ02 -1.593e+00 6.529e-03 -243.931 < 2e-16 ***
ORIGIN_SZDTSZ03 -2.794e+00 1.303e-02 -214.461 < 2e-16 ***
ORIGIN_SZDTSZ04 -4.246e+00 1.187e-01 -35.767 < 2e-16 ***
ORIGIN_SZDTSZ05 -3.123e+00 2.223e-02 -140.470 < 2e-16 ***
ORIGIN_SZDTSZ06 -2.942e+00 1.720e-02 -171.002 < 2e-16 ***
ORIGIN_SZDTSZ07 -2.082e+00 1.988e-02 -104.738 < 2e-16 ***
ORIGIN_SZDTSZ08 -2.397e+00 9.915e-03 -241.762 < 2e-16 ***
ORIGIN_SZDTSZ09 -3.023e+00 2.058e-02 -146.870 < 2e-16 ***
ORIGIN_SZDTSZ10 -2.239e+00 1.075e-02 -208.330 < 2e-16 ***
ORIGIN_SZDTSZ11 -2.373e+00 1.080e-02 -219.608 < 2e-16 ***
ORIGIN_SZDTSZ12 -3.467e+00 2.369e-02 -146.326 < 2e-16 ***
ORIGIN_SZDTSZ13 -2.423e+00 1.246e-02 -194.391 < 2e-16 ***
ORIGIN_SZGLSZ01 -1.381e+00 7.253e-03 -190.392 < 2e-16 ***
ORIGIN_SZGLSZ02 1.921e-01 4.087e-03 47.015 < 2e-16 ***
ORIGIN_SZGLSZ03 1.909e-01 4.069e-03 46.924 < 2e-16 ***
ORIGIN_SZGLSZ04 9.414e-01 3.434e-03 274.148 < 2e-16 ***
ORIGIN_SZGLSZ05 6.514e-01 3.621e-03 179.913 < 2e-16 ***
ORIGIN_SZHGSZ01 2.314e-01 4.013e-03 57.649 < 2e-16 ***
ORIGIN_SZHGSZ02 6.052e-01 3.859e-03 156.850 < 2e-16 ***
ORIGIN_SZHGSZ03 2.649e-01 4.231e-03 62.603 < 2e-16 ***
ORIGIN_SZHGSZ04 9.243e-01 3.596e-03 257.013 < 2e-16 ***
ORIGIN_SZHGSZ05 1.247e+00 3.536e-03 352.717 < 2e-16 ***
ORIGIN_SZHGSZ06 7.758e-02 4.308e-03 18.008 < 2e-16 ***
ORIGIN_SZHGSZ07 7.093e-01 3.730e-03 190.183 < 2e-16 ***
ORIGIN_SZHGSZ08 1.983e-01 4.310e-03 46.003 < 2e-16 ***
ORIGIN_SZHGSZ09 -5.470e-01 5.820e-03 -93.993 < 2e-16 ***
ORIGIN_SZHGSZ10 -3.655e+00 3.904e-02 -93.627 < 2e-16 ***
ORIGIN_SZJESZ01 3.619e-01 4.152e-03 87.174 < 2e-16 ***
ORIGIN_SZJESZ02 2.474e-01 4.139e-03 59.780 < 2e-16 ***
ORIGIN_SZJESZ03 2.392e-01 4.372e-03 54.712 < 2e-16 ***
ORIGIN_SZJESZ04 -1.194e+00 7.807e-03 -152.922 < 2e-16 ***
ORIGIN_SZJESZ05 -2.042e+00 1.203e-02 -169.679 < 2e-16 ***
ORIGIN_SZJESZ06 3.185e-01 4.061e-03 78.441 < 2e-16 ***
ORIGIN_SZJESZ07 -1.784e+00 9.420e-03 -189.428 < 2e-16 ***
ORIGIN_SZJESZ08 -6.383e-01 8.460e-03 -75.453 < 2e-16 ***
ORIGIN_SZJESZ09 4.333e-01 4.225e-03 102.547 < 2e-16 ***
ORIGIN_SZJESZ10 -2.043e+00 1.706e-02 -119.703 < 2e-16 ***
ORIGIN_SZJESZ11 -2.129e+00 1.705e-02 -124.901 < 2e-16 ***
ORIGIN_SZJWSZ01 2.500e-01 5.421e-03 46.110 < 2e-16 ***
ORIGIN_SZJWSZ02 8.518e-01 3.806e-03 223.792 < 2e-16 ***
ORIGIN_SZJWSZ03 1.247e+00 3.534e-03 352.929 < 2e-16 ***
ORIGIN_SZJWSZ04 1.331e+00 3.587e-03 371.005 < 2e-16 ***
ORIGIN_SZJWSZ05 -1.336e+00 1.064e-02 -125.601 < 2e-16 ***
ORIGIN_SZJWSZ06 -1.082e+00 8.856e-03 -122.201 < 2e-16 ***
ORIGIN_SZJWSZ07 -2.722e+00 2.313e-02 -117.725 < 2e-16 ***
ORIGIN_SZJWSZ08 1.938e+00 3.444e-03 562.704 < 2e-16 ***
ORIGIN_SZJWSZ09 1.849e+00 3.274e-03 564.679 < 2e-16 ***
ORIGIN_SZKLSZ01 1.724e-01 3.934e-03 43.818 < 2e-16 ***
ORIGIN_SZKLSZ02 -4.691e-01 4.916e-03 -95.421 < 2e-16 ***
ORIGIN_SZKLSZ03 -4.563e-01 4.913e-03 -92.879 < 2e-16 ***
ORIGIN_SZKLSZ04 -1.701e+00 7.133e-03 -238.440 < 2e-16 ***
ORIGIN_SZKLSZ05 -9.591e-01 6.692e-03 -143.325 < 2e-16 ***
ORIGIN_SZKLSZ06 -5.925e-01 4.621e-03 -128.235 < 2e-16 ***
ORIGIN_SZKLSZ07 -9.390e-01 6.052e-03 -155.149 < 2e-16 ***
ORIGIN_SZKLSZ08 -8.177e-01 5.286e-03 -154.692 < 2e-16 ***
ORIGIN_SZKLSZ09 -1.555e+00 6.666e-03 -233.195 < 2e-16 ***
ORIGIN_SZLKSZ01 -3.144e+00 3.138e-02 -100.181 < 2e-16 ***
ORIGIN_SZMDSZ01 -2.390e+00 2.165e-02 -110.385 < 2e-16 ***
ORIGIN_SZMDSZ02 -1.105e+00 9.805e-03 -112.726 < 2e-16 ***
ORIGIN_SZMDSZ03 -1.814e+00 1.333e-02 -136.054 < 2e-16 ***
ORIGIN_SZMPSZ01 -1.043e+00 6.616e-03 -157.671 < 2e-16 ***
ORIGIN_SZMPSZ02 -5.017e-01 5.396e-03 -92.981 < 2e-16 ***
ORIGIN_SZMPSZ03 9.798e-02 4.286e-03 22.859 < 2e-16 ***
ORIGIN_SZMSSZ01 -8.790e+00 5.774e-01 -15.225 < 2e-16 ***
ORIGIN_SZMUSZ01 -1.278e+00 5.864e-03 -217.983 < 2e-16 ***
ORIGIN_SZMUSZ02 -3.199e+00 1.508e-02 -212.130 < 2e-16 ***
ORIGIN_SZMUSZ03 -1.859e+00 6.896e-03 -269.625 < 2e-16 ***
ORIGIN_SZNTSZ01 -2.371e+00 2.539e-02 -93.365 < 2e-16 ***
ORIGIN_SZNTSZ02 -2.678e+00 1.430e-02 -187.274 < 2e-16 ***
ORIGIN_SZNTSZ03 -8.217e-01 5.910e-03 -139.032 < 2e-16 ***
ORIGIN_SZNTSZ05 -3.086e+00 3.681e-02 -83.854 < 2e-16 ***
ORIGIN_SZNTSZ06 -3.480e+00 4.117e-02 -84.526 < 2e-16 ***
ORIGIN_SZNVSZ01 6.222e-01 3.576e-03 174.005 < 2e-16 ***
ORIGIN_SZNVSZ02 -4.850e-01 4.908e-03 -98.813 < 2e-16 ***
ORIGIN_SZNVSZ03 -1.138e+00 6.105e-03 -186.411 < 2e-16 ***
ORIGIN_SZNVSZ04 -1.341e+00 7.156e-03 -187.457 < 2e-16 ***
ORIGIN_SZNVSZ05 -2.592e+00 1.333e-02 -194.516 < 2e-16 ***
ORIGIN_SZORSZ01 -2.926e+00 2.704e-02 -108.241 < 2e-16 ***
ORIGIN_SZORSZ02 -1.247e+00 5.849e-03 -213.260 < 2e-16 ***
ORIGIN_SZORSZ03 -1.701e+00 6.994e-03 -243.285 < 2e-16 ***
ORIGIN_SZOTSZ01 -1.707e+00 7.271e-03 -234.799 < 2e-16 ***
ORIGIN_SZOTSZ02 -1.752e+00 8.125e-03 -215.586 < 2e-16 ***
ORIGIN_SZOTSZ03 -8.029e-01 5.301e-03 -151.463 < 2e-16 ***
ORIGIN_SZOTSZ04 -7.437e-01 8.307e-03 -89.520 < 2e-16 ***
ORIGIN_SZPGSZ01 -7.216e-01 9.296e-03 -77.619 < 2e-16 ***
ORIGIN_SZPGSZ02 -3.412e-01 5.675e-03 -60.112 < 2e-16 ***
ORIGIN_SZPGSZ03 1.139e+00 3.636e-03 313.334 < 2e-16 ***
ORIGIN_SZPGSZ04 1.229e+00 3.635e-03 338.242 < 2e-16 ***
ORIGIN_SZPGSZ05 4.388e-01 4.535e-03 96.768 < 2e-16 ***
ORIGIN_SZPLSZ01 -5.164e-01 7.882e-03 -65.514 < 2e-16 ***
ORIGIN_SZPLSZ02 -1.422e+00 1.136e-02 -125.152 < 2e-16 ***
ORIGIN_SZPLSZ03 -2.962e+00 3.150e-02 -94.020 < 2e-16 ***
ORIGIN_SZPLSZ04 -2.561e+00 3.713e-02 -68.962 < 2e-16 ***
ORIGIN_SZPLSZ05 -2.215e+00 1.717e-02 -128.966 < 2e-16 ***
ORIGIN_SZPNSZ01 1.514e+00 3.830e-03 395.240 < 2e-16 ***
ORIGIN_SZPNSZ02 -5.651e-01 9.775e-03 -57.807 < 2e-16 ***
ORIGIN_SZPNSZ03 -1.941e+00 1.725e-02 -112.539 < 2e-16 ***
ORIGIN_SZPNSZ04 -2.605e+00 2.486e-02 -104.808 < 2e-16 ***
ORIGIN_SZPNSZ05 -1.791e+00 1.781e-02 -100.565 < 2e-16 ***
ORIGIN_SZPRSZ01 -7.687e-01 9.568e-03 -80.340 < 2e-16 ***
ORIGIN_SZPRSZ02 1.057e+00 3.815e-03 277.009 < 2e-16 ***
ORIGIN_SZPRSZ03 8.494e-01 3.812e-03 222.816 < 2e-16 ***
ORIGIN_SZPRSZ04 -2.751e-01 6.179e-03 -44.527 < 2e-16 ***
ORIGIN_SZPRSZ05 1.279e+00 3.636e-03 351.679 < 2e-16 ***
ORIGIN_SZPRSZ06 -5.141e-01 6.870e-03 -74.835 < 2e-16 ***
ORIGIN_SZPRSZ07 -2.565e+00 1.709e-02 -150.036 < 2e-16 ***
ORIGIN_SZPRSZ08 1.155e-01 5.028e-03 22.967 < 2e-16 ***
ORIGIN_SZQTSZ01 -4.434e-01 5.419e-03 -81.814 < 2e-16 ***
ORIGIN_SZQTSZ02 -7.177e-01 5.071e-03 -141.517 < 2e-16 ***
ORIGIN_SZQTSZ03 -2.192e-01 4.605e-03 -47.590 < 2e-16 ***
ORIGIN_SZQTSZ04 -1.216e+00 6.380e-03 -190.643 < 2e-16 ***
ORIGIN_SZQTSZ05 -2.460e-01 4.629e-03 -53.133 < 2e-16 ***
ORIGIN_SZQTSZ06 -6.079e-01 5.311e-03 -114.465 < 2e-16 ***
ORIGIN_SZQTSZ07 -1.609e+00 7.926e-03 -202.991 < 2e-16 ***
ORIGIN_SZQTSZ08 -3.407e-01 4.864e-03 -70.045 < 2e-16 ***
ORIGIN_SZQTSZ09 -5.850e-01 5.418e-03 -107.967 < 2e-16 ***
ORIGIN_SZQTSZ10 -4.832e-01 5.414e-03 -89.248 < 2e-16 ***
ORIGIN_SZQTSZ11 -1.363e+00 7.880e-03 -172.991 < 2e-16 ***
ORIGIN_SZQTSZ12 -8.072e-01 6.660e-03 -121.197 < 2e-16 ***
ORIGIN_SZQTSZ13 -1.359e-01 4.889e-03 -27.794 < 2e-16 ***
ORIGIN_SZQTSZ14 -1.378e+00 7.243e-03 -190.293 < 2e-16 ***
ORIGIN_SZQTSZ15 -7.662e-01 8.389e-03 -91.335 < 2e-16 ***
ORIGIN_SZRCSZ01 -6.511e-01 5.320e-03 -122.393 < 2e-16 ***
ORIGIN_SZRCSZ02 -2.302e+00 1.519e-02 -151.558 < 2e-16 ***
ORIGIN_SZRCSZ03 -1.487e+00 7.473e-03 -198.987 < 2e-16 ***
ORIGIN_SZRCSZ04 -2.264e+00 1.142e-02 -198.187 < 2e-16 ***
ORIGIN_SZRCSZ05 -2.852e+00 1.395e-02 -204.500 < 2e-16 ***
ORIGIN_SZRCSZ06 -4.847e-01 6.998e-03 -69.267 < 2e-16 ***
ORIGIN_SZRCSZ08 -2.701e+00 1.696e-02 -159.243 < 2e-16 ***
ORIGIN_SZRCSZ09 -1.979e+00 1.218e-02 -162.421 < 2e-16 ***
ORIGIN_SZRCSZ10 -1.814e+00 6.961e-03 -260.529 < 2e-16 ***
ORIGIN_SZRVSZ01 -2.924e+00 1.381e-02 -211.798 < 2e-16 ***
ORIGIN_SZRVSZ02 -1.247e+00 6.825e-03 -182.764 < 2e-16 ***
ORIGIN_SZRVSZ03 -1.986e+00 1.024e-02 -194.024 < 2e-16 ***
ORIGIN_SZRVSZ04 -1.956e+00 1.446e-02 -135.317 < 2e-16 ***
ORIGIN_SZRVSZ05 -2.334e+00 1.301e-02 -179.333 < 2e-16 ***
ORIGIN_SZSBSZ01 8.094e-01 4.404e-03 183.770 < 2e-16 ***
ORIGIN_SZSBSZ02 -5.981e-01 6.511e-03 -91.864 < 2e-16 ***
ORIGIN_SZSBSZ03 9.692e-01 3.866e-03 250.676 < 2e-16 ***
ORIGIN_SZSBSZ04 8.023e-01 4.374e-03 183.402 < 2e-16 ***
ORIGIN_SZSBSZ05 -1.164e-01 5.535e-03 -21.025 < 2e-16 ***
ORIGIN_SZSBSZ06 -1.658e+00 1.366e-02 -121.375 < 2e-16 ***
ORIGIN_SZSBSZ07 -8.111e-01 9.288e-03 -87.327 < 2e-16 ***
ORIGIN_SZSBSZ08 -1.018e+00 9.606e-03 -105.978 < 2e-16 ***
ORIGIN_SZSBSZ09 -4.647e-01 7.091e-03 -65.531 < 2e-16 ***
ORIGIN_SZSESZ02 1.129e+00 3.612e-03 312.605 < 2e-16 ***
ORIGIN_SZSESZ03 1.266e+00 3.468e-03 364.995 < 2e-16 ***
ORIGIN_SZSESZ04 9.950e-01 3.944e-03 252.294 < 2e-16 ***
ORIGIN_SZSESZ05 -1.600e-01 4.861e-03 -32.917 < 2e-16 ***
ORIGIN_SZSESZ06 9.714e-01 3.773e-03 257.433 < 2e-16 ***
ORIGIN_SZSESZ07 -2.192e+00 1.401e-02 -156.422 < 2e-16 ***
ORIGIN_SZSGSZ01 -9.081e-01 6.942e-03 -130.802 < 2e-16 ***
ORIGIN_SZSGSZ02 -1.121e+00 8.242e-03 -136.001 < 2e-16 ***
ORIGIN_SZSGSZ03 2.957e-01 4.310e-03 68.610 < 2e-16 ***
ORIGIN_SZSGSZ04 3.137e-01 3.974e-03 78.934 < 2e-16 ***
ORIGIN_SZSGSZ05 -1.706e+00 8.608e-03 -198.143 < 2e-16 ***
ORIGIN_SZSGSZ06 4.110e-01 3.784e-03 108.625 < 2e-16 ***
ORIGIN_SZSGSZ07 -5.318e-01 5.001e-03 -106.342 < 2e-16 ***
ORIGIN_SZSKSZ01 -1.875e-01 6.551e-03 -28.626 < 2e-16 ***
ORIGIN_SZSKSZ02 2.968e-01 4.819e-03 61.597 < 2e-16 ***
ORIGIN_SZSKSZ03 -4.481e-01 6.175e-03 -72.566 < 2e-16 ***
ORIGIN_SZSKSZ04 -2.421e+00 2.210e-02 -109.528 < 2e-16 ***
ORIGIN_SZSKSZ05 -1.178e+00 1.220e-02 -96.532 < 2e-16 ***
ORIGIN_SZSLSZ01 -3.191e+00 2.594e-02 -123.005 < 2e-16 ***
ORIGIN_SZSLSZ04 -3.569e-01 5.839e-03 -61.120 < 2e-16 ***
ORIGIN_SZSRSZ01 -1.609e+00 7.335e-03 -219.281 < 2e-16 ***
ORIGIN_SZSRSZ02 -1.802e+00 7.294e-03 -247.088 < 2e-16 ***
ORIGIN_SZSRSZ03 -2.852e+00 1.575e-02 -181.121 < 2e-16 ***
ORIGIN_SZSVSZ01 -2.793e+00 2.818e-02 -99.109 < 2e-16 ***
ORIGIN_SZTHSZ01 -3.163e+00 5.294e-02 -59.752 < 2e-16 ***
ORIGIN_SZTHSZ03 -1.531e+00 1.250e-02 -122.477 < 2e-16 ***
ORIGIN_SZTHSZ04 -2.909e+00 2.509e-02 -115.956 < 2e-16 ***
ORIGIN_SZTHSZ06 -1.759e+00 1.152e-02 -152.645 < 2e-16 ***
ORIGIN_SZTMSZ01 9.814e-01 4.123e-03 238.018 < 2e-16 ***
ORIGIN_SZTMSZ02 2.210e+00 3.182e-03 694.536 < 2e-16 ***
ORIGIN_SZTMSZ03 1.473e+00 3.403e-03 432.821 < 2e-16 ***
ORIGIN_SZTMSZ04 9.293e-01 3.912e-03 237.539 < 2e-16 ***
ORIGIN_SZTMSZ05 -3.070e-01 6.257e-03 -49.070 < 2e-16 ***
ORIGIN_SZTNSZ01 -1.213e+00 6.224e-03 -194.863 < 2e-16 ***
ORIGIN_SZTNSZ02 -1.086e+00 5.885e-03 -184.612 < 2e-16 ***
ORIGIN_SZTNSZ03 -1.618e+00 7.881e-03 -205.292 < 2e-16 ***
ORIGIN_SZTNSZ04 -7.424e-01 5.668e-03 -130.978 < 2e-16 ***
ORIGIN_SZTPSZ01 -6.483e-01 5.130e-03 -126.371 < 2e-16 ***
ORIGIN_SZTPSZ02 4.440e-01 3.626e-03 122.443 < 2e-16 ***
ORIGIN_SZTPSZ03 -5.077e-01 5.183e-03 -97.963 < 2e-16 ***
ORIGIN_SZTPSZ04 -2.901e-01 4.842e-03 -59.906 < 2e-16 ***
ORIGIN_SZTPSZ05 -2.393e-01 5.031e-03 -47.566 < 2e-16 ***
ORIGIN_SZTPSZ06 2.111e-01 5.031e-03 41.961 < 2e-16 ***
ORIGIN_SZTPSZ07 -1.719e-01 5.057e-03 -33.998 < 2e-16 ***
ORIGIN_SZTPSZ08 -7.438e-01 6.660e-03 -111.688 < 2e-16 ***
ORIGIN_SZTPSZ09 -5.087e-01 5.378e-03 -94.589 < 2e-16 ***
ORIGIN_SZTPSZ10 -3.941e-01 5.422e-03 -72.683 < 2e-16 ***
ORIGIN_SZTPSZ11 2.500e-01 4.233e-03 59.058 < 2e-16 ***
ORIGIN_SZTPSZ12 -5.277e-01 5.332e-03 -98.959 < 2e-16 ***
ORIGIN_SZTSSZ01 -3.376e+00 3.881e-02 -86.987 < 2e-16 ***
ORIGIN_SZTSSZ02 5.017e-01 5.840e-03 85.906 < 2e-16 ***
ORIGIN_SZTSSZ03 5.258e-01 5.880e-03 89.425 < 2e-16 ***
ORIGIN_SZTSSZ04 5.202e-01 6.017e-03 86.453 < 2e-16 ***
ORIGIN_SZTSSZ05 -9.979e-01 1.152e-02 -86.635 < 2e-16 ***
ORIGIN_SZTSSZ06 -1.205e+00 1.400e-02 -86.054 < 2e-16 ***
ORIGIN_SZWCSZ01 2.115e-01 5.900e-03 35.847 < 2e-16 ***
ORIGIN_SZWCSZ02 -2.919e+00 2.780e-02 -105.000 < 2e-16 ***
ORIGIN_SZWCSZ03 -4.669e+00 1.401e-01 -33.336 < 2e-16 ***
ORIGIN_SZWDSZ01 1.450e+00 3.462e-03 419.015 < 2e-16 ***
ORIGIN_SZWDSZ02 1.094e+00 3.941e-03 277.563 < 2e-16 ***
ORIGIN_SZWDSZ03 2.218e+00 3.368e-03 658.639 < 2e-16 ***
ORIGIN_SZWDSZ04 1.165e+00 4.118e-03 282.874 < 2e-16 ***
ORIGIN_SZWDSZ05 5.627e-01 4.172e-03 134.882 < 2e-16 ***
ORIGIN_SZWDSZ06 1.265e+00 3.822e-03 330.935 < 2e-16 ***
ORIGIN_SZWDSZ07 4.847e-02 5.877e-03 8.247 < 2e-16 ***
ORIGIN_SZWDSZ08 -3.533e-01 6.566e-03 -53.810 < 2e-16 ***
ORIGIN_SZWDSZ09 1.737e+00 3.541e-03 490.615 < 2e-16 ***
ORIGIN_SZYSSZ01 -7.839e-02 4.671e-03 -16.784 < 2e-16 ***
ORIGIN_SZYSSZ02 9.726e-01 4.154e-03 234.118 < 2e-16 ***
ORIGIN_SZYSSZ03 1.971e+00 3.468e-03 568.377 < 2e-16 ***
ORIGIN_SZYSSZ04 9.231e-01 3.740e-03 246.851 < 2e-16 ***
ORIGIN_SZYSSZ05 1.667e-01 4.618e-03 36.094 < 2e-16 ***
ORIGIN_SZYSSZ06 -7.774e-01 7.463e-03 -104.159 < 2e-16 ***
ORIGIN_SZYSSZ07 -5.435e-01 7.184e-03 -75.656 < 2e-16 ***
ORIGIN_SZYSSZ08 -2.803e-02 5.144e-03 -5.449 5.06e-08 ***
ORIGIN_SZYSSZ09 1.388e+00 3.573e-03 388.457 < 2e-16 ***
log(DESTIN_AGE25_64) 2.695e-02 6.755e-05 398.928 < 2e-16 ***
log(dist) -7.131e-01 1.031e-04 -6917.129 < 2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
(Dispersion parameter for poisson family taken to be 1)
Null deviance: 102725286 on 20848 degrees of freedom
Residual deviance: 44624132 on 20538 degrees of freedom
AIC: 44756335
Number of Fisher Scoring iterations: 7Let’s check the R-square values of origin constrained SIM model this time.
CalcRSquared(orcSIM$data$TRIPS, orcSIM$fitted.values)[1] 0.3886462Notably, R-squared improves from the unconstrained SIM model.
7.5 Destination Constrained Spatial Interaction Model
In this section, we will calibrate an destination constrained SIM. For destination constrained SIM, only explanatory variables which represent how propulsive the origins are will be used. This is because such models emphasize the demand or attractiveness of the destinations rather than the limitations or capacities of the origins. The demand or attractiveness of the destination sites determines the potential for generating interactions or flows.
decSIM <- glm(formula = TRIPS ~
DESTIN_SZ +
log(ORIGIN_AGE25_64) +
log(dist),
family = poisson(link = "log"),
data = SIM_data,
na.action = na.exclude)
summary(decSIM)
Call:
glm(formula = TRIPS ~ DESTIN_SZ + log(ORIGIN_AGE25_64) + log(dist),
family = poisson(link = "log"), data = SIM_data, na.action = na.exclude)
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 11.1867022 0.0027974 3998.965 < 2e-16 ***
DESTIN_SZAMSZ02 0.2157520 0.0035355 61.024 < 2e-16 ***
DESTIN_SZAMSZ03 0.2435289 0.0034884 69.811 < 2e-16 ***
DESTIN_SZAMSZ04 -0.9023616 0.0051831 -174.097 < 2e-16 ***
DESTIN_SZAMSZ05 -0.9484973 0.0048931 -193.844 < 2e-16 ***
DESTIN_SZAMSZ06 -0.7571980 0.0047785 -158.458 < 2e-16 ***
DESTIN_SZAMSZ07 -1.8444727 0.0085349 -216.110 < 2e-16 ***
DESTIN_SZAMSZ08 -0.9299210 0.0054401 -170.937 < 2e-16 ***
DESTIN_SZAMSZ09 -0.9908230 0.0052101 -190.174 < 2e-16 ***
DESTIN_SZAMSZ10 0.4562387 0.0035574 128.250 < 2e-16 ***
DESTIN_SZAMSZ11 0.1019878 0.0061208 16.662 < 2e-16 ***
DESTIN_SZAMSZ12 -0.1309597 0.0044830 -29.212 < 2e-16 ***
DESTIN_SZBDSZ01 0.5260496 0.0032088 163.940 < 2e-16 ***
DESTIN_SZBDSZ02 -0.1623807 0.0040830 -39.770 < 2e-16 ***
DESTIN_SZBDSZ03 0.0142683 0.0036127 3.949 7.83e-05 ***
DESTIN_SZBDSZ04 1.0129906 0.0029404 344.511 < 2e-16 ***
DESTIN_SZBDSZ05 0.4240820 0.0032733 129.560 < 2e-16 ***
DESTIN_SZBDSZ06 0.1356215 0.0036775 36.879 < 2e-16 ***
DESTIN_SZBDSZ07 -0.9548822 0.0074517 -128.142 < 2e-16 ***
DESTIN_SZBDSZ08 -1.4063187 0.0075549 -186.146 < 2e-16 ***
DESTIN_SZBKSZ01 -1.0781080 0.0053701 -200.761 < 2e-16 ***
DESTIN_SZBKSZ02 -0.1825131 0.0043597 -41.864 < 2e-16 ***
DESTIN_SZBKSZ03 -0.4468506 0.0045296 -98.652 < 2e-16 ***
DESTIN_SZBKSZ04 0.0702745 0.0039971 17.582 < 2e-16 ***
DESTIN_SZBKSZ05 -0.6713685 0.0046264 -145.116 < 2e-16 ***
DESTIN_SZBKSZ06 -0.9025531 0.0052372 -172.334 < 2e-16 ***
DESTIN_SZBKSZ07 0.2745051 0.0033978 80.789 < 2e-16 ***
DESTIN_SZBKSZ08 -1.0320480 0.0057424 -179.724 < 2e-16 ***
DESTIN_SZBKSZ09 -0.2317999 0.0040906 -56.667 < 2e-16 ***
DESTIN_SZBLSZ01 -0.3307294 0.0058398 -56.634 < 2e-16 ***
DESTIN_SZBLSZ02 0.6534505 0.0056360 115.943 < 2e-16 ***
DESTIN_SZBLSZ03 1.4781734 0.0060542 244.157 < 2e-16 ***
DESTIN_SZBLSZ04 0.0790456 0.0107855 7.329 2.32e-13 ***
DESTIN_SZBMSZ01 -0.1559248 0.0037381 -41.712 < 2e-16 ***
DESTIN_SZBMSZ02 -0.4971169 0.0039272 -126.584 < 2e-16 ***
DESTIN_SZBMSZ03 -0.9102203 0.0048043 -189.461 < 2e-16 ***
DESTIN_SZBMSZ04 -0.5705684 0.0041013 -139.118 < 2e-16 ***
DESTIN_SZBMSZ05 -0.5949305 0.0049705 -119.692 < 2e-16 ***
DESTIN_SZBMSZ06 -1.6558842 0.0084636 -195.648 < 2e-16 ***
DESTIN_SZBMSZ07 -0.1676643 0.0036334 -46.145 < 2e-16 ***
DESTIN_SZBMSZ08 -1.1095628 0.0048276 -229.837 < 2e-16 ***
DESTIN_SZBMSZ09 -2.0402834 0.0076816 -265.606 < 2e-16 ***
DESTIN_SZBMSZ10 -1.5562070 0.0060101 -258.932 < 2e-16 ***
DESTIN_SZBMSZ11 -1.5936693 0.0060162 -264.895 < 2e-16 ***
DESTIN_SZBMSZ12 -1.2308134 0.0066734 -184.435 < 2e-16 ***
DESTIN_SZBMSZ13 -0.2571982 0.0038941 -66.048 < 2e-16 ***
DESTIN_SZBMSZ14 -1.0457304 0.0061118 -171.101 < 2e-16 ***
DESTIN_SZBMSZ15 -1.2507466 0.0057444 -217.732 < 2e-16 ***
DESTIN_SZBMSZ16 -1.6149520 0.0061458 -262.773 < 2e-16 ***
DESTIN_SZBMSZ17 -1.5270738 0.0072589 -210.373 < 2e-16 ***
DESTIN_SZBPSZ01 -0.4907416 0.0044484 -110.319 < 2e-16 ***
DESTIN_SZBPSZ02 -1.4328876 0.0071285 -201.009 < 2e-16 ***
DESTIN_SZBPSZ03 -1.1435451 0.0066159 -172.849 < 2e-16 ***
DESTIN_SZBPSZ04 -0.6297222 0.0049102 -128.247 < 2e-16 ***
DESTIN_SZBPSZ05 0.5236032 0.0032571 160.759 < 2e-16 ***
DESTIN_SZBPSZ06 -0.9980571 0.0067305 -148.289 < 2e-16 ***
DESTIN_SZBPSZ07 -0.5826701 0.0063629 -91.574 < 2e-16 ***
DESTIN_SZBSSZ01 -0.1697808 0.0037485 -45.293 < 2e-16 ***
DESTIN_SZBSSZ02 -0.7762707 0.0043558 -178.214 < 2e-16 ***
DESTIN_SZBSSZ03 0.3075821 0.0031998 96.126 < 2e-16 ***
DESTIN_SZBTSZ01 0.2760882 0.0033675 81.986 < 2e-16 ***
DESTIN_SZBTSZ02 -0.7722114 0.0056050 -137.772 < 2e-16 ***
DESTIN_SZBTSZ03 0.0059834 0.0038328 1.561 0.11850
DESTIN_SZBTSZ04 -1.6335853 0.0080414 -203.148 < 2e-16 ***
DESTIN_SZBTSZ05 -0.7492179 0.0056614 -132.339 < 2e-16 ***
DESTIN_SZBTSZ06 -0.8945211 0.0050858 -175.887 < 2e-16 ***
DESTIN_SZBTSZ07 -2.0193686 0.0082480 -244.831 < 2e-16 ***
DESTIN_SZBTSZ08 -1.3129105 0.0069456 -189.027 < 2e-16 ***
DESTIN_SZCCSZ01 -0.6226949 0.0058387 -106.650 < 2e-16 ***
DESTIN_SZCHSZ01 -1.0906540 0.0075296 -144.848 < 2e-16 ***
DESTIN_SZCHSZ02 -0.0128257 0.0045923 -2.793 0.00522 **
DESTIN_SZCHSZ03 1.6641822 0.0032262 515.832 < 2e-16 ***
DESTIN_SZCKSZ01 -0.1486951 0.0040634 -36.594 < 2e-16 ***
DESTIN_SZCKSZ02 -0.4306479 0.0044366 -97.067 < 2e-16 ***
DESTIN_SZCKSZ03 0.7344433 0.0032448 226.348 < 2e-16 ***
DESTIN_SZCKSZ04 -0.6586195 0.0051131 -128.809 < 2e-16 ***
DESTIN_SZCKSZ05 -0.3488802 0.0054369 -64.170 < 2e-16 ***
DESTIN_SZCKSZ06 0.7880887 0.0037632 209.419 < 2e-16 ***
DESTIN_SZCLSZ01 0.4778561 0.0039409 121.257 < 2e-16 ***
DESTIN_SZCLSZ02 -2.3268839 0.0112332 -207.144 < 2e-16 ***
DESTIN_SZCLSZ03 -0.9012108 0.0060507 -148.943 < 2e-16 ***
DESTIN_SZCLSZ04 0.1294757 0.0036022 35.944 < 2e-16 ***
DESTIN_SZCLSZ05 -1.2726360 0.0070980 -179.295 < 2e-16 ***
DESTIN_SZCLSZ06 0.1795746 0.0034057 52.728 < 2e-16 ***
DESTIN_SZCLSZ07 -0.6190253 0.0044971 -137.649 < 2e-16 ***
DESTIN_SZCLSZ08 -0.4619840 0.0048835 -94.600 < 2e-16 ***
DESTIN_SZCLSZ09 0.3280972 0.0054454 60.253 < 2e-16 ***
DESTIN_SZDTSZ01 -0.8845669 0.0045026 -196.458 < 2e-16 ***
DESTIN_SZDTSZ02 -0.8568300 0.0043877 -195.279 < 2e-16 ***
DESTIN_SZDTSZ03 -1.0874333 0.0052439 -207.372 < 2e-16 ***
DESTIN_SZDTSZ04 -1.1494437 0.0112537 -102.140 < 2e-16 ***
DESTIN_SZDTSZ05 -1.0362291 0.0087132 -118.927 < 2e-16 ***
DESTIN_SZDTSZ06 -1.1614611 0.0058684 -197.917 < 2e-16 ***
DESTIN_SZDTSZ07 -1.8984868 0.0179356 -105.850 < 2e-16 ***
DESTIN_SZDTSZ08 -0.7077496 0.0042454 -166.711 < 2e-16 ***
DESTIN_SZDTSZ09 -1.5752269 0.0096965 -162.453 < 2e-16 ***
DESTIN_SZDTSZ10 -1.3537512 0.0075965 -178.207 < 2e-16 ***
DESTIN_SZDTSZ11 -0.8090332 0.0045013 -179.731 < 2e-16 ***
DESTIN_SZDTSZ12 -2.4347627 0.0148595 -163.853 < 2e-16 ***
DESTIN_SZDTSZ13 -2.0504803 0.0094558 -216.850 < 2e-16 ***
DESTIN_SZGLSZ01 0.0598692 0.0041609 14.388 < 2e-16 ***
DESTIN_SZGLSZ02 -0.3192364 0.0038784 -82.312 < 2e-16 ***
DESTIN_SZGLSZ03 0.4035376 0.0032534 124.037 < 2e-16 ***
DESTIN_SZGLSZ04 0.3500292 0.0031921 109.653 < 2e-16 ***
DESTIN_SZGLSZ05 0.2092624 0.0032941 63.527 < 2e-16 ***
DESTIN_SZHGSZ01 0.4285880 0.0032318 132.617 < 2e-16 ***
DESTIN_SZHGSZ02 -0.5794121 0.0044398 -130.505 < 2e-16 ***
DESTIN_SZHGSZ03 -1.0762078 0.0053359 -201.692 < 2e-16 ***
DESTIN_SZHGSZ04 -0.2662674 0.0037779 -70.481 < 2e-16 ***
DESTIN_SZHGSZ05 -0.1466547 0.0037801 -38.796 < 2e-16 ***
DESTIN_SZHGSZ06 -0.6782939 0.0044485 -152.479 < 2e-16 ***
DESTIN_SZHGSZ07 0.2513178 0.0034156 73.580 < 2e-16 ***
DESTIN_SZHGSZ08 -0.2891581 0.0040130 -72.055 < 2e-16 ***
DESTIN_SZHGSZ09 0.3032494 0.0041747 72.640 < 2e-16 ***
DESTIN_SZHGSZ10 -3.3530740 0.0292581 -114.603 < 2e-16 ***
DESTIN_SZJESZ01 -0.1716230 0.0042420 -40.458 < 2e-16 ***
DESTIN_SZJESZ02 -0.4806391 0.0043000 -111.776 < 2e-16 ***
DESTIN_SZJESZ03 -0.6281332 0.0047290 -132.826 < 2e-16 ***
DESTIN_SZJESZ04 -0.1385887 0.0051022 -27.162 < 2e-16 ***
DESTIN_SZJESZ05 -0.8062739 0.0075357 -106.994 < 2e-16 ***
DESTIN_SZJESZ06 0.4153823 0.0034451 120.572 < 2e-16 ***
DESTIN_SZJESZ07 -0.8640741 0.0062388 -138.499 < 2e-16 ***
DESTIN_SZJESZ08 -0.5818943 0.0064926 -89.625 < 2e-16 ***
DESTIN_SZJESZ09 -0.4955426 0.0047394 -104.559 < 2e-16 ***
DESTIN_SZJESZ10 0.7060401 0.0062127 113.645 < 2e-16 ***
DESTIN_SZJESZ11 1.0253229 0.0055311 185.375 < 2e-16 ***
DESTIN_SZJWSZ01 -0.3670256 0.0053684 -68.368 < 2e-16 ***
DESTIN_SZJWSZ02 -0.3526519 0.0045263 -77.912 < 2e-16 ***
DESTIN_SZJWSZ03 0.6441760 0.0033769 190.758 < 2e-16 ***
DESTIN_SZJWSZ04 1.0326753 0.0031443 328.427 < 2e-16 ***
DESTIN_SZJWSZ05 -0.1502422 0.0050433 -29.790 < 2e-16 ***
DESTIN_SZJWSZ06 0.3959462 0.0046070 85.944 < 2e-16 ***
DESTIN_SZJWSZ07 -0.7665967 0.0191890 -39.950 < 2e-16 ***
DESTIN_SZJWSZ08 0.3818886 0.0039317 97.130 < 2e-16 ***
DESTIN_SZJWSZ09 1.4839379 0.0028767 515.856 < 2e-16 ***
DESTIN_SZKLSZ01 -0.6023948 0.0040944 -147.128 < 2e-16 ***
DESTIN_SZKLSZ02 -0.7999697 0.0046892 -170.598 < 2e-16 ***
DESTIN_SZKLSZ03 -1.2547986 0.0051793 -242.274 < 2e-16 ***
DESTIN_SZKLSZ04 -1.7610618 0.0067522 -260.815 < 2e-16 ***
DESTIN_SZKLSZ05 -1.1618041 0.0068174 -170.417 < 2e-16 ***
DESTIN_SZKLSZ06 -0.9580560 0.0045104 -212.409 < 2e-16 ***
DESTIN_SZKLSZ07 -1.0710104 0.0050836 -210.680 < 2e-16 ***
DESTIN_SZKLSZ08 -0.1243401 0.0036596 -33.976 < 2e-16 ***
DESTIN_SZKLSZ09 -1.8312633 0.0066208 -276.592 < 2e-16 ***
DESTIN_SZLKSZ01 -1.6497161 0.0190096 -86.783 < 2e-16 ***
DESTIN_SZMDSZ01 -1.3105988 0.0158618 -82.626 < 2e-16 ***
DESTIN_SZMDSZ02 -1.0992364 0.0093242 -117.890 < 2e-16 ***
DESTIN_SZMDSZ03 -2.8656212 0.0210282 -136.275 < 2e-16 ***
DESTIN_SZMPSZ01 -1.3100389 0.0066619 -196.645 < 2e-16 ***
DESTIN_SZMPSZ02 -0.8719388 0.0048810 -178.641 < 2e-16 ***
DESTIN_SZMPSZ03 -0.2223471 0.0039632 -56.103 < 2e-16 ***
DESTIN_SZMSSZ01 -3.8589249 0.0737669 -52.312 < 2e-16 ***
DESTIN_SZMUSZ01 -1.1915341 0.0049498 -240.724 < 2e-16 ***
DESTIN_SZMUSZ02 -1.4516013 0.0071123 -204.096 < 2e-16 ***
DESTIN_SZMUSZ03 -1.2439342 0.0049899 -249.292 < 2e-16 ***
DESTIN_SZNTSZ01 -2.6749916 0.0221410 -120.816 < 2e-16 ***
DESTIN_SZNTSZ02 -2.0703748 0.0090174 -229.597 < 2e-16 ***
DESTIN_SZNTSZ03 -1.2352944 0.0061847 -199.734 < 2e-16 ***
DESTIN_SZNTSZ05 -1.7770413 0.0165151 -107.601 < 2e-16 ***
DESTIN_SZNTSZ06 -2.9854061 0.0275374 -108.413 < 2e-16 ***
DESTIN_SZNVSZ01 -0.2807211 0.0035950 -78.087 < 2e-16 ***
DESTIN_SZNVSZ02 -0.4938925 0.0041125 -120.095 < 2e-16 ***
DESTIN_SZNVSZ03 -0.5998929 0.0043582 -137.646 < 2e-16 ***
DESTIN_SZNVSZ04 -2.1388131 0.0087051 -245.698 < 2e-16 ***
DESTIN_SZNVSZ05 -1.8866286 0.0075631 -249.452 < 2e-16 ***
DESTIN_SZORSZ01 -1.9908872 0.0173765 -114.573 < 2e-16 ***
DESTIN_SZORSZ02 -0.1725616 0.0037208 -46.378 < 2e-16 ***
DESTIN_SZORSZ03 -0.9527747 0.0048150 -197.874 < 2e-16 ***
DESTIN_SZOTSZ01 -1.5938496 0.0061762 -258.062 < 2e-16 ***
DESTIN_SZOTSZ02 -0.8247973 0.0053824 -153.241 < 2e-16 ***
DESTIN_SZOTSZ03 -1.4953393 0.0057915 -258.196 < 2e-16 ***
DESTIN_SZOTSZ04 -1.5827262 0.0083269 -190.074 < 2e-16 ***
DESTIN_SZPGSZ01 -2.1519057 0.0143339 -150.127 < 2e-16 ***
DESTIN_SZPGSZ02 -0.7694243 0.0053698 -143.288 < 2e-16 ***
DESTIN_SZPGSZ03 0.5762599 0.0032828 175.538 < 2e-16 ***
DESTIN_SZPGSZ04 0.1405343 0.0036966 38.017 < 2e-16 ***
DESTIN_SZPGSZ05 -0.9141066 0.0061275 -149.182 < 2e-16 ***
DESTIN_SZPLSZ01 -0.1188682 0.0057893 -20.532 < 2e-16 ***
DESTIN_SZPLSZ02 -1.3085587 0.0105564 -123.959 < 2e-16 ***
DESTIN_SZPLSZ03 -0.1078859 0.0083997 -12.844 < 2e-16 ***
DESTIN_SZPLSZ04 -0.0097287 0.0080244 -1.212 0.22536
DESTIN_SZPLSZ05 -0.6868459 0.0098218 -69.931 < 2e-16 ***
DESTIN_SZPNSZ01 1.1092381 0.0042567 260.585 < 2e-16 ***
DESTIN_SZPNSZ02 1.6409198 0.0055758 294.291 < 2e-16 ***
DESTIN_SZPNSZ03 0.9437528 0.0063421 148.809 < 2e-16 ***
DESTIN_SZPNSZ04 1.7505749 0.0063198 276.998 < 2e-16 ***
DESTIN_SZPNSZ05 0.9235380 0.0091335 101.116 < 2e-16 ***
DESTIN_SZPRSZ01 -0.6866781 0.0062158 -110.474 < 2e-16 ***
DESTIN_SZPRSZ02 -0.1412976 0.0041929 -33.699 < 2e-16 ***
DESTIN_SZPRSZ03 0.7883502 0.0031559 249.801 < 2e-16 ***
DESTIN_SZPRSZ04 -0.7862784 0.0071030 -110.696 < 2e-16 ***
DESTIN_SZPRSZ05 -0.0580155 0.0039437 -14.711 < 2e-16 ***
DESTIN_SZPRSZ06 0.4698562 0.0041687 112.712 < 2e-16 ***
DESTIN_SZPRSZ07 -1.4826676 0.0099226 -149.424 < 2e-16 ***
DESTIN_SZPRSZ08 -0.8120437 0.0055879 -145.322 < 2e-16 ***
DESTIN_SZQTSZ01 -1.6130987 0.0084135 -191.727 < 2e-16 ***
DESTIN_SZQTSZ02 -1.4973414 0.0060894 -245.895 < 2e-16 ***
DESTIN_SZQTSZ03 -0.8872366 0.0053600 -165.529 < 2e-16 ***
DESTIN_SZQTSZ04 -1.1556653 0.0055741 -207.328 < 2e-16 ***
DESTIN_SZQTSZ05 -0.9777188 0.0049348 -198.127 < 2e-16 ***
DESTIN_SZQTSZ06 -1.2069541 0.0052575 -229.566 < 2e-16 ***
DESTIN_SZQTSZ07 -1.7714823 0.0089330 -198.308 < 2e-16 ***
DESTIN_SZQTSZ08 0.0692739 0.0038135 18.165 < 2e-16 ***
DESTIN_SZQTSZ09 -0.3608995 0.0045591 -79.161 < 2e-16 ***
DESTIN_SZQTSZ10 -0.6856482 0.0047309 -144.930 < 2e-16 ***
DESTIN_SZQTSZ11 -0.1083675 0.0044482 -24.362 < 2e-16 ***
DESTIN_SZQTSZ12 -0.3829882 0.0053038 -72.209 < 2e-16 ***
DESTIN_SZQTSZ13 0.0493821 0.0040402 12.223 < 2e-16 ***
DESTIN_SZQTSZ14 -0.1559213 0.0044786 -34.815 < 2e-16 ***
DESTIN_SZQTSZ15 0.0229944 0.0054617 4.210 2.55e-05 ***
DESTIN_SZRCSZ01 -1.0723876 0.0052719 -203.415 < 2e-16 ***
DESTIN_SZRCSZ02 -2.3479621 0.0146221 -160.577 < 2e-16 ***
DESTIN_SZRCSZ03 -1.1278286 0.0070968 -158.922 < 2e-16 ***
DESTIN_SZRCSZ04 -2.4707507 0.0103338 -239.093 < 2e-16 ***
DESTIN_SZRCSZ05 -2.3692312 0.0097386 -243.283 < 2e-16 ***
DESTIN_SZRCSZ06 -2.1999028 0.0126499 -173.907 < 2e-16 ***
DESTIN_SZRCSZ08 -2.1638375 0.0106204 -203.744 < 2e-16 ***
DESTIN_SZRCSZ09 -1.6944317 0.0100154 -169.182 < 2e-16 ***
DESTIN_SZRCSZ10 -1.1881633 0.0051502 -230.703 < 2e-16 ***
DESTIN_SZRVSZ01 -2.1807893 0.0088319 -246.921 < 2e-16 ***
DESTIN_SZRVSZ02 -2.3688550 0.0118593 -199.747 < 2e-16 ***
DESTIN_SZRVSZ03 -2.4412943 0.0101240 -241.140 < 2e-16 ***
DESTIN_SZRVSZ04 -1.9470769 0.0126953 -153.370 < 2e-16 ***
DESTIN_SZRVSZ05 -2.0894753 0.0115159 -181.443 < 2e-16 ***
DESTIN_SZSBSZ01 -0.0373022 0.0047058 -7.927 2.25e-15 ***
DESTIN_SZSBSZ02 -0.9827428 0.0062625 -156.925 < 2e-16 ***
DESTIN_SZSBSZ03 0.7163226 0.0034531 207.446 < 2e-16 ***
DESTIN_SZSBSZ04 0.1235750 0.0044279 27.908 < 2e-16 ***
DESTIN_SZSBSZ05 -0.7382915 0.0057162 -129.158 < 2e-16 ***
DESTIN_SZSBSZ06 -2.5927009 0.0213720 -121.313 < 2e-16 ***
DESTIN_SZSBSZ07 -0.7182676 0.0152514 -47.095 < 2e-16 ***
DESTIN_SZSBSZ08 1.5018362 0.0042113 356.623 < 2e-16 ***
DESTIN_SZSBSZ09 0.8844992 0.0041582 212.713 < 2e-16 ***
DESTIN_SZSESZ02 -0.1318437 0.0038825 -33.959 < 2e-16 ***
DESTIN_SZSESZ03 0.6859754 0.0030859 222.293 < 2e-16 ***
DESTIN_SZSESZ04 -0.5494408 0.0044621 -123.136 < 2e-16 ***
DESTIN_SZSESZ05 -0.1377203 0.0038190 -36.062 < 2e-16 ***
DESTIN_SZSESZ06 -0.5091891 0.0047354 -107.528 < 2e-16 ***
DESTIN_SZSESZ07 -2.8997136 0.0196781 -147.358 < 2e-16 ***
DESTIN_SZSGSZ01 -0.4409464 0.0048826 -90.309 < 2e-16 ***
DESTIN_SZSGSZ02 0.0247217 0.0043341 5.704 1.17e-08 ***
DESTIN_SZSGSZ03 -0.3655516 0.0040393 -90.500 < 2e-16 ***
DESTIN_SZSGSZ04 -0.3761233 0.0040275 -93.388 < 2e-16 ***
DESTIN_SZSGSZ05 -2.2288402 0.0082919 -268.799 < 2e-16 ***
DESTIN_SZSGSZ06 0.4483146 0.0031507 142.291 < 2e-16 ***
DESTIN_SZSGSZ07 -0.4451934 0.0041120 -108.267 < 2e-16 ***
DESTIN_SZSISZ01 -1.0137317 0.0138601 -73.140 < 2e-16 ***
DESTIN_SZSKSZ01 -0.1076271 0.0059259 -18.162 < 2e-16 ***
DESTIN_SZSKSZ02 0.6887225 0.0042226 163.103 < 2e-16 ***
DESTIN_SZSKSZ03 -0.0129496 0.0050096 -2.585 0.00974 **
DESTIN_SZSKSZ04 -0.6781705 0.0128267 -52.872 < 2e-16 ***
DESTIN_SZSKSZ05 0.0433668 0.0095630 4.535 5.77e-06 ***
DESTIN_SZSLSZ01 -0.3657900 0.0065940 -55.473 < 2e-16 ***
DESTIN_SZSLSZ04 -0.4759952 0.0053388 -89.158 < 2e-16 ***
DESTIN_SZSRSZ01 -1.6889598 0.0064558 -261.619 < 2e-16 ***
DESTIN_SZSRSZ02 -1.6488407 0.0076359 -215.932 < 2e-16 ***
DESTIN_SZSRSZ03 -1.6281126 0.0069575 -234.008 < 2e-16 ***
DESTIN_SZSVSZ01 -2.1556160 0.0283356 -76.074 < 2e-16 ***
DESTIN_SZTHSZ01 -2.9255328 0.0338052 -86.541 < 2e-16 ***
DESTIN_SZTHSZ03 -1.9163698 0.0175137 -109.421 < 2e-16 ***
DESTIN_SZTHSZ04 -2.1355042 0.0175317 -121.808 < 2e-16 ***
DESTIN_SZTHSZ06 -1.5210235 0.0118415 -128.449 < 2e-16 ***
DESTIN_SZTMSZ01 -0.0005737 0.0043358 -0.132 0.89473
DESTIN_SZTMSZ02 1.7103810 0.0027607 619.542 < 2e-16 ***
DESTIN_SZTMSZ03 0.7838316 0.0031108 251.974 < 2e-16 ***
DESTIN_SZTMSZ04 0.8054132 0.0032377 248.758 < 2e-16 ***
DESTIN_SZTMSZ05 0.5909874 0.0040592 145.591 < 2e-16 ***
DESTIN_SZTNSZ01 -0.7528564 0.0045661 -164.880 < 2e-16 ***
DESTIN_SZTNSZ02 -1.5227198 0.0059542 -255.739 < 2e-16 ***
DESTIN_SZTNSZ03 -1.5568524 0.0073448 -211.966 < 2e-16 ***
DESTIN_SZTNSZ04 -1.1567363 0.0056275 -205.550 < 2e-16 ***
DESTIN_SZTPSZ01 -0.5816183 0.0045480 -127.883 < 2e-16 ***
DESTIN_SZTPSZ02 0.1851227 0.0031499 58.771 < 2e-16 ***
DESTIN_SZTPSZ03 -0.5136588 0.0045773 -112.218 < 2e-16 ***
DESTIN_SZTPSZ04 -1.5613506 0.0062201 -251.019 < 2e-16 ***
DESTIN_SZTPSZ05 -0.9004356 0.0048203 -186.800 < 2e-16 ***
DESTIN_SZTPSZ06 -0.2972237 0.0054240 -54.798 < 2e-16 ***
DESTIN_SZTPSZ07 -1.7714502 0.0089233 -198.520 < 2e-16 ***
DESTIN_SZTPSZ08 -1.3378522 0.0065895 -203.029 < 2e-16 ***
DESTIN_SZTPSZ09 -0.6317225 0.0049380 -127.931 < 2e-16 ***
DESTIN_SZTPSZ10 -0.6875757 0.0060740 -113.199 < 2e-16 ***
DESTIN_SZTPSZ11 -0.3770786 0.0040938 -92.110 < 2e-16 ***
DESTIN_SZTPSZ12 -0.8361943 0.0051187 -163.362 < 2e-16 ***
DESTIN_SZTSSZ01 -0.5639514 0.0187948 -30.006 < 2e-16 ***
DESTIN_SZTSSZ02 1.0132592 0.0074397 136.196 < 2e-16 ***
DESTIN_SZTSSZ03 1.5035121 0.0055094 272.901 < 2e-16 ***
DESTIN_SZTSSZ04 1.6312507 0.0055972 291.440 < 2e-16 ***
DESTIN_SZTSSZ05 1.8458811 0.0059379 310.864 < 2e-16 ***
DESTIN_SZTSSZ06 0.9211767 0.0095065 96.900 < 2e-16 ***
DESTIN_SZWCSZ01 1.2914958 0.0040273 320.689 < 2e-16 ***
DESTIN_SZWCSZ02 0.1942770 0.0078781 24.660 < 2e-16 ***
DESTIN_SZWCSZ03 -1.3867410 0.0247189 -56.100 < 2e-16 ***
DESTIN_SZWDSZ01 1.6801486 0.0029093 577.519 < 2e-16 ***
DESTIN_SZWDSZ02 -0.2711460 0.0047331 -57.287 < 2e-16 ***
DESTIN_SZWDSZ03 1.3042810 0.0030644 425.620 < 2e-16 ***
DESTIN_SZWDSZ04 0.1532090 0.0044491 34.436 < 2e-16 ***
DESTIN_SZWDSZ05 0.2225277 0.0042637 52.191 < 2e-16 ***
DESTIN_SZWDSZ06 0.5839541 0.0034028 171.609 < 2e-16 ***
DESTIN_SZWDSZ07 1.0576879 0.0045226 233.866 < 2e-16 ***
DESTIN_SZWDSZ08 0.7419707 0.0051173 144.994 < 2e-16 ***
DESTIN_SZWDSZ09 0.7615138 0.0036993 205.854 < 2e-16 ***
DESTIN_SZYSSZ01 1.2905530 0.0031461 410.210 < 2e-16 ***
DESTIN_SZYSSZ02 0.2413510 0.0041445 58.234 < 2e-16 ***
DESTIN_SZYSSZ03 -0.0572067 0.0043515 -13.147 < 2e-16 ***
DESTIN_SZYSSZ04 0.0286800 0.0041318 6.941 3.88e-12 ***
DESTIN_SZYSSZ05 -1.5645361 0.0086049 -181.819 < 2e-16 ***
DESTIN_SZYSSZ06 -1.2542840 0.0065130 -192.582 < 2e-16 ***
DESTIN_SZYSSZ07 -0.7731698 0.0077014 -100.394 < 2e-16 ***
DESTIN_SZYSSZ08 0.7327456 0.0032481 225.595 < 2e-16 ***
DESTIN_SZYSSZ09 0.3957525 0.0033560 117.923 < 2e-16 ***
log(ORIGIN_AGE25_64) 0.2264127 0.0001051 2154.705 < 2e-16 ***
log(dist) -0.7145351 0.0001033 -6918.842 < 2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
(Dispersion parameter for poisson family taken to be 1)
Null deviance: 102725286 on 20848 degrees of freedom
Residual deviance: 43341354 on 20537 degrees of freedom
AIC: 43473559
Number of Fisher Scoring iterations: 7Let’s check the R-square values of destination constrained SIM model this time.
CalcRSquared(decSIM$data$TRIPS, decSIM$fitted.values)[1] 0.49284967.6 Doubly Constrained Spatial Interaction Model
In this section, we will calibrate a doubly constrained SIM. For doubly constrained SIM, both the attractiveness at the destinations and the propulsiveness at the origins are considered. The model is typically expressed in the form of a distance function between the origin and destination.
dbcSIM <- glm(formula = TRIPS ~
ORIGIN_SZ +
DESTIN_SZ +
log(dist),
family = poisson(link = "log"),
data = SIM_data,
na.action = na.exclude)
summary(dbcSIM)
Call:
glm(formula = TRIPS ~ ORIGIN_SZ + DESTIN_SZ + log(dist), family = poisson(link = "log"),
data = SIM_data, na.action = na.exclude)
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 12.8286466 0.0037011 3466.136 < 2e-16 ***
ORIGIN_SZAMSZ02 1.0080358 0.0038100 264.576 < 2e-16 ***
ORIGIN_SZAMSZ03 0.5715510 0.0038787 147.356 < 2e-16 ***
ORIGIN_SZAMSZ04 0.1374926 0.0043172 31.847 < 2e-16 ***
ORIGIN_SZAMSZ05 0.0142679 0.0049202 2.900 0.00373 **
ORIGIN_SZAMSZ06 0.4955344 0.0045099 109.876 < 2e-16 ***
ORIGIN_SZAMSZ07 -0.8472976 0.0075058 -112.886 < 2e-16 ***
ORIGIN_SZAMSZ08 -0.6618457 0.0070454 -93.940 < 2e-16 ***
ORIGIN_SZAMSZ09 0.4153695 0.0046580 89.173 < 2e-16 ***
ORIGIN_SZAMSZ10 0.3732637 0.0041066 90.893 < 2e-16 ***
ORIGIN_SZAMSZ11 -1.5503320 0.0093025 -166.658 < 2e-16 ***
ORIGIN_SZAMSZ12 -1.7090680 0.0092239 -185.288 < 2e-16 ***
ORIGIN_SZBDSZ01 0.7695777 0.0037988 202.583 < 2e-16 ***
ORIGIN_SZBDSZ02 0.4022140 0.0043709 92.021 < 2e-16 ***
ORIGIN_SZBDSZ03 0.7807270 0.0038631 202.099 < 2e-16 ***
ORIGIN_SZBDSZ04 1.4078084 0.0033953 414.636 < 2e-16 ***
ORIGIN_SZBDSZ05 0.5084625 0.0038593 131.750 < 2e-16 ***
ORIGIN_SZBDSZ06 0.8299127 0.0039217 211.623 < 2e-16 ***
ORIGIN_SZBDSZ07 -0.8665200 0.0072092 -120.196 < 2e-16 ***
ORIGIN_SZBDSZ08 -0.6900134 0.0068115 -101.301 < 2e-16 ***
ORIGIN_SZBKSZ01 -0.1824104 0.0055044 -33.139 < 2e-16 ***
ORIGIN_SZBKSZ02 0.5411432 0.0044973 120.326 < 2e-16 ***
ORIGIN_SZBKSZ03 0.9599842 0.0041657 230.451 < 2e-16 ***
ORIGIN_SZBKSZ04 -0.1069206 0.0051199 -20.883 < 2e-16 ***
ORIGIN_SZBKSZ05 0.0733120 0.0048543 15.103 < 2e-16 ***
ORIGIN_SZBKSZ06 0.2512732 0.0050615 49.644 < 2e-16 ***
ORIGIN_SZBKSZ07 0.7846140 0.0038612 203.203 < 2e-16 ***
ORIGIN_SZBKSZ08 0.2624991 0.0045044 58.276 < 2e-16 ***
ORIGIN_SZBKSZ09 0.0705537 0.0047531 14.844 < 2e-16 ***
ORIGIN_SZBLSZ01 -1.7588408 0.0107643 -163.396 < 2e-16 ***
ORIGIN_SZBLSZ02 -2.9437818 0.0161582 -182.185 < 2e-16 ***
ORIGIN_SZBLSZ03 -5.1413715 0.0346478 -148.390 < 2e-16 ***
ORIGIN_SZBLSZ04 -2.3630948 0.0170632 -138.490 < 2e-16 ***
ORIGIN_SZBMSZ01 0.2236721 0.0041224 54.258 < 2e-16 ***
ORIGIN_SZBMSZ02 -1.0115464 0.0058265 -173.611 < 2e-16 ***
ORIGIN_SZBMSZ03 -0.0972568 0.0046902 -20.736 < 2e-16 ***
ORIGIN_SZBMSZ04 0.3084163 0.0041625 74.094 < 2e-16 ***
ORIGIN_SZBMSZ05 -1.0599524 0.0062863 -168.613 < 2e-16 ***
ORIGIN_SZBMSZ06 -1.3267571 0.0099021 -133.988 < 2e-16 ***
ORIGIN_SZBMSZ07 -0.3268714 0.0046780 -69.874 < 2e-16 ***
ORIGIN_SZBMSZ08 -0.1929591 0.0046241 -41.729 < 2e-16 ***
ORIGIN_SZBMSZ09 -0.7386238 0.0059190 -124.790 < 2e-16 ***
ORIGIN_SZBMSZ10 -0.8340636 0.0063085 -132.212 < 2e-16 ***
ORIGIN_SZBMSZ11 -0.5115401 0.0054974 -93.051 < 2e-16 ***
ORIGIN_SZBMSZ12 -0.8496475 0.0074682 -113.768 < 2e-16 ***
ORIGIN_SZBMSZ13 0.0000855 0.0046329 0.018 0.98528
ORIGIN_SZBMSZ14 -0.2930598 0.0055118 -53.169 < 2e-16 ***
ORIGIN_SZBMSZ15 0.0564604 0.0049621 11.378 < 2e-16 ***
ORIGIN_SZBMSZ16 -0.9920353 0.0062767 -158.051 < 2e-16 ***
ORIGIN_SZBMSZ17 -1.4210546 0.0096204 -147.712 < 2e-16 ***
ORIGIN_SZBPSZ01 0.5103191 0.0046953 108.687 < 2e-16 ***
ORIGIN_SZBPSZ02 0.7412988 0.0052841 140.289 < 2e-16 ***
ORIGIN_SZBPSZ03 0.9746610 0.0049525 196.802 < 2e-16 ***
ORIGIN_SZBPSZ04 0.7126562 0.0042822 166.425 < 2e-16 ***
ORIGIN_SZBPSZ05 0.5441855 0.0039826 136.639 < 2e-16 ***
ORIGIN_SZBPSZ06 -0.7842357 0.0068811 -113.970 < 2e-16 ***
ORIGIN_SZBPSZ07 -0.7499548 0.0069139 -108.471 < 2e-16 ***
ORIGIN_SZBSSZ01 0.0731246 0.0044884 16.292 < 2e-16 ***
ORIGIN_SZBSSZ02 0.5352792 0.0040103 133.476 < 2e-16 ***
ORIGIN_SZBSSZ03 0.2261247 0.0039791 56.829 < 2e-16 ***
ORIGIN_SZBTSZ01 -0.0477017 0.0044633 -10.687 < 2e-16 ***
ORIGIN_SZBTSZ02 -0.7958437 0.0062268 -127.809 < 2e-16 ***
ORIGIN_SZBTSZ03 -0.1054728 0.0047370 -22.266 < 2e-16 ***
ORIGIN_SZBTSZ04 -0.4790439 0.0082831 -57.834 < 2e-16 ***
ORIGIN_SZBTSZ05 -1.4056079 0.0088385 -159.032 < 2e-16 ***
ORIGIN_SZBTSZ06 -0.5347751 0.0058567 -91.309 < 2e-16 ***
ORIGIN_SZBTSZ07 -1.6335294 0.0092343 -176.898 < 2e-16 ***
ORIGIN_SZBTSZ08 -0.9386709 0.0072661 -129.185 < 2e-16 ***
ORIGIN_SZCCSZ01 -1.6704687 0.0114550 -145.829 < 2e-16 ***
ORIGIN_SZCHSZ01 -1.1822296 0.0095951 -123.212 < 2e-16 ***
ORIGIN_SZCHSZ02 -0.9169778 0.0069920 -131.148 < 2e-16 ***
ORIGIN_SZCHSZ03 -0.4786555 0.0049759 -96.195 < 2e-16 ***
ORIGIN_SZCKSZ01 0.4725580 0.0042831 110.330 < 2e-16 ***
ORIGIN_SZCKSZ02 1.1845462 0.0043794 270.482 < 2e-16 ***
ORIGIN_SZCKSZ03 0.8569898 0.0039849 215.060 < 2e-16 ***
ORIGIN_SZCKSZ04 1.8643409 0.0041375 450.597 < 2e-16 ***
ORIGIN_SZCKSZ05 1.5314798 0.0051465 297.577 < 2e-16 ***
ORIGIN_SZCKSZ06 1.1678178 0.0054526 214.174 < 2e-16 ***
ORIGIN_SZCLSZ01 -0.6635639 0.0060139 -110.338 < 2e-16 ***
ORIGIN_SZCLSZ02 -1.2339220 0.0107523 -114.759 < 2e-16 ***
ORIGIN_SZCLSZ03 -0.2975066 0.0055782 -53.334 < 2e-16 ***
ORIGIN_SZCLSZ04 0.8892746 0.0037698 235.891 < 2e-16 ***
ORIGIN_SZCLSZ05 -1.5469666 0.0107228 -144.268 < 2e-16 ***
ORIGIN_SZCLSZ06 0.8472355 0.0036524 231.965 < 2e-16 ***
ORIGIN_SZCLSZ07 -0.1036529 0.0046093 -22.488 < 2e-16 ***
ORIGIN_SZCLSZ08 0.2787814 0.0052161 53.446 < 2e-16 ***
ORIGIN_SZCLSZ09 -2.2763627 0.0143377 -158.767 < 2e-16 ***
ORIGIN_SZDTSZ01 -1.3542815 0.0070310 -192.617 < 2e-16 ***
ORIGIN_SZDTSZ02 -1.2946921 0.0065791 -196.790 < 2e-16 ***
ORIGIN_SZDTSZ03 -2.4887167 0.0130587 -190.579 < 2e-16 ***
ORIGIN_SZDTSZ04 -3.9155534 0.1187448 -32.975 < 2e-16 ***
ORIGIN_SZDTSZ05 -2.8513019 0.0222503 -128.147 < 2e-16 ***
ORIGIN_SZDTSZ06 -2.7709355 0.0172199 -160.915 < 2e-16 ***
ORIGIN_SZDTSZ07 -1.5728569 0.0199524 -78.831 < 2e-16 ***
ORIGIN_SZDTSZ08 -2.1236255 0.0099750 -212.894 < 2e-16 ***
ORIGIN_SZDTSZ09 -2.6310751 0.0206638 -127.328 < 2e-16 ***
ORIGIN_SZDTSZ10 -1.8811873 0.0108120 -173.991 < 2e-16 ***
ORIGIN_SZDTSZ11 -2.0080252 0.0109002 -184.219 < 2e-16 ***
ORIGIN_SZDTSZ12 -3.0383727 0.0237245 -128.069 < 2e-16 ***
ORIGIN_SZDTSZ13 -2.0325313 0.0125101 -162.471 < 2e-16 ***
ORIGIN_SZGLSZ01 -1.4410978 0.0073489 -196.098 < 2e-16 ***
ORIGIN_SZGLSZ02 0.2299975 0.0041793 55.032 < 2e-16 ***
ORIGIN_SZGLSZ03 0.0413994 0.0041648 9.940 < 2e-16 ***
ORIGIN_SZGLSZ04 0.9505227 0.0035341 268.960 < 2e-16 ***
ORIGIN_SZGLSZ05 0.5792325 0.0037121 156.041 < 2e-16 ***
ORIGIN_SZHGSZ01 0.1587759 0.0041046 38.682 < 2e-16 ***
ORIGIN_SZHGSZ02 0.6983976 0.0039539 176.637 < 2e-16 ***
ORIGIN_SZHGSZ03 0.4526194 0.0043124 104.959 < 2e-16 ***
ORIGIN_SZHGSZ04 0.9544562 0.0036843 259.059 < 2e-16 ***
ORIGIN_SZHGSZ05 1.2893079 0.0036401 354.192 < 2e-16 ***
ORIGIN_SZHGSZ06 0.1989697 0.0043796 45.431 < 2e-16 ***
ORIGIN_SZHGSZ07 0.6315046 0.0038285 164.948 < 2e-16 ***
ORIGIN_SZHGSZ08 0.1231067 0.0044066 27.937 < 2e-16 ***
ORIGIN_SZHGSZ09 -0.7501850 0.0060031 -124.965 < 2e-16 ***
ORIGIN_SZHGSZ10 -3.1265291 0.0391147 -79.932 < 2e-16 ***
ORIGIN_SZJESZ01 0.4288625 0.0043546 98.485 < 2e-16 ***
ORIGIN_SZJESZ02 0.3474167 0.0043021 80.756 < 2e-16 ***
ORIGIN_SZJESZ03 0.3648711 0.0045394 80.378 < 2e-16 ***
ORIGIN_SZJESZ04 -1.2601204 0.0079938 -157.638 < 2e-16 ***
ORIGIN_SZJESZ05 -1.9448605 0.0121589 -159.954 < 2e-16 ***
ORIGIN_SZJESZ06 0.2098534 0.0042352 49.550 < 2e-16 ***
ORIGIN_SZJESZ07 -1.7159895 0.0095137 -180.371 < 2e-16 ***
ORIGIN_SZJESZ08 -0.7169182 0.0087853 -81.605 < 2e-16 ***
ORIGIN_SZJESZ09 0.4193497 0.0044663 93.893 < 2e-16 ***
ORIGIN_SZJESZ10 -2.8678926 0.0176560 -162.431 < 2e-16 ***
ORIGIN_SZJESZ11 -3.2358276 0.0175961 -183.894 < 2e-16 ***
ORIGIN_SZJWSZ01 0.5101264 0.0058424 87.315 < 2e-16 ***
ORIGIN_SZJWSZ02 0.9610102 0.0040039 240.019 < 2e-16 ***
ORIGIN_SZJWSZ03 1.1521401 0.0037948 303.610 < 2e-16 ***
ORIGIN_SZJWSZ04 0.9109277 0.0039031 233.386 < 2e-16 ***
ORIGIN_SZJWSZ05 -1.6940137 0.0108284 -156.442 < 2e-16 ***
ORIGIN_SZJWSZ06 -1.3507828 0.0090717 -148.900 < 2e-16 ***
ORIGIN_SZJWSZ07 -2.4328195 0.0236166 -103.013 < 2e-16 ***
ORIGIN_SZJWSZ08 2.0095909 0.0038202 526.038 < 2e-16 ***
ORIGIN_SZJWSZ09 1.3978884 0.0035538 393.348 < 2e-16 ***
ORIGIN_SZKLSZ01 0.3121000 0.0040190 77.657 < 2e-16 ***
ORIGIN_SZKLSZ02 -0.3029470 0.0050001 -60.588 < 2e-16 ***
ORIGIN_SZKLSZ03 -0.2983065 0.0049845 -59.846 < 2e-16 ***
ORIGIN_SZKLSZ04 -1.3230496 0.0071793 -184.287 < 2e-16 ***
ORIGIN_SZKLSZ05 -0.5929324 0.0068531 -86.520 < 2e-16 ***
ORIGIN_SZKLSZ06 -0.3022173 0.0046926 -64.404 < 2e-16 ***
ORIGIN_SZKLSZ07 -0.6579762 0.0061338 -107.271 < 2e-16 ***
ORIGIN_SZKLSZ08 -0.7332902 0.0053679 -136.607 < 2e-16 ***
ORIGIN_SZKLSZ09 -1.2143607 0.0067143 -180.863 < 2e-16 ***
ORIGIN_SZLKSZ01 -2.4850760 0.0318602 -77.999 < 2e-16 ***
ORIGIN_SZMDSZ01 -1.5811172 0.0232731 -67.938 < 2e-16 ***
ORIGIN_SZMDSZ02 -1.0951370 0.0100671 -108.784 < 2e-16 ***
ORIGIN_SZMDSZ03 -1.3536856 0.0136002 -99.535 < 2e-16 ***
ORIGIN_SZMPSZ01 -0.8673314 0.0067115 -129.230 < 2e-16 ***
ORIGIN_SZMPSZ02 -0.4562093 0.0054836 -83.195 < 2e-16 ***
ORIGIN_SZMPSZ03 0.0916973 0.0043973 20.853 < 2e-16 ***
ORIGIN_SZMSSZ01 -7.9746474 0.5801853 -13.745 < 2e-16 ***
ORIGIN_SZMUSZ01 -0.9378756 0.0059225 -158.358 < 2e-16 ***
ORIGIN_SZMUSZ02 -2.7559471 0.0151242 -182.221 < 2e-16 ***
ORIGIN_SZMUSZ03 -1.5194920 0.0069417 -218.894 < 2e-16 ***
ORIGIN_SZNTSZ01 -2.2763222 0.0254083 -89.590 < 2e-16 ***
ORIGIN_SZNTSZ02 -2.2479701 0.0143513 -156.639 < 2e-16 ***
ORIGIN_SZNTSZ03 -0.4887519 0.0060138 -81.272 < 2e-16 ***
ORIGIN_SZNTSZ05 -2.6582077 0.0368387 -72.158 < 2e-16 ***
ORIGIN_SZNTSZ06 -3.0437598 0.0412123 -73.856 < 2e-16 ***
ORIGIN_SZNVSZ01 0.8054664 0.0036721 219.346 < 2e-16 ***
ORIGIN_SZNVSZ02 -0.3042935 0.0050050 -60.798 < 2e-16 ***
ORIGIN_SZNVSZ03 -1.0012104 0.0061807 -161.989 < 2e-16 ***
ORIGIN_SZNVSZ04 -1.0785055 0.0072052 -149.685 < 2e-16 ***
ORIGIN_SZNVSZ05 -2.3436312 0.0133544 -175.495 < 2e-16 ***
ORIGIN_SZORSZ01 -2.5128809 0.0270712 -92.825 < 2e-16 ***
ORIGIN_SZORSZ02 -1.0481000 0.0059300 -176.744 < 2e-16 ***
ORIGIN_SZORSZ03 -1.3795670 0.0070565 -195.503 < 2e-16 ***
ORIGIN_SZOTSZ01 -1.3155959 0.0073440 -179.138 < 2e-16 ***
ORIGIN_SZOTSZ02 -1.4783711 0.0082265 -179.708 < 2e-16 ***
ORIGIN_SZOTSZ03 -0.4297072 0.0053938 -79.666 < 2e-16 ***
ORIGIN_SZOTSZ04 -0.5129011 0.0083879 -61.147 < 2e-16 ***
ORIGIN_SZPGSZ01 0.6623215 0.0108073 61.285 < 2e-16 ***
ORIGIN_SZPGSZ02 -0.2388834 0.0058034 -41.163 < 2e-16 ***
ORIGIN_SZPGSZ03 1.0047350 0.0037882 265.226 < 2e-16 ***
ORIGIN_SZPGSZ04 1.1945429 0.0037758 316.368 < 2e-16 ***
ORIGIN_SZPGSZ05 0.6371460 0.0047143 135.151 < 2e-16 ***
ORIGIN_SZPLSZ01 -0.5935397 0.0082098 -72.297 < 2e-16 ***
ORIGIN_SZPLSZ02 -1.1081991 0.0115462 -95.979 < 2e-16 ***
ORIGIN_SZPLSZ03 -3.2681291 0.0321596 -101.622 < 2e-16 ***
ORIGIN_SZPLSZ04 -2.8842133 0.0371538 -77.629 < 2e-16 ***
ORIGIN_SZPLSZ05 -2.2182972 0.0175021 -126.745 < 2e-16 ***
ORIGIN_SZPNSZ01 1.0210921 0.0047269 216.018 < 2e-16 ***
ORIGIN_SZPNSZ02 -1.7788939 0.0113206 -157.138 < 2e-16 ***
ORIGIN_SZPNSZ03 -2.6923480 0.0177421 -151.749 < 2e-16 ***
ORIGIN_SZPNSZ04 -4.6135977 0.0256462 -179.894 < 2e-16 ***
ORIGIN_SZPNSZ05 -3.1475469 0.0200624 -156.888 < 2e-16 ***
ORIGIN_SZPRSZ01 -0.6516517 0.0098332 -66.270 < 2e-16 ***
ORIGIN_SZPRSZ02 0.9968282 0.0039774 250.620 < 2e-16 ***
ORIGIN_SZPRSZ03 0.4647176 0.0039697 117.067 < 2e-16 ***
ORIGIN_SZPRSZ04 -0.4244243 0.0064856 -65.441 < 2e-16 ***
ORIGIN_SZPRSZ05 1.1207043 0.0037897 295.726 < 2e-16 ***
ORIGIN_SZPRSZ06 -0.9757300 0.0070733 -137.946 < 2e-16 ***
ORIGIN_SZPRSZ07 -2.6051739 0.0172091 -151.383 < 2e-16 ***
ORIGIN_SZPRSZ08 0.0319438 0.0051812 6.165 7.03e-10 ***
ORIGIN_SZQTSZ01 0.0177811 0.0055780 3.188 0.00143 **
ORIGIN_SZQTSZ02 -0.4111178 0.0051582 -79.702 < 2e-16 ***
ORIGIN_SZQTSZ03 0.1146756 0.0047544 24.120 < 2e-16 ***
ORIGIN_SZQTSZ04 -1.0392568 0.0064735 -160.539 < 2e-16 ***
ORIGIN_SZQTSZ05 0.0815554 0.0047666 17.110 < 2e-16 ***
ORIGIN_SZQTSZ06 -0.2795372 0.0054275 -51.504 < 2e-16 ***
ORIGIN_SZQTSZ07 -1.2937812 0.0080004 -161.715 < 2e-16 ***
ORIGIN_SZQTSZ08 -0.3396025 0.0050344 -67.457 < 2e-16 ***
ORIGIN_SZQTSZ09 -0.4809761 0.0055455 -86.733 < 2e-16 ***
ORIGIN_SZQTSZ10 -0.3470876 0.0055552 -62.479 < 2e-16 ***
ORIGIN_SZQTSZ11 -1.4293687 0.0080162 -178.309 < 2e-16 ***
ORIGIN_SZQTSZ12 -0.7675547 0.0069079 -111.112 < 2e-16 ***
ORIGIN_SZQTSZ13 -0.2201302 0.0050975 -43.184 < 2e-16 ***
ORIGIN_SZQTSZ14 -1.3580373 0.0073862 -183.861 < 2e-16 ***
ORIGIN_SZQTSZ15 -1.0918424 0.0089246 -122.341 < 2e-16 ***
ORIGIN_SZRCSZ01 -0.3868138 0.0054130 -71.461 < 2e-16 ***
ORIGIN_SZRCSZ02 -1.9458168 0.0152168 -127.873 < 2e-16 ***
ORIGIN_SZRCSZ03 -0.9424591 0.0076035 -123.951 < 2e-16 ***
ORIGIN_SZRCSZ04 -2.0113838 0.0114492 -175.678 < 2e-16 ***
ORIGIN_SZRCSZ05 -2.3747795 0.0139802 -169.867 < 2e-16 ***
ORIGIN_SZRCSZ06 -0.1848037 0.0070599 -26.177 < 2e-16 ***
ORIGIN_SZRCSZ08 -2.4004375 0.0170293 -140.959 < 2e-16 ***
ORIGIN_SZRCSZ09 -1.6632104 0.0122261 -136.037 < 2e-16 ***
ORIGIN_SZRCSZ10 -1.4594717 0.0070136 -208.092 < 2e-16 ***
ORIGIN_SZRVSZ01 -2.5350492 0.0138487 -183.053 < 2e-16 ***
ORIGIN_SZRVSZ02 -0.7017935 0.0069268 -101.316 < 2e-16 ***
ORIGIN_SZRVSZ03 -1.3619629 0.0103197 -131.977 < 2e-16 ***
ORIGIN_SZRVSZ04 -1.5363160 0.0144993 -105.958 < 2e-16 ***
ORIGIN_SZRVSZ05 -1.7329389 0.0132236 -131.049 < 2e-16 ***
ORIGIN_SZSBSZ01 0.8699884 0.0049782 174.760 < 2e-16 ***
ORIGIN_SZSBSZ02 -0.6640952 0.0066767 -99.465 < 2e-16 ***
ORIGIN_SZSBSZ03 0.5219493 0.0041925 124.497 < 2e-16 ***
ORIGIN_SZSBSZ04 0.4026841 0.0048325 83.327 < 2e-16 ***
ORIGIN_SZSBSZ05 -0.0942157 0.0057793 -16.302 < 2e-16 ***
ORIGIN_SZSBSZ06 -1.0296651 0.0143389 -71.809 < 2e-16 ***
ORIGIN_SZSBSZ07 -0.1720061 0.0100076 -17.188 < 2e-16 ***
ORIGIN_SZSBSZ08 -2.0697634 0.0100444 -206.061 < 2e-16 ***
ORIGIN_SZSBSZ09 -1.1670954 0.0074298 -157.083 < 2e-16 ***
ORIGIN_SZSESZ02 1.1934115 0.0037168 321.087 < 2e-16 ***
ORIGIN_SZSESZ03 1.1018771 0.0035787 307.902 < 2e-16 ***
ORIGIN_SZSESZ04 1.1550026 0.0041073 281.209 < 2e-16 ***
ORIGIN_SZSESZ05 -0.1715327 0.0049537 -34.627 < 2e-16 ***
ORIGIN_SZSESZ06 1.1579692 0.0038919 297.534 < 2e-16 ***
ORIGIN_SZSESZ07 -1.9221669 0.0140370 -136.936 < 2e-16 ***
ORIGIN_SZSGSZ01 -0.8655946 0.0070989 -121.933 < 2e-16 ***
ORIGIN_SZSGSZ02 -1.2008708 0.0083503 -143.811 < 2e-16 ***
ORIGIN_SZSGSZ03 0.3270329 0.0044015 74.300 < 2e-16 ***
ORIGIN_SZSGSZ04 0.4073009 0.0040513 100.537 < 2e-16 ***
ORIGIN_SZSGSZ05 -1.5450851 0.0086361 -178.910 < 2e-16 ***
ORIGIN_SZSGSZ06 0.3411563 0.0038657 88.253 < 2e-16 ***
ORIGIN_SZSGSZ07 -0.4503131 0.0050675 -88.863 < 2e-16 ***
ORIGIN_SZSKSZ01 -0.2306384 0.0069419 -33.224 < 2e-16 ***
ORIGIN_SZSKSZ02 0.0951760 0.0052541 18.115 < 2e-16 ***
ORIGIN_SZSKSZ03 -0.4150873 0.0064099 -64.758 < 2e-16 ***
ORIGIN_SZSKSZ04 -2.0984590 0.0235127 -89.248 < 2e-16 ***
ORIGIN_SZSKSZ05 -1.1052254 0.0146423 -75.482 < 2e-16 ***
ORIGIN_SZSLSZ01 -2.9910706 0.0260625 -114.765 < 2e-16 ***
ORIGIN_SZSLSZ04 -0.3299003 0.0059429 -55.511 < 2e-16 ***
ORIGIN_SZSRSZ01 -1.1501834 0.0074256 -154.894 < 2e-16 ***
ORIGIN_SZSRSZ02 -1.4204906 0.0073588 -193.032 < 2e-16 ***
ORIGIN_SZSRSZ03 -2.3720718 0.0158444 -149.710 < 2e-16 ***
ORIGIN_SZSVSZ01 -2.2152636 0.0418137 -52.979 < 2e-16 ***
ORIGIN_SZTHSZ01 -3.3790533 0.0529657 -63.797 < 2e-16 ***
ORIGIN_SZTHSZ03 -0.8874917 0.0127598 -69.554 < 2e-16 ***
ORIGIN_SZTHSZ04 -2.0903014 0.0252899 -82.654 < 2e-16 ***
ORIGIN_SZTHSZ06 -1.3866058 0.0116481 -119.042 < 2e-16 ***
ORIGIN_SZTMSZ01 0.7630424 0.0043488 175.460 < 2e-16 ***
ORIGIN_SZTMSZ02 1.5946931 0.0033379 477.751 < 2e-16 ***
ORIGIN_SZTMSZ03 1.1501736 0.0035485 324.130 < 2e-16 ***
ORIGIN_SZTMSZ04 0.4101586 0.0041431 98.999 < 2e-16 ***
ORIGIN_SZTMSZ05 -0.8927171 0.0065734 -135.807 < 2e-16 ***
ORIGIN_SZTNSZ01 -0.9507277 0.0063118 -150.627 < 2e-16 ***
ORIGIN_SZTNSZ02 -0.7253285 0.0059708 -121.480 < 2e-16 ***
ORIGIN_SZTNSZ03 -1.2663473 0.0079725 -158.840 < 2e-16 ***
ORIGIN_SZTNSZ04 -0.3659710 0.0057831 -63.283 < 2e-16 ***
ORIGIN_SZTPSZ01 -0.4339873 0.0052290 -82.995 < 2e-16 ***
ORIGIN_SZTPSZ02 0.5099532 0.0037079 137.533 < 2e-16 ***
ORIGIN_SZTPSZ03 -0.4838407 0.0052883 -91.493 < 2e-16 ***
ORIGIN_SZTPSZ04 -0.0873479 0.0049189 -17.758 < 2e-16 ***
ORIGIN_SZTPSZ05 0.0550444 0.0051471 10.694 < 2e-16 ***
ORIGIN_SZTPSZ06 0.5842031 0.0057890 100.916 < 2e-16 ***
ORIGIN_SZTPSZ07 0.0566139 0.0051614 10.969 < 2e-16 ***
ORIGIN_SZTPSZ08 -0.4241616 0.0067731 -62.624 < 2e-16 ***
ORIGIN_SZTPSZ09 -0.4978890 0.0054959 -90.593 < 2e-16 ***
ORIGIN_SZTPSZ10 -0.1701460 0.0055506 -30.654 < 2e-16 ***
ORIGIN_SZTPSZ11 0.3031065 0.0043540 69.615 < 2e-16 ***
ORIGIN_SZTPSZ12 -0.4717124 0.0054169 -87.082 < 2e-16 ***
ORIGIN_SZTSSZ01 -3.3300478 0.0404286 -82.369 < 2e-16 ***
ORIGIN_SZTSSZ02 0.2307557 0.0078216 29.502 < 2e-16 ***
ORIGIN_SZTSSZ03 0.0754197 0.0076139 9.906 < 2e-16 ***
ORIGIN_SZTSSZ04 -0.4032132 0.0082678 -48.769 < 2e-16 ***
ORIGIN_SZTSSZ05 -2.6194676 0.0132954 -197.021 < 2e-16 ***
ORIGIN_SZTSSZ06 -3.1081134 0.0189244 -164.238 < 2e-16 ***
ORIGIN_SZWCSZ01 -0.6246864 0.0068479 -91.223 < 2e-16 ***
ORIGIN_SZWCSZ02 -2.9493271 0.0281276 -104.855 < 2e-16 ***
ORIGIN_SZWCSZ03 -4.8725330 0.1400733 -34.786 < 2e-16 ***
ORIGIN_SZWDSZ01 0.8601421 0.0036482 235.770 < 2e-16 ***
ORIGIN_SZWDSZ02 0.9493368 0.0041702 227.646 < 2e-16 ***
ORIGIN_SZWDSZ03 1.6444988 0.0037658 436.695 < 2e-16 ***
ORIGIN_SZWDSZ04 1.1133839 0.0045495 244.728 < 2e-16 ***
ORIGIN_SZWDSZ05 0.4227436 0.0044126 95.803 < 2e-16 ***
ORIGIN_SZWDSZ06 0.9470220 0.0041164 230.060 < 2e-16 ***
ORIGIN_SZWDSZ07 -0.4126556 0.0062927 -65.576 < 2e-16 ***
ORIGIN_SZWDSZ08 -0.8643052 0.0069038 -125.193 < 2e-16 ***
ORIGIN_SZWDSZ09 1.5889438 0.0039038 407.025 < 2e-16 ***
ORIGIN_SZYSSZ01 -0.5072205 0.0048488 -104.607 < 2e-16 ***
ORIGIN_SZYSSZ02 0.9660373 0.0045272 213.383 < 2e-16 ***
ORIGIN_SZYSSZ03 2.2755249 0.0037466 607.360 < 2e-16 ***
ORIGIN_SZYSSZ04 0.9319198 0.0038859 239.823 < 2e-16 ***
ORIGIN_SZYSSZ05 0.4672048 0.0047129 99.133 < 2e-16 ***
ORIGIN_SZYSSZ06 -0.5368483 0.0075816 -70.809 < 2e-16 ***
ORIGIN_SZYSSZ07 -0.4486462 0.0077892 -57.598 < 2e-16 ***
ORIGIN_SZYSSZ08 -0.4125001 0.0053438 -77.192 < 2e-16 ***
ORIGIN_SZYSSZ09 1.1701985 0.0037680 310.565 < 2e-16 ***
DESTIN_SZAMSZ02 -0.0319791 0.0036566 -8.746 < 2e-16 ***
DESTIN_SZAMSZ03 0.1016514 0.0035758 28.428 < 2e-16 ***
DESTIN_SZAMSZ04 -0.9147795 0.0052425 -174.492 < 2e-16 ***
DESTIN_SZAMSZ05 -0.8526911 0.0049581 -171.981 < 2e-16 ***
DESTIN_SZAMSZ06 -0.6910791 0.0049061 -140.861 < 2e-16 ***
DESTIN_SZAMSZ07 -1.6280332 0.0086015 -189.272 < 2e-16 ***
DESTIN_SZAMSZ08 -0.8158992 0.0055327 -147.469 < 2e-16 ***
DESTIN_SZAMSZ09 -1.0414618 0.0053076 -196.221 < 2e-16 ***
DESTIN_SZAMSZ10 0.0835468 0.0036462 22.913 < 2e-16 ***
DESTIN_SZAMSZ11 -0.0249127 0.0063331 -3.934 8.36e-05 ***
DESTIN_SZAMSZ12 0.0550569 0.0045492 12.103 < 2e-16 ***
DESTIN_SZBDSZ01 0.3940684 0.0033449 117.813 < 2e-16 ***
DESTIN_SZBDSZ02 -0.3124403 0.0042357 -73.764 < 2e-16 ***
DESTIN_SZBDSZ03 -0.1153709 0.0037433 -30.820 < 2e-16 ***
DESTIN_SZBDSZ04 0.6957570 0.0030808 225.837 < 2e-16 ***
DESTIN_SZBDSZ05 0.4175773 0.0033945 123.017 < 2e-16 ***
DESTIN_SZBDSZ06 -0.0382765 0.0038547 -9.930 < 2e-16 ***
DESTIN_SZBDSZ07 -0.6175696 0.0077178 -80.019 < 2e-16 ***
DESTIN_SZBDSZ08 -1.4147210 0.0076463 -185.020 < 2e-16 ***
DESTIN_SZBKSZ01 -1.2162782 0.0055018 -221.070 < 2e-16 ***
DESTIN_SZBKSZ02 -0.3822412 0.0046338 -82.489 < 2e-16 ***
DESTIN_SZBKSZ03 -0.8325920 0.0047788 -174.228 < 2e-16 ***
DESTIN_SZBKSZ04 -0.0661815 0.0042274 -15.655 < 2e-16 ***
DESTIN_SZBKSZ05 -0.6456351 0.0047870 -134.873 < 2e-16 ***
DESTIN_SZBKSZ06 -1.0680079 0.0054193 -197.076 < 2e-16 ***
DESTIN_SZBKSZ07 0.0951491 0.0035956 26.463 < 2e-16 ***
DESTIN_SZBKSZ08 -1.0880412 0.0059037 -184.299 < 2e-16 ***
DESTIN_SZBKSZ09 -0.1975852 0.0042692 -46.282 < 2e-16 ***
DESTIN_SZBLSZ01 -0.6962720 0.0059797 -116.440 < 2e-16 ***
DESTIN_SZBLSZ02 0.4264182 0.0058655 72.700 < 2e-16 ***
DESTIN_SZBLSZ03 1.7024805 0.0063957 266.193 < 2e-16 ***
DESTIN_SZBLSZ04 -0.3596634 0.0113939 -31.566 < 2e-16 ***
DESTIN_SZBMSZ01 -0.0495675 0.0038547 -12.859 < 2e-16 ***
DESTIN_SZBMSZ02 -0.2008501 0.0040124 -50.057 < 2e-16 ***
DESTIN_SZBMSZ03 -0.5693816 0.0049240 -115.633 < 2e-16 ***
DESTIN_SZBMSZ04 -0.2923587 0.0042404 -68.946 < 2e-16 ***
DESTIN_SZBMSZ05 -0.3968469 0.0050687 -78.294 < 2e-16 ***
DESTIN_SZBMSZ06 -1.1050662 0.0086155 -128.265 < 2e-16 ***
DESTIN_SZBMSZ07 0.1429083 0.0037470 38.139 < 2e-16 ***
DESTIN_SZBMSZ08 -0.7470911 0.0049200 -151.848 < 2e-16 ***
DESTIN_SZBMSZ09 -1.5595446 0.0077563 -201.069 < 2e-16 ***
DESTIN_SZBMSZ10 -1.0747502 0.0060958 -176.310 < 2e-16 ***
DESTIN_SZBMSZ11 -1.2214197 0.0060970 -200.332 < 2e-16 ***
DESTIN_SZBMSZ12 -0.7485205 0.0068683 -108.981 < 2e-16 ***
DESTIN_SZBMSZ13 0.0065838 0.0040488 1.626 0.10393
DESTIN_SZBMSZ14 -0.6236215 0.0063521 -98.176 < 2e-16 ***
DESTIN_SZBMSZ15 -0.9548592 0.0058947 -161.986 < 2e-16 ***
DESTIN_SZBMSZ16 -1.2784708 0.0062068 -205.979 < 2e-16 ***
DESTIN_SZBMSZ17 -1.3618668 0.0073263 -185.886 < 2e-16 ***
DESTIN_SZBPSZ01 -0.7964555 0.0046774 -170.276 < 2e-16 ***
DESTIN_SZBPSZ02 -1.6684606 0.0074592 -223.680 < 2e-16 ***
DESTIN_SZBPSZ03 -1.5247773 0.0070197 -217.213 < 2e-16 ***
DESTIN_SZBPSZ04 -0.8658276 0.0051348 -168.620 < 2e-16 ***
DESTIN_SZBPSZ05 0.2883371 0.0034728 83.027 < 2e-16 ***
DESTIN_SZBPSZ06 -0.9232838 0.0069064 -133.685 < 2e-16 ***
DESTIN_SZBPSZ07 -0.5162287 0.0066217 -77.960 < 2e-16 ***
DESTIN_SZBSSZ01 0.0294154 0.0038429 7.655 1.94e-14 ***
DESTIN_SZBSSZ02 -0.7819237 0.0044253 -176.693 < 2e-16 ***
DESTIN_SZBSSZ03 0.3757079 0.0032714 114.848 < 2e-16 ***
DESTIN_SZBTSZ01 0.3352039 0.0035069 95.583 < 2e-16 ***
DESTIN_SZBTSZ02 -0.6894708 0.0057279 -120.370 < 2e-16 ***
DESTIN_SZBTSZ03 0.1822522 0.0039813 45.777 < 2e-16 ***
DESTIN_SZBTSZ04 -1.2274421 0.0083292 -147.366 < 2e-16 ***
DESTIN_SZBTSZ05 -0.3918730 0.0058027 -67.532 < 2e-16 ***
DESTIN_SZBTSZ06 -0.6453936 0.0051819 -124.548 < 2e-16 ***
DESTIN_SZBTSZ07 -1.6158463 0.0083014 -194.648 < 2e-16 ***
DESTIN_SZBTSZ08 -0.7876561 0.0071096 -110.788 < 2e-16 ***
DESTIN_SZCCSZ01 -0.5225479 0.0059473 -87.863 < 2e-16 ***
DESTIN_SZCHSZ01 -1.1640940 0.0076489 -152.190 < 2e-16 ***
DESTIN_SZCHSZ02 -0.0413568 0.0047475 -8.711 < 2e-16 ***
DESTIN_SZCHSZ03 1.2674388 0.0034117 371.493 < 2e-16 ***
DESTIN_SZCKSZ01 -0.4162431 0.0043038 -96.715 < 2e-16 ***
DESTIN_SZCKSZ02 -0.9357724 0.0047503 -196.992 < 2e-16 ***
DESTIN_SZCKSZ03 0.2640587 0.0035530 74.320 < 2e-16 ***
DESTIN_SZCKSZ04 -1.5691123 0.0055248 -284.010 < 2e-16 ***
DESTIN_SZCKSZ05 -1.2783168 0.0063768 -200.463 < 2e-16 ***
DESTIN_SZCKSZ06 0.1064760 0.0050942 20.901 < 2e-16 ***
DESTIN_SZCLSZ01 0.2288951 0.0041424 55.256 < 2e-16 ***
DESTIN_SZCLSZ02 -2.1703685 0.0113247 -191.649 < 2e-16 ***
DESTIN_SZCLSZ03 -0.7991076 0.0062298 -128.272 < 2e-16 ***
DESTIN_SZCLSZ04 -0.0559249 0.0037892 -14.759 < 2e-16 ***
DESTIN_SZCLSZ05 -0.9646551 0.0072269 -133.481 < 2e-16 ***
DESTIN_SZCLSZ06 0.0650755 0.0035542 18.309 < 2e-16 ***
DESTIN_SZCLSZ07 -0.4914307 0.0046185 -106.405 < 2e-16 ***
DESTIN_SZCLSZ08 -0.4277559 0.0052051 -82.180 < 2e-16 ***
DESTIN_SZCLSZ09 0.3590441 0.0056904 63.097 < 2e-16 ***
DESTIN_SZDTSZ01 -0.5937033 0.0045634 -130.101 < 2e-16 ***
DESTIN_SZDTSZ02 -0.8038994 0.0044412 -181.011 < 2e-16 ***
DESTIN_SZDTSZ03 -1.0521104 0.0052893 -198.912 < 2e-16 ***
DESTIN_SZDTSZ04 -0.5710812 0.0113649 -50.249 < 2e-16 ***
DESTIN_SZDTSZ05 -0.6460195 0.0087850 -73.537 < 2e-16 ***
DESTIN_SZDTSZ06 -1.1079338 0.0059052 -187.620 < 2e-16 ***
DESTIN_SZDTSZ07 -1.9259046 0.0179834 -107.093 < 2e-16 ***
DESTIN_SZDTSZ08 -0.5847820 0.0043016 -135.944 < 2e-16 ***
DESTIN_SZDTSZ09 -1.7043172 0.0097281 -175.194 < 2e-16 ***
DESTIN_SZDTSZ10 -1.1931992 0.0076743 -155.480 < 2e-16 ***
DESTIN_SZDTSZ11 -0.6995185 0.0045569 -153.508 < 2e-16 ***
DESTIN_SZDTSZ12 -2.3706613 0.0148907 -159.204 < 2e-16 ***
DESTIN_SZDTSZ13 -1.7773676 0.0095016 -187.061 < 2e-16 ***
DESTIN_SZGLSZ01 0.1024918 0.0042575 24.073 < 2e-16 ***
DESTIN_SZGLSZ02 -0.2417847 0.0039771 -60.794 < 2e-16 ***
DESTIN_SZGLSZ03 0.3842995 0.0033403 115.051 < 2e-16 ***
DESTIN_SZGLSZ04 0.3044211 0.0033045 92.123 < 2e-16 ***
DESTIN_SZGLSZ05 0.1933056 0.0033871 57.071 < 2e-16 ***
DESTIN_SZHGSZ01 0.2764496 0.0033162 83.364 < 2e-16 ***
DESTIN_SZHGSZ02 -0.7709951 0.0045450 -169.637 < 2e-16 ***
DESTIN_SZHGSZ03 -1.2473364 0.0054264 -229.865 < 2e-16 ***
DESTIN_SZHGSZ04 -0.5027846 0.0038772 -129.679 < 2e-16 ***
DESTIN_SZHGSZ05 -0.5550290 0.0039144 -141.792 < 2e-16 ***
DESTIN_SZHGSZ06 -0.7505517 0.0045232 -165.934 < 2e-16 ***
DESTIN_SZHGSZ07 0.1050218 0.0035192 29.843 < 2e-16 ***
DESTIN_SZHGSZ08 -0.3181972 0.0041026 -77.560 < 2e-16 ***
DESTIN_SZHGSZ09 -0.0295870 0.0042751 -6.921 4.49e-12 ***
DESTIN_SZHGSZ10 -3.4766582 0.0292962 -118.673 < 2e-16 ***
DESTIN_SZJESZ01 -0.3714994 0.0044685 -83.137 < 2e-16 ***
DESTIN_SZJESZ02 -0.6269424 0.0044625 -140.490 < 2e-16 ***
DESTIN_SZJESZ03 -0.7949010 0.0049051 -162.056 < 2e-16 ***
DESTIN_SZJESZ04 -0.2192921 0.0053043 -41.343 < 2e-16 ***
DESTIN_SZJESZ05 -0.9005872 0.0076649 -117.495 < 2e-16 ***
DESTIN_SZJESZ06 0.1989816 0.0036183 54.993 < 2e-16 ***
DESTIN_SZJESZ07 -0.9521763 0.0063393 -150.202 < 2e-16 ***
DESTIN_SZJESZ08 -0.9160649 0.0067272 -136.174 < 2e-16 ***
DESTIN_SZJESZ09 -0.5341769 0.0049586 -107.727 < 2e-16 ***
DESTIN_SZJESZ10 0.6130739 0.0065635 93.406 < 2e-16 ***
DESTIN_SZJESZ11 1.0791907 0.0059048 182.765 < 2e-16 ***
DESTIN_SZJWSZ01 -0.8762600 0.0058197 -150.567 < 2e-16 ***
DESTIN_SZJWSZ02 -0.8067463 0.0047679 -169.204 < 2e-16 ***
DESTIN_SZJWSZ03 0.2259927 0.0036892 61.257 < 2e-16 ***
DESTIN_SZJWSZ04 0.7349551 0.0034603 212.398 < 2e-16 ***
DESTIN_SZJWSZ05 -0.4431286 0.0051966 -85.274 < 2e-16 ***
DESTIN_SZJWSZ06 -0.2577411 0.0048291 -53.373 < 2e-16 ***
DESTIN_SZJWSZ07 -1.5867481 0.0197637 -80.286 < 2e-16 ***
DESTIN_SZJWSZ08 -0.7365217 0.0044349 -166.074 < 2e-16 ***
DESTIN_SZJWSZ09 0.9154590 0.0031699 288.796 < 2e-16 ***
DESTIN_SZKLSZ01 -0.4509526 0.0041850 -107.755 < 2e-16 ***
DESTIN_SZKLSZ02 -0.6645796 0.0047718 -139.273 < 2e-16 ***
DESTIN_SZKLSZ03 -1.0829651 0.0052448 -206.484 < 2e-16 ***
DESTIN_SZKLSZ04 -1.5785854 0.0067942 -232.343 < 2e-16 ***
DESTIN_SZKLSZ05 -0.9395983 0.0069813 -134.588 < 2e-16 ***
DESTIN_SZKLSZ06 -0.6884537 0.0045843 -150.177 < 2e-16 ***
DESTIN_SZKLSZ07 -0.8218312 0.0051648 -159.123 < 2e-16 ***
DESTIN_SZKLSZ08 0.0160975 0.0037377 4.307 1.66e-05 ***
DESTIN_SZKLSZ09 -1.4894050 0.0066663 -223.422 < 2e-16 ***
DESTIN_SZLKSZ01 -2.0028692 0.0193335 -103.596 < 2e-16 ***
DESTIN_SZMDSZ01 -1.5414873 0.0170599 -90.357 < 2e-16 ***
DESTIN_SZMDSZ02 -1.2259043 0.0095505 -128.360 < 2e-16 ***
DESTIN_SZMDSZ03 -2.2796718 0.0213844 -106.604 < 2e-16 ***
DESTIN_SZMPSZ01 -0.9703954 0.0067630 -143.486 < 2e-16 ***
DESTIN_SZMPSZ02 -0.7245757 0.0049650 -145.935 < 2e-16 ***
DESTIN_SZMPSZ03 -0.1246202 0.0040772 -30.565 < 2e-16 ***
DESTIN_SZMSSZ01 -1.3125115 0.0744552 -17.628 < 2e-16 ***
DESTIN_SZMUSZ01 -1.0455081 0.0050085 -208.748 < 2e-16 ***
DESTIN_SZMUSZ02 -1.0149063 0.0071783 -141.386 < 2e-16 ***
DESTIN_SZMUSZ03 -1.1182935 0.0050393 -221.916 < 2e-16 ***
DESTIN_SZNTSZ01 -2.4064130 0.0221736 -108.526 < 2e-16 ***
DESTIN_SZNTSZ02 -1.7009221 0.0090669 -187.597 < 2e-16 ***
DESTIN_SZNTSZ03 -1.0877002 0.0062883 -172.972 < 2e-16 ***
DESTIN_SZNTSZ05 -1.7151117 0.0165765 -103.466 < 2e-16 ***
DESTIN_SZNTSZ06 -3.0200774 0.0275709 -109.538 < 2e-16 ***
DESTIN_SZNVSZ01 -0.2454668 0.0036909 -66.506 < 2e-16 ***
DESTIN_SZNVSZ02 -0.2903825 0.0042079 -69.009 < 2e-16 ***
DESTIN_SZNVSZ03 -0.3395519 0.0044344 -76.573 < 2e-16 ***
DESTIN_SZNVSZ04 -1.8639054 0.0087534 -212.936 < 2e-16 ***
DESTIN_SZNVSZ05 -1.5570749 0.0076006 -204.862 < 2e-16 ***
DESTIN_SZORSZ01 -1.3851704 0.0174520 -79.370 < 2e-16 ***
DESTIN_SZORSZ02 0.0387143 0.0038109 10.159 < 2e-16 ***
DESTIN_SZORSZ03 -0.7512487 0.0048856 -153.769 < 2e-16 ***
DESTIN_SZOTSZ01 -1.1063264 0.0062492 -177.035 < 2e-16 ***
DESTIN_SZOTSZ02 -0.4423678 0.0054952 -80.500 < 2e-16 ***
DESTIN_SZOTSZ03 -1.2143759 0.0058722 -206.799 < 2e-16 ***
DESTIN_SZOTSZ04 -1.4630078 0.0083582 -175.038 < 2e-16 ***
DESTIN_SZPGSZ01 -2.4794604 0.0166041 -149.328 < 2e-16 ***
DESTIN_SZPGSZ02 -0.8478807 0.0055376 -153.114 < 2e-16 ***
DESTIN_SZPGSZ03 0.2296091 0.0034558 66.442 < 2e-16 ***
DESTIN_SZPGSZ04 -0.2862819 0.0038761 -73.859 < 2e-16 ***
DESTIN_SZPGSZ05 -1.1176630 0.0064187 -174.125 < 2e-16 ***
DESTIN_SZPLSZ01 -0.6107283 0.0060482 -100.977 < 2e-16 ***
DESTIN_SZPLSZ02 -1.7030140 0.0107340 -158.656 < 2e-16 ***
DESTIN_SZPLSZ03 -0.1742585 0.0086479 -20.150 < 2e-16 ***
DESTIN_SZPLSZ04 -0.0891321 0.0082045 -10.864 < 2e-16 ***
DESTIN_SZPLSZ05 -0.8206606 0.0100449 -81.699 < 2e-16 ***
DESTIN_SZPNSZ01 -0.3076205 0.0049172 -62.560 < 2e-16 ***
DESTIN_SZPNSZ02 0.9275884 0.0066573 139.334 < 2e-16 ***
DESTIN_SZPNSZ03 0.0950709 0.0067179 14.152 < 2e-16 ***
DESTIN_SZPNSZ04 1.7869438 0.0073358 243.593 < 2e-16 ***
DESTIN_SZPNSZ05 0.9944933 0.0106714 93.193 < 2e-16 ***
DESTIN_SZPRSZ01 -1.0883324 0.0063979 -170.109 < 2e-16 ***
DESTIN_SZPRSZ02 -0.4975984 0.0044093 -112.853 < 2e-16 ***
DESTIN_SZPRSZ03 0.6272557 0.0033042 189.835 < 2e-16 ***
DESTIN_SZPRSZ04 -0.6356411 0.0074834 -84.940 < 2e-16 ***
DESTIN_SZPRSZ05 -0.3763500 0.0041476 -90.740 < 2e-16 ***
DESTIN_SZPRSZ06 0.1841643 0.0043211 42.620 < 2e-16 ***
DESTIN_SZPRSZ07 -1.0653637 0.0100317 -106.200 < 2e-16 ***
DESTIN_SZPRSZ08 -0.8634218 0.0057522 -150.102 < 2e-16 ***
DESTIN_SZQTSZ01 -1.5137864 0.0086259 -175.493 < 2e-16 ***
DESTIN_SZQTSZ02 -1.2508941 0.0061698 -202.743 < 2e-16 ***
DESTIN_SZQTSZ03 -0.6998385 0.0055441 -126.232 < 2e-16 ***
DESTIN_SZQTSZ04 -0.7328309 0.0056708 -129.228 < 2e-16 ***
DESTIN_SZQTSZ05 -0.6554222 0.0050774 -129.087 < 2e-16 ***
DESTIN_SZQTSZ06 -0.8980212 0.0053673 -167.315 < 2e-16 ***
DESTIN_SZQTSZ07 -1.5404086 0.0090090 -170.986 < 2e-16 ***
DESTIN_SZQTSZ08 0.2336931 0.0039637 58.958 < 2e-16 ***
DESTIN_SZQTSZ09 -0.3835458 0.0046844 -81.877 < 2e-16 ***
DESTIN_SZQTSZ10 -0.4984970 0.0048580 -102.614 < 2e-16 ***
DESTIN_SZQTSZ11 -0.0144119 0.0045726 -3.152 0.00162 **
DESTIN_SZQTSZ12 -0.2262282 0.0055135 -41.032 < 2e-16 ***
DESTIN_SZQTSZ13 0.1474976 0.0042318 34.855 < 2e-16 ***
DESTIN_SZQTSZ14 0.1112488 0.0046067 24.149 < 2e-16 ***
DESTIN_SZQTSZ15 0.0964562 0.0058069 16.611 < 2e-16 ***
DESTIN_SZRCSZ01 -0.8259567 0.0053674 -153.885 < 2e-16 ***
DESTIN_SZRCSZ02 -2.2641402 0.0146502 -154.547 < 2e-16 ***
DESTIN_SZRCSZ03 -1.0174766 0.0072116 -141.089 < 2e-16 ***
DESTIN_SZRCSZ04 -2.1520712 0.0103677 -207.574 < 2e-16 ***
DESTIN_SZRCSZ05 -2.1798021 0.0097692 -223.130 < 2e-16 ***
DESTIN_SZRCSZ06 -1.8927620 0.0126836 -149.229 < 2e-16 ***
DESTIN_SZRCSZ08 -1.5944519 0.0107021 -148.985 < 2e-16 ***
DESTIN_SZRCSZ09 -1.3664499 0.0100596 -135.835 < 2e-16 ***
DESTIN_SZRCSZ10 -0.8106750 0.0052088 -155.635 < 2e-16 ***
DESTIN_SZRVSZ01 -1.6488292 0.0088804 -185.671 < 2e-16 ***
DESTIN_SZRVSZ02 -2.2122648 0.0119376 -185.319 < 2e-16 ***
DESTIN_SZRVSZ03 -1.9433384 0.0101965 -190.588 < 2e-16 ***
DESTIN_SZRVSZ04 -1.5638031 0.0127438 -122.711 < 2e-16 ***
DESTIN_SZRVSZ05 -1.3421459 0.0117115 -114.601 < 2e-16 ***
DESTIN_SZSBSZ01 -0.3826176 0.0054100 -70.725 < 2e-16 ***
DESTIN_SZSBSZ02 -0.9281934 0.0064270 -144.421 < 2e-16 ***
DESTIN_SZSBSZ03 0.6003815 0.0038044 157.811 < 2e-16 ***
DESTIN_SZSBSZ04 0.0383456 0.0049441 7.756 8.78e-15 ***
DESTIN_SZSBSZ05 -0.7642864 0.0059839 -127.725 < 2e-16 ***
DESTIN_SZSBSZ06 -2.1547719 0.0223953 -96.215 < 2e-16 ***
DESTIN_SZSBSZ07 -2.0184111 0.0163921 -123.133 < 2e-16 ***
DESTIN_SZSBSZ08 1.0585725 0.0045531 232.495 < 2e-16 ***
DESTIN_SZSBSZ09 0.5416971 0.0044614 121.420 < 2e-16 ***
DESTIN_SZSESZ02 -0.5331643 0.0040259 -132.434 < 2e-16 ***
DESTIN_SZSESZ03 0.3651744 0.0032038 113.983 < 2e-16 ***
DESTIN_SZSESZ04 -0.8717163 0.0046738 -186.511 < 2e-16 ***
DESTIN_SZSESZ05 -0.2586572 0.0039102 -66.150 < 2e-16 ***
DESTIN_SZSESZ06 -0.8930919 0.0049280 -181.227 < 2e-16 ***
DESTIN_SZSESZ07 -3.2801732 0.0197048 -166.466 < 2e-16 ***
DESTIN_SZSGSZ01 -0.2307557 0.0049986 -46.164 < 2e-16 ***
DESTIN_SZSGSZ02 -0.2503996 0.0044074 -56.814 < 2e-16 ***
DESTIN_SZSGSZ03 -0.4642544 0.0041304 -112.398 < 2e-16 ***
DESTIN_SZSGSZ04 -0.3866295 0.0041114 -94.038 < 2e-16 ***
DESTIN_SZSGSZ05 -2.0698393 0.0083196 -248.790 < 2e-16 ***
DESTIN_SZSGSZ06 0.3958431 0.0032265 122.683 < 2e-16 ***
DESTIN_SZSGSZ07 -0.4343179 0.0041796 -103.914 < 2e-16 ***
DESTIN_SZSISZ01 -0.5844619 0.0139428 -41.919 < 2e-16 ***
DESTIN_SZSKSZ01 -0.6708219 0.0062908 -106.635 < 2e-16 ***
DESTIN_SZSKSZ02 0.2019226 0.0046474 43.449 < 2e-16 ***
DESTIN_SZSKSZ03 -0.5039124 0.0052688 -95.640 < 2e-16 ***
DESTIN_SZSKSZ04 -0.9159282 0.0137135 -66.790 < 2e-16 ***
DESTIN_SZSKSZ05 -0.4787141 0.0115189 -41.559 < 2e-16 ***
DESTIN_SZSLSZ01 -0.7672697 0.0066932 -114.635 < 2e-16 ***
DESTIN_SZSLSZ04 -0.9177537 0.0054430 -168.611 < 2e-16 ***
DESTIN_SZSRSZ01 -1.1818764 0.0065352 -180.849 < 2e-16 ***
DESTIN_SZSRSZ02 -1.3030182 0.0077031 -169.155 < 2e-16 ***
DESTIN_SZSRSZ03 -1.4193324 0.0070032 -202.668 < 2e-16 ***
DESTIN_SZSVSZ01 -1.0638836 0.0421548 -25.238 < 2e-16 ***
DESTIN_SZTHSZ01 -3.7186433 0.0338624 -109.816 < 2e-16 ***
DESTIN_SZTHSZ03 -2.4807672 0.0178452 -139.016 < 2e-16 ***
DESTIN_SZTHSZ04 -2.4477557 0.0177022 -138.274 < 2e-16 ***
DESTIN_SZTHSZ06 -1.9277088 0.0119892 -160.788 < 2e-16 ***
DESTIN_SZTMSZ01 -0.3957891 0.0046099 -85.856 < 2e-16 ***
DESTIN_SZTMSZ02 1.2843697 0.0029280 438.650 < 2e-16 ***
DESTIN_SZTMSZ03 0.5255595 0.0032808 160.192 < 2e-16 ***
DESTIN_SZTMSZ04 0.7492490 0.0034407 217.759 < 2e-16 ***
DESTIN_SZTMSZ05 0.5931173 0.0042822 138.507 < 2e-16 ***
DESTIN_SZTNSZ01 -0.4535588 0.0046495 -97.549 < 2e-16 ***
DESTIN_SZTNSZ02 -1.0934837 0.0060338 -181.225 < 2e-16 ***
DESTIN_SZTNSZ03 -1.1159434 0.0074351 -150.092 < 2e-16 ***
DESTIN_SZTNSZ04 -0.9995060 0.0057250 -174.588 < 2e-16 ***
DESTIN_SZTPSZ01 -0.3696675 0.0046378 -79.707 < 2e-16 ***
DESTIN_SZTPSZ02 0.2525639 0.0032305 78.181 < 2e-16 ***
DESTIN_SZTPSZ03 -0.2528258 0.0047075 -53.707 < 2e-16 ***
DESTIN_SZTPSZ04 -1.5279995 0.0062841 -243.154 < 2e-16 ***
DESTIN_SZTPSZ05 -0.9404446 0.0049204 -191.131 < 2e-16 ***
DESTIN_SZTPSZ06 -0.4598113 0.0063548 -72.357 < 2e-16 ***
DESTIN_SZTPSZ07 -1.7673404 0.0090777 -194.690 < 2e-16 ***
DESTIN_SZTPSZ08 -1.2383224 0.0067084 -184.593 < 2e-16 ***
DESTIN_SZTPSZ09 -0.3746430 0.0050668 -73.941 < 2e-16 ***
DESTIN_SZTPSZ10 -0.9533305 0.0062146 -153.401 < 2e-16 ***
DESTIN_SZTPSZ11 -0.3056595 0.0042168 -72.487 < 2e-16 ***
DESTIN_SZTPSZ12 -0.6867893 0.0052094 -131.836 < 2e-16 ***
DESTIN_SZTSSZ01 -0.8762131 0.0198461 -44.150 < 2e-16 ***
DESTIN_SZTSSZ02 -0.6843938 0.0094213 -72.643 < 2e-16 ***
DESTIN_SZTSSZ03 -0.0904515 0.0070799 -12.776 < 2e-16 ***
DESTIN_SZTSSZ04 0.3743395 0.0075914 49.311 < 2e-16 ***
DESTIN_SZTSSZ05 1.3853172 0.0075725 182.941 < 2e-16 ***
DESTIN_SZTSSZ06 1.5296750 0.0134882 113.409 < 2e-16 ***
DESTIN_SZWCSZ01 0.9091963 0.0047097 193.047 < 2e-16 ***
DESTIN_SZWCSZ02 -0.6439147 0.0081845 -78.675 < 2e-16 ***
DESTIN_SZWCSZ03 -2.1459650 0.0247981 -86.538 < 2e-16 ***
DESTIN_SZWDSZ01 0.9549661 0.0031046 307.594 < 2e-16 ***
DESTIN_SZWDSZ02 -0.7933981 0.0050239 -157.924 < 2e-16 ***
DESTIN_SZWDSZ03 0.5997619 0.0035370 169.568 < 2e-16 ***
DESTIN_SZWDSZ04 -0.4695093 0.0050382 -93.190 < 2e-16 ***
DESTIN_SZWDSZ05 -0.1335226 0.0045530 -29.327 < 2e-16 ***
DESTIN_SZWDSZ06 0.1507121 0.0036921 40.820 < 2e-16 ***
DESTIN_SZWDSZ07 0.0983013 0.0049505 19.857 < 2e-16 ***
DESTIN_SZWDSZ08 0.0429259 0.0054738 7.842 4.43e-15 ***
DESTIN_SZWDSZ09 -0.2510980 0.0041672 -60.256 < 2e-16 ***
DESTIN_SZYSSZ01 0.9397554 0.0033699 278.866 < 2e-16 ***
DESTIN_SZYSSZ02 -0.3479625 0.0045725 -76.100 < 2e-16 ***
DESTIN_SZYSSZ03 -1.2752017 0.0047773 -266.929 < 2e-16 ***
DESTIN_SZYSSZ04 -0.4025264 0.0043523 -92.486 < 2e-16 ***
DESTIN_SZYSSZ05 -1.9068813 0.0087608 -217.660 < 2e-16 ***
DESTIN_SZYSSZ06 -1.3395696 0.0066157 -202.482 < 2e-16 ***
DESTIN_SZYSSZ07 -0.6224806 0.0083748 -74.328 < 2e-16 ***
DESTIN_SZYSSZ08 0.6438403 0.0033832 190.307 < 2e-16 ***
DESTIN_SZYSSZ09 -0.0171702 0.0035489 -4.838 1.31e-06 ***
log(dist) -0.6926466 0.0001089 -6358.780 < 2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
(Dispersion parameter for poisson family taken to be 1)
Null deviance: 102725286 on 20848 degrees of freedom
Residual deviance: 34895044 on 20230 degrees of freedom
AIC: 35027862
Number of Fisher Scoring iterations: 7Again, let’s check the R-square values of doubly constrained SIM model this time.
CalcRSquared(dbcSIM$data$TRIPS, dbcSIM$fitted.values)[1] 0.5676511Notice that there is a relatively larger improvement in the R-squared value than destination constrained SIM.
8 Model Comparison
Another useful model performance measure for continuous dependent variable is Root Mean Squared Error. In this sub-section, you will learn how to use compare_performance() of performance package.
First of all, let us create a list called model_list by using the code chunk below.
model_list <- list(unconstrained=uncSIM,
originConstrained=orcSIM,
destinationConstrained=decSIM,
doublyConstrained=dbcSIM)Next, we will compute the RMSE of all the models in model_list file by using the code chunk below.
compare_performance(model_list,
metrics = "RMSE")# Comparison of Model Performance Indices
Name | Model | RMSE
-----------------------------------------
unconstrained | glm | 5371.511
originConstrained | glm | 4735.612
destinationConstrained | glm | 4290.835
doublyConstrained | glm | 4114.493The print above reveals that doubly constrained SIM is the best model among all the four SIMs because it has the smallest RMSE value of 4114.493
9 Visualising fitted values
Firstly we will extract the fitted values from each model by using the code chunk below.
df_unc <- as.data.frame(uncSIM$fitted.values) %>%
round(digits = 0)Next, we will join the values to SIM_data data frame.
SIM_data <- SIM_data %>%
cbind(df_unc) %>%
rename(uncTRIPS = "uncSIM$fitted.values")We repeat the same step by for Origin Constrained SIM (i.e. orcSIM)
df_orc <- as.data.frame(orcSIM$fitted.values) %>%
round(digits = 0)SIM_data <- SIM_data %>%
cbind(df_orc) %>%
rename(orcTRIPS = "orcSIM$fitted.values")We repeat the same step by for Destination Constrained SIM (i.e. decSIM)
df_dec <- as.data.frame(decSIM$fitted.values) %>%
round(digits = 0)SIM_data <- SIM_data %>%
cbind(df_dec) %>%
rename(decTRIPS = "decSIM$fitted.values")We repeat the same step by for Doubly Constrained SIM (i.e. dbcSIM)
df_dbc <- as.data.frame(dbcSIM$fitted.values) %>%
round(digits = 0)SIM_data <- SIM_data %>%
cbind(df_dbc) %>%
rename(dbcTRIPS = "dbcSIM$fitted.values")Next, we will create plots for each model.
Show the code
unc_p <- ggplot(data = SIM_data,
aes(x = uncTRIPS,
y = TRIPS)) +
geom_point() +
geom_smooth(method = lm) +
theme(
plot.title = element_text(face= 'bold', size = 10),
panel.grid.major = element_line(colour = "#ede5de", linetype = 1, linewidth = 0.5),
panel.grid.minor = element_line(colour = "#ede5de", linetype = 1, linewidth= 0.5),
plot.background = element_rect(fill="#E4D5C9",colour="#E4D5C9"),
panel.border = element_blank(),
panel.background = element_blank(),
axis.ticks = element_blank(),
axis.title.y = element_text(hjust=1, angle=0, size = 8),
axis.title.x = element_text(size = 8),
strip.text = element_text(face= 'bold'),
strip.background = element_rect(color="#E4D5C9", fill="#E4D5C9")
)
orc_p <- ggplot(data = SIM_data,
aes(x = orcTRIPS,
y = TRIPS)) +
geom_point() +
geom_smooth(method = lm) +
theme(
plot.title = element_text(face= 'bold', size = 10),
panel.grid.major = element_line(colour = "#ede5de", linetype = 1, linewidth = 0.5),
panel.grid.minor = element_line(colour = "#ede5de", linetype = 1, linewidth= 0.5),
plot.background = element_rect(fill="#E4D5C9",colour="#E4D5C9"),
panel.border = element_blank(),
panel.background = element_blank(),
axis.ticks = element_blank(),
axis.title.y = element_text(hjust=1, angle=0, size = 8),
axis.title.x = element_text(size = 8),
strip.text = element_text(face= 'bold'),
strip.background = element_rect(color="#E4D5C9", fill="#E4D5C9")
)
dec_p <- ggplot(data = SIM_data,
aes(x = decTRIPS,
y = TRIPS)) +
geom_point() +
geom_smooth(method = lm) +
theme(
plot.title = element_text(face= 'bold', size = 10),
panel.grid.major = element_line(colour = "#ede5de", linetype = 1, linewidth = 0.5),
panel.grid.minor = element_line(colour = "#ede5de", linetype = 1, linewidth= 0.5),
plot.background = element_rect(fill="#E4D5C9",colour="#E4D5C9"),
panel.border = element_blank(),
panel.background = element_blank(),
axis.ticks = element_blank(),
axis.title.y = element_text(hjust=1, angle=0, size = 8),
axis.title.x = element_text(size = 8),
strip.text = element_text(face= 'bold'),
strip.background = element_rect(color="#E4D5C9", fill="#E4D5C9")
)
dbc_p <- ggplot(data = SIM_data,
aes(x = dbcTRIPS,
y = TRIPS)) +
geom_point() +
geom_smooth(method = lm) +
theme(
plot.title = element_text(face= 'bold', size = 10),
panel.grid.major = element_line(colour = "#ede5de", linetype = 1, linewidth = 0.5),
panel.grid.minor = element_line(colour = "#ede5de", linetype = 1, linewidth= 0.5),
plot.background = element_rect(fill="#E4D5C9",colour="#E4D5C9"),
panel.border = element_blank(),
panel.background = element_blank(),
axis.ticks = element_blank(),
axis.title.y = element_text(hjust=1, angle=0, size = 8),
axis.title.x = element_text(size = 8),
strip.text = element_text(face= 'bold'),
strip.background = element_rect(color="#E4D5C9", fill="#E4D5C9")
)Now, we will put all the graphs into a single visual for better comparison by using the code chunk below.
ggarrange(unc_p, orc_p, dec_p, dbc_p,
ncol = 2,
nrow = 2)
10 Reference
T. S., Kam. (2023). Calibrating Spatial Interaction Models with R. R4GDSA. https://r4gdsa.netlify.app/chap16