Qingyang Li | Qingyang.Li@Childmind.org
By far, we have generated 6 Smoothness estimation reports, That are:
Smoothness Estimation, part 5: Part 5 and part 4 are all the same except that when estimate the smoothness, -detrend option was on for part 5. The -detrend option did not change anything for methed 1). but for method 2), the smoothness values are smaller than the targets.
Smoothness Estimation, part 6: Part 6 and part 5 are all the same except that when estimate the smoothness for each subject, instead of using 3dFWHMx on single files, I concatenate the drivetive files into 4D as input and use -out flag to output their FWHM.
The target smoothness for each derivative is listed below
Note: the code below is just for plot and sumarroize the estimated smoothness. The code which I used to estimate the smoothness is not here (but here)
df_Smoothed_sub = read.csv("/home2/data/Projects/NKI_TRT_paper/NewAnalysisWithCPAC/Smoothestimation/SmoothnessofDerivatives_Smoothed_estimated_byAFNI_detrend_sub.csv")
# here we only show the data for one of the strategy
df_Smoothed_sub = subset(df_Smoothed_sub, df_Smoothed_sub$pipeline == "pipeline_NewStrawn" &
(df_Smoothed_sub$scan == "mx_1400_First5Min" | df_Smoothed_sub$scan == "mx_1400_First5MinDS" |
df_Smoothed_sub$scan == "mx_1400_First5MinDSNF" | df_Smoothed_sub$scan ==
"mx_645_First5Min" | df_Smoothed_sub$scan == "mx_645_First5MinDS" |
df_Smoothed_sub$scan == "mx_645_First5MinDSNF" | df_Smoothed_sub$scan ==
"std_2500"))
# df_Smoothed_grp =
# read.csv('/home2/data/Projects/NKI_TRT_paper/NewAnalysisWithCPAC/Smoothestimation/SmoothnessofDerivatives_Smoothed_estimated_byAFNI_detrend_grp.csv')
# here we only show the data for one of the strategy df_Smoothed_grp =
# subset(df_Smoothed_grp, df_Smoothed_grp$pipeline == 'pipeline_NewStrawn'
# & (df_Smoothed_grp$scan == 'mx_1400_First5Min' | df_Smoothed_grp$scan ==
# 'mx_1400_First5MinDS' | df_Smoothed_grp$scan == 'mx_1400_First5MinDSNF'
# |df_Smoothed_grp$scan == 'mx_645_First5Min' | df_Smoothed_grp$scan ==
# 'mx_645_First5MinDS' | df_Smoothed_grp$scan == 'mx_645_First5MinDSNF'
# |df_Smoothed_grp$scan == 'std_2500') )
# df_Smoothed_grp$session = as.factor(df_Smoothed_grp$session)
df_Smoothed_sub$session = as.factor(df_Smoothed_sub$session)
df_Smoothed_sub$subid = as.factor(df_Smoothed_sub$subid)
summary(df_Smoothed_sub)
## pipeline scan
## pipeline_NewStrawn :9548 mx_1400_First5Min :1364
## pipeline_WoodworthTown: 0 mx_1400_First5MinDS :1364
## mx_1400_First5MinDSNF:1364
## mx_645_First5Min :1364
## mx_645_First5MinDS :1364
## mx_645_First5MinDSNF :1364
## (Other) :1364
## derivative session subid
## alff : 308 1:4774 21002 : 434
## degree_centrality_binarize : 308 2:4774 21006 : 434
## degree_centrality_weighted : 308 21018 : 434
## falff : 308 21024 : 434
## reho : 308 1427581: 434
## sca_roi_PNAS_smith10_plus5_1: 308 1793622: 434
## (Other) :7700 (Other):6944
## FWHMx FWHMy FWHMz
## Min. : 0.00 Min. : 3.82 Min. : 0.00
## 1st Qu.: 6.03 1st Qu.: 5.85 1st Qu.: 5.85
## Median : 6.48 Median : 6.28 Median : 6.37
## Mean : 6.86 Mean : 6.53 Mean : 6.62
## 3rd Qu.: 7.03 3rd Qu.: 6.83 3rd Qu.: 6.97
## Max. :18.37 Max. :13.96 Max. :98.97
##
# summary(df_Smoothed_grp)
for (i in derivs) {
print(plotDeriv_sub(df_Smoothed_sub, i, "x", outSummary = "on"))
# print(plotDeriv_grp(df_Smoothed_grp,i, 'x', outSummary = 'on'))
print(plotDeriv_sub(df_Smoothed_sub, i, "y"))
# print(plotDeriv_grp(df_Smoothed_grp,i, 'y'))
print(plotDeriv_sub(df_Smoothed_sub, i, "z"))
# print(plotDeriv_grp(df_Smoothed_grp,i, 'z'))
}
## pipeline scan
## pipeline_NewStrawn :308 mx_1400_First5Min :44
## pipeline_WoodworthTown: 0 mx_1400_First5MinDS :44
## mx_1400_First5MinDSNF:44
## mx_645_First5Min :44
## mx_645_First5MinDS :44
## mx_645_First5MinDSNF :44
## (Other) :44
## derivative session subid FWHMx
## alff :308 1:154 21002 : 14 Min. :4.49
## degree_centrality_binarize : 0 2:154 21006 : 14 1st Qu.:5.80
## degree_centrality_weighted : 0 21018 : 14 Median :6.23
## falff : 0 21024 : 14 Mean :6.25
## reho : 0 1427581: 14 3rd Qu.:6.72
## sca_roi_PNAS_smith10_plus5_1: 0 1793622: 14 Max. :8.28
## (Other) : 0 (Other):224
## FWHMy FWHMz
## Min. :3.82 Min. :4.26
## 1st Qu.:5.03 1st Qu.:5.58
## Median :5.48 Median :5.97
## Mean :5.47 Mean :5.91
## 3rd Qu.:5.94 3rd Qu.:6.29
## Max. :6.76 Max. :6.96
##
## pipeline scan
## pipeline_NewStrawn :308 mx_1400_First5Min :44
## pipeline_WoodworthTown: 0 mx_1400_First5MinDS :44
## mx_1400_First5MinDSNF:44
## mx_645_First5Min :44
## mx_645_First5MinDS :44
## mx_645_First5MinDSNF :44
## (Other) :44
## derivative session subid FWHMx
## degree_centrality_binarize :308 1:154 21002 : 14 Min. : 7.37
## alff : 0 2:154 21006 : 14 1st Qu.:10.05
## degree_centrality_weighted : 0 21018 : 14 Median :10.52
## falff : 0 21024 : 14 Mean :10.68
## reho : 0 1427581: 14 3rd Qu.:11.28
## sca_roi_PNAS_smith10_plus5_1: 0 1793622: 14 Max. :15.84
## (Other) : 0 (Other):224
## FWHMy FWHMz
## Min. : 7.20 Min. : 7.07
## 1st Qu.: 8.80 1st Qu.: 8.96
## Median : 9.19 Median : 9.63
## Mean : 9.21 Mean : 9.77
## 3rd Qu.: 9.62 3rd Qu.:10.47
## Max. :11.27 Max. :24.79
##
## pipeline scan
## pipeline_NewStrawn :308 mx_1400_First5Min :44
## pipeline_WoodworthTown: 0 mx_1400_First5MinDS :44
## mx_1400_First5MinDSNF:44
## mx_645_First5Min :44
## mx_645_First5MinDS :44
## mx_645_First5MinDSNF :44
## (Other) :44
## derivative session subid FWHMx
## degree_centrality_weighted :308 1:154 21002 : 14 Min. : 7.22
## alff : 0 2:154 21006 : 14 1st Qu.: 9.96
## degree_centrality_binarize : 0 21018 : 14 Median :10.54
## falff : 0 21024 : 14 Mean :10.65
## reho : 0 1427581: 14 3rd Qu.:11.15
## sca_roi_PNAS_smith10_plus5_1: 0 1793622: 14 Max. :16.79
## (Other) : 0 (Other):224
## FWHMy FWHMz
## Min. : 7.11 Min. : 6.90
## 1st Qu.: 8.80 1st Qu.: 8.83
## Median : 9.19 Median : 9.48
## Mean : 9.21 Mean : 9.54
## 3rd Qu.: 9.63 3rd Qu.:10.15
## Max. :10.91 Max. :13.34
##
## pipeline scan
## pipeline_NewStrawn :308 mx_1400_First5Min :44
## pipeline_WoodworthTown: 0 mx_1400_First5MinDS :44
## mx_1400_First5MinDSNF:44
## mx_645_First5Min :44
## mx_645_First5MinDS :44
## mx_645_First5MinDSNF :44
## (Other) :44
## derivative session subid FWHMx
## falff :308 1:154 21002 : 14 Min. :5.00
## alff : 0 2:154 21006 : 14 1st Qu.:5.47
## degree_centrality_binarize : 0 21018 : 14 Median :5.74
## degree_centrality_weighted : 0 21024 : 14 Mean :5.82
## reho : 0 1427581: 14 3rd Qu.:6.08
## sca_roi_PNAS_smith10_plus5_1: 0 1793622: 14 Max. :8.59
## (Other) : 0 (Other):224
## FWHMy FWHMz
## Min. :4.44 Min. :4.66
## 1st Qu.:5.15 1st Qu.:5.22
## Median :5.40 Median :5.48
## Mean :5.47 Mean :5.56
## 3rd Qu.:5.70 3rd Qu.:5.79
## Max. :8.51 Max. :9.22
##
## pipeline scan
## pipeline_NewStrawn :308 mx_1400_First5Min :44
## pipeline_WoodworthTown: 0 mx_1400_First5MinDS :44
## mx_1400_First5MinDSNF:44
## mx_645_First5Min :44
## mx_645_First5MinDS :44
## mx_645_First5MinDSNF :44
## (Other) :44
## derivative session subid FWHMx
## reho :308 1:154 21002 : 14 Min. : 7.84
## alff : 0 2:154 21006 : 14 1st Qu.: 9.81
## degree_centrality_binarize : 0 21018 : 14 Median :10.62
## degree_centrality_weighted : 0 21024 : 14 Mean :10.74
## falff : 0 1427581: 14 3rd Qu.:11.47
## sca_roi_PNAS_smith10_plus5_1: 0 1793622: 14 Max. :18.37
## (Other) : 0 (Other):224
## FWHMy FWHMz
## Min. : 6.58 Min. : 6.32
## 1st Qu.: 8.21 1st Qu.: 7.48
## Median : 8.84 Median : 7.90
## Mean : 8.90 Mean : 8.22
## 3rd Qu.: 9.49 3rd Qu.: 8.57
## Max. :12.80 Max. :14.20
##
## pipeline scan
## pipeline_NewStrawn :308 mx_1400_First5Min :44
## pipeline_WoodworthTown: 0 mx_1400_First5MinDS :44
## mx_1400_First5MinDSNF:44
## mx_645_First5Min :44
## mx_645_First5MinDS :44
## mx_645_First5MinDSNF :44
## (Other) :44
## derivative session subid FWHMx
## vmhc :308 1:154 21002 : 14 Min. : 0.00
## alff : 0 2:154 21006 : 14 1st Qu.: 6.96
## degree_centrality_binarize: 0 21018 : 14 Median : 7.26
## degree_centrality_weighted: 0 21024 : 14 Mean : 7.39
## falff : 0 1427581: 14 3rd Qu.: 7.63
## reho : 0 1793622: 14 Max. :16.03
## (Other) : 0 (Other):224
## FWHMy FWHMz
## Min. :5.71 Min. :5.80
## 1st Qu.:6.43 1st Qu.:6.42
## Median :6.63 Median :6.68
## Mean :6.64 Mean :6.69
## 3rd Qu.:6.86 3rd Qu.:6.93
## Max. :7.76 Max. :8.06
##
## pipeline scan
## pipeline_NewStrawn :308 mx_1400_First5Min :44
## pipeline_WoodworthTown: 0 mx_1400_First5MinDS :44
## mx_1400_First5MinDSNF:44
## mx_645_First5Min :44
## mx_645_First5MinDS :44
## mx_645_First5MinDSNF :44
## (Other) :44
## derivative session
## spatial_regression_spatial_map_PNAS_Smith09_rsn10_0000:308 1:154
## alff : 0 2:154
## degree_centrality_binarize : 0
## degree_centrality_weighted : 0
## falff : 0
## reho : 0
## (Other) : 0
## subid FWHMx FWHMy FWHMz
## 21002 : 14 Min. :5.03 Min. :4.87 Min. : 4.83
## 21006 : 14 1st Qu.:6.06 1st Qu.:5.93 1st Qu.: 5.99
## 21018 : 14 Median :6.29 Median :6.18 Median : 6.26
## 21024 : 14 Mean :6.33 Mean :6.24 Mean : 6.32
## 1427581: 14 3rd Qu.:6.62 3rd Qu.:6.49 3rd Qu.: 6.59
## 1793622: 14 Max. :8.03 Max. :8.84 Max. :11.16
## (Other):224
## pipeline scan
## pipeline_NewStrawn :308 mx_1400_First5Min :44
## pipeline_WoodworthTown: 0 mx_1400_First5MinDS :44
## mx_1400_First5MinDSNF:44
## mx_645_First5Min :44
## mx_645_First5MinDS :44
## mx_645_First5MinDSNF :44
## (Other) :44
## derivative session
## spatial_regression_spatial_map_PNAS_Smith09_rsn10_0001:308 1:154
## alff : 0 2:154
## degree_centrality_binarize : 0
## degree_centrality_weighted : 0
## falff : 0
## reho : 0
## (Other) : 0
## subid FWHMx FWHMy FWHMz
## 21002 : 14 Min. :5.49 Min. : 6.04 Min. : 0.00
## 21006 : 14 1st Qu.:6.42 1st Qu.: 6.69 1st Qu.: 6.62
## 21018 : 14 Median :6.80 Median : 6.98 Median : 7.04
## 21024 : 14 Mean :6.88 Mean : 7.12 Mean : 7.26
## 1427581: 14 3rd Qu.:7.17 3rd Qu.: 7.36 3rd Qu.: 7.52
## 1793622: 14 Max. :9.78 Max. :11.01 Max. :48.62
## (Other):224
## pipeline scan
## pipeline_NewStrawn :308 mx_1400_First5Min :44
## pipeline_WoodworthTown: 0 mx_1400_First5MinDS :44
## mx_1400_First5MinDSNF:44
## mx_645_First5Min :44
## mx_645_First5MinDS :44
## mx_645_First5MinDSNF :44
## (Other) :44
## derivative session
## spatial_regression_spatial_map_PNAS_Smith09_rsn10_0002:308 1:154
## alff : 0 2:154
## degree_centrality_binarize : 0
## degree_centrality_weighted : 0
## falff : 0
## reho : 0
## (Other) : 0
## subid FWHMx FWHMy FWHMz
## 21002 : 14 Min. : 5.72 Min. :5.34 Min. : 0.00
## 21006 : 14 1st Qu.: 6.34 1st Qu.:6.33 1st Qu.: 6.26
## 21018 : 14 Median : 6.64 Median :6.62 Median : 6.59
## 21024 : 14 Mean : 6.75 Mean :6.70 Mean : 6.67
## 1427581: 14 3rd Qu.: 7.00 3rd Qu.:6.93 3rd Qu.: 6.94
## 1793622: 14 Max. :10.30 Max. :9.77 Max. :18.94
## (Other):224
## pipeline scan
## pipeline_NewStrawn :308 mx_1400_First5Min :44
## pipeline_WoodworthTown: 0 mx_1400_First5MinDS :44
## mx_1400_First5MinDSNF:44
## mx_645_First5Min :44
## mx_645_First5MinDS :44
## mx_645_First5MinDSNF :44
## (Other) :44
## derivative session
## spatial_regression_spatial_map_PNAS_Smith09_rsn10_0003:308 1:154
## alff : 0 2:154
## degree_centrality_binarize : 0
## degree_centrality_weighted : 0
## falff : 0
## reho : 0
## (Other) : 0
## subid FWHMx FWHMy FWHMz
## 21002 : 14 Min. :5.84 Min. :5.44 Min. : 5.39
## 21006 : 14 1st Qu.:6.39 1st Qu.:6.05 1st Qu.: 6.13
## 21018 : 14 Median :6.62 Median :6.30 Median : 6.39
## 21024 : 14 Mean :6.67 Mean :6.34 Mean : 6.44
## 1427581: 14 3rd Qu.:6.88 3rd Qu.:6.55 3rd Qu.: 6.62
## 1793622: 14 Max. :8.29 Max. :8.58 Max. :14.77
## (Other):224
## pipeline scan
## pipeline_NewStrawn :308 mx_1400_First5Min :44
## pipeline_WoodworthTown: 0 mx_1400_First5MinDS :44
## mx_1400_First5MinDSNF:44
## mx_645_First5Min :44
## mx_645_First5MinDS :44
## mx_645_First5MinDSNF :44
## (Other) :44
## derivative session
## spatial_regression_spatial_map_PNAS_Smith09_rsn10_0004:308 1:154
## alff : 0 2:154
## degree_centrality_binarize : 0
## degree_centrality_weighted : 0
## falff : 0
## reho : 0
## (Other) : 0
## subid FWHMx FWHMy FWHMz
## 21002 : 14 Min. : 5.17 Min. : 5.77 Min. : 0.00
## 21006 : 14 1st Qu.: 6.58 1st Qu.: 6.64 1st Qu.: 6.79
## 21018 : 14 Median : 6.94 Median : 7.04 Median : 7.18
## 21024 : 14 Mean : 7.05 Mean : 7.19 Mean : 7.21
## 1427581: 14 3rd Qu.: 7.39 3rd Qu.: 7.48 3rd Qu.: 7.65
## 1793622: 14 Max. :10.86 Max. :13.96 Max. :18.12
## (Other):224
## pipeline scan
## pipeline_NewStrawn :308 mx_1400_First5Min :44
## pipeline_WoodworthTown: 0 mx_1400_First5MinDS :44
## mx_1400_First5MinDSNF:44
## mx_645_First5Min :44
## mx_645_First5MinDS :44
## mx_645_First5MinDSNF :44
## (Other) :44
## derivative session
## spatial_regression_spatial_map_PNAS_Smith09_rsn10_0005:308 1:154
## alff : 0 2:154
## degree_centrality_binarize : 0
## degree_centrality_weighted : 0
## falff : 0
## reho : 0
## (Other) : 0
## subid FWHMx FWHMy FWHMz
## 21002 : 14 Min. :5.40 Min. : 5.41 Min. : 0.00
## 21006 : 14 1st Qu.:6.34 1st Qu.: 6.27 1st Qu.: 6.40
## 21018 : 14 Median :6.58 Median : 6.59 Median : 6.66
## 21024 : 14 Mean :6.64 Mean : 6.64 Mean : 6.78
## 1427581: 14 3rd Qu.:6.86 3rd Qu.: 6.87 3rd Qu.: 7.00
## 1793622: 14 Max. :8.98 Max. :10.09 Max. :27.31
## (Other):224
## pipeline scan
## pipeline_NewStrawn :308 mx_1400_First5Min :44
## pipeline_WoodworthTown: 0 mx_1400_First5MinDS :44
## mx_1400_First5MinDSNF:44
## mx_645_First5Min :44
## mx_645_First5MinDS :44
## mx_645_First5MinDSNF :44
## (Other) :44
## derivative session
## spatial_regression_spatial_map_PNAS_Smith09_rsn10_0006:308 1:154
## alff : 0 2:154
## degree_centrality_binarize : 0
## degree_centrality_weighted : 0
## falff : 0
## reho : 0
## (Other) : 0
## subid FWHMx FWHMy FWHMz
## 21002 : 14 Min. :5.78 Min. :5.58 Min. : 0.00
## 21006 : 14 1st Qu.:6.53 1st Qu.:6.33 1st Qu.: 6.61
## 21018 : 14 Median :6.80 Median :6.59 Median : 6.92
## 21024 : 14 Mean :6.85 Mean :6.64 Mean : 7.04
## 1427581: 14 3rd Qu.:7.12 3rd Qu.:6.91 3rd Qu.: 7.28
## 1793622: 14 Max. :8.99 Max. :8.92 Max. :17.57
## (Other):224
## pipeline scan
## pipeline_NewStrawn :308 mx_1400_First5Min :44
## pipeline_WoodworthTown: 0 mx_1400_First5MinDS :44
## mx_1400_First5MinDSNF:44
## mx_645_First5Min :44
## mx_645_First5MinDS :44
## mx_645_First5MinDSNF :44
## (Other) :44
## derivative session
## spatial_regression_spatial_map_PNAS_Smith09_rsn10_0007:308 1:154
## alff : 0 2:154
## degree_centrality_binarize : 0
## degree_centrality_weighted : 0
## falff : 0
## reho : 0
## (Other) : 0
## subid FWHMx FWHMy FWHMz
## 21002 : 14 Min. : 5.67 Min. : 5.79 Min. : 0.00
## 21006 : 14 1st Qu.: 6.62 1st Qu.: 6.60 1st Qu.: 6.72
## 21018 : 14 Median : 6.96 Median : 6.84 Median : 7.02
## 21024 : 14 Mean : 7.03 Mean : 6.95 Mean : 7.09
## 1427581: 14 3rd Qu.: 7.34 3rd Qu.: 7.18 3rd Qu.: 7.34
## 1793622: 14 Max. :10.23 Max. :10.64 Max. :13.55
## (Other):224
## pipeline scan
## pipeline_NewStrawn :308 mx_1400_First5Min :44
## pipeline_WoodworthTown: 0 mx_1400_First5MinDS :44
## mx_1400_First5MinDSNF:44
## mx_645_First5Min :44
## mx_645_First5MinDS :44
## mx_645_First5MinDSNF :44
## (Other) :44
## derivative session
## spatial_regression_spatial_map_PNAS_Smith09_rsn10_0008:308 1:154
## alff : 0 2:154
## degree_centrality_binarize : 0
## degree_centrality_weighted : 0
## falff : 0
## reho : 0
## (Other) : 0
## subid FWHMx FWHMy FWHMz
## 21002 : 14 Min. :6.09 Min. :5.91 Min. : 0.00
## 21006 : 14 1st Qu.:6.63 1st Qu.:6.39 1st Qu.: 6.51
## 21018 : 14 Median :6.94 Median :6.70 Median : 6.81
## 21024 : 14 Mean :6.98 Mean :6.74 Mean : 6.92
## 1427581: 14 3rd Qu.:7.22 3rd Qu.:6.97 3rd Qu.: 7.19
## 1793622: 14 Max. :9.00 Max. :8.50 Max. :23.52
## (Other):224
## pipeline scan
## pipeline_NewStrawn :308 mx_1400_First5Min :44
## pipeline_WoodworthTown: 0 mx_1400_First5MinDS :44
## mx_1400_First5MinDSNF:44
## mx_645_First5Min :44
## mx_645_First5MinDS :44
## mx_645_First5MinDSNF :44
## (Other) :44
## derivative session
## spatial_regression_spatial_map_PNAS_Smith09_rsn10_0009:308 1:154
## alff : 0 2:154
## degree_centrality_binarize : 0
## degree_centrality_weighted : 0
## falff : 0
## reho : 0
## (Other) : 0
## subid FWHMx FWHMy FWHMz
## 21002 : 14 Min. :5.78 Min. :5.88 Min. : 0.00
## 21006 : 14 1st Qu.:6.62 1st Qu.:6.37 1st Qu.: 6.52
## 21018 : 14 Median :6.84 Median :6.63 Median : 6.79
## 21024 : 14 Mean :6.91 Mean :6.71 Mean : 6.87
## 1427581: 14 3rd Qu.:7.14 3rd Qu.:6.96 3rd Qu.: 7.11
## 1793622: 14 Max. :8.65 Max. :8.90 Max. :28.02
## (Other):224
showing correlations of mean smoothness with the number of voxels with ICC > 0.5, for each statistics.
## pipeline scan
## pipeline_NewStrawn :224 mx_1400_First5Min :32
## pipeline_WoodworthTown: 0 mx_1400_First5MinDS :32
## mx_1400_First5MinDSNF:32
## mx_645_First5Min :32
## mx_645_First5MinDS :32
## mx_645_First5MinDSNF :32
## (Other) :32
## derivative session FWHMx FWHMy
## alff : 14 1:112 Min. :5.41 Min. :5.27
## degree_centrality_binarize: 14 2:112 1st Qu.:6.12 1st Qu.:5.97
## degree_centrality_weighted: 14 Median :6.32 Median :6.23
## falff : 14 Mean :6.65 Mean :6.48
## reho : 14 3rd Qu.:6.51 3rd Qu.:6.47
## spatial_regression_0000 : 14 Max. :9.56 Max. :8.99
## (Other) :140
## FWHMz nVox FWHMm
## Min. :5.21 Min. : 4313 Min. :5.29
## 1st Qu.:5.93 1st Qu.: 10458 1st Qu.:6.04
## Median :6.20 Median : 14062 Median :6.25
## Mean :6.44 Mean : 33233 Mean :6.52
## 3rd Qu.:6.46 3rd Qu.: 28317 3rd Qu.:6.46
## Max. :9.24 Max. :143992 Max. :9.26
##
