R Markdown

##This is now both phases ## For now - used last observation carried forward for missing data

##IGF1 ##First looking at Baseline levels of IGF-1 to look for carryover effects ## Then IGF1 response (z score) across arms

clin_labs_igf1z1
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 0.1262963 1.464228 29
B Baseline -0.0007143 1.397548 29
clin_labs_igf1z1
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 0.1262963 1.464228 29
A Week_12 0.0006897 1.290329 29
B Baseline -0.0007143 1.397548 29
B Week_12 1.9628571 1.640325 29

##Based on this, it doesn’t look like overall there was carryover of IGF1

##IGF1 model ##1) look for carryover

## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: clin_labs_igf1z1 ~ EVENT + DRUG + DRUG * EVENT + GROUP * DRUG +  
##     (1 | record_id)
##    Data: ghin_person
## 
## REML criterion at convergence: 361.3
## 
## Scaled residuals: 
##      Min       1Q   Median       3Q      Max 
## -2.30595 -0.48916 -0.07007  0.39400  2.36009 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 1.1179   1.057   
##  Residual              0.9643   0.982   
## Number of obs: 112, groups:  record_id, 29
## 
## Fixed effects:
##                    Estimate Std. Error      df t value Pr(>|t|)    
## (Intercept)          0.4113     0.3523 48.1690   1.168    0.249    
## EVENTWeek_12        -0.0893     0.2640 79.8319  -0.338    0.736    
## DRUGB               -0.3403     0.3207 80.4153  -1.061    0.292    
## GROUPB              -0.6656     0.4737 37.9042  -1.405    0.168    
## EVENTWeek_12:DRUGB   2.0529     0.3723 79.6416   5.515  4.2e-07 ***
## DRUGB:GROUPB         0.5164     0.3748 80.3126   1.378    0.172    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB  GROUPB EVENTW_12:
## EVENTWek_12 -0.375                                    
## DRUGB       -0.439  0.412                             
## GROUPB      -0.629 -0.026      0.201                  
## EVENTW_12:D  0.266 -0.709     -0.580  0.019           
## DRUGB:GROUP  0.231  0.033     -0.554 -0.394 -0.024

##here what we care about is group*drug interaction which is looking for carryover ## term not significant, so just on IGF1 it seems no carryover which is good

##MSC Bands 1-4, P60-N75 Difference - Long

vep_long_msc_band_1
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 0.3869310 0.1595297 29
A Week_12 0.3899655 0.1751706 29
B Baseline 0.3750000 0.1495586 26
B Week_12 0.3893846 0.1576759 26

vep_long_msc_band_2
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 0.3574828 0.1866616 29
A Week_12 0.3755517 0.2064568 29
B Baseline 0.3456154 0.1851177 26
B Week_12 0.3769231 0.1633859 26

vep_long_msc_band_3
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 0.2172414 0.1538758 29
A Week_12 0.1824138 0.1555469 29
B Baseline 0.1960385 0.1397242 26
B Week_12 0.2070000 0.1302097 26 

## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: formula
##    Data: data
## 
## REML criterion at convergence: -148.2
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -1.9441 -0.5711 -0.1579  0.4586  2.7582 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 0.014626 0.12094 
##  Residual              0.007158 0.08461 
## Number of obs: 110, groups:  record_id, 30
## 
## Fixed effects:
##                    Estimate Std. Error       df t value Pr(>|t|)    
## (Intercept)         0.22580    0.02716 48.28116   8.313 7.11e-11 ***
## EVENTWeek_12       -0.03483    0.02222 76.52746  -1.567   0.1211    
## DRUGB              -0.03928    0.02328 77.88610  -1.687   0.0956 .  
## EVENTWeek_12:DRUGB  0.04579    0.03232 76.52746   1.417   0.1606    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.409                  
## DRUGB       -0.397  0.477           
## EVENTW_12:D  0.281 -0.688     -0.694
vep_long_msc_band_4
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 0.1196552 0.0504446 29
A Week_12 0.1465517 0.0765980 29
B Baseline 0.1279231 0.0711303 26
B Week_12 0.1526154 0.0735011 26

PNAMP
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 4.157586 3.487043 29
A Week_12 4.316207 4.076654 29
B Baseline 3.729231 2.768107 26
B Week_12 4.730385 4.665066 26

iASD and PMS separate

Modeling for VEP, only in people with BL and W12

## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: vep_long_msc_band_2 ~ EVENT + DRUG + DRUG * EVENT + (1 | record_id)
##    Data: filtered_vep_long
## 
## REML criterion at convergence: -102.6
## 
## Scaled residuals: 
##      Min       1Q   Median       3Q      Max 
## -1.80611 -0.59700 -0.07776  0.59873  2.14426 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 0.02383  0.1544  
##  Residual              0.01081  0.1040  
## Number of obs: 110, groups:  record_id, 30
## 
## Fixed effects:
##                    Estimate Std. Error       df t value Pr(>|t|)    
## (Intercept)         0.36222    0.03424 47.69081  10.579 4.21e-14 ***
## EVENTWeek_12        0.01807    0.02730 77.10254   0.662    0.510    
## DRUGB              -0.03042    0.02862 78.37206  -1.063    0.291    
## EVENTWeek_12:DRUGB  0.01324    0.03971 77.10254   0.333    0.740    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.399                  
## DRUGB       -0.387  0.477           
## EVENTW_12:D  0.274 -0.688     -0.694
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: PNAMP ~ EVENT + DRUG + DRUG * EVENT + GROUP * DRUG + (1 | record_id)
##    Data: filtered_vep_long
## 
## REML criterion at convergence: 542.5
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -2.9456 -0.4546 -0.1376  0.3238  4.1422 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 8.930    2.988   
##  Residual              5.368    2.317   
## Number of obs: 110, groups:  record_id, 30
## 
## Fixed effects:
##                    Estimate Std. Error      df t value Pr(>|t|)    
## (Intercept)          3.9040     0.9310 42.9199   4.193 0.000135 ***
## EVENTWeek_12         0.1586     0.6084 76.6392   0.261 0.795022    
## DRUGB               -0.9895     0.7761 78.0688  -1.275 0.206105    
## GROUPB               0.5712     1.2535 35.7452   0.456 0.651376    
## EVENTWeek_12:DRUGB   0.8425     0.8849 76.6392   0.952 0.344047    
## DRUGB:GROUPB         0.6569     0.9174 79.5925   0.716 0.476074    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB  GROUPB EVENTW_12:
## EVENTWek_12 -0.327                                    
## DRUGB       -0.376  0.392                             
## GROUPB      -0.663  0.000      0.184                  
## EVENTW_12:D  0.225 -0.688     -0.570  0.000           
## DRUGB:GROUP  0.209  0.000     -0.571 -0.333  0.000
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: PNAMP ~ EVENT + DRUG + DRUG * EVENT + (1 | record_id)
##    Data: filtered_vep_long
## 
## REML criterion at convergence: 547.4
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -3.0051 -0.4408 -0.1187  0.3758  4.2588 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 8.841    2.973   
##  Residual              5.322    2.307   
## Number of obs: 110, groups:  record_id, 30
## 
## Fixed effects:
##                    Estimate Std. Error      df t value Pr(>|t|)    
## (Intercept)          4.1961     0.6932 53.6200   6.053 1.44e-07 ***
## EVENTWeek_12         0.1586     0.6058 77.6637   0.262    0.794    
## DRUGB               -0.6682     0.6344 79.1930  -1.053    0.295    
## EVENTWeek_12:DRUGB   0.8425     0.8811 77.6637   0.956    0.342    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.437                  
## DRUGB       -0.424  0.478           
## EVENTW_12:D  0.300 -0.688     -0.695
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: PNAMP ~ DRUG + EVENT + DRUG * EVENT + diag_dx___1 + (1 | record_id)
##    Data: filtered_vep_long
## 
## REML criterion at convergence: 543.5
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -3.0311 -0.4615 -0.0846  0.4143  4.2279 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 8.548    2.924   
##  Residual              5.323    2.307   
## Number of obs: 110, groups:  record_id, 30
## 
## Fixed effects:
##                    Estimate Std. Error      df t value Pr(>|t|)    
## (Intercept)          4.9721     0.8952 39.9246   5.554 2.01e-06 ***
## DRUGB               -0.6590     0.6344 79.1908  -1.039    0.302    
## EVENTWeek_12         0.1586     0.6059 77.6354   0.262    0.794    
## diag_dx___1PMS      -1.5657     1.1593 28.4886  -1.351    0.187    
## DRUGB:EVENTWeek_12   0.8425     0.8812 77.6354   0.956    0.342    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) DRUGB  EVENTW d____1
## DRUGB       -0.323                     
## EVENTWek_12 -0.338  0.478              
## dg_d___1PMS -0.642 -0.008  0.000       
## DRUGB:EVENT  0.233 -0.695 -0.688  0.000
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: PNAMP ~ EVENT + DRUG + DRUG * EVENT + (1 | record_id)
##    Data: ghin_person_pms
## 
## REML criterion at convergence: 555.3
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -3.0424 -0.4458 -0.1173  0.3664  4.2956 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 8.797    2.966   
##  Residual              5.212    2.283   
## Number of obs: 112, groups:  record_id, 30
## 
## Fixed effects:
##                    Estimate Std. Error      df t value Pr(>|t|)    
## (Intercept)          4.1959     0.6894 53.6920   6.086 1.26e-07 ***
## EVENTWeek_12         0.1586     0.5995 79.5640   0.265    0.792    
## DRUGB               -0.6503     0.6112 80.4549  -1.064    0.291    
## EVENTWeek_12:DRUGB   0.8292     0.8660 79.7246   0.957    0.341    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.435                  
## DRUGB       -0.433  0.490           
## EVENTW_12:D  0.301 -0.692     -0.691
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: vep_long_msc_band_2 ~ EVENT + DRUG + DRUG * EVENT + (1 | record_id)
##    Data: ghin_person_pms
## 
## REML criterion at convergence: -107.1
## 
## Scaled residuals: 
##      Min       1Q   Median       3Q      Max 
## -1.81734 -0.54208 -0.09023  0.60292  2.15939 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 0.02401  0.1549  
##  Residual              0.01054  0.1026  
## Number of obs: 112, groups:  record_id, 30
## 
## Fixed effects:
##                    Estimate Std. Error       df t value Pr(>|t|)    
## (Intercept)         0.36229    0.03418 47.48189  10.598 4.18e-14 ***
## EVENTWeek_12        0.01807    0.02696 79.06923   0.670    0.505    
## DRUGB              -0.03261    0.02750 79.79833  -1.186    0.239    
## EVENTWeek_12:DRUGB  0.01468    0.03894 79.19570   0.377    0.707    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.394                  
## DRUGB       -0.393  0.490           
## EVENTW_12:D  0.273 -0.692     -0.690

##Looking at VEP Responders in Group B Only - by amplitude difference - anything 0+

EVENT VEP_Responder MEAN_PNAMPl SD_PNAMP N_PNAMPl
Baseline 0 4.022500 2.831091 12
Week_12 0 2.887500 2.713787 12
Baseline 1 3.477857 2.793920 14
Week_12 1 6.310000 5.458746 14

EVENT VEP_Responder MEAN_sand_total SD_sand_total N_sand_total
Baseline 0 36.16667 12.939252 12
Week_12 0 31.25000 12.785823 12
Baseline 1 37.78571 9.641087 14
Week_12 1 34.21429 8.988699 14

EVENT VEP_Responder MEAN_sand_hyper SD_sand_hyper N_sand_hyper
Baseline 0 2.750000 3.744693 12
Week_12 0 1.333333 2.640018 12
Baseline 1 3.928571 3.647218 14
Week_12 1 3.285714 2.998168 14

EVENT VEP_Responder MEAN_sand_hypo SD_sand_hypo N_sand_hypo
Baseline 0 16.16667 7.334022 12
Week_12 0 14.00000 7.639610 12
Baseline 1 15.57143 6.477348 14
Week_12 1 13.28571 7.897697 14

EVENT VEP_Responder MEAN_sand_seek SD_sand_seek N_sand_seek
Baseline 0 17.25000 8.058028 12
Week_12 0 15.91667 6.430868 12
Baseline 1 18.28571 7.205004 14
Week_12 1 17.64286 6.096748 14 

## 
## Call:
## lm(formula = ghin_personDRUGW12$PNAMPDIFFSHORT ~ ghin_personDRUGW12$SANDTOTDIFF)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -8.8153 -2.9426 -1.2460  0.5828 27.1840 
## 
## Coefficients:
##                                Estimate Std. Error t value Pr(>|t|)
## (Intercept)                      2.8890     2.0666   1.398    0.176
## ghin_personDRUGW12$SANDTOTDIFF   0.2377     0.3283   0.724    0.477
## 
## Residual standard error: 7.01 on 22 degrees of freedom
##   (2 observations deleted due to missingness)
## Multiple R-squared:  0.02327,    Adjusted R-squared:  -0.02113 
## F-statistic: 0.5241 on 1 and 22 DF,  p-value: 0.4767

###NOW LOOK AT SHORT VEP

PNAMPshort
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 7.040000 5.342126 26
A Week_12 7.619615 6.473463 26
B Baseline 6.510800 3.575760 25
B Week_12 8.346400 7.352989 25 

## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: formula
##    Data: data
## 
## REML criterion at convergence: 616.3
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -2.1162 -0.4550 -0.1337  0.3631  4.1517 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 14.59    3.819   
##  Residual              19.28    4.391   
## Number of obs: 102, groups:  record_id, 28
## 
## Fixed effects:
##                    Estimate Std. Error      df t value Pr(>|t|)    
## (Intercept)          6.9134     1.1308 69.4301   6.114 5.05e-08 ***
## EVENTWeek_12         0.5796     1.2179 72.0315   0.476    0.636    
## DRUGB               -0.3139     1.2492 74.4845  -0.251    0.802    
## EVENTWeek_12:DRUGB   1.2560     1.7396 72.0315   0.722    0.473    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.539                  
## DRUGB       -0.539  0.488           
## EVENTW_12:D  0.377 -0.700     -0.696
vep_short_msc_band_1
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 0.3081538 0.1128478 26
A Week_12 0.3017308 0.1279195 26
B Baseline 0.2432400 0.1321238 25
B Week_12 0.3006800 0.1269313 25 

## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: formula
##    Data: data
## 
## REML criterion at convergence: -155.6
## 
## Scaled residuals: 
##      Min       1Q   Median       3Q      Max 
## -1.97113 -0.46657  0.03764  0.45023  2.73097 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 0.009077 0.09527 
##  Residual              0.006394 0.07996 
## Number of obs: 102, groups:  record_id, 28
## 
## Fixed effects:
##                     Estimate Std. Error        df t value Pr(>|t|)    
## (Intercept)         0.307119   0.024014 54.177804  12.789  < 2e-16 ***
## EVENTWeek_12       -0.006423   0.022177 71.374934  -0.290  0.77294    
## DRUGB              -0.064152   0.022830 73.159989  -2.810  0.00635 ** 
## EVENTWeek_12:DRUGB  0.063863   0.031676 71.374934   2.016  0.04755 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.462                  
## DRUGB       -0.463  0.486           
## EVENTW_12:D  0.323 -0.700     -0.694
vep_short_msc_band_2
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 0.3165385 0.1215984 26
A Week_12 0.3054615 0.1531201 26
B Baseline 0.2904000 0.1067891 25
B Week_12 0.3198800 0.1437594 25 

## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: formula
##    Data: data
## 
## REML criterion at convergence: -149.8
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -1.8769 -0.5436 -0.1025  0.5895  2.5027 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 0.01102  0.10496 
##  Residual              0.00651  0.08069 
## Number of obs: 102, groups:  record_id, 28
## 
## Fixed effects:
##                    Estimate Std. Error       df t value Pr(>|t|)    
## (Intercept)         0.31213    0.02551 50.34272  12.235   <2e-16 ***
## EVENTWeek_12       -0.01108    0.02238 71.14406  -0.495    0.622    
## DRUGB              -0.01883    0.02306 72.74789  -0.817    0.417    
## EVENTWeek_12:DRUGB  0.04056    0.03196 71.14406   1.269    0.209    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.439                  
## DRUGB       -0.440  0.485           
## EVENTW_12:D  0.307 -0.700     -0.693
vep_short_msc_band_3
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 0.2107692 0.1112968 26
A Week_12 0.2227692 0.1450455 26
B Baseline 0.2009200 0.1007724 25
B Week_12 0.2150400 0.1502951 25 

## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: formula
##    Data: data
## 
## REML criterion at convergence: -128.9
## 
## Scaled residuals: 
##      Min       1Q   Median       3Q      Max 
## -1.64474 -0.56054 -0.08382  0.39612  3.03860 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 0.006519 0.08074 
##  Residual              0.009865 0.09933 
## Number of obs: 102, groups:  record_id, 28
## 
## Fixed effects:
##                     Estimate Std. Error        df t value Pr(>|t|)    
## (Intercept)         0.209249   0.024901 72.522218   8.403 2.57e-12 ***
## EVENTWeek_12        0.012000   0.027548 71.961462   0.436    0.664    
## DRUGB              -0.009336   0.028229 74.562500  -0.331    0.742    
## EVENTWeek_12:DRUGB  0.002120   0.039346 71.961462   0.054    0.957    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.553                  
## DRUGB       -0.553  0.488           
## EVENTW_12:D  0.387 -0.700     -0.697

#look at correlations 

## 
## Call:
## lm(formula = ghin_person_filtered_vep_short_W12DRUG$PNMSC1DIFFSHORT ~ 
##     ghin_person_filtered_vep_short_W12DRUG$SANDTOTDIFFSHORT)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.178423 -0.059311 -0.003867  0.079300  0.210913 
## 
## Coefficients:
##                                                          Estimate Std. Error
## (Intercept)                                             -0.002245   0.023143
## ghin_person_filtered_vep_short_W12DRUG$SANDTOTDIFFSHORT -0.004389   0.003870
##                                                         t value Pr(>|t|)
## (Intercept)                                              -0.097    0.924
## ghin_person_filtered_vep_short_W12DRUG$SANDTOTDIFFSHORT  -1.134    0.268
## 
## Residual standard error: 0.1106 on 23 degrees of freedom
## Multiple R-squared:  0.05296,    Adjusted R-squared:  0.01178 
## F-statistic: 1.286 on 1 and 23 DF,  p-value: 0.2684

## 
## Call:
## lm(formula = ghin_person_filtered_vep_short_W12DRUG$PNMSC1DIFFSHORT ~ 
##     ghin_person_filtered_vep_short_W12DRUG$ABCSWDIFFSHORT)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.160866 -0.084722  0.003061  0.093848  0.237419 
## 
## Coefficients:
##                                                       Estimate Std. Error
## (Intercept)                                           0.006794   0.023255
## ghin_person_filtered_vep_short_W12DRUG$ABCSWDIFFSHORT 0.001642   0.006343
##                                                       t value Pr(>|t|)
## (Intercept)                                             0.292    0.773
## ghin_person_filtered_vep_short_W12DRUG$ABCSWDIFFSHORT   0.259    0.798
## 
## Residual standard error: 0.1135 on 23 degrees of freedom
## Multiple R-squared:  0.002907,   Adjusted R-squared:  -0.04045 
## F-statistic: 0.06705 on 1 and 23 DF,  p-value: 0.798

## 
## Call:
## lm(formula = ghin_person_filtered_vep_short_W12DRUG$PNMSC1DIFFSHORT ~ 
##     ghin_person_filtered_vep_short_W12DRUG$CGIIDIFFSHORT)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -0.16377 -0.07955 -0.00832  0.09245  0.23300 
## 
## Coefficients:
##                                                       Estimate Std. Error
## (Intercept)                                           0.006545   0.024761
## ghin_person_filtered_vep_short_W12DRUG$CGIIDIFFSHORT -0.001775   0.016398
##                                                      t value Pr(>|t|)
## (Intercept)                                            0.264    0.794
## ghin_person_filtered_vep_short_W12DRUG$CGIIDIFFSHORT  -0.108    0.915
## 
## Residual standard error: 0.1136 on 23 degrees of freedom
## Multiple R-squared:  0.0005092,  Adjusted R-squared:  -0.04295 
## F-statistic: 0.01172 on 1 and 23 DF,  p-value: 0.9147

#NEED TO FINISH

PNAMP
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 2.872143 3.533136 14
A Week_12 3.725714 4.004147 14
B Baseline 3.430769 2.758757 13
B Week_12 3.627692 3.085820 13 

## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: formula
##    Data: data
## 
## REML criterion at convergence: 227.3
## 
## Scaled residuals: 
##      Min       1Q   Median       3Q      Max 
## -2.00059 -0.51343 -0.01047  0.38931  2.74783 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 8.901    2.983   
##  Residual              2.065    1.437   
## Number of obs: 54, groups:  record_id, 15
## 
## Fixed effects:
##                    Estimate Std. Error      df t value Pr(>|t|)   
## (Intercept)          3.0125     0.8633 19.2425   3.489  0.00242 **
## EVENTWeek_12         0.8536     0.5431 36.1728   1.572  0.12475   
## DRUGB                0.1768     0.5684 36.6261   0.311  0.75749   
## EVENTWeek_12:DRUGB  -0.6566     0.7827 36.1728  -0.839  0.40701   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.315                  
## DRUGB       -0.312  0.478           
## EVENTW_12:D  0.218 -0.694     -0.689
PNAMP
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 5.357333 3.087550 15
A Week_12 4.867333 4.204300 15
B Baseline 4.027692 2.856505 13
B Week_12 5.833077 5.760421 13 

## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: formula
##    Data: data
## 
## REML criterion at convergence: 290.2
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -2.3165 -0.4477 -0.1338  0.4247  3.4652 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 8.208    2.865   
##  Residual              8.425    2.903   
## Number of obs: 56, groups:  record_id, 15
## 
## Fixed effects:
##                    Estimate Std. Error     df t value Pr(>|t|)    
## (Intercept)           5.357      1.053 31.377   5.088 1.62e-05 ***
## EVENTWeek_12         -0.490      1.060 38.372  -0.462    0.646    
## DRUGB                -1.498      1.113 39.124  -1.346    0.186    
## EVENTWeek_12:DRUGB    2.295      1.556 38.372   1.476    0.148    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.503                  
## DRUGB       -0.479  0.476           
## EVENTW_12:D  0.343 -0.681     -0.699

vep_short_msc_band_1
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 0.3130909 0.1209276 11
A Week_12 0.3145455 0.1258303 11
B Baseline 0.2637143 0.1444099 14
B Week_12 0.3027500 0.1483118 12 

## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: formula
##    Data: data
## 
## REML criterion at convergence: -54
## 
## Scaled residuals: 
##      Min       1Q   Median       3Q      Max 
## -1.65713 -0.42658 -0.09334  0.39448  2.27724 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 0.008844 0.09404 
##  Residual              0.008959 0.09465 
## Number of obs: 48, groups:  record_id, 14
## 
## Fixed effects:
##                     Estimate Std. Error        df t value Pr(>|t|)    
## (Intercept)         0.307928   0.038891 32.766408   7.918 4.15e-09 ***
## EVENTWeek_12        0.001455   0.040361 32.133379   0.036    0.971    
## DRUGB              -0.044214   0.038997 33.497059  -1.134    0.265    
## EVENTWeek_12:DRUGB  0.031701   0.055268 32.544233   0.574    0.570    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.519                  
## DRUGB       -0.581  0.517           
## EVENTW_12:D  0.388 -0.730     -0.684
vep_short_msc_band_1
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 0.3045333 0.1107332 15
A Week_12 0.2923333 0.1329934 15
B Baseline 0.2246923 0.1090110 13
B Week_12 0.2987692 0.1097794 13 

## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: formula
##    Data: data
## 
## REML criterion at convergence: -89.9
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -2.0911 -0.4697  0.1387  0.4677  1.8141 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 0.008839 0.09402 
##  Residual              0.004926 0.07019 
## Number of obs: 56, groups:  record_id, 15
## 
## Fixed effects:
##                    Estimate Std. Error       df t value Pr(>|t|)    
## (Intercept)         0.30453    0.03029 24.13637  10.053 4.18e-10 ***
## EVENTWeek_12       -0.01220    0.02563 37.88537  -0.476   0.6368    
## DRUGB              -0.07264    0.02697 38.40942  -2.694   0.0104 *  
## EVENTWeek_12:DRUGB  0.08628    0.03761 37.88537   2.294   0.0274 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.423                  
## DRUGB       -0.402  0.475           
## EVENTW_12:D  0.288 -0.681     -0.697

##EEG

ASSR20hz
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 0.2020436 0.0679459 29
A Week_12 0.1937031 0.0787733 29
B Baseline 0.1960980 0.1151193 29
B Week_12 0.2238967 0.1269851 29 

## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: formula
##    Data: data
## 
## REML criterion at convergence: -184
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -1.4444 -0.6002 -0.1934  0.2910  4.4693 
## 
## Random effects:
##  Groups    Name        Variance  Std.Dev.
##  record_id (Intercept) 0.0002036 0.01427 
##  Residual              0.0098463 0.09923 
## Number of obs: 116, groups:  record_id, 30
## 
## Fixed effects:
##                      Estimate Std. Error         df t value Pr(>|t|)    
## (Intercept)          0.201917   0.018616 111.878145  10.847   <2e-16 ***
## EVENTWeek_12        -0.008341   0.026059  83.553720  -0.320    0.750    
## DRUGB               -0.005792   0.026068  85.015849  -0.222    0.825    
## EVENTWeek_12:DRUGB   0.036139   0.036852  83.553720   0.981    0.330    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.700                  
## DRUGB       -0.700  0.500           
## EVENTW_12:D  0.495 -0.707     -0.707
ASSR40hz
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 0.2094240 0.0708492 29
A Week_12 0.2118860 0.0813338 29
B Baseline 0.2089573 0.0896188 29
B Week_12 0.1953981 0.0698684 29 

## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: formula
##    Data: data
## 
## REML criterion at convergence: -243.4
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -1.4263 -0.6737 -0.1092  0.4607  2.6197 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 0.001107 0.03327 
##  Residual              0.005044 0.07102 
## Number of obs: 116, groups:  record_id, 30
## 
## Fixed effects:
##                      Estimate Std. Error         df t value Pr(>|t|)    
## (Intercept)          0.208262   0.014551 102.636790  14.313   <2e-16 ***
## EVENTWeek_12         0.002462   0.018652  82.044760   0.132    0.895    
## DRUGB                0.001293   0.018701  83.425361   0.069    0.945    
## EVENTWeek_12:DRUGB  -0.016021   0.026378  82.044760  -0.607    0.545    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.641                  
## DRUGB       -0.643  0.499           
## EVENTW_12:D  0.453 -0.707     -0.705
Alpha
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline -0.1787149 0.1788329 29
A Week_12 -0.2008053 0.1580765 29
B Baseline -0.2231300 0.2349026 29
B Week_12 -0.2359144 0.2666709 29

Beta
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline -0.5216432 0.1697310 29
A Week_12 -0.5544043 0.1522185 29
B Baseline -0.5742507 0.2290404 29
B Week_12 -0.5629814 0.2618475 29

Delta
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 0.1196562 0.2088669 29
A Week_12 0.0864912 0.1576466 29
B Baseline 0.0814892 0.2559082 29
B Week_12 0.0515846 0.2859633 29

Gamma
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline -1.001249 0.1887615 29
A Week_12 -1.064516 0.1584065 29
B Baseline -1.075087 0.2550124 29
B Week_12 -1.081301 0.2781952 29

Theta
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline -0.0223727 0.1878619 29
A Week_12 -0.0492523 0.1605966 29
B Baseline -0.0579275 0.2359005 29
B Week_12 -0.0875284 0.2659000 29

Total
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline -0.6359379 0.1743302 29
A Week_12 -0.6769025 0.1445330 29
B Baseline -0.6909275 0.2369968 29
B Week_12 -0.7061221 0.2684389 29

##SAND Results

sand_total
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 37.00000 10.67708 29
A Week_12 35.06897 10.58277 29
B Baseline 36.62069 11.73339 29
B Week_12 32.51724 11.45993 29
sand_hyperreactivity
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 3.448276 4.289407 29
A Week_12 2.620690 3.288658 29
B Baseline 3.448276 3.679928 29
B Week_12 2.517241 3.019232 29
sand_hyporeactivity
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 16.51724 8.279151 29
A Week_12 15.17241 6.990141 29
B Baseline 15.62069 7.153490 29
B Week_12 13.34483 7.719529 29
sand_seeking
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 17.03448 5.615557 29
A Week_12 17.27586 6.250123 29
B Baseline 17.55172 7.174290 29
B Week_12 16.65517 6.096110 29

#stats

## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: sand_hyporeactivity ~ EVENT + DRUG + DRUG * EVENT + GROUP * DRUG +  
##     (1 | record_id)
##    Data: ghin_person
## 
## REML criterion at convergence: 658
## 
## Scaled residuals: 
##      Min       1Q   Median       3Q      Max 
## -2.31862 -0.53638  0.03826  0.47352  2.11892 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 50.72    7.122   
##  Residual               8.80    2.967   
## Number of obs: 116, groups:  record_id, 30
## 
## Fixed effects:
##                    Estimate Std. Error      df t value Pr(>|t|)    
## (Intercept)         16.4724     1.9562 32.6256   8.421 1.09e-09 ***
## EVENTWeek_12        -1.3448     0.7790 81.9551  -1.726   0.0881 .  
## DRUGB               -0.6185     0.9646 82.2393  -0.641   0.5232    
## GROUPB              -0.6564     2.7184 30.4535  -0.241   0.8108    
## EVENTWeek_12:DRUGB  -0.9310     1.1017 81.9551  -0.845   0.4005    
## DRUGB:GROUPB         0.3404     1.1197 82.3769   0.304   0.7619    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB  GROUPB EVENTW_12:
## EVENTWek_12 -0.199                                    
## DRUGB       -0.236  0.404                             
## GROUPB      -0.691  0.000      0.112                  
## EVENTW_12:D  0.141 -0.707     -0.571  0.000           
## DRUGB:GROUP  0.134  0.000     -0.580 -0.206  0.000
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: formula
##    Data: data
## 
## REML criterion at convergence: 737.6
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -2.6036 -0.5406 -0.0356  0.4705  2.1507 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 118.38   10.88   
##  Residual               15.61    3.95   
## Number of obs: 116, groups:  record_id, 30
## 
## Fixed effects:
##                    Estimate Std. Error      df t value Pr(>|t|)    
## (Intercept)         36.0451     2.1196 34.6025  17.006   <2e-16 ***
## EVENTWeek_12        -1.9310     1.0374 82.6890  -1.861   0.0662 .  
## DRUGB                0.5245     1.0461 82.8709   0.501   0.6174    
## EVENTWeek_12:DRUGB  -2.1724     1.4671 82.6890  -1.481   0.1425    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.245                  
## DRUGB       -0.247  0.496           
## EVENTW_12:D  0.173 -0.707     -0.701
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: formula
##    Data: data
## 
## REML criterion at convergence: 516.1
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -1.8186 -0.5663 -0.0194  0.3350  3.3607 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 10.298   3.209   
##  Residual               2.513   1.585   
## Number of obs: 116, groups:  record_id, 30
## 
## Fixed effects:
##                    Estimate Std. Error      df t value Pr(>|t|)    
## (Intercept)          3.3551     0.6567 39.9882   5.109 8.36e-06 ***
## EVENTWeek_12        -0.8276     0.4163 83.1022  -1.988   0.0501 .  
## DRUGB                0.1410     0.4196 83.4111   0.336   0.7376    
## EVENTWeek_12:DRUGB  -0.1034     0.5887 83.1022  -0.176   0.8609    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.317                  
## DRUGB       -0.319  0.496           
## EVENTW_12:D  0.224 -0.707     -0.702
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: formula
##    Data: data
## 
## REML criterion at convergence: 663.9
## 
## Scaled residuals: 
##      Min       1Q   Median       3Q      Max 
## -2.36443 -0.56527  0.06095  0.50252  2.15989 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 48.912   6.994   
##  Residual               8.708   2.951   
## Number of obs: 116, groups:  record_id, 30
## 
## Fixed effects:
##                    Estimate Std. Error      df t value Pr(>|t|)    
## (Intercept)         16.1511     1.3911 36.9297  11.610 6.93e-14 ***
## EVENTWeek_12        -1.3448     0.7749 82.9435  -1.735   0.0864 .  
## DRUGB               -0.4506     0.7813 83.1800  -0.577   0.5657    
## EVENTWeek_12:DRUGB  -0.9310     1.0959 82.9435  -0.850   0.3980    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.279                  
## DRUGB       -0.281  0.496           
## EVENTW_12:D  0.197 -0.707     -0.701
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: formula
##    Data: data
## 
## REML criterion at convergence: 631.1
## 
## Scaled residuals: 
##      Min       1Q   Median       3Q      Max 
## -1.86079 -0.62820  0.05043  0.60648  1.75844 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 36.071   6.006   
##  Residual               6.518   2.553   
## Number of obs: 116, groups:  record_id, 30
## 
## Fixed effects:
##                    Estimate Std. Error      df t value Pr(>|t|)    
## (Intercept)         16.5625     1.1961 36.6916  13.847 3.82e-16 ***
## EVENTWeek_12         0.2414     0.6704 82.6337   0.360    0.720    
## DRUGB                0.8095     0.6759 82.8748   1.198    0.234    
## EVENTWeek_12:DRUGB  -1.1379     0.9482 82.6337  -1.200    0.234    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.280                  
## DRUGB       -0.283  0.496           
## EVENTW_12:D  0.198 -0.707     -0.701

##CGI Improvement next up

cgii_statistics
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 0.0000000 0.0000000 29
A Week_4 0.7241379 0.7971625 29
A Week_8 0.6206897 0.9788402 29
A Week_12 0.8620690 0.9533452 29
B Baseline 0.0000000 0.0000000 29
B Week_4 0.5172414 0.9863602 29
B Week_8 0.9655172 1.0516232 29
B Week_12 1.1379310 1.1251710 29 

## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: cgii_statistics ~ EVENT + DRUG + DRUG * EVENT + GROUP * DRUG +  
##     (1 | record_id)
##    Data: ghin_person
## 
## REML criterion at convergence: 264.2
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -4.0476 -0.1812  0.0121  0.2475  1.6259 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 0.04046  0.2012  
##  Residual              0.50912  0.7135  
## Number of obs: 116, groups:  record_id, 30
## 
## Fixed effects:
##                    Estimate Std. Error       df t value Pr(>|t|)    
## (Intercept)         0.06897    0.16866 97.64861   0.409    0.684    
## EVENTWeek_12        0.86207    0.18738 80.93740   4.601 1.54e-05 ***
## DRUGB              -0.13793    0.22995 82.88868  -0.600    0.550    
## GROUPB             -0.15396    0.20178 69.02047  -0.763    0.448    
## EVENTWeek_12:DRUGB  0.27586    0.26500 80.93740   1.041    0.301    
## DRUGB:GROUPB        0.28729    0.26578 83.85848   1.081    0.283    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB  GROUPB EVENTW_12:
## EVENTWek_12 -0.556                                    
## DRUGB       -0.664  0.407                             
## GROUPB      -0.578  0.000      0.366                  
## EVENTW_12:D  0.393 -0.707     -0.576  0.000           
## DRUGB:GROUP  0.379  0.000     -0.578 -0.659  0.000
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: cgii_statistics ~ EVENT + DRUG + DRUG * EVENT + (1 | record_id)
##    Data: ghin_person
## 
## REML criterion at convergence: 262.6
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -4.1784 -0.1317 -0.0009  0.1977  1.5260 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 0.03122  0.1767  
##  Residual              0.51263  0.7160  
## Number of obs: 116, groups:  record_id, 30
## 
## Fixed effects:
##                      Estimate Std. Error         df t value Pr(>|t|)    
## (Intercept)         -0.004003   0.136929 110.973121  -0.029    0.977    
## EVENTWeek_12         0.862069   0.188026  82.251801   4.585 1.61e-05 ***
## DRUGB                0.004277   0.188202  83.744485   0.023    0.982    
## EVENTWeek_12:DRUGB   0.275862   0.265908  82.251801   1.037    0.303    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.687                  
## DRUGB       -0.687  0.500           
## EVENTW_12:D  0.485 -0.707     -0.706
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: cgii_statistics ~ EVENT + DRUG + DRUG * EVENT + (1 | record_id)
##    Data: gh
## 
## REML criterion at convergence: 575.6
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -3.7910 -0.4663  0.0453  0.5183  2.1239 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 0.1471   0.3836  
##  Residual              0.5905   0.7685  
## Number of obs: 232, groups:  record_id, 30
## 
## Fixed effects:
##                     Estimate Std. Error        df t value Pr(>|t|)    
## (Intercept)         -0.02425    0.15922 175.96192  -0.152 0.879118    
## EVENTWeek_4          0.72414    0.20181 192.97476   3.588 0.000422 ***
## EVENTWeek_8          0.62069    0.20181 192.97476   3.076 0.002405 ** 
## EVENTWeek_12         0.86207    0.20181 192.97476   4.272 3.04e-05 ***
## DRUGB                0.02808    0.20225 194.36814   0.139 0.889717    
## EVENTWeek_4:DRUGB   -0.20690    0.28540 192.97476  -0.725 0.469368    
## EVENTWeek_8:DRUGB    0.34483    0.28540 192.97476   1.208 0.228438    
## EVENTWeek_12:DRUGB   0.27586    0.28540 192.97476   0.967 0.334962    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_4 EVENTWk_8 EVENTWk_12 DRUGB  EVENTW_4: EVENTW_8:
## EVENTWeek_4 -0.634                                                          
## EVENTWeek_8 -0.634  0.500                                                   
## EVENTWek_12 -0.634  0.500     0.500                                         
## DRUGB       -0.635  0.499     0.499     0.499                               
## EVENTW_4:DR  0.448 -0.707    -0.354    -0.354     -0.706                    
## EVENTW_8:DR  0.448 -0.354    -0.707    -0.354     -0.706  0.500             
## EVENTW_12:D  0.448 -0.354    -0.354    -0.707     -0.706  0.500     0.500

##no significant carry over ##CGI-I was different BL to W12 before accounting for drop out, once LOCF it is no longer significant ##CGI-I not significantly different when we look at all timepoints

##ABC

abc_subscale1_irritability
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 10.758621 9.660875 29
A Week_4 9.068965 8.280489 29
A Week_8 9.000000 7.101308 29
A Week_12 8.344828 7.162780 29
B Baseline 10.896552 9.347822 29
B Week_4 8.344828 8.844998 29
B Week_8 7.827586 7.755500 29
B Week_12 7.793103 7.451650 29
abc_subscale2_lethargy
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 13.65517 8.469592 29
A Week_4 12.37931 7.997383 29
A Week_8 11.17241 7.942062 29
A Week_12 10.65517 7.560593 29
B Baseline 14.68966 10.295989 29
B Week_4 11.44828 8.862385 29
B Week_8 10.24138 7.404898 29
B Week_12 10.31034 7.714103 29
abc_subscale3_stereotypy
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 7.172414 5.251319 29
A Week_4 5.551724 4.428691 29
A Week_8 5.724138 4.207956 29
A Week_12 5.413793 5.348680 29
B Baseline 6.965517 5.421404 29
B Week_4 5.620690 4.828589 29
B Week_8 4.793103 4.091623 29
B Week_12 5.000000 4.590363 29
abc_subscale4_hyperact
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 23.44828 11.296606 29
A Week_4 20.20690 10.366082 29
A Week_8 18.89655 9.828329 29
A Week_12 18.82759 9.914164 29
B Baseline 24.96552 9.578558 29
B Week_4 18.96552 9.469812 29
B Week_8 18.89655 9.611544 29
B Week_12 18.68966 10.423543 29
abc_subscale5_in_speech
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 3.034483 3.877432 29
A Week_4 2.379310 3.004922 29
A Week_8 2.379310 2.704676 29
A Week_12 2.172414 2.633154 29
B Baseline 2.793103 3.233522 29
B Week_4 2.103448 2.743106 29
B Week_8 2.103448 2.730055 29
B Week_12 1.931034 3.161499 29

#ABC graphs

##MCDI

cdi_vocab_understood
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 196.3929 141.4982 29
A Week_12 222.9655 138.8899 29
B Baseline 213.8966 135.3615 29
B Week_12 232.6207 147.7223 29

cdi_vocab_produced
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 70.57143 113.2066 29
A Week_12 87.89655 129.8299 29
B Baseline 72.17241 119.6143 29
B Week_12 110.00000 145.6422 29

cdi_vocab_produced_amended_score
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 75.35714 114.7536 29
A Week_12 95.31034 132.4403 29
B Baseline 96.24138 133.7865 29
B Week_12 126.13793 158.8379 29 

## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: formula
##    Data: data
## 
## REML criterion at convergence: 1281.8
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -5.5921 -0.3811 -0.0512  0.3814  2.5877 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 15530    124.6   
##  Residual               2276     47.7   
## Number of obs: 115, groups:  record_id, 30
## 
## Fixed effects:
##                    Estimate Std. Error      df t value Pr(>|t|)    
## (Intercept)          89.205     24.516  35.600   3.639 0.000862 ***
## EVENTWeek_12          7.974     12.674  81.698   0.629 0.530994    
## DRUGB               -20.026     12.781  81.902  -1.567 0.121009    
## EVENTWeek_12:DRUGB   29.853     17.821  81.676   1.675 0.097717 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.264                  
## DRUGB       -0.267  0.508           
## EVENTW_12:D  0.188 -0.711     -0.706

##RBS-R

stereo_behav_score
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 5.551724 4.005231 29
A Week_4 4.310345 3.285286 29
A Week_8 3.931034 3.283411 29
A Week_12 3.448276 3.344358 29
B Baseline 5.068966 4.471585 29
B Week_4 4.517241 3.728469 29
B Week_8 3.896552 3.199831 29
B Week_12 3.965517 3.489604 29

self_injur_score
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 2.827586 3.742316 29
A Week_4 2.517241 3.323302 29
A Week_8 2.068966 2.814897 29
A Week_12 2.068966 3.022901 29
B Baseline 2.724138 3.981175 29
B Week_4 1.931034 3.453063 29
B Week_8 1.827586 3.117565 29
B Week_12 1.931034 3.442705 29

comp_behav_score
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 3.413793 3.914417 29
A Week_4 2.827586 3.129000 29
A Week_8 1.931034 2.576418 29
A Week_12 2.241379 3.043203 29
B Baseline 2.241379 3.521923 29
B Week_4 2.137931 3.430162 29
B Week_8 2.482759 2.785430 29
B Week_12 2.034483 2.665691 29

ritu_behav_score
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 2.827586 3.306584 29
A Week_4 2.206897 2.677676 29
A Week_8 2.310345 3.263472 29
A Week_12 2.034483 2.337318 29
B Baseline 2.448276 3.611353 29
B Week_4 2.448276 3.123485 29
B Week_8 2.241379 3.019640 29
B Week_12 2.000000 2.591194 29

same_behav_score
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 5.517241 6.127343 29
A Week_4 4.862069 4.348433 29
A Week_8 4.862069 4.381162 29
A Week_12 4.068966 4.308886 29
B Baseline 4.931034 4.832413 29
B Week_4 4.586207 4.982234 29
B Week_8 4.206897 4.776275 29
B Week_12 4.068966 5.063391 29

restric_behav_score
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 3.689655 2.916743 29
A Week_4 3.827586 2.928541 29
A Week_8 2.965517 2.337318 29
A Week_12 3.172414 2.592144 29
B Baseline 3.931034 2.939036 29
B Week_4 3.724138 2.914208 29
B Week_8 3.241379 2.262350 29
B Week_12 3.275862 2.589768 29

overall_score
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 23.82759 17.01275 29
A Week_4 20.55172 14.77108 29
A Week_8 18.06897 13.51964 29
A Week_12 17.03448 13.68389 29
B Baseline 21.34483 18.06590 29
B Week_4 19.34483 16.41053 29
B Week_8 17.89655 14.19594 29
B Week_12 17.27586 14.45411 29

##ADAMS

adams_manic_score
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 8.620690 3.363818 29
A Week_4 7.896552 3.885365 29
A Week_8 6.827586 3.864388 29
A Week_12 7.551724 4.330838 29
B Baseline 8.551724 4.014138 29
B Week_4 7.758621 3.621909 29
B Week_8 7.137931 3.681267 29
B Week_12 6.931034 3.453063 29

adams_depressed_score
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 3.896552 3.233141 29
A Week_4 2.344828 2.224472 29
A Week_8 2.172414 1.712746 29
A Week_12 2.206897 2.144417 29
B Baseline 2.965517 3.448424 29
B Week_4 2.758621 3.180916 29
B Week_8 2.482759 2.293168 29
B Week_12 2.413793 2.260171 29

adams_soc_avoid_score
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 8.620690 4.813773 29
A Week_4 7.655172 5.293364 29
A Week_8 7.241379 5.068496 29
A Week_12 7.413793 4.939037 29
B Baseline 8.413793 5.130561 29
B Week_4 7.344828 4.915543 29
B Week_8 6.620690 3.839771 29
B Week_12 6.551724 4.917797 29

adams_obsessive_score
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 2.827586 2.315614 29
A Week_4 2.000000 2.104417 29
A Week_8 1.965517 1.917639 29
A Week_12 1.827586 1.852970 29
B Baseline 2.344828 2.125960 29
B Week_4 1.931034 1.998768 29
B Week_8 2.172414 1.909917 29
B Week_12 2.413793 1.973345 29

adams_gen_anxiety_score
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 5.344828 4.064444 29
A Week_4 4.172414 3.882828 29
A Week_8 3.551724 3.322931 29
A Week_12 3.655172 3.096955 29
B Baseline 4.896552 4.237122 29
B Week_4 4.103448 3.608965 29
B Week_8 3.517241 3.112425 29
B Week_12 3.931034 3.575023 29

##VAS table and graph

caregiver_vas_total
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 230.9310 70.15184 29
A Week_4 211.5517 75.60356 29
A Week_8 208.4138 74.61162 29
A Week_12 208.2759 74.80541 29
B Baseline 225.9310 67.95898 29
B Week_4 222.1379 68.39576 29
B Week_8 205.2759 79.89810 29
B Week_12 189.1034 87.67609 29

##Clinician VAS

clin_vas_1_speech_severe_pms_fu
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 86.55172 14.82658 29
A Week_4 82.41379 16.50869 29
A Week_8 78.44828 21.38522 29
A Week_12 75.68966 19.76307 29
B Baseline 85.51724 19.79109 29
B Week_4 79.93103 21.81404 29
B Week_8 75.96552 20.85406 29
B Week_12 70.17241 23.62213 29
clin_vas_2_think_severe_pms_fu
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 81.20690 16.18284 29
A Week_4 75.17241 15.95097 29
A Week_8 71.72414 18.09084 29
A Week_12 66.72414 19.83305 29
B Baseline 81.03448 17.74720 29
B Week_4 75.17241 18.77892 29
B Week_8 71.72414 17.99186 29
B Week_12 64.65517 19.81752 29
clin_vas_3_seizures_severe_pms_fu
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 7.241379 16.88121 29
A Week_4 4.827586 14.54566 29
A Week_8 3.103448 11.68132 29
A Week_12 3.793103 12.07530 29
B Baseline 8.620690 19.40621 29
B Week_4 8.620690 19.40621 29
B Week_8 8.620690 19.40621 29
B Week_12 8.620690 19.40621 29
clin_vas_4_gross_severe_pms_fu
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 38.96552 22.84964 29
A Week_4 37.24138 21.53014 29
A Week_8 35.68966 21.11720 29
A Week_12 35.00000 24.53132 29
B Baseline 38.96552 26.53769 29
B Week_4 36.65517 23.92291 29
B Week_8 34.65517 22.59626 29
B Week_12 32.06897 21.52728 29
clin_vas_5_restrict_severe_pms_fu
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 71.55172 28.31908 29
A Week_4 68.44828 29.00442 29
A Week_8 66.89655 29.56270 29
A Week_12 67.58621 30.66585 29
B Baseline 71.20690 28.64983 29
B Week_4 68.62069 29.54707 29
B Week_8 64.48276 30.21781 29
B Week_12 61.20690 29.56999 29
clin_vas_6_soc_comm_severe_pms_fu
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 73.10345 23.42975 29
A Week_4 70.00000 23.52810 29
A Week_8 68.62069 24.60025 29
A Week_12 65.75862 25.67469 29
B Baseline 76.03448 23.99328 29
B Week_4 70.00000 22.63846 29
B Week_8 65.51724 23.88009 29
B Week_12 63.44828 23.53202 29
clin_vas_7_sensory_severe_pms_fu
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 65.24138 30.59371 29
A Week_4 66.72414 30.33373 29
A Week_8 63.96552 31.00632 29
A Week_12 63.79310 31.47151 29
B Baseline 67.24138 32.41985 29
B Week_4 65.51724 33.81528 29
B Week_8 63.44828 33.30335 29
B Week_12 60.17241 32.96326 29
clin_vas_8_daily_act_severe_pms_fu
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 71.58621 28.14696 29
A Week_4 68.48276 28.46630 29
A Week_8 66.72414 28.79347 29
A Week_12 63.62069 28.65735 29
B Baseline 70.44828 27.94329 29
B Week_4 63.10345 31.00731 29
B Week_8 63.62069 30.49873 29
B Week_12 59.48276 30.97651 29
clin_vas_9_sleep_severe_pms_fu
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 35.17241 30.13312 29
A Week_4 31.03448 30.13005 29
A Week_8 33.10345 31.03609 29
A Week_12 30.86207 29.12623 29
B Baseline 31.20690 28.71209 29
B Week_4 31.55172 25.18161 29
B Week_8 28.44828 25.46368 29
B Week_12 26.37931 24.60025 29

##Sensory Profile

sp_care_quad_1_low_reg
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 46.17241 9.878074 29
A Week_4 48.65517 10.958610 29
A Week_8 51.20690 8.821157 29
A Week_12 48.93103 11.423444 29
B Baseline 48.03448 11.718370 29
B Week_4 52.13793 12.322233 29
B Week_8 51.28571 10.309496 29
B Week_12 50.93103 9.938852 29

sp_care_quad_2_sen_seek
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 91.82759 13.98486 29
A Week_4 95.37931 15.58941 29
A Week_8 96.79310 13.36675 29
A Week_12 93.68966 16.01585 29
B Baseline 94.10345 18.24896 29
B Week_4 97.37931 15.35395 29
B Week_8 97.21429 16.46176 29
B Week_12 97.55172 15.97674 29

sp_care_quad_3_sen_sen
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 77.72414 11.62785 29
A Week_4 80.03448 11.59582 29
A Week_8 81.68966 10.79717 29
A Week_12 79.41379 11.51557 29
B Baseline 80.06897 11.87113 29
B Week_4 80.24138 12.36832 29
B Week_8 80.32143 10.47414 29
B Week_12 81.24138 11.96619 29

sp_care_quad_4_sen_avoid
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 112.5517 16.47202 29
A Week_4 115.1379 14.94448 29
A Week_8 115.7586 14.86043 29
A Week_12 113.6207 16.24899 29
B Baseline 114.8276 14.25550 29
B Week_4 116.1034 16.47584 29
B Week_8 116.0714 15.53712 29
B Week_12 118.1034 13.76939 29

sp_care_total
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 453.4483 53.84341 29
A Week_4 469.1034 54.24241 29
A Week_8 475.8276 51.91894 29
A Week_12 462.9655 59.61153 29
B Baseline 465.3103 61.77847 29
B Week_4 476.6552 62.15032 29
B Week_8 475.7143 60.87805 29
B Week_12 478.4483 56.98282 29

#vineland GSVs

vin3_recept_growth_scale
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 78.89655 24.04511 29
A Week_12 84.58621 19.49087 29
B Baseline 82.58621 22.85125 29
B Week_12 85.65517 20.70072 29

vin3_express_growth_scale
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 92.82759 28.07524 29
A Week_12 95.03448 29.53507 29
B Baseline 94.55172 29.89217 29
B Week_12 96.41379 30.28026 29

vin3_written_growth_scale
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 36.51724 26.74966 29
A Week_12 39.51724 26.71359 29
B Baseline 39.89286 27.06007 29
B Week_12 41.17857 27.03092 29

vin3_personal_growth_scale
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 92.10345 16.53426 29
A Week_12 94.13793 17.12918 29
B Baseline 92.06897 17.60871 29
B Week_12 94.51724 17.15813 29

vin3_domestic_growth_scale
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 18.68966 12.80654 29
A Week_12 18.58621 13.10539 29
B Baseline 18.21429 12.85923 29
B Week_12 19.85714 13.39904 29

vin3_commun_growth_scale
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 25.55172 16.09034 29
A Week_12 27.20690 16.17401 29
B Baseline 27.07143 15.15580 29
B Week_12 29.82143 15.18784 29

vin3_interper_growth_scale
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 67.44828 11.88512 29
A Week_12 67.96552 11.98357 29
B Baseline 67.75862 12.51472 29
B Week_12 69.51724 12.30163 29

vin3_playleis_growth_scale
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 60.44828 11.77948 29
A Week_12 59.68966 14.92626 29
B Baseline 58.24138 13.71822 29
B Week_12 60.44828 17.12889 29

vin3_coping_growth_scale
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 45.93103 5.762756 29
A Week_12 47.24138 6.378844 29
B Baseline 47.72414 6.192718 29
B Week_12 48.00000 5.700877 29

vin3_grossmot_growth_scale
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 135.75 10.75015 29
A Week_12 138.15 12.49115 29
B Baseline 136.15 12.28296 29
B Week_12 138.45 11.17080 29

vin3_finemot_growth_scale
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 101.45 19.36757 29
A Week_12 102.85 20.30109 29
B Baseline 103.20 17.44646 29
B Week_12 103.80 19.82449 29

#SAND

sand_total
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 37.00000 10.67708 29
A Week_12 35.06897 10.58277 29
B Baseline 36.62069 11.73339 29
B Week_12 32.51724 11.45993 29
sand_hyperreactivity
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 3.448276 4.289407 29
A Week_12 2.620690 3.288658 29
B Baseline 3.448276 3.679928 29
B Week_12 2.517241 3.019232 29
sand_hyporeactivity
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 16.51724 8.279151 29
A Week_12 15.17241 6.990141 29
B Baseline 15.62069 7.153490 29
B Week_12 13.34483 7.719529 29
sand_seeking
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 17.03448 5.615557 29
A Week_12 17.27586 6.250123 29
B Baseline 17.55172 7.174290 29
B Week_12 16.65517 6.096110 29

#Mullen

mull_grossmot_raw_score
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 29.51724 4.469106 29
A Week_12 29.41379 4.395755 29
B Baseline 29.27586 4.905511 29
B Week_12 29.62069 4.655364 29

mull_vis_recept_raw_score
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 27.89655 12.00103 29
A Week_12 28.96552 11.87579 29
B Baseline 29.31034 11.49877 29
B Week_12 29.48276 12.15560 29

mull_finemot_raw_score
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 25.62069 10.75843 29
A Week_12 25.68966 11.33296 29
B Baseline 26.17241 11.04235 29
B Week_12 26.41379 11.29323 29

mull_recept_lang_raw_score
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 21.68966 11.18386 29
A Week_12 21.41379 11.59285 29
B Baseline 22.27586 11.12557 29
B Week_12 22.89655 11.74231 29

mull_expres_lang_raw_score
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 16.93103 11.04190 29
A Week_12 16.65517 10.96513 29
B Baseline 17.62069 11.12081 29
B Week_12 17.89655 11.83320 29

mull_expres_lang_raw_score
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 16.93103 11.04190 29
A Week_12 16.65517 10.96513 29
B Baseline 17.62069 11.12081 29
B Week_12 17.89655 11.83320 29

#PPVT EVT

ppvt5_raw
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 35.80952 38.91223 29
A Week_12 38.72727 40.21156 29
B Baseline 34.60870 36.59948 29
B Week_12 43.40909 40.42001 29

ppvt5_gsv
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 438.6667 24.14608 29
A Week_12 438.0455 26.29227 29
B Baseline 439.0435 23.43943 29
B Week_12 444.0000 24.14243 29 

## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: formula
##    Data: data
## 
## REML criterion at convergence: 680.1
## 
## Scaled residuals: 
##      Min       1Q   Median       3Q      Max 
## -2.67318 -0.43958  0.05167  0.49670  1.56963 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 549.13   23.434  
##  Residual               60.28    7.764  
## Number of obs: 88, groups:  record_id, 26
## 
## Fixed effects:
##                    Estimate Std. Error       df t value Pr(>|t|)    
## (Intercept)        438.1703     4.9252  30.4921  88.966   <2e-16 ***
## EVENTWeek_12        -0.6888     2.4231  59.2734  -0.284    0.777    
## DRUGB               -1.2310     2.4147  59.3912  -0.510    0.612    
## EVENTWeek_12:DRUGB   3.0132     3.3847  59.3441   0.890    0.377    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.253                  
## DRUGB       -0.259  0.518           
## EVENTW_12:D  0.177 -0.715     -0.697
evt3_raw
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 25.57143 29.51066 29
A Week_12 28.92857 33.06279 29
B Baseline 32.92308 29.13892 29
B Week_12 32.35714 30.60319 29

evt3_gsv
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 434.2857 28.33163 29
A Week_12 435.1429 31.90336 29
B Baseline 442.6154 26.54411 29
B Week_12 442.1538 29.68683 29

#Both Caregiver Strain Index

cgsq_os
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 2.711599 1.149167 29
A Week_4 2.633229 1.094011 29
A Week_8 2.520376 1.127863 29
A Week_12 2.661442 1.027341 29
B Baseline 2.840125 1.127015 29
B Week_4 2.664577 1.108384 29
B Week_8 2.652038 1.131566 29
B Week_12 2.721003 1.112325 29

cgsq_ses
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 1.939655 0.9651406 29
A Week_4 2.103448 0.9601833 29
A Week_8 1.991379 1.0077440 29
A Week_12 2.034483 0.9653001 29
B Baseline 2.215517 0.9489753 29
B Week_4 2.120690 0.8676237 29
B Week_8 2.163793 0.9240726 29
B Week_12 2.051724 0.9691199 29

cgsq_sis
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 3.741379 0.9683342 29
A Week_4 3.442529 1.0401426 29
A Week_8 3.362069 1.1038074 29
A Week_12 3.517241 1.0730694 29
B Baseline 3.660919 0.9589623 29
B Week_4 3.465517 1.0795532 29
B Week_8 3.390805 1.1929654 29
B Week_12 3.459770 1.0261029 29

cgsq_global
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 8.392633 2.667625 29
A Week_4 8.179206 2.664118 29
A Week_8 7.873824 2.880652 29
A Week_12 8.213166 2.719040 29
B Baseline 8.716562 2.675804 29
B Week_4 8.250784 2.705640 29
B Week_8 8.206635 2.860244 29
B Week_12 8.232497 2.758848 29

csi_total
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 7.310345 3.230474 29
A Week_4 6.535714 3.756300 29
A Week_8 6.724138 3.788107 29
A Week_12 7.103448 3.735399 29
B Baseline 7.827586 3.665174 29
B Week_4 7.137931 3.925413 29
B Week_8 6.758621 3.942943 29
B Week_12 7.310345 3.606576 29

Children Sleep Habit

summarize_variable (gh, "cshq_total_score")
## `summarise()` has grouped output by 'DRUG'. You can override using the
## `.groups` argument.
cshq_total_score
DRUG EVENT Mean Standard Deviation Sample Size
A Baseline 2.026820 0.4979280 29
A Week_4 2.019157 0.4884639 29
A Week_8 2.038314 0.4772384 29
A Week_12 1.948276 0.4416906 29
B Baseline 2.086207 0.4355727 29
B Week_4 2.053640 0.5279830 29
B Week_8 1.946360 0.4606382 29
B Week_12 1.896552 0.4704036 29 
plot_variable_by_event (gh, cshq_total_score)+ ggtitle(" CSHQ Total") +ylim (0,3)

#starting the hunt #this looks at people who reported properly on the ABCSW - were there differences in these people? 

## 
## Call:
## lm(formula = abc_exact_age_months ~ SW_right, data = ghin_personBL1)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -43.46 -27.32 -10.39  25.80  67.51 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    99.19       7.86  12.620  4.5e-13 ***
## SW_right      -23.67      12.98  -1.823    0.079 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 34.26 on 28 degrees of freedom
## Multiple R-squared:  0.1061, Adjusted R-squared:  0.0742 
## F-statistic: 3.324 on 1 and 28 DF,  p-value: 0.07896

## 
## Call:
## lm(formula = cognitive_fsiq_fsdq_r01_gh ~ SW_right, data = ghin_personBL1)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -26.182  -9.304  -1.232  12.620  34.918 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   26.079      3.460   7.537 3.29e-08 ***
## SW_right1     12.003      5.714   2.100   0.0448 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 15.08 on 28 degrees of freedom
## Multiple R-squared:  0.1361, Adjusted R-squared:  0.1053 
## F-statistic: 4.412 on 1 and 28 DF,  p-value: 0.04482

## 
## Call:
## lm(formula = cdi_vocab_produced ~ SW_right, data = ghin_personBL1)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -103.73  -52.37  -48.87   46.52  312.63 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)  
## (Intercept)    52.37      26.18   2.001   0.0552 .
## SW_right1      51.36      43.23   1.188   0.2448  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 114.1 on 28 degrees of freedom
## Multiple R-squared:  0.048,  Adjusted R-squared:  0.014 
## F-statistic: 1.412 on 1 and 28 DF,  p-value: 0.2448
## 
## Call:
## lm(formula = PNAMP ~ SW_right, data = ghin_personBL1)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -4.394 -2.333 -0.635  1.387  9.650 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   3.2900     0.7592   4.334 0.000171 ***
## SW_right1     2.0036     1.2537   1.598 0.121242    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.309 on 28 degrees of freedom
## Multiple R-squared:  0.08359,    Adjusted R-squared:  0.05086 
## F-statistic: 2.554 on 1 and 28 DF,  p-value: 0.1212

## 
## Call:
## lm(formula = PNAMPshort ~ SW_right, data = ghin_personBL1)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -5.414 -3.291 -1.657  2.679 11.433 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    5.517      1.037   5.321 1.63e-05 ***
## SW_right1      2.077      1.704   1.219    0.234    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 4.275 on 25 degrees of freedom
##   (3 observations deleted due to missingness)
## Multiple R-squared:  0.0561, Adjusted R-squared:  0.01834 
## F-statistic: 1.486 on 1 and 25 DF,  p-value: 0.2342

## 
## Call:
## lm(formula = vep_short_msc_band_1 ~ SW_right, data = ghin_personBL1)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -0.21682 -0.04471  0.00540  0.02918  0.24018 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  0.23882    0.02309  10.343 1.62e-10 ***
## SW_right1    0.11378    0.03794   2.999  0.00606 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.09521 on 25 degrees of freedom
##   (3 observations deleted due to missingness)
## Multiple R-squared:  0.2645, Adjusted R-squared:  0.2351 
## F-statistic: 8.992 on 1 and 25 DF,  p-value: 0.006058

## 
## Call:
## lm(formula = clin_labs_igf1z1 ~ SW_right, data = ghin_personW12DRUG)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -4.1178 -1.0965  0.1322  1.0476  3.1522 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   1.6178     0.3774   4.287  0.00022 ***
## SW_right1     0.9662     0.6315   1.530  0.13806    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.601 on 26 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.08261,    Adjusted R-squared:  0.04733 
## F-statistic: 2.341 on 1 and 26 DF,  p-value: 0.1381

## 
## Call:
## lm(formula = clin_labs_igf1z1 ~ SW_right, data = ghin_personBL1)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.7068 -0.9160  0.0932  0.3932  5.3840 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.09316    0.33010  -0.282    0.780
## SW_right1    0.30916    0.56215   0.550    0.587
## 
## Residual standard error: 1.439 on 27 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.01108,    Adjusted R-squared:  -0.02555 
## F-statistic: 0.3025 on 1 and 27 DF,  p-value: 0.5869

## 
## Call:
## lm(formula = abc_subscale2_lethargy_2 ~ SW_right, data = ghin_personBL1)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -16.0000  -6.9605  -0.6053   5.9474  28.0000 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   15.211      2.389   6.366 6.88e-07 ***
## SW_right1      2.789      3.946   0.707    0.485    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 10.41 on 28 degrees of freedom
## Multiple R-squared:  0.01754,    Adjusted R-squared:  -0.01755 
## F-statistic: 0.4998 on 1 and 28 DF,  p-value: 0.4854

##    
##     iASD PMS
##   0    9  10
##   1    6   5
## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  tab
## X-squared = 0, df = 1, p-value = 1
##    
##      1  2
##   0  3 16
##   1  5  6
## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  tab2
## X-squared = 1.8016, df = 1, p-value = 0.1795

#what about people who repsonded on the VEP, were there differences in these people? 

## 
## Call:
## lm(formula = abc_exact_age_months ~ VEP_Responder, data = ghin_personDRUGBL)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -54.70 -23.92  -4.56  24.40  62.34 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)      86.764      9.818   8.837 5.19e-09 ***
## VEP_Responder1    7.701     13.380   0.576     0.57    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 34.01 on 24 degrees of freedom
## Multiple R-squared:  0.01362,    Adjusted R-squared:  -0.02748 
## F-statistic: 0.3313 on 1 and 24 DF,  p-value: 0.5703

## 
## Call:
## lm(formula = cognitive_fsiq_fsdq_r01_gh ~ VEP_Responder, data = ghin_personDRUGBL)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -13.85 -11.65  -3.10  11.25  21.25 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)      31.200      5.743   5.432 0.000206 ***
## VEP_Responder1   -5.450      7.321  -0.744 0.472238    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 12.84 on 11 degrees of freedom
##   (13 observations deleted due to missingness)
## Multiple R-squared:  0.04796,    Adjusted R-squared:  -0.03859 
## F-statistic: 0.5541 on 1 and 11 DF,  p-value: 0.4722

## 
## Call:
## lm(formula = cdi_vocab_produced ~ VEP_Responder, data = ghin_personDRUGBL)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -82.79 -79.28 -75.77  50.75 299.25 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)  
## (Intercept)      77.750     36.476   2.132   0.0435 *
## VEP_Responder1    5.036     49.708   0.101   0.9201  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 126.4 on 24 degrees of freedom
## Multiple R-squared:  0.0004274,  Adjusted R-squared:  -0.04122 
## F-statistic: 0.01026 on 1 and 24 DF,  p-value: 0.9201

## 
## Call:
## lm(formula = PNAMP ~ VEP_Responder, data = ghin_personDRUGBL)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.7025 -2.0690 -0.3079  0.9387  6.7421 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)      4.0225     0.8115   4.957 4.64e-05 ***
## VEP_Responder1  -0.5446     1.1058  -0.493    0.627    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 2.811 on 24 degrees of freedom
## Multiple R-squared:  0.01001,    Adjusted R-squared:  -0.03124 
## F-statistic: 0.2426 on 1 and 24 DF,  p-value: 0.6268

## 
## Call:
## lm(formula = vep_short_msc_band_1 ~ VEP_Responder, data = ghin_personDRUGBL)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -0.24423 -0.04523 -0.00950  0.02477  0.39977 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)     0.21750    0.03828   5.682 8.75e-06 ***
## VEP_Responder1  0.04873    0.05309   0.918    0.368    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.1326 on 23 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.03534,    Adjusted R-squared:  -0.006602 
## F-statistic: 0.8426 on 1 and 23 DF,  p-value: 0.3682

## 
## Call:
## lm(formula = clin_labs_igf1z1 ~ VEP_Responder, data = ghin_personDRUGBL)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.4343 -0.9343 -0.1118  0.7657  4.0882 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)
## (Intercept)      0.1118     0.4134   0.270    0.789
## VEP_Responder1  -0.4775     0.5524  -0.864    0.396
## 
## Residual standard error: 1.371 on 23 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.03146,    Adjusted R-squared:  -0.01065 
## F-statistic: 0.7472 on 1 and 23 DF,  p-value: 0.3963

##    
##     iASD PMS
##   0    5   7
##   1    8   6
## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  tab11
## X-squared = 0.15476, df = 1, p-value = 0.694
##    
##      1  2
##   0  4  8
##   1  3 11
## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  tab21
## X-squared = 0.057018, df = 1, p-value = 0.8113