Source: https://github.com/markziemann/asd_meth/cohort_characterisation.Rmd

library("kableExtra")
library("magrittr")

Intro

Clinical characterisation of the cohort, shown in table 1 of the manuscript.

There are 3 samplesheets provided. Separate ones for blood and guthrie and also a combined one where the

19 ASD_blood_summarysheet.csv 40 ASD_guthrie_blood_combined_summarysheet.csv 23 ASD_guthrie_summarysheet.csv

datg <- read.table("ASD_guthrie_summarysheet.csv",sep=",",header=TRUE)

datb <- read.table("ASD_blood_summarysheet.csv",sep=",",header=TRUE)

dat <- datg

dat$Guthrie_sample <- TRUE

dat$Blood_sample <- datg$SampleName %in% datb$Sample_Name

dat %>%
  kbl(caption="Full sample sheet of the anonymous participant data") %>%
  kable_paper("hover", full_width = FALSE)
Full sample sheet of the anonymous participant data
SampleID ChipID Sentrix_ID Sentrix_Position SampleName Sex Age_months Age_years Birth_mode Birth_complications gestational_age Birth_weight.g. Birth_order Help_breathing_start Motor_skills Social_interaction_on_ADOS Diagnosis SRS_social_scores language IQ Family_ID Guthrie_sample Blood_sample
S3 202047860135_R03C01 202047860135 R03C01 ASD001.1 1 163 17 1 0 35.0 2112.000 1 1 2 3 0 103 53 62 1 TRUE TRUE
S4 202047860135_R04C01 202047860135 R04C01 ASD001.2 1 163 17 1 0 35.0 2680.000 2 1 2 12 2 121 61 68 1 TRUE TRUE
S31 202010290151_R07C01 202010290151 R07C01 ASD002.1 2 95 11 1 0 36.0 2240.000 1 0 2 9 1 84 113 98 2 TRUE TRUE
S32 202010290151_R08C01 202010290151 R08C01 ASD002.2 2 95 11 1 0 36.0 3395.000 2 0 0 13 2 47 113 105 2 TRUE TRUE
S11 202013790140_R03C01 202013790140 R03C01 ASD003.1 1 44 7 1 0 36.0 2220.000 2 1 1 15 2 74 102 110 3 TRUE TRUE
S12 202013790140_R04C01 202013790140 R04C01 ASD003.2 1 44 7 1 0 36.0 2020.000 1 1 1 12 1 94 108 103 3 TRUE TRUE
S13 202013790140_R05C01 202013790140 R05C01 ASD004.1 2 80 10 1 0 36.7 2400.000 1 0 0 1 0 35 99 94 4 TRUE TRUE
S14 202013790140_R06C01 202013790140 R06C01 ASD004.2 2 80 10 1 0 36.7 1900.000 2 0 0 8 1 61 66 67 4 TRUE TRUE
S17 202010290146_R01C01 202010290146 R01C01 ASD005.1 1 73 8 2 0 37.5 2296.000 1 0 0 16 2 95 50 0 5 TRUE FALSE
S18 202010290146_R02C01 202010290146 R02C01 ASD005.2 1 73 8 2 3 37.5 3120.000 2 1 0 5 0 79 68 102 5 TRUE FALSE
S21 202010290146_R05C01 202010290146 R05C01 ASD006.1 1 37 5 1 0 36.4 2400.000 1 0 0 16 2 125 50 53 6 TRUE TRUE
S22 202010290146_R06C01 202010290146 R06C01 ASD006.2 1 37 5 1 0 36.4 2.500 2 0 0 11 1 73 81 86 6 TRUE TRUE
S27 202010290151_R03C01 202010290151 R03C01 ASD007.1 2 131 14 0 3 34.0 1927.700 1 0 2 10 2 78 66 71 7 TRUE TRUE
S28 202010290151_R04C01 202010290151 R04C01 ASD007.2 2 10 14 3 3 34.0 1814.000 2 0 2 18 2 100 50 62 7 TRUE TRUE
S7 202047860135_R07C01 202047860135 R07C01 ASD008.1 2 69 8 1 0 37.0 2610.000 1 0 2 18 2 111 50 60 8 TRUE TRUE
S8 202047860135_R08C01 202047860135 R08C01 ASD008.2 2 69 8 1 0 37.0 2684.000 2 0 2 20 2 104 50 64 8 TRUE TRUE
S23 202010290146_R07C01 202010290146 R07C01 ASD009.1 2 41 6 2 1 31.0 1817.000 1 1 2 7 0 102 95 100 9 TRUE FALSE
S24 202010290146_R08C01 202010290146 R08C01 ASD009.2 2 41 6 2 1 31.0 2093.000 2 1 2 9 1 127 87 105 9 TRUE FALSE
S39 202010290176_R07C01 202010290176 R07C01 ASD012.1 2 136 13 1 0 38.0 2336.001 1 0 0 17 2 45 99 88 12 TRUE TRUE
S40 202010290176_R08C01 202010290176 R08C01 ASD012.2 2 136 13 1 0 38.0 2630.840 2 0 0 12 1 52 79 80 12 TRUE TRUE
S35 202010290176_R03C01 202010290176 R03C01 ASD013.1 2 12 13 2 2 31.0 1551.000 1 1 0 1 0 4 120 106 13 TRUE TRUE
S36 202010290176_R04C01 202010290176 R04C01 ASD013.2 2 146 13 2 2 31.0 1055.000 2 1 0 13 2 34 41 60 13 TRUE TRUE

Guthrie card arm summary stats

Provide some summary statistics for the Guthrie arm of the project.

message("summary stats")
## summary stats
summary(datg)
##    SampleID            ChipID            Sentrix_ID       Sentrix_Position  
##  Length:22          Length:22          Min.   :2.02e+11   Length:22         
##  Class :character   Class :character   1st Qu.:2.02e+11   Class :character  
##  Mode  :character   Mode  :character   Median :2.02e+11   Mode  :character  
##                                        Mean   :2.02e+11                     
##                                        3rd Qu.:2.02e+11                     
##                                        Max.   :2.02e+11                     
##   SampleName             Sex          Age_months       Age_years    
##  Length:22          Min.   :1.000   Min.   : 10.00   Min.   : 5.00  
##  Class :character   1st Qu.:1.000   1st Qu.: 41.75   1st Qu.: 7.25  
##  Mode  :character   Median :2.000   Median : 73.00   Median :10.00  
##                     Mean   :1.636   Mean   : 80.68   Mean   :10.18  
##                     3rd Qu.:2.000   3rd Qu.:122.00   3rd Qu.:13.00  
##                     Max.   :2.000   Max.   :163.00   Max.   :17.00  
##    Birth_mode    Birth_complications gestational_age Birth_weight.g. 
##  Min.   :0.000   Min.   :0.0000      Min.   :31.00   Min.   :   2.5  
##  1st Qu.:1.000   1st Qu.:0.0000      1st Qu.:34.25   1st Qu.:1906.9  
##  Median :1.000   Median :0.0000      Median :36.00   Median :2230.0  
##  Mean   :1.318   Mean   :0.6818      Mean   :35.33   Mean   :2150.2  
##  3rd Qu.:2.000   3rd Qu.:1.0000      3rd Qu.:36.92   3rd Qu.:2557.5  
##  Max.   :3.000   Max.   :3.0000      Max.   :38.00   Max.   :3395.0  
##   Birth_order  Help_breathing_start  Motor_skills    Social_interaction_on_ADOS
##  Min.   :1.0   Min.   :0.0000       Min.   :0.0000   Min.   : 1.00             
##  1st Qu.:1.0   1st Qu.:0.0000       1st Qu.:0.0000   1st Qu.: 8.25             
##  Median :1.5   Median :0.0000       Median :0.5000   Median :12.00             
##  Mean   :1.5   Mean   :0.4091       Mean   :0.9091   Mean   :11.18             
##  3rd Qu.:2.0   3rd Qu.:1.0000       3rd Qu.:2.0000   3rd Qu.:15.75             
##  Max.   :2.0   Max.   :1.0000       Max.   :2.0000   Max.   :20.00             
##    Diagnosis     SRS_social_scores    language            IQ        
##  Min.   :0.000   Min.   :  4.00    Min.   : 41.00   Min.   :  0.00  
##  1st Qu.:1.000   1st Qu.: 54.25    1st Qu.: 50.75   1st Qu.: 62.50  
##  Median :1.500   Median : 81.50    Median : 73.50   Median : 83.00  
##  Mean   :1.273   Mean   : 79.45    Mean   : 77.32   Mean   : 79.27  
##  3rd Qu.:2.000   3rd Qu.:102.75    3rd Qu.: 99.00   3rd Qu.:101.50  
##  Max.   :2.000   Max.   :127.00    Max.   :120.00   Max.   :110.00  
##    Family_ID     
##  Min.   : 1.000  
##  1st Qu.: 3.250  
##  Median : 6.000  
##  Mean   : 6.364  
##  3rd Qu.: 8.750  
##  Max.   :13.000
message("standard deviation")
## standard deviation
mysd <- apply(datg,2,sd)
## Warning in var(if (is.vector(x) || is.factor(x)) x else as.double(x), na.rm =
## na.rm): NAs introduced by coercion

## Warning in var(if (is.vector(x) || is.factor(x)) x else as.double(x), na.rm =
## na.rm): NAs introduced by coercion

## Warning in var(if (is.vector(x) || is.factor(x)) x else as.double(x), na.rm =
## na.rm): NAs introduced by coercion

## Warning in var(if (is.vector(x) || is.factor(x)) x else as.double(x), na.rm =
## na.rm): NAs introduced by coercion
format(mysd,scientific=FALSE)
##                   SampleID                     ChipID 
##         "              NA"         "              NA" 
##                 Sentrix_ID           Sentrix_Position 
##         "14586836.0205855"         "              NA" 
##                 SampleName                        Sex 
##         "              NA"         "       0.4923660" 
##                 Age_months                  Age_years 
##         "      47.1283823"         "       3.6727616" 
##                 Birth_mode        Birth_complications 
##         "       0.6463350"         "       1.1291111" 
##            gestational_age            Birth_weight.g. 
##         "       2.3486940"         "     697.6346866" 
##                Birth_order       Help_breathing_start 
##         "       0.5117663"         "       0.5032363" 
##               Motor_skills Social_interaction_on_ADOS 
##         "       0.9714540"         "       5.4566374" 
##                  Diagnosis          SRS_social_scores 
##         "       0.8270325"         "      32.9223144" 
##                   language                         IQ 
##         "      25.3853846"         "      25.8883718" 
##                  Family_ID 
##         "       3.8115795"
message("categorical counts")
## categorical counts
mytab <- apply(datg,2,table)
mytab
## $SampleID
## 
## S11 S12 S13 S14 S17 S18 S21 S22 S23 S24 S27 S28  S3 S31 S32 S35 S36 S39  S4 S40 
##   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1 
##  S7  S8 
##   1   1 
## 
## $ChipID
## 
## 202010290146_R01C01 202010290146_R02C01 202010290146_R05C01 202010290146_R06C01 
##                   1                   1                   1                   1 
## 202010290146_R07C01 202010290146_R08C01 202010290151_R03C01 202010290151_R04C01 
##                   1                   1                   1                   1 
## 202010290151_R07C01 202010290151_R08C01 202010290176_R03C01 202010290176_R04C01 
##                   1                   1                   1                   1 
## 202010290176_R07C01 202010290176_R08C01 202013790140_R03C01 202013790140_R04C01 
##                   1                   1                   1                   1 
## 202013790140_R05C01 202013790140_R06C01 202047860135_R03C01 202047860135_R04C01 
##                   1                   1                   1                   1 
## 202047860135_R07C01 202047860135_R08C01 
##                   1                   1 
## 
## $Sentrix_ID
## 
## 202010290146 202010290151 202010290176 202013790140 202047860135 
##            6            4            4            4            4 
## 
## $Sentrix_Position
## 
## R01C01 R02C01 R03C01 R04C01 R05C01 R06C01 R07C01 R08C01 
##      1      1      4      4      2      2      4      4 
## 
## $SampleName
## 
## ASD001.1 ASD001.2 ASD002.1 ASD002.2 ASD003.1 ASD003.2 ASD004.1 ASD004.2 
##        1        1        1        1        1        1        1        1 
## ASD005.1 ASD005.2 ASD006.1 ASD006.2 ASD007.1 ASD007.2 ASD008.1 ASD008.2 
##        1        1        1        1        1        1        1        1 
## ASD009.1 ASD009.2 ASD012.1 ASD012.2 ASD013.1 ASD013.2 
##        1        1        1        1        1        1 
## 
## $Sex
## 
##  1  2 
##  8 14 
## 
## $Age_months
## 
##  10  12  37  41  44  69  73  80  95 131 136 146 163 
##   1   1   2   2   2   2   2   2   2   1   2   1   2 
## 
## $Age_years
## 
##  5  6  7  8 10 11 13 14 17 
##  2  2  2  4  2  2  4  2  2 
## 
## $Birth_mode
## 
##  0  1  2  3 
##  1 14  6  1 
## 
## $Birth_complications
## 
##  0  1  2  3 
## 15  2  2  3 
## 
## $gestational_age
## 
## 31.0 34.0 35.0 36.0 36.4 36.7 37.0 37.5 38.0 
##    4    2    2    4    2    2    2    2    2 
## 
## $Birth_weight.g.
## 
##    2.500 1055.000 1551.000 1814.000 1817.000 1900.000 1927.700 2020.000 
##        1        1        1        1        1        1        1        1 
## 2093.000 2112.000 2220.000 2240.000 2296.000 2336.001 2400.000 2610.000 
##        1        1        1        1        1        1        2        1 
## 2630.840 2680.000 2684.000 3120.000 3395.000 
##        1        1        1        1        1 
## 
## $Birth_order
## 
##  1  2 
## 11 11 
## 
## $Help_breathing_start
## 
##  0  1 
## 13  9 
## 
## $Motor_skills
## 
##  0  1  2 
## 11  2  9 
## 
## $Social_interaction_on_ADOS
## 
##  1  3  5  7  8  9 10 11 12 13 15 16 17 18 20 
##  2  1  1  1  1  2  1  1  3  2  1  2  1  2  1 
## 
## $Diagnosis
## 
##  0  1  2 
##  5  6 11 
## 
## $SRS_social_scores
## 
##   4  34  35  45  47  52  61  73  74  78  79  84  94  95 100 102 103 104 111 121 
##   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1 
## 125 127 
##   1   1 
## 
## $language
## 
##  41  50  53  61  66  68  79  81  87  95  99 102 108 113 120 
##   1   5   1   1   2   1   1   1   1   1   2   1   1   2   1 
## 
## $IQ
## 
##   0  53  60  62  64  67  68  71  80  86  88  94  98 100 102 103 105 106 110 
##   1   1   2   2   1   1   1   1   1   1   1   1   1   1   1   1   2   1   1 
## 
## $Family_ID
## 
##  1  2  3  4  5  6  7  8  9 12 13 
##  2  2  2  2  2  2  2  2  2  2  2

Blood sample arm summary stats

Now for the follow-up blood samples.

message("summary stats")
## summary stats
summary(datb)
##    SampleID            ChipID            Sentrix_ID       Sentrix_Position  
##  Length:18          Length:18          Min.   :2.02e+11   Length:18         
##  Class :character   Class :character   1st Qu.:2.02e+11   Class :character  
##  Mode  :character   Mode  :character   Median :2.02e+11   Mode  :character  
##                                        Mean   :2.02e+11                     
##                                        3rd Qu.:2.02e+11                     
##                                        Max.   :2.02e+11                     
##  Sample_Name             Sex          Age_months       Age_years    
##  Length:18          Min.   :1.000   Min.   : 10.00   Min.   : 5.00  
##  Class :character   1st Qu.:1.000   1st Qu.: 44.00   1st Qu.: 8.00  
##  Mode  :character   Median :2.000   Median : 80.00   Median :11.00  
##                     Mean   :1.667   Mean   : 85.94   Mean   :10.89  
##                     3rd Qu.:2.000   3rd Qu.:134.75   3rd Qu.:13.00  
##                     Max.   :2.000   Max.   :163.00   Max.   :17.00  
##    Birth_mode    Birth_complications gestational_age Birth_weight.g. 
##  Min.   :0.000   Min.   :0.0000      Min.   :31.00   Min.   :   2.5  
##  1st Qu.:1.000   1st Qu.:0.0000      1st Qu.:35.00   1st Qu.:1906.9  
##  Median :1.000   Median :0.0000      Median :36.00   Median :2230.0  
##  Mean   :1.167   Mean   :0.5556      Mean   :35.57   Mean   :2109.9  
##  3rd Qu.:1.000   3rd Qu.:0.0000      3rd Qu.:36.70   3rd Qu.:2557.5  
##  Max.   :3.000   Max.   :3.0000      Max.   :38.00   Max.   :3395.0  
##   Birth_order  Help_breathing_start  Motor_skills    Social_interaction_on_ADOS
##  Min.   :1.0   Min.   :0.0000       Min.   :0.0000   Min.   : 1.00             
##  1st Qu.:1.0   1st Qu.:0.0000       1st Qu.:0.0000   1st Qu.: 9.25             
##  Median :1.5   Median :0.0000       Median :0.5000   Median :12.00             
##  Mean   :1.5   Mean   :0.3333       Mean   :0.8889   Mean   :11.61             
##  3rd Qu.:2.0   3rd Qu.:1.0000       3rd Qu.:2.0000   3rd Qu.:15.75             
##  Max.   :2.0   Max.   :1.0000       Max.   :2.0000   Max.   :20.00             
##    Diagnosis     SRS_social_scores    language            IQ        
##  Min.   :0.000   Min.   :  4.00    Min.   : 41.00   Min.   : 53.00  
##  1st Qu.:1.000   1st Qu.: 48.25    1st Qu.: 50.75   1st Qu.: 62.50  
##  Median :2.000   Median : 76.00    Median : 72.50   Median : 75.50  
##  Mean   :1.389   Mean   : 74.72    Mean   : 77.83   Mean   : 79.83  
##  3rd Qu.:2.000   3rd Qu.:102.25    3rd Qu.:101.25   3rd Qu.: 97.00  
##  Max.   :2.000   Max.   :125.00    Max.   :120.00   Max.   :110.00  
##    Family_ID     
##  Min.   : 1.000  
##  1st Qu.: 3.000  
##  Median : 6.000  
##  Mean   : 6.222  
##  3rd Qu.: 8.000  
##  Max.   :13.000
message("standard deviation")
## standard deviation
mysd <- apply(datb,2,sd)
## Warning in var(if (is.vector(x) || is.factor(x)) x else as.double(x), na.rm =
## na.rm): NAs introduced by coercion

## Warning in var(if (is.vector(x) || is.factor(x)) x else as.double(x), na.rm =
## na.rm): NAs introduced by coercion

## Warning in var(if (is.vector(x) || is.factor(x)) x else as.double(x), na.rm =
## na.rm): NAs introduced by coercion

## Warning in var(if (is.vector(x) || is.factor(x)) x else as.double(x), na.rm =
## na.rm): NAs introduced by coercion
format(mysd,scientific=FALSE)
##                   SampleID                     ChipID 
##         "              NA"         "              NA" 
##                 Sentrix_ID           Sentrix_Position 
##         "15710030.5380391"         "              NA" 
##                Sample_Name                        Sex 
##         "              NA"         "       0.4850713" 
##                 Age_months                  Age_years 
##         "      50.2212425"         "       3.6764575" 
##                 Birth_mode        Birth_complications 
##         "       0.6183469"         "       1.0966378" 
##            gestational_age            Birth_weight.g. 
##         "       1.9988232"         "     732.2514469" 
##                Birth_order       Help_breathing_start 
##         "       0.5144958"         "       0.4850713" 
##               Motor_skills Social_interaction_on_ADOS 
##         "       0.9633818"         "       5.6270240" 
##                  Diagnosis          SRS_social_scores 
##         "       0.7775443"         "      33.7356602" 
##                   language                         IQ 
##         "      26.8859246"         "      19.0610474" 
##                  Family_ID 
##         "       4.1096093"
message("categorical counts")
## categorical counts
mytab <- apply(datb,2,table)
mytab
## $SampleID
## 
##  S1 S10 S15 S16 S19  S2 S20 S25 S26 S29 S30 S33 S34 S37 S38  S5  S6  S9 
##   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1 
## 
## $ChipID
## 
## 202010290146_R03C01 202010290146_R04C01 202010290151_R01C01 202010290151_R02C01 
##                   1                   1                   1                   1 
## 202010290151_R05C01 202010290151_R06C01 202010290176_R01C01 202010290176_R02C01 
##                   1                   1                   1                   1 
## 202010290176_R05C01 202010290176_R06C01 202013790140_R01C01 202013790140_R02C01 
##                   1                   1                   1                   1 
## 202013790140_R07C01 202013790140_R08C01 202047860135_R01C01 202047860135_R02C01 
##                   1                   1                   1                   1 
## 202047860135_R05C01 202047860135_R06C01 
##                   1                   1 
## 
## $Sentrix_ID
## 
## 202010290146 202010290151 202010290176 202013790140 202047860135 
##            2            4            4            4            4 
## 
## $Sentrix_Position
## 
## R01C01 R02C01 R03C01 R04C01 R05C01 R06C01 R07C01 R08C01 
##      4      4      1      1      3      3      1      1 
## 
## $Sample_Name
## 
## ASD001.1 ASD001.2 ASD002.1 ASD002.2 ASD003.1 ASD003.2 ASD004.1 ASD004.2 
##        1        1        1        1        1        1        1        1 
## ASD006.1 ASD006.2 ASD007.1 ASD007.2 ASD008.1 ASD008.2 ASD012.1 ASD012.2 
##        1        1        1        1        1        1        1        1 
## ASD013.1 ASD013.2 
##        1        1 
## 
## $Sex
## 
##  1  2 
##  6 12 
## 
## $Age_months
## 
##  10  12  37  44  69  80  95 131 136 146 163 
##   1   1   2   2   2   2   2   1   2   1   2 
## 
## $Age_years
## 
##  5  7  8 10 11 13 14 17 
##  2  2  2  2  2  4  2  2 
## 
## $Birth_mode
## 
##  0  1  2  3 
##  1 14  2  1 
## 
## $Birth_complications
## 
##  0  2  3 
## 14  2  2 
## 
## $gestational_age
## 
## 31.0 34.0 35.0 36.0 36.4 36.7 37.0 38.0 
##    2    2    2    4    2    2    2    2 
## 
## $Birth_weight.g.
## 
##    2.500 1055.000 1551.000 1814.000 1900.000 1927.700 2020.000 2112.000 
##        1        1        1        1        1        1        1        1 
## 2220.000 2240.000 2336.001 2400.000 2610.000 2630.840 2680.000 2684.000 
##        1        1        1        2        1        1        1        1 
## 3395.000 
##        1 
## 
## $Birth_order
## 
## 1 2 
## 9 9 
## 
## $Help_breathing_start
## 
##  0  1 
## 12  6 
## 
## $Motor_skills
## 
## 0 1 2 
## 9 2 7 
## 
## $Social_interaction_on_ADOS
## 
##  1  3  8  9 10 11 12 13 15 16 17 18 20 
##  2  1  1  1  1  1  3  2  1  1  1  2  1 
## 
## $Diagnosis
## 
##  0  1  2 
##  3  5 10 
## 
## $SRS_social_scores
## 
##   4  34  35  45  47  52  61  73  74  78  84  94 100 103 104 111 121 125 
##   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1 
## 
## $language
## 
##  41  50  53  61  66  79  81  99 102 108 113 120 
##   1   4   1   1   2   1   1   2   1   1   2   1 
## 
## $IQ
## 
##  53  60  62  64  67  68  71  80  86  88  94  98 103 105 106 110 
##   1   2   2   1   1   1   1   1   1   1   1   1   1   1   1   1 
## 
## $Family_ID
## 
##  1  2  3  4  6  7  8 12 13 
##  2  2  2  2  2  2  2  2  2

Session information

sessionInfo()
## R version 4.2.2 Patched (2022-11-10 r83330)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 22.04.1 LTS
## 
## Matrix products: default
## BLAS:   /usr/lib/x86_64-linux-gnu/blas/libblas.so.3.10.0
## LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.10.0
## 
## locale:
##  [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
##  [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
##  [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
##  [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
##  [9] LC_ADDRESS=C               LC_TELEPHONE=C            
## [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       
## 
## attached base packages:
## [1] stats     graphics  grDevices utils     datasets  methods   base     
## 
## other attached packages:
## [1] magrittr_2.0.3   kableExtra_1.3.4
## 
## loaded via a namespace (and not attached):
##  [1] rstudioapi_0.14   knitr_1.41        xml2_1.3.3        munsell_0.5.0    
##  [5] rvest_1.0.3       viridisLite_0.4.1 colorspace_2.0-3  R6_2.5.1         
##  [9] rlang_1.0.6       fastmap_1.1.0     highr_0.9         stringr_1.5.0    
## [13] httr_1.4.4        tools_4.2.2       webshot_0.5.4     xfun_0.35        
## [17] cli_3.4.1         jquerylib_0.1.4   systemfonts_1.0.4 htmltools_0.5.4  
## [21] yaml_2.3.6        digest_0.6.31     lifecycle_1.0.3   sass_0.4.4       
## [25] vctrs_0.5.1       glue_1.6.2        cachem_1.0.6      evaluate_0.19    
## [29] rmarkdown_2.19    stringi_1.7.8     compiler_4.2.2    bslib_0.4.2      
## [33] scales_1.2.1      svglite_2.1.0     jsonlite_1.8.4