R Markdown

##This is now both phases

##So far we have: MSC Band 2 for Long VEP, SAND Results, CGI-I Results, ABS, RBS-R, ADAMS, Caregiver VAS, Sensory Profile

For now - used last observation carried forward for missing data

##IGF1 ##This is Baseline levels only to look for carryover of IGF1

IGF1
GROUP redcap_event_name DRUG MEAN_igf1 SD_igf1 N_igf1
A baseline_phase_1_arm_1 A 135.6667 114.19574 15
A baseline_phase_2_arm_1 B 143.0667 100.32196 15
B baseline_phase_1_arm_1 B 119.4286 63.57759 15
B baseline_phase_2_arm_1 A 123.4167 66.92527 14

##Based on this, it doesn’t look like overall there was carryover of IGF1 ##mean values at baseline 1 and baseline 2 are similar - i think we should z score this but not sure the population levels that we used before

##IGF1 model ##1) look for carryover

## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: IGF1 ~ EVENT + DRUG + DRUG * EVENT + GROUP * DRUG + (1 | record_id)
##    Data: ghin_person
## 
## REML criterion at convergence: 1207.2
## 
## Scaled residuals: 
##      Min       1Q   Median       3Q      Max 
## -1.98331 -0.49626 -0.07543  0.40580  2.58199 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 5965     77.24   
##  Residual              2082     45.63   
## Number of obs: 113, groups:  record_id, 29
## 
## Fixed effects:
##                    Estimate Std. Error      df t value Pr(>|t|)    
## (Intercept)         133.936     22.467  36.373   5.962 7.53e-07 ***
## EVENTWeek_12         -1.139     12.275  80.290  -0.093    0.926    
## DRUGB                -6.043     14.597  80.150  -0.414    0.680    
## GROUPB               -1.294     31.228  31.848  -0.041    0.967    
## EVENTWeek_12:DRUGB   97.839     17.251  80.244   5.671 2.16e-07 ***
## DRUGB:GROUPB          9.086     17.270  80.247   0.526    0.600    
## ---
## 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.273                                    
## DRUGB       -0.326  0.420                             
## GROUPB      -0.661 -0.019      0.144                  
## EVENTW_12:D  0.194 -0.712     -0.578  0.014           
## DRUGB:GROUP  0.169  0.035     -0.542 -0.281 -0.036

##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-2, P60-N75 Difference

EVENT DRUG MEAN_msc_band_1 SD_msc_band_1 NlongMSC1
Baseline A 0.3869310 0.1595297 29
Week_12 A 0.3899655 0.1751706 29
Baseline B 0.3598929 0.1495443 30
Week_12 B 0.3799200 0.1605194 29
EVENT DRUG MEAN_msc_band_2 SD_msc_band_2 N
Baseline A 0.3574828 0.1866616 29
Week_12 A 0.3755517 0.2064568 29
Baseline B 0.3283571 0.1882139 30
Week_12 B 0.3754400 0.1589497 29
EVENT DRUG MEAN_msc_band_3 SD_msc_band_3 NlongMSC3
Baseline A 0.2172414 0.1538758 29
Week_12 A 0.1824138 0.1555469 29
Baseline B 0.1956429 0.1340266 30
Week_12 B 0.2046400 0.1275967 29
EVENT DRUG MEAN_msc_band_4 SD_msc_band_4 NlongMSC4
Baseline A 0.1196552 0.0504446 29
Week_12 A 0.1465517 0.0765980 29
Baseline B 0.1283214 0.0680626 30
Week_12 B 0.1464800 0.0723776 29
EVENT DRUG MEAN_P60_N75 SD_P60_N75 NP60_N75
Baseline A 4.157586 3.487043 29
Week_12 A 4.317931 4.075812 29
Baseline B 3.590714 2.722051 30
Week_12 B 4.924000 4.659851 29

Plots for VEP

All patients together

iASD and PMS separate

Modeling for VEP

## 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
## 
## REML criterion at convergence: -100.9
## 
## Scaled residuals: 
##      Min       1Q   Median       3Q      Max 
## -1.72815 -0.49735 -0.07145  0.64216  2.08483 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 0.02330  0.1527  
##  Residual              0.01129  0.1062  
## Number of obs: 111, groups:  record_id, 30
## 
## Fixed effects:
##                    Estimate Std. Error       df t value Pr(>|t|)    
## (Intercept)         0.36233    0.03422 49.13318  10.587 2.83e-14 ***
## EVENTWeek_12        0.01807    0.02790 78.12721   0.648    0.519    
## DRUGB              -0.03934    0.02846 78.90107  -1.382    0.171    
## EVENTWeek_12:DRUGB  0.01908    0.04058 78.33697   0.470    0.640    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.408                  
## DRUGB       -0.406  0.490           
## EVENTW_12:D  0.280 -0.687     -0.685
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: vep_long_msc_band_2 ~ EVENT + DRUG + DRUG * EVENT + GROUP + (1 |  
##     record_id)
##    Data: ghin_person
## 
## REML criterion at convergence: -97.1
## 
## Scaled residuals: 
##      Min       1Q   Median       3Q      Max 
## -1.72947 -0.49515 -0.05799  0.64552  2.06891 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 0.02419  0.1555  
##  Residual              0.01129  0.1062  
## Number of obs: 111, groups:  record_id, 30
## 
## Fixed effects:
##                    Estimate Std. Error       df t value Pr(>|t|)    
## (Intercept)         0.35428    0.04584 36.94267   7.729 3.13e-09 ***
## EVENTWeek_12        0.01807    0.02790 78.14316   0.648    0.519    
## DRUGB              -0.03958    0.02846 78.83357  -1.391    0.168    
## GROUPB              0.01628    0.06042 28.10272   0.269    0.790    
## EVENTWeek_12:DRUGB  0.01906    0.04058 78.34507   0.470    0.640    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB  GROUPB
## EVENTWek_12 -0.304                         
## DRUGB       -0.288  0.490                  
## GROUPB      -0.655  0.000     -0.022       
## EVENTW_12:D  0.207 -0.687     -0.685  0.003
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: vep_long_msc_band_2 ~ DRUG + EVENT + DRUG * EVENT + diag_dx___1 +  
##     (1 | record_id)
##    Data: ghin_person
## 
## REML criterion at convergence: -97.1
## 
## Scaled residuals: 
##      Min       1Q   Median       3Q      Max 
## -1.73979 -0.49360 -0.06645  0.63904  2.07831 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 0.02425  0.1557  
##  Residual              0.01129  0.1062  
## Number of obs: 111, groups:  record_id, 30
## 
## Fixed effects:
##                     Estimate Std. Error        df t value Pr(>|t|)    
## (Intercept)         0.360955   0.045880 36.899960   7.867 2.09e-09 ***
## DRUGB              -0.039448   0.028466 78.811279  -1.386    0.170    
## EVENTWeek_12        0.018069   0.027900 78.121652   0.648    0.519    
## diag_dx___1PMS      0.002838   0.060488 28.081288   0.047    0.963    
## DRUGB:EVENTWeek_12  0.019029   0.040583 78.323282   0.469    0.640    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) DRUGB  EVENTW d____1
## DRUGB       -0.288                     
## EVENTWek_12 -0.304  0.490              
## dg_d___1PMS -0.655 -0.022  0.000       
## DRUGB:EVENT  0.207 -0.685 -0.687  0.003
## 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: ghin_person
## 
## REML criterion at convergence: 545.7
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -2.9786 -0.4760 -0.1578  0.3201  4.1551 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 8.836    2.973   
##  Residual              5.268    2.295   
## Number of obs: 111, groups:  record_id, 30
## 
## Fixed effects:
##                    Estimate Std. Error      df t value Pr(>|t|)    
## (Intercept)          3.9048     0.9249 43.1214   4.222 0.000123 ***
## EVENTWeek_12         0.1603     0.6028 77.6271   0.266 0.790934    
## DRUGB               -0.9392     0.7691 79.1378  -1.221 0.225631    
## GROUPB               0.5701     1.2456 35.9552   0.458 0.649949    
## EVENTWeek_12:DRUGB   0.9524     0.8767 77.8712   1.086 0.280662    
## DRUGB:GROUPB         0.5426     0.8947 79.4856   0.606 0.545964    
## ---
## 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.326                                    
## DRUGB       -0.375  0.392                             
## GROUPB      -0.664  0.000      0.184                  
## EVENTW_12:D  0.224 -0.688     -0.555 -0.001           
## DRUGB:GROUP  0.212  0.000     -0.601 -0.336  0.011
## 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
## 
## REML criterion at convergence: 550.4
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -3.0187 -0.4648 -0.1220  0.3472  4.2583 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 8.674    2.945   
##  Residual              5.225    2.286   
## Number of obs: 111, groups:  record_id, 30
## 
## Fixed effects:
##                    Estimate Std. Error      df t value Pr(>|t|)    
## (Intercept)          4.1942     0.6867 53.9258   6.107 1.15e-07 ***
## EVENTWeek_12         0.1603     0.6003 78.6101   0.267    0.790    
## DRUGB               -0.6483     0.6119 79.5006  -1.059    0.293    
## EVENTWeek_12:DRUGB   0.9449     0.8730 78.8585   1.082    0.282    
## ---
## 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.435  0.490           
## EVENTW_12:D  0.300 -0.688     -0.686
## 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: ghin_person
## 
## REML criterion at convergence: 546.3
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -3.0384 -0.4657 -0.1156  0.3706  4.2261 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 8.321    2.885   
##  Residual              5.230    2.287   
## Number of obs: 111, groups:  record_id, 30
## 
## Fixed effects:
##                    Estimate Std. Error      df t value Pr(>|t|)    
## (Intercept)          4.9863     0.8837 40.0184   5.642  1.5e-06 ***
## DRUGB               -0.6265     0.6124 79.3103  -1.023    0.309    
## EVENTWeek_12         0.1603     0.6006 78.4636   0.267    0.790    
## diag_dx___1PMS      -1.5988     1.1429 28.3945  -1.399    0.173    
## DRUGB:EVENTWeek_12   0.9427     0.8734 78.7197   1.079    0.284    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) DRUGB  EVENTW d____1
## DRUGB       -0.322                     
## EVENTWek_12 -0.340  0.490              
## dg_d___1PMS -0.641 -0.024  0.000       
## DRUGB:EVENT  0.231 -0.686 -0.688  0.003

##SAND Results

EVENT DRUG MEAN_sand_total SD_sand_total N_sand_total
Baseline A 36.58621 10.601608 29
Week_12 A 36.25926 9.448616 29
Baseline B 35.81481 11.662026 30
Week_12 B 32.10714 11.548094 29
EVENT DRUG MEAN_sand_hyper SD_sand_hyper N_sand_hyper
Baseline A 3.379310 4.091924 29
Week_12 A 2.740741 3.181589 29
Baseline B 3.555556 3.734798 30
Week_12 B 2.517241 3.019232 29
EVENT DRUG MEAN_sand_hypo SD_sand_hypo N_sand_hypo
Baseline A 16.31034 8.237993 29
Week_12 A 15.62963 6.575877 29
Baseline B 15.11111 7.078316 30
Week_12 B 13.03571 7.676236 29
EVENT DRUG MEAN_sand_seek SD_sand_seek N_sand_seek
Baseline A 16.89655 5.690499 29
Week_12 A 17.88889 6.027714 29
Baseline B 17.14815 7.193500 30
Week_12 B 16.53571 6.221124 29

##CGI Improvement next up

EVENT DRUG MEAN_cgii SD_cgii Ncgii
Baseline A 0.0000000 0.0000000 29
Week_4 A 0.7241379 0.7971625 29
Week_8 A 0.6206897 0.9788402 29
Week_12 A 0.8620690 0.9533452 29
Baseline B 0.0000000 0.0000000 30
Week_4 B 0.5172414 0.9863602 29
Week_8 B 1.0357143 0.9993384 29
Week_12 B 1.3333333 0.8770580 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: 230.6
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -4.8725 -0.2488  0.0250  0.2493  1.8758 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 0.03551  0.1884  
##  Residual              0.37800  0.6148  
## Number of obs: 115, groups:  record_id, 30
## 
## Fixed effects:
##                    Estimate Std. Error       df t value Pr(>|t|)    
## (Intercept)         0.06897    0.14658 96.63274   0.471   0.6391    
## EVENTWeek_12        0.86207    0.16146 82.57645   5.339 8.05e-07 ***
## DRUGB              -0.11559    0.19590 82.88805  -0.590   0.5568    
## GROUPB             -0.14314    0.17604 69.56975  -0.813   0.4189    
## EVENTWeek_12:DRUGB  0.46964    0.22975 83.58663   2.044   0.0441 *  
## DRUGB:GROUPB        0.23639    0.22967 83.57470   1.029   0.3063    
## ---
## 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.551                                    
## DRUGB       -0.666  0.412                             
## GROUPB      -0.580  0.000      0.364                  
## EVENTW_12:D  0.387 -0.703     -0.564  0.002           
## DRUGB:GROUP  0.374  0.000     -0.575 -0.648 -0.011
## 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: 228.4
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -4.9990 -0.1818  0.0198  0.2214  1.7695 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 0.03116  0.1765  
##  Residual              0.37811  0.6149  
## Number of obs: 115, groups:  record_id, 30
## 
## Fixed effects:
##                      Estimate Std. Error         df t value Pr(>|t|)    
## (Intercept)        -1.905e-16  1.188e-01  1.093e+02   0.000    1.000    
## EVENTWeek_12        8.621e-01  1.615e-01  8.361e+01   5.338  7.9e-07 ***
## DRUGB               4.146e-17  1.602e-01  8.440e+01   0.000    1.000    
## EVENTWeek_12:DRUGB  4.719e-01  2.297e-01  8.463e+01   2.054    0.043 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.680                  
## DRUGB       -0.687  0.504           
## EVENTW_12:D  0.479 -0.703     -0.697
## 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: 545.8
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -4.1772 -0.5033  0.0167  0.5245  2.2461 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 0.1148   0.3389  
##  Residual              0.5346   0.7311  
## Number of obs: 230, groups:  record_id, 30
## 
## Fixed effects:
##                      Estimate Std. Error         df t value Pr(>|t|)    
## (Intercept)         -0.007877   0.149510 184.862319  -0.053 0.958041    
## EVENTWeek_4          0.724138   0.192005 193.379430   3.771 0.000216 ***
## EVENTWeek_8          0.620690   0.192005 193.379430   3.233 0.001442 ** 
## EVENTWeek_12         0.862069   0.192005 193.379430   4.490 1.22e-05 ***
## DRUGB                0.007877   0.190641 194.475474   0.041 0.967085    
## EVENTWeek_4:DRUGB   -0.203951   0.270520 193.639546  -0.754 0.451811    
## EVENTWeek_8:DRUGB    0.409917   0.271911 194.186760   1.508 0.133298    
## EVENTWeek_12:DRUGB   0.466838   0.273319 194.377783   1.708 0.089227 .  
## ---
## 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.642                                                          
## EVENTWeek_8 -0.642  0.500                                                   
## EVENTWek_12 -0.642  0.500     0.500                                         
## DRUGB       -0.650  0.504     0.504     0.504                               
## EVENTW_4:DR  0.456 -0.710    -0.355    -0.355     -0.703                    
## EVENTW_8:DR  0.456 -0.353    -0.706    -0.353     -0.701  0.494             
## EVENTW_12:D  0.453 -0.351    -0.351    -0.702     -0.697  0.491     0.490

##no significant carry over ##CGI-I is different from BL to W12 in drug and placebo *NOT accounting for dropout yet though ##CGI-I not significantly different when we look at all timepoints

##ABC

ABC Irritability
EVENT DRUG MEAN_abc_irrit SD_abc_irrit N_abc_irrit
Baseline A 10.758621 9.660875 29
Week_4 A 9.068965 8.280489 29
Week_8 A 8.821429 6.869436 29
Week_12 A 8.344828 7.162780 29
Baseline B 11.100000 9.252586 30
Week_4 B 8.344828 8.844998 29
Week_8 B 7.689655 7.860859 29
Week_12 B 7.793103 7.451650 29
ABC Social Withdrawal
EVENT DRUG MEAN_abc_SW SD_abc_SW N_abc_SW
Baseline A 13.65517 8.469592 29
Week_4 A 12.37931 7.997383 29
Week_8 A 11.20690 7.902890 29
Week_12 A 10.65517 7.560593 29
Baseline B 14.76667 10.125704 30
Week_4 B 11.44828 8.862385 29
Week_8 B 10.20690 7.442059 29
Week_12 B 10.31034 7.714103 29
ABC Stereotypy
EVENT DRUG MEAN_abc_Stereo SD_abc_Stereo N_abc_Stereo
Baseline A 7.172414 5.251319 29
Week_4 A 5.551724 4.428691 29
Week_8 A 5.750000 4.282825 29
Week_12 A 5.413793 5.348680 29
Baseline B 6.866667 5.354556 30
Week_4 B 5.620690 4.828589 29
Week_8 B 4.793103 4.091623 29
Week_12 B 5.000000 4.590363 29
ABC Hyperactivity
EVENT DRUG MEAN_abc_Hyp SD_abc_Hyp N_abc_Hyp
Baseline A 23.44828 11.296606 29
Week_4 A 20.20690 10.366082 29
Week_8 A 19.06897 9.819554 29
Week_12 A 18.82759 9.914164 29
Baseline B 24.70000 9.523655 30
Week_4 B 18.96552 9.469812 29
Week_8 B 18.72414 9.617308 29
Week_12 B 18.68966 10.423543 29
ABC Inappropriate Speech
EVENT DRUG MEAN_abc_InSpeech SD_abc_InSpeech N_abc_InSpeech
Baseline A 3.034483 3.877432 29
Week_4 A 2.379310 3.004922 29
Week_8 A 2.448276 2.667077 29
Week_12 A 2.172414 2.633154 29
Baseline B 2.733333 3.194104 30
Week_4 B 2.103448 2.743106 29
Week_8 B 2.034483 2.757882 29
Week_12 B 1.931034 3.161499 29

#ABC graphs

##MCDI

MCDI Amended
EVENT DRUG MEAN_mcdi_amend SD_mcdi_amend N_mcdi_amend
Baseline A 75.35714 114.7536 29
Week_12 A 84.57143 121.3335 29
Baseline B 93.06667 132.6047 30
Week_12 B 117.17857 154.1099 29

##RBS-R

RBS-R Stereotyped Behavior
EVENT DRUG MEAN_rbs_stereo SD_rbs_stereo N__rbs_stereo
Baseline A 5.551724 4.005231 29
Week_4 A 4.310345 3.285286 29
Week_8 A 3.931034 3.283411 29
Week_12 A 3.571429 3.338092 29
Baseline B 4.933333 4.456172 30
Week_4 B 4.250000 3.502645 29
Week_8 B 3.896552 3.199831 29
Week_12 B 3.928571 3.547859 29
RBS-R Self Injury
EVENT DRUG MEAN_rbs_self SD_rbs_self N__rbs_self
Baseline A 2.827586 3.742316 29
Week_4 A 2.517241 3.323302 29
Week_8 A 2.068966 2.814897 29
Week_12 A 2.142857 3.051585 29
Baseline B 2.733333 3.912256 30
Week_4 B 1.500000 2.603417 29
Week_8 B 1.827586 3.117565 29
Week_12 B 2.000000 3.485419 29
RBS-R Compulsive Behavior
EVENT DRUG MEAN_rbs_comp SD_rbs_comp N__rbs_comp
Baseline A 3.413793 3.914417 29
Week_4 A 2.827586 3.129000 29
Week_8 A 1.965517 2.556263 29
Week_12 A 2.321429 3.067797 29
Baseline B 2.166667 3.484778 30
Week_4 B 1.607143 1.931019 29
Week_8 B 2.448276 2.810519 29
Week_12 B 2.035714 2.714599 29
RBS-R Ritualistic Behavior
EVENT DRUG MEAN_rbs_rit SD_rbs_rit N_rbs_rit
Baseline A 2.827586 3.306584 29
Week_4 A 2.206897 2.677676 29
Week_8 A 2.310345 3.263472 29
Week_12 A 2.107143 2.346617 29
Baseline B 2.400000 3.558380 30
Week_4 B 2.142857 2.704102 29
Week_8 B 2.241379 3.019640 29
Week_12 B 2.000000 2.638743 29
RBS-R Sameness Behavior
EVENT DRUG MEAN_rbs_same SD_rbs_same N_rbs_same
Baseline A 5.517241 6.127343 29
Week_4 A 4.862069 4.348433 29
Week_8 A 4.862069 4.381162 29
Week_12 A 4.214286 4.314979 29
Baseline B 4.933333 4.748381 30
Week_4 B 4.321429 4.861440 29
Week_8 B 4.206897 4.776275 29
Week_12 B 4.142857 5.140358 29
RBS-R Restricted Behavior
EVENT DRUG MEAN_rbs_resitrc SD_rbs_resitrc N_rbs_resitrc
Baseline A 3.689655 2.916743 29
Week_4 A 3.827586 2.928541 29
Week_8 A 3.000000 2.299068 29
Week_12 A 3.285714 2.565543 29
Baseline B 3.833333 2.937080 30
Week_4 B 3.535714 2.781981 29
Week_8 B 3.206897 2.304953 29
Week_12 B 3.321429 2.625425 29
RBS-R Total Score
EVENT DRUG MEAN_rbs_total SD_rbs_total N_rbs_total
Baseline A 23.82759 17.01275 29
Week_4 A 20.55172 14.77108 29
Week_8 A 18.13793 13.43430 29
Week_12 A 17.64286 13.52970 29
Baseline B 21.00000 17.85188 30
Week_4 B 17.35714 12.66729 29
Week_8 B 17.82759 14.27553 29
Week_12 B 17.42857 14.69550 29

##ADAMS

ADAMS Manic Score
EVENT DRUG MEAN_adams_manic SD_adams_manic N_adams_manic
Baseline A 8.620690 3.363818 29
Week_4 A 7.896552 3.885365 29
Week_8 A 6.965517 3.904967 29
Week_12 A 7.551724 4.330838 29
Baseline B 8.433333 3.997269 30
Week_4 B 7.758621 3.621909 29
Week_8 B 7.000000 3.644957 29
Week_12 B 6.827586 3.433750 29
ADAMS Depressed Score
EVENT DRUG MEAN_adams_dep SD_adams_dep N_adams_dep
Baseline A 3.896552 3.233141 29
Week_4 A 2.344828 2.224472 29
Week_8 A 2.172414 1.712746 29
Week_12 A 2.206897 2.144417 29
Baseline B 3.033333 3.408745 30
Week_4 B 2.758621 3.180916 29
Week_8 B 2.482759 2.293168 29
Week_12 B 2.275862 2.281862 29
ADAMS Avoid Score
EVENT DRUG MEAN_adams_avoid SD_adams_avoid N_adams_avoid
Baseline A 8.620690 4.813773 29
Week_4 A 7.655172 5.293364 29
Week_8 A 7.241379 5.068496 29
Week_12 A 7.413793 4.939037 29
Baseline B 8.500000 5.063391 30
Week_4 B 7.344828 4.915543 29
Week_8 B 6.620690 3.839771 29
Week_12 B 6.517241 4.954224 29
ADAMS Obsessive Score
EVENT DRUG MEAN_adams_pbs SD_adamsobs N_adams_obs
Baseline A 2.827586 2.315614 29
Week_4 A 2.000000 2.104417 29
Week_8 A 1.965517 1.917639 29
Week_12 A 1.896552 1.819435 29
Baseline B 2.333333 2.089932 30
Week_4 B 1.931034 1.998768 29
Week_8 B 2.172414 1.909917 29
Week_12 B 2.310345 2.019608 29
ADAMS Anxiety Score
EVENT DRUG MEAN_adams_anx SD_adams_anx N_adams_anx
Baseline A 5.344828 4.064444 29
Week_4 A 4.172414 3.882828 29
Week_8 A 3.586207 3.300619 29
Week_12 A 3.586207 3.156821 29
Baseline B 4.933333 4.168298 30
Week_4 B 4.103448 3.608965 29
Week_8 B 3.482759 3.135291 29
Week_12 B 3.655172 3.477231 29

#graphs

##VAS table and graph

Caregiver VAS Total
EVENT DRUG MEAN_vas SD_vas N_vas
Baseline A 239.1786 55.29788 29
Week_4 A 217.5000 69.73893 29
Week_8 A 216.1786 63.62393 29
Week_12 A 215.7143 64.33584 29
Baseline B 234.2414 53.85672 30
Week_4 B 222.5714 69.61025 29
Week_8 B 212.2500 71.24534 29
Week_12 B 187.3214 88.74852 29

Caregiver VAS stats 1) look for carryover

## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: vas_tot ~ EVENT + DRUG + DRUG * EVENT + GROUP * DRUG + (1 | record_id)
##    Data: ghin_person
## 
## REML criterion at convergence: 1178
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -3.3534 -0.5371  0.0603  0.5822  1.7068 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 2531     50.31   
##  Residual              1875     43.30   
## Number of obs: 113, groups:  record_id, 29
## 
## Fixed effects:
##                    Estimate Std. Error     df t value Pr(>|t|)    
## (Intercept)          241.62      16.77  44.93  14.408   <2e-16 ***
## EVENTWeek_12         -23.46      11.57  80.34  -2.028   0.0459 *  
## DRUGB                -24.93      14.16  80.37  -1.761   0.0821 .  
## GROUPB                -5.16      22.06  36.65  -0.234   0.8164    
## EVENTWeek_12:DRUGB   -24.51      16.31  80.42  -1.502   0.1370    
## DRUGB:GROUPB          39.09      16.43  81.05   2.380   0.0197 *  
## ---
## 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.345                                    
## DRUGB       -0.423  0.409                             
## GROUPB      -0.670  0.000      0.214                  
## EVENTW_12:D  0.245 -0.709     -0.563  0.000           
## DRUGB:GROUP  0.243  0.000     -0.583 -0.378 -0.011

##Does look like carryover effects here - would need to sort out before analyzing (either phase 1 only or something more complicated)

##Sensory Profile

Sensory Profile Low Registration
EVENT DRUG MEAN_sp1 SD_sp1 N_sp1
Baseline A 46.17241 9.878074 29
Week_4 A 48.65517 10.958610 29
Week_8 A 50.96552 8.837755 29
Week_12 A 48.93103 11.423444 29
Baseline B 48.13333 11.527279 30
Week_4 B 51.03448 12.918654 29
Week_8 B 51.24138 10.397991 29
Week_12 B 51.00000 10.053429 29
Sensory Profile Sensory Seeking
EVENT DRUG MEAN_sp2 SD_sp2 N_sp2
Baseline A 91.82759 13.98486 29
Week_4 A 95.37931 15.58941 29
Week_8 A 96.17241 13.58063 29
Week_12 A 93.68966 16.01585 29
Baseline B 93.96667 17.94721 30
Week_4 B 97.37931 15.35395 29
Week_8 B 97.10345 16.45198 29
Week_12 B 97.55172 15.97674 29
Sensory Profile Sensory Sensitivity
EVENT DRUG MEAN_sp3 SD_sp3 N_sp3
Baseline A 77.72414 11.62785 29
Week_4 A 80.03448 11.59582 29
Week_8 A 81.48276 10.66912 29
Week_12 A 79.41379 11.51557 29
Baseline B 79.46667 12.12218 30
Week_4 B 80.03448 12.62787 29
Week_8 B 80.41379 10.27659 29
Week_12 B 81.03448 11.82456 29
Sensory Profile Sensory Avoidance
EVENT DRUG MEAN_sp4 SD_sp4 N_sp4
Baseline A 112.5517 16.47202 29
Week_4 A 115.1379 14.94448 29
Week_8 A 115.5862 14.64464 29
Week_12 A 113.6207 16.24899 29
Baseline B 113.8333 15.02890 30
Week_4 B 115.9655 16.71715 29
Week_8 B 116.1379 15.66325 29
Week_12 B 118.2414 13.91211 29
Sensory Profile Total
EVENT DRUG MEAN_sptot SD_sptot N_sptot
Baseline A 453.4483 53.84341 29
Week_4 A 469.1034 54.24241 29
Week_8 A 474.3793 51.28242 29
Week_12 A 462.9655 59.61153 29
Baseline B 463.9000 61.19350 30
Week_4 B 474.5517 64.62838 29
Week_8 B 475.4483 61.52211 29
Week_12 B 478.6552 57.29203 29

##Sensory Profile Graphs

##Clinician VAS

Clinician VAS - Speech
EVENT DRUG MEAN_vas1 SD_vas1 N_vas1
Baseline A 86.55172 14.82658 29
Week_4 A 82.41379 16.50869 29
Week_8 A 79.65517 18.46478 29
Week_12 A 75.68966 19.76307 29
Baseline B 84.83333 19.80436 30
Week_4 B 80.64286 21.86866 29
Week_8 B 74.75862 23.28190 29
Week_12 B 70.17241 23.62213 29
Clinician VAS - Think
EVENT DRUG MEAN_vas2 SD_vas2 N_vas2
Baseline A 81.20690 16.18284 29
Week_4 A 75.17241 15.95097 29
Week_8 A 72.41379 16.98848 29
Week_12 A 66.72414 19.83305 29
Baseline B 79.50000 19.35824 30
Week_4 B 75.89286 18.71094 29
Week_8 B 71.03448 19.01030 29
Week_12 B 64.65517 19.81752 29
Clinician VAS - Seizure
EVENT DRUG MEAN_vas3 SD_vas3 N_vas3
Baseline A 7.241379 16.88121 29
Week_4 A 4.827586 14.54566 29
Week_8 A 3.103448 11.68132 29
Week_12 A 3.793103 12.07530 29
Baseline B 9.000000 19.18153 30
Week_4 B 8.928571 19.69006 29
Week_8 B 8.620690 19.40621 29
Week_12 B 8.620690 19.40621 29
Clinician VAS - Gross Motor
EVENT DRUG MEAN_vas4 SD_vas4 N_vas4
Baseline A 38.96552 22.84964 29
Week_4 A 37.24138 21.53014 29
Week_8 A 36.20690 21.02954 29
Week_12 A 35.00000 24.53132 29
Baseline B 38.83333 26.08617 30
Week_4 B 37.60714 23.79595 29
Week_8 B 34.13793 22.64118 29
Week_12 B 32.06897 21.52728 29
Clinician VAS - Restricted Behavior
EVENT DRUG MEAN_vas5 SD_vas5 N_vas5
Baseline A 71.55172 28.31908 29
Week_4 A 68.44828 29.00442 29
Week_8 A 67.24138 28.80309 29
Week_12 A 67.58621 30.66585 29
Baseline B 70.33333 28.55525 30
Week_4 B 68.92857 30.04186 29
Week_8 B 64.13793 30.91084 29
Week_12 B 61.20690 29.56999 29
Clinician VAS - Social Communication
EVENT DRUG MEAN_vas6 SD_vas6 N_vas6
Baseline A 73.10345 23.42975 29
Week_4 A 70.00000 23.52810 29
Week_8 A 68.62069 24.60025 29
Week_12 A 65.75862 25.67469 29
Baseline B 76.66667 23.82889 30
Week_4 B 69.82143 23.03307 29
Week_8 B 65.51724 23.88009 29
Week_12 B 63.44828 23.53202 29
Clinician VAS - Sensory
EVENT DRUG MEAN_vas7 SD_vas7 N_vas7
Baseline A 65.24138 30.59371 29
Week_4 A 66.72414 30.33373 29
Week_8 A 64.65517 30.67488 29
Week_12 A 63.79310 31.47151 29
Baseline B 67.33333 31.85997 30
Week_4 B 67.32143 32.98378 29
Week_8 B 62.75862 33.58322 29
Week_12 B 60.17241 32.96326 29
Clinician VAS - Daily Activites
EVENT DRUG MEAN_vas8 SD_vas8 N_vas8
Baseline A 71.58621 28.14696 29
Week_4 A 68.48276 28.46630 29
Week_8 A 66.72414 28.79347 29
Week_12 A 63.62069 28.65735 29
Baseline B 69.93333 27.60177 30
Week_4 B 65.00000 29.81424 29
Week_8 B 63.62069 30.49873 29
Week_12 B 59.48276 30.97651 29
Clinician VAS - Sleep
EVENT DRUG MEAN_vas9 SD_vas9 N_vas9
Baseline A 35.17241 30.13312 29
Week_4 A 31.03448 30.13005 29
Week_8 A 33.10345 31.03609 29
Week_12 A 30.86207 29.12623 29
Baseline B 30.33333 28.61557 30
Week_4 B 31.96429 25.54369 29
Week_8 B 28.44828 25.46368 29
Week_12 B 26.37931 24.60025 29

#Clinician VAS Graphs

#Clinician VAS social communication stats

vassoci <- lmer (clin_vas_6_soc_comm_severe_pms ~  EVENT + DRUG  + DRUG*EVENT + (1|record_id), data = ghin_person)
summary(vassoci)
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: clin_vas_6_soc_comm_severe_pms ~ EVENT + DRUG + DRUG * EVENT +  
##     (1 | record_id)
##    Data: ghin_person
## 
## REML criterion at convergence: 932
## 
## Scaled residuals: 
##      Min       1Q   Median       3Q      Max 
## -2.29200 -0.55426  0.01607  0.54275  2.32985 
## 
## Random effects:
##  Groups    Name        Variance Std.Dev.
##  record_id (Intercept) 486.85   22.065  
##  Residual               89.83    9.478  
## Number of obs: 117, groups:  record_id, 30
## 
## Fixed effects:
##                    Estimate Std. Error     df t value Pr(>|t|)    
## (Intercept)          73.405      4.401 37.424  16.677   <2e-16 ***
## EVENTWeek_12         -7.345      2.489 84.063  -2.951   0.0041 ** 
## DRUGB                 3.262      2.478 84.178   1.317   0.1916    
## EVENTWeek_12:DRUGB   -5.074      3.510 84.090  -1.446   0.1520    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) EVENTWk_12 DRUGB 
## EVENTWek_12 -0.283                  
## DRUGB       -0.288  0.502           
## EVENTW_12:D  0.200 -0.709     -0.700