Supplemental material 5: Prediction using bifactor-(S-1) model and matrix regression

knitr::opts_chunk$set(echo = TRUE, warning = TRUE, message = FALSE)

---

Preparation

load packages

options(scipen = 99999)

library(psych)
library(knitr)
library(dplyr)
library(tidyverse)
#library(QuantPsyc)
library(datarobot)
library(xtable)
library(reshape2)
library(naniar)
library(pander)
library(lavaan)
library(data.table)
library(corrplot)
library(GPArotation)
library(writexl)
library(car)
library(semTable)
library(ggplot2)
library(kableExtra)

Load correlation matrix

load("scale_cor_pred.rda")

Prediction of PSF and average affective well-being with S-1 Model with N8 as reference and orthogonal S factors

Bifactor-(S-1) model with PF (N8) as reference

s1wn8asref = "
# general factor
N8_PF_Problems =~  PF_affect_comm + PF_affect_diff + PF_affect_tolrnce +  PF_forming_relshps + PF_holistic_percptn + PF_identity + PF_impulse_reg +PF_intern_relshp_mod + PF_self_reflection + PF_self_est_reg + PF_anticipation +
SAD_soc_inter + SAD_avoidance + SAD_performing + SAD_devaluation + SAD_being_observ + SAD_eating + SAD_attractive_pers  + Spec_phob_avoid + Agoraphobia + DEP_hopelessn + DEP_depr_mood + DEP_dimin_selfconf + DEP_worthlessn + DEP_suicid_intent + DEP_dimin_interest + DEP_guilt + DEP_crying + DEP_decision_probl  + Panic_doctor_visits + Panic_attacks + Panic_phys_symp + Panic_dysf_thought + Panic_fear_recurr + Panic_avoid_behv + Somatic_symptoms + Hypochondriasis + GAD_free_float_anx + GAD_nervousness + GAD_irritability + GAD_excess_worry + DEP_rumination + DEP_psycmot_agit  + OCD_sympt + DI_perfectionism + EATA_restricted + EATA_omit_meal + EATA_time_frames + EATA_count_calories + EATA_purging_behv + EATA_exces_exercise + EATB_loss_ctrl + EATB_cope_negaff + EATB_overeating + EATB_intrus_thoughts + DEP_incr_appetite + EATB_soc_disguised + PS_unusual_beliefs + PS_percept_dysreg + Psychotic_symptoms + Dissociation + PS_eccentricity  + Manic_sympt    + NA_sep_insecur + DT_suspiciousness   + NA_perseveration  + DEP_psycmot_inhib + DEP_reduc_energy + DEP_hypersomnia + DEP_circad_mood_fluct + DEP_concentration + DEP_loss_appetite + DEP_loss_libido + DT_intimcy_avoid + DT_restr_affectivity + DT_withdrawal + DEP_emot_numbn + DEP_crying_inab + DT_anhedonia + AN_callousness + AN_hostility + AN_deceitfulness + AN_manipulativns + AN_entitlement + AN_grandiosity  + DI_impulsivity + Spec_phob_fear + Spec_phob_animals + Spec_phob_blood + EATA_body_selfest + EATB_worries + EATA_fear_weightg + EATA_fear_weightl + Subst_use_drug + DI_normviolation + Subst_use_alc + DI_risk_taking + DI_irresponsibility


#N9
N9_Somatic_Dep =~ DEP_psycmot_inhib + DEP_reduc_energy + DEP_hypersomnia + DEP_circad_mood_fluct + DEP_concentration + DEP_loss_appetite + DEP_loss_libido + DT_anhedonia + DEP_depr_mood + DEP_dimin_selfconf + DEP_dimin_interest + DEP_decision_probl + Somatic_symptoms


#N2
N2_Cognitive_Dep =~  DEP_hopelessn + DEP_depr_mood + DEP_dimin_selfconf + DEP_worthlessn + DEP_suicid_intent + DEP_dimin_interest + DEP_guilt + DEP_crying + DEP_decision_probl + PF_self_reflection + DEP_psycmot_agit + DEP_reduc_energy + DEP_concentration + DT_anhedonia


#N4
N4_GAD_OCD =~ GAD_free_float_anx + GAD_nervousness + GAD_irritability + GAD_excess_worry + DEP_psycmot_agit + OCD_sympt + DI_perfectionism + DEP_concentration + DEP_decision_probl + PF_self_reflection

#N1
N1_SAD_Agora =~ SAD_soc_inter + SAD_avoidance + SAD_performing + SAD_devaluation + SAD_being_observ + SAD_eating + SAD_attractive_pers + Spec_phob_avoid + Agoraphobia + DT_restr_affectivity + DT_withdrawal

#N3
N3_Panic =~ Panic_doctor_visits + Panic_attacks + Panic_phys_symp + Panic_dysf_thought + Panic_fear_recurr + Panic_avoid_behv + Somatic_symptoms + Hypochondriasis

#N12
N12_Spec_Phob =~ Spec_phob_fear + Spec_phob_animals + Spec_phob_blood + SAD_performing + Spec_phob_avoid

#N13
N13_ED_thoughts =~  EATA_body_selfest + EATB_worries + EATA_fear_weightg + EATA_fear_weightl + EATA_exces_exercise + EATB_soc_disguised

#N6
N6_ED_Bulimic =~ EATB_loss_ctrl + EATB_cope_negaff + EATB_overeating + EATB_intrus_thoughts + DEP_incr_appetite + EATB_soc_disguised + EATA_body_selfest + EATB_worries + EATA_fear_weightg

#N5
N5_ED_Anorectic =~ EATA_restricted + EATA_omit_meal + EATA_time_frames + EATA_count_calories + EATA_purging_behv + EATA_exces_exercise + EATB_soc_disguised + EATA_body_selfest + EATB_worries + EATA_fear_weightg

#N10
N10_Detachment =~ DT_intimcy_avoid + DT_restr_affectivity + DT_withdrawal + DEP_emot_numbn + DEP_crying_inab + DT_anhedonia + DEP_dimin_interest

#N7
N7_Thought_Disorder =~ PS_unusual_beliefs + PS_percept_dysreg + Psychotic_symptoms + Dissociation + Manic_sympt  + DEP_crying_inab + PF_self_reflection + OCD_sympt

#N11
N11_Antagonism =~  AN_callousness + AN_hostility + AN_deceitfulness + AN_manipulativns + AN_entitlement + AN_grandiosity + DI_irresponsibility


#N14
N14_Disinhibition =~  Subst_use_drug + DI_normviolation + Subst_use_alc + DI_risk_taking + DI_irresponsibility + Manic_sympt



#P Functioning as S factor
N8_PF_Problems ~~ 0*N5_ED_Anorectic
N8_PF_Problems ~~ 0*N1_SAD_Agora
N8_PF_Problems ~~ 0*N2_Cognitive_Dep
N8_PF_Problems ~~ 0*N3_Panic
N8_PF_Problems ~~ 0*N4_GAD_OCD
N8_PF_Problems ~~ 0*N6_ED_Bulimic
N8_PF_Problems ~~ 0*N7_Thought_Disorder
N8_PF_Problems ~~ 0*N9_Somatic_Dep
N8_PF_Problems ~~ 0*N10_Detachment
N8_PF_Problems ~~ 0*N11_Antagonism
N8_PF_Problems ~~ 0*N12_Spec_Phob
N8_PF_Problems ~~ 0*N13_ED_thoughts
N8_PF_Problems ~~ 0*N14_Disinhibition
"


# fit model
s1model =  cfa(model = s1wn8asref, sample.cov = scale_cov, sample.nobs = 27173)



# calculate unbiased SRMR and highest standardized residual covariance and average R2 for manifest variables according to Shi et al. 2018
test = lavResiduals(s1model)
save(test, file="resids1.rda")
uSRMR = data.frame(test$summary)
uSRMR$value = uSRMR$cov
names(uSRMR) = c("fitindex", "value")
uSRMR$fitindex = rownames(uSRMR)
uSRMR = uSRMR[5:8,]

kable(uSRMR)

	fitindex	value
usrmr	usrmr	0.0438474
usrmr.se	usrmr.se	0.0001920
usrmr.ci.lower	usrmr.ci.lower	0.0435316
usrmr.ci.upper	usrmr.ci.upper	0.0441633

kable(data.frame(max_std_cov = max(test$cov)))

max_std_cov
0.3745343

kable(data.frame(average_R2_manifest_vars_x_0.10 =mean(lavInspect(s1model, what = "rsquare")) * .1))

average_R2_manifest_vars_x_0.10
0.0555514

summary(s1model, std=T)

## lavaan 0.6.16 ended normally after 206 iterations
##
##   Estimator                                         ML
##   Optimization method                           NLMINB
##   Number of model parameters                       388
##
##   Number of observations                         27173
##
## Model Test User Model:
##
##   Test statistic                              426130.451
##   Degrees of freedom                                4463
##   P-value (Chi-square)                             0.000
##
## Parameter Estimates:
##
##   Standard errors                             Standard
##   Information                                 Expected
##   Information saturated (h1) model          Structured
##
## Latent Variables:
##                          Estimate  Std.Err  z-value  P(>|z|)   Std.lv  Std.all
##   N8_PF_Problems =~
##     PF_affect_comm          1.000                               0.906    0.749
##     PF_affect_diff          1.199    0.008  149.168    0.000    1.087    0.845
##     PF_affct_tlrnc          1.107    0.008  142.565    0.000    1.004    0.813
##     PF_frmng_rlshp          0.710    0.009   81.604    0.000    0.644    0.489
##     PF_hlstc_prcpt          0.723    0.007  102.046    0.000    0.655    0.603
##     PF_identity             1.052    0.007  157.995    0.000    0.954    0.887
##     PF_impulse_reg          0.736    0.007  103.385    0.000    0.667    0.610
##     PF_ntrn_rlshp_          0.923    0.008  121.996    0.000    0.836    0.709
##     PF_self_rflctn          1.173    0.008  144.717    0.000    1.063    0.825
##     PF_self_est_rg          0.826    0.007  120.067    0.000    0.748    0.699
##     PF_anticipatin          0.911    0.008  111.676    0.000    0.826    0.655
##     SAD_soc_inter           0.957    0.010   98.449    0.000    0.867    0.585
##     SAD_avoidance           1.046    0.010  107.321    0.000    0.948    0.633
##     SAD_performing          0.780    0.008   98.136    0.000    0.707    0.583
##     SAD_devaluatin          1.259    0.010  121.733    0.000    1.142    0.709
##     SAD_being_bsrv          0.914    0.009  102.362    0.000    0.829    0.606
##     SAD_eating              0.858    0.008  103.857    0.000    0.778    0.614
##     SAD_ttrctv_prs          0.780    0.010   78.868    0.000    0.707    0.474
##     Spec_phob_avod          0.949    0.010   94.504    0.000    0.860    0.563
##     Agoraphobia             0.709    0.008   93.625    0.000    0.643    0.557
##     DEP_hopelessn           0.862    0.008  105.814    0.000    0.781    0.626
##     DEP_depr_mood           0.857    0.007  121.145    0.000    0.777    0.707
##     DEP_dmn_slfcnf          0.874    0.008  113.870    0.000    0.792    0.669
##     DEP_worthlessn          1.035    0.009  117.696    0.000    0.938    0.689
##     DEP_suicd_ntnt          0.734    0.008   97.334    0.000    0.666    0.578
##     DEP_dimn_ntrst          0.900    0.008  114.940    0.000    0.816    0.674
##     DEP_guilt               1.034    0.009  119.525    0.000    0.937    0.698
##     DEP_crying              0.606    0.008   77.500    0.000    0.549    0.466
##     DEP_decsn_prbl          1.026    0.008  121.041    0.000    0.930    0.705
##     Panc_dctr_vsts          0.203    0.005   40.367    0.000    0.184    0.248
##     Panic_attacks           0.663    0.008   79.647    0.000    0.601    0.480
##     Panc_phys_symp          0.800    0.011   75.678    0.000    0.725    0.457
##     Pnc_dysf_thght          0.754    0.008   92.834    0.000    0.684    0.554
##     Panic_fer_rcrr          0.860    0.010   88.771    0.000    0.779    0.531
##     Panic_avod_bhv          0.879    0.009   96.258    0.000    0.796    0.572
##     Somatc_symptms          0.660    0.006  109.461    0.000    0.598    0.644
##     Hypochondriass          0.645    0.008   78.486    0.000    0.585    0.472
##     GAD_fre_flt_nx          1.080    0.009  116.374    0.000    0.979    0.682
##     GAD_nervousnss          1.137    0.010  116.835    0.000    1.030    0.685
##     GAD_irritablty          1.136    0.010  108.888    0.000    1.030    0.642
##     GAD_excss_wrry          1.048    0.008  125.276    0.000    0.950    0.728
##     DEP_rumination          1.053    0.009  123.127    0.000    0.955    0.715
##     DEP_psycmot_gt          0.977    0.008  119.446    0.000    0.885    0.697
##     OCD_sympt               0.686    0.007  103.852    0.000    0.622    0.614
##     DI_perfectinsm          0.585    0.008   75.682    0.000    0.530    0.456
##     EATA_restrictd          0.351    0.009   41.198    0.000    0.318    0.253
##     EATA_omit_meal          0.352    0.006   56.163    0.000    0.319    0.343
##     EATA_time_frms          0.332    0.006   54.169    0.000    0.301    0.331
##     EATA_cont_clrs          0.228    0.006   37.085    0.000    0.206    0.228
##     EATA_prgng_bhv          0.195    0.004   51.429    0.000    0.177    0.314
##     EATA_excs_xrcs          0.241    0.006   43.195    0.000    0.219    0.264
##     EATB_loss_ctrl          0.404    0.010   40.621    0.000    0.366    0.249
##     EATB_cope_ngff          0.415    0.010   40.185    0.000    0.376    0.247
##     EATB_overeatng          0.441    0.011   40.503    0.000    0.400    0.249
##     EATB_ntrs_thgh          0.470    0.007   65.435    0.000    0.426    0.397
##     DEP_incr_apptt          0.371    0.007   52.754    0.000    0.336    0.322
##     EATB_soc_dsgsd          0.584    0.009   65.301    0.000    0.529    0.396
##     PS_unusul_blfs          0.463    0.005   92.007    0.000    0.419    0.550
##     PS_prcpt_dysrg          0.614    0.006  105.784    0.000    0.557    0.624
##     Psychtc_symptm          0.520    0.005   97.688    0.000    0.471    0.581
##     Dissociation            0.747    0.006  120.283    0.000    0.677    0.701
##     PS_eccentricty          0.975    0.008  126.209    0.000    0.884    0.731
##     Manic_sympt             0.258    0.004   64.959    0.000    0.234    0.394
##     NA_sep_insecur          0.822    0.008  103.510    0.000    0.745    0.611
##     DT_suspicisnss          0.952    0.009  106.684    0.000    0.862    0.628
##     NA_perseveratn          0.765    0.007  111.372    0.000    0.694    0.653
##     DEP_psycmt_nhb          0.601    0.006   93.425    0.000    0.544    0.556
##     DEP_reduc_nrgy          0.778    0.007  113.121    0.000    0.705    0.665
##     DEP_hypersomni          0.428    0.007   63.127    0.000    0.388    0.383
##     DEP_crcd_md_fl          0.520    0.008   66.743    0.000    0.471    0.404
##     DEP_concentrtn          0.931    0.008  115.515    0.000    0.844    0.677
##     DEP_loss_apptt          0.444    0.006   78.705    0.000    0.402    0.473
##     DEP_loss_libid          0.313    0.008   39.856    0.000    0.284    0.244
##     DT_intimcy_avd          0.525    0.007   77.192    0.000    0.476    0.465
##     DT_rstr_ffctvt          0.786    0.008  100.832    0.000    0.712    0.598
##     DT_withdrawal           0.628    0.008   81.513    0.000    0.569    0.490
##     DEP_emot_numbn          0.830    0.008  102.495    0.000    0.752    0.606
##     DEP_crying_inb          0.791    0.008   95.457    0.000    0.717    0.567
##     DT_anhedonia            0.788    0.008   97.780    0.000    0.714    0.581
##     AN_callousness          0.229    0.004   53.572    0.000    0.207    0.327
##     AN_hostility            0.533    0.006   84.401    0.000    0.483    0.506
##     AN_deceitflnss          0.488    0.006   80.684    0.000    0.442    0.485
##     AN_manipultvns          0.345    0.006   60.344    0.000    0.313    0.367
##     AN_entitlement          0.180    0.005   35.629    0.000    0.163    0.219
##     AN_grandiosity          0.528    0.007   73.401    0.000    0.478    0.443
##     DI_impulsivity          0.798    0.007  109.609    0.000    0.723    0.644
##     Spec_phob_fear          0.682    0.009   79.933    0.000    0.618    0.481
##     Spec_phob_nmls          0.496    0.009   57.597    0.000    0.450    0.351
##     Spec_phob_blod          0.325    0.007   48.520    0.000    0.295    0.296
##     EATA_bdy_slfst          0.721    0.010   72.163    0.000    0.654    0.437
##     EATB_worries            0.738    0.009   79.707    0.000    0.669    0.480
##     EATA_fer_wghtg          0.644    0.010   62.411    0.000    0.584    0.380
##     EATA_fer_wghtl          0.275    0.007   39.559    0.000    0.250    0.242
##     Subst_use_drug          0.189    0.004   45.106    0.000    0.171    0.276
##     DI_normviolatn          0.382    0.007   56.559    0.000    0.347    0.345
##     Subst_use_alc           0.019    0.005    3.826    0.000    0.017    0.024
##     DI_risk_taking          0.353    0.005   64.726    0.000    0.320    0.392
##     DI_irrspnsblty          0.351    0.005   69.567    0.000    0.318    0.420
##   N9_Somatic_Dep =~
##     DEP_psycmt_nhb          1.000                               0.333    0.340
##     DEP_reduc_nrgy          1.636    0.030   55.041    0.000    0.545    0.514
##     DEP_hypersomni          0.863    0.024   36.257    0.000    0.287    0.284
##     DEP_crcd_md_fl          0.923    0.027   34.635    0.000    0.307    0.263
##     DEP_concentrtn          1.244    0.027   46.218    0.000    0.414    0.332
##     DEP_loss_apptt          0.442    0.017   25.456    0.000    0.147    0.173
##     DEP_loss_libid          1.028    0.029   35.916    0.000    0.342    0.294
##     DT_anhedonia            0.767    0.023   33.238    0.000    0.256    0.208
##     DEP_depr_mood           0.722    0.018   39.054    0.000    0.241    0.219
##     DEP_dmn_slfcnf          1.065    0.023   46.138    0.000    0.355    0.299
##     DEP_dimn_ntrst          1.391    0.027   51.068    0.000    0.463    0.383
##     DEP_decsn_prbl          0.806    0.023   35.227    0.000    0.269    0.204
##     Somatc_symptms          0.549    0.016   33.420    0.000    0.183    0.197
##   N2_Cognitive_Dep =~
##     DEP_hopelessn           1.000                               0.641    0.513
##     DEP_depr_mood           0.730    0.009   77.762    0.000    0.468    0.425
##     DEP_dmn_slfcnf          0.863    0.011   79.415    0.000    0.553    0.467
##     DEP_worthlessn          1.005    0.012   85.092    0.000    0.644    0.473
##     DEP_suicd_ntnt          0.452    0.010   43.734    0.000    0.290    0.252
##     DEP_dimn_ntrst          0.599    0.011   57.014    0.000    0.384    0.318
##     DEP_guilt               0.890    0.011   78.901    0.000    0.570    0.425
##     DEP_crying              0.411    0.011   36.166    0.000    0.263    0.223
##     DEP_decsn_prbl          0.527    0.011   49.456    0.000    0.338    0.256
##     PF_self_rflctn          0.042    0.011    3.864    0.000    0.027    0.021
##     DEP_psycmot_gt          0.321    0.010   33.004    0.000    0.206    0.162
##     DEP_reduc_nrgy          0.548    0.010   53.808    0.000    0.351    0.331
##     DEP_concentrtn          0.412    0.011   38.176    0.000    0.264    0.211
##     DT_anhedonia            0.376    0.010   35.913    0.000    0.241    0.196
##   N4_GAD_OCD =~
##     GAD_fre_flt_nx          1.000                               0.922    0.643
##     GAD_nervousnss          1.000    0.006  166.045    0.000    0.922    0.613
##     GAD_irritablty          0.871    0.008  116.157    0.000    0.803    0.501
##     GAD_excss_wrry          0.775    0.005  154.980    0.000    0.715    0.548
##     DEP_psycmot_gt          0.214    0.006   34.164    0.000    0.198    0.156
##     OCD_sympt               0.160    0.006   27.133    0.000    0.148    0.146
##     DI_perfectinsm          0.149    0.007   20.432    0.000    0.138    0.118
##     DEP_concentrtn          0.044    0.006    7.772    0.000    0.041    0.033
##     DEP_decsn_prbl          0.079    0.006   13.277    0.000    0.073    0.055
##     PF_self_rflctn         -0.067    0.006  -11.461    0.000   -0.062   -0.048
##   N1_SAD_Agora =~
##     SAD_soc_inter           1.000                               1.017    0.686
##     SAD_avoidance           0.986    0.006  175.999    0.000    1.003    0.670
##     SAD_performing          0.724    0.005  132.131    0.000    0.736    0.607
##     SAD_devaluatin          0.918    0.006  162.214    0.000    0.933    0.580
##     SAD_being_bsrv          0.749    0.006  128.040    0.000    0.761    0.556
##     SAD_eating              0.649    0.005  118.335    0.000    0.660    0.521
##     SAD_ttrctv_prs          0.663    0.008   86.758    0.000    0.674    0.452
##     Spec_phob_avod          0.461    0.008   61.385    0.000    0.469    0.307
##     Agoraphobia             0.433    0.006   76.376    0.000    0.441    0.382
##     DT_rstr_ffctvt          0.196    0.006   34.825    0.000    0.200    0.168
##     DT_withdrawal           0.309    0.006   51.509    0.000    0.314    0.270
##   N3_Panic =~
##     Panc_dctr_vsts          1.000                               0.557    0.749
##     Panic_attacks           1.606    0.012  139.128    0.000    0.894    0.714
##     Panc_phys_symp          2.093    0.015  141.523    0.000    1.165    0.734
##     Pnc_dysf_thght          1.499    0.011  138.656    0.000    0.834    0.676
##     Panic_fer_rcrr          1.672    0.013  126.121    0.000    0.931    0.634
##     Panic_avod_bhv          1.101    0.013   86.695    0.000    0.613    0.440
##     Somatc_symptms          0.452    0.008   56.996    0.000    0.251    0.271
##     Hypochondriass          0.543    0.013   43.413    0.000    0.302    0.244
##   N12_Spec_Phob =~
##     Spec_phob_fear          1.000                               0.960    0.748
##     Spec_phob_nmls          0.792    0.009   88.536    0.000    0.761    0.593
##     Spec_phob_blod          0.557    0.007   80.547    0.000    0.535    0.538
##     SAD_performing          0.044    0.005    8.325    0.000    0.042    0.035
##     Spec_phob_avod          0.481    0.009   54.568    0.000    0.461    0.302
##   N13_ED_thoughts =~
##     EATA_bdy_slfst          1.000                               1.237    0.826
##     EATB_worries            0.854    0.007  126.586    0.000    1.056    0.758
##     EATA_fer_wghtg          0.954    0.008  123.926    0.000    1.180    0.767
##     EATA_fer_wghtl          0.279    0.007   39.750    0.000    0.345    0.334
##     EATA_excs_xrcs          0.180    0.005   37.848    0.000    0.222    0.269
##     EATB_soc_dsgsd          0.492    0.006   81.401    0.000    0.609    0.456
##   N6_ED_Bulimic =~
##     EATB_loss_ctrl          1.000                               1.367    0.931
##     EATB_cope_ngff          0.953    0.004  254.442    0.000    1.302    0.854
##     EATB_overeatng          1.037    0.004  281.933    0.000    1.417    0.881
##     EATB_ntrs_thgh          0.517    0.003  156.385    0.000    0.707    0.658
##     DEP_incr_apptt          0.408    0.004  106.792    0.000    0.557    0.532
##     EATB_soc_dsgsd          0.364    0.007   52.327    0.000    0.497    0.372
##     EATA_bdy_slfst          0.163    0.013   12.969    0.000    0.222    0.148
##     EATB_worries            0.200    0.011   18.673    0.000    0.273    0.196
##     EATA_fer_wghtg          0.249    0.012   20.736    0.000    0.341    0.222
##   N5_ED_Anorectic =~
##     EATA_restrictd          1.000                               1.097    0.871
##     EATA_omit_meal          0.760    0.003  266.152    0.000    0.833    0.896
##     EATA_time_frms          0.679    0.003  210.553    0.000    0.745    0.817
##     EATA_cont_clrs          0.661    0.004  187.903    0.000    0.724    0.799
##     EATA_prgng_bhv          0.290    0.003  111.659    0.000    0.318    0.565
##     EATA_excs_xrcs          0.281    0.007   40.155    0.000    0.308    0.372
##     EATB_soc_dsgsd          0.112    0.014    7.808    0.000    0.123    0.092
##     EATA_bdy_slfst         -0.247    0.027   -9.039    0.000   -0.271   -0.181
##     EATB_worries           -0.153    0.023   -6.542    0.000   -0.167   -0.120
##     EATA_fer_wghtg         -0.141    0.026   -5.392    0.000   -0.155   -0.100
##   N10_Detachment =~
##     DT_intimcy_avd          1.000                               0.399    0.389
##     DT_rstr_ffctvt          1.407    0.026   54.217    0.000    0.561    0.471
##     DT_withdrawal           1.111    0.023   47.835    0.000    0.443    0.381
##     DEP_emot_numbn          1.388    0.026   52.785    0.000    0.553    0.446
##     DEP_crying_inb          0.853    0.022   38.455    0.000    0.340    0.269
##     DT_anhedonia            1.232    0.024   51.034    0.000    0.491    0.400
##     DEP_dimn_ntrst          0.552    0.014   38.580    0.000    0.220    0.182
##   N7_Thought_Disorder =~
##     PS_unusul_blfs          1.000                               0.482    0.632
##     PS_prcpt_dysrg          0.708    0.011   67.043    0.000    0.341    0.382
##     Psychtc_symptm          0.930    0.011   86.025    0.000    0.448    0.553
##     Dissociation            0.534    0.010   51.869    0.000    0.257    0.266
##     Manic_sympt             0.272    0.008   32.463    0.000    0.131    0.220
##     DEP_crying_inb          0.102    0.015    6.869    0.000    0.049    0.039
##     PF_self_rflctn         -0.233    0.015  -15.042    0.000   -0.112   -0.087
##     OCD_sympt               0.386    0.012   31.250    0.000    0.186    0.183
##   N11_Antagonism =~
##     AN_callousness          1.000                               0.363    0.573
##     AN_hostility            1.048    0.018   58.288    0.000    0.381    0.399
##     AN_deceitflnss          1.264    0.018   68.529    0.000    0.460    0.504
##     AN_manipultvns          1.315    0.019   70.422    0.000    0.478    0.560
##     AN_entitlement          0.931    0.016   58.601    0.000    0.338    0.454
##     AN_grandiosity          1.199    0.021   57.122    0.000    0.436    0.403
##     DI_irrspnsblty          0.300    0.016   18.341    0.000    0.109    0.144
##   N14_Disinhibition =~
##     Subst_use_drug          1.000                               0.274    0.441
##     DI_normviolatn          2.700    0.048   56.002    0.000    0.741    0.737
##     Subst_use_alc           0.862    0.022   39.487    0.000    0.237    0.326
##     DI_risk_taking          1.127    0.024   46.458    0.000    0.309    0.379
##     DI_irrspnsblty          0.382    0.022   17.237    0.000    0.105    0.138
##     Manic_sympt             0.538    0.016   32.863    0.000    0.148    0.248
##
## Covariances:
##                          Estimate  Std.Err  z-value  P(>|z|)   Std.lv  Std.all
##   N8_PF_Problems ~~
##     N5_ED_Anorectc          0.000                               0.000    0.000
##     N1_SAD_Agora            0.000                               0.000    0.000
##     N2_Cognitiv_Dp          0.000                               0.000    0.000
##     N3_Panic                0.000                               0.000    0.000
##     N4_GAD_OCD              0.000                               0.000    0.000
##     N6_ED_Bulimic           0.000                               0.000    0.000
##     N7_Thght_Dsrdr          0.000                               0.000    0.000
##     N9_Somatic_Dep          0.000                               0.000    0.000
##     N10_Detachment          0.000                               0.000    0.000
##     N11_Antagonism          0.000                               0.000    0.000
##     N12_Spec_Phob           0.000                               0.000    0.000
##     N13_ED_thoghts          0.000                               0.000    0.000
##     N14_Disinhibtn          0.000                               0.000    0.000
##   N9_Somatic_Dep ~~
##     N2_Cognitiv_Dp          0.016    0.003    6.231    0.000    0.075    0.075
##     N4_GAD_OCD              0.050    0.003   19.114    0.000    0.164    0.164
##     N1_SAD_Agora            0.003    0.003    1.201    0.230    0.009    0.009
##     N3_Panic                0.015    0.001   10.095    0.000    0.081    0.081
##     N12_Spec_Phob           0.007    0.003    2.795    0.005    0.023    0.023
##     N13_ED_thoghts         -0.025    0.003   -7.345    0.000   -0.060   -0.060
##     N6_ED_Bulimic           0.044    0.003   12.600    0.000    0.096    0.096
##     N5_ED_Anorectc         -0.027    0.003   -9.893    0.000   -0.075   -0.075
##     N10_Detachment          0.025    0.001   18.025    0.000    0.188    0.188
##     N7_Thght_Dsrdr         -0.012    0.002   -7.692    0.000   -0.073   -0.073
##     N11_Antagonism         -0.006    0.001   -5.524    0.000   -0.051   -0.051
##     N14_Disinhibtn         -0.003    0.001   -3.919    0.000   -0.035   -0.035
##   N2_Cognitive_Dep ~~
##     N4_GAD_OCD              0.075    0.005   16.023    0.000    0.127    0.127
##     N1_SAD_Agora            0.067    0.005   13.390    0.000    0.103    0.103
##     N3_Panic               -0.047    0.003  -16.876    0.000   -0.130   -0.130
##     N12_Spec_Phob          -0.053    0.005  -10.683    0.000   -0.086   -0.086
##     N13_ED_thoghts          0.141    0.006   22.824    0.000    0.178    0.178
##     N6_ED_Bulimic           0.035    0.006    5.522    0.000    0.040    0.040
##     N5_ED_Anorectc          0.014    0.005    2.743    0.006    0.020    0.020
##     N10_Detachment          0.031    0.002   12.807    0.000    0.121    0.121
##     N7_Thght_Dsrdr         -0.150    0.003  -49.119    0.000   -0.485   -0.485
##     N11_Antagonism         -0.082    0.002  -36.809    0.000   -0.354   -0.354
##     N14_Disinhibtn         -0.048    0.002  -28.050    0.000   -0.275   -0.275
##   N4_GAD_OCD ~~
##     N1_SAD_Agora            0.158    0.007   23.583    0.000    0.169    0.169
##     N3_Panic                0.136    0.004   36.071    0.000    0.265    0.265
##     N12_Spec_Phob           0.158    0.007   23.679    0.000    0.179    0.179
##     N13_ED_thoghts          0.011    0.008    1.332    0.183    0.009    0.009
##     N6_ED_Bulimic          -0.049    0.008   -5.822    0.000   -0.039   -0.039
##     N5_ED_Anorectc         -0.064    0.007   -9.321    0.000   -0.063   -0.063
##     N10_Detachment         -0.088    0.003  -26.984    0.000   -0.240   -0.240
##     N7_Thght_Dsrdr         -0.113    0.004  -31.427    0.000   -0.255   -0.255
##     N11_Antagonism         -0.079    0.003  -28.672    0.000   -0.237   -0.237
##     N14_Disinhibtn         -0.056    0.002  -25.393    0.000   -0.220   -0.220
##   N1_SAD_Agora ~~
##     N3_Panic                0.071    0.004   18.109    0.000    0.126    0.126
##     N12_Spec_Phob           0.382    0.008   49.229    0.000    0.391    0.391
##     N13_ED_thoghts          0.143    0.009   15.604    0.000    0.114    0.114
##     N6_ED_Bulimic           0.024    0.009    2.595    0.009    0.017    0.017
##     N5_ED_Anorectc          0.116    0.007   15.665    0.000    0.104    0.104
##     N10_Detachment          0.035    0.004    9.942    0.000    0.087    0.087
##     N7_Thght_Dsrdr         -0.002    0.004   -0.483    0.629   -0.004   -0.004
##     N11_Antagonism         -0.034    0.003  -11.761    0.000   -0.091   -0.091
##     N14_Disinhibtn         -0.051    0.002  -21.910    0.000   -0.181   -0.181
##   N3_Panic ~~
##     N12_Spec_Phob           0.175    0.004   42.089    0.000    0.326    0.326
##     N13_ED_thoghts         -0.008    0.005   -1.786    0.074   -0.012   -0.012
##     N6_ED_Bulimic          -0.046    0.005   -9.128    0.000   -0.060   -0.060
##     N5_ED_Anorectc          0.003    0.004    0.846    0.398    0.006    0.006
##     N10_Detachment         -0.057    0.002  -28.889    0.000   -0.258   -0.258
##     N7_Thght_Dsrdr          0.031    0.002   14.847    0.000    0.115    0.115
##     N11_Antagonism         -0.028    0.002  -17.329    0.000   -0.136   -0.136
##     N14_Disinhibtn         -0.003    0.001   -2.194    0.028   -0.017   -0.017
##   N12_Spec_Phob ~~
##     N13_ED_thoghts          0.028    0.009    3.252    0.001    0.023    0.023
##     N6_ED_Bulimic          -0.002    0.009   -0.260    0.795   -0.002   -0.002
##     N5_ED_Anorectc         -0.016    0.007   -2.157    0.031   -0.015   -0.015
##     N10_Detachment         -0.076    0.003  -22.030    0.000   -0.198   -0.198
##     N7_Thght_Dsrdr          0.068    0.004   17.924    0.000    0.147    0.147
##     N11_Antagonism         -0.013    0.003   -4.655    0.000   -0.038   -0.038
##     N14_Disinhibtn         -0.050    0.002  -21.859    0.000   -0.192   -0.192
##   N13_ED_thoughts ~~
##     N6_ED_Bulimic           0.589    0.028   21.345    0.000    0.349    0.349
##     N5_ED_Anorectc          0.898    0.034   26.187    0.000    0.662    0.662
##     N10_Detachment         -0.011    0.004   -2.776    0.006   -0.022   -0.022
##     N7_Thght_Dsrdr          0.010    0.005    2.111    0.035    0.017    0.017
##     N11_Antagonism         -0.025    0.004   -6.995    0.000   -0.055   -0.055
##     N14_Disinhibtn          0.006    0.003    2.171    0.030    0.018    0.018
##   N6_ED_Bulimic ~~
##     N5_ED_Anorectc          0.470    0.010   46.593    0.000    0.314    0.314
##     N10_Detachment          0.010    0.004    2.360    0.018    0.018    0.018
##     N7_Thght_Dsrdr          0.008    0.005    1.750    0.080    0.013    0.013
##     N11_Antagonism          0.044    0.004   11.767    0.000    0.088    0.088
##     N14_Disinhibtn          0.002    0.003    0.722    0.470    0.005    0.005
##   N5_ED_Anorectic ~~
##     N10_Detachment          0.030    0.003    8.831    0.000    0.068    0.068
##     N7_Thght_Dsrdr          0.073    0.004   18.571    0.000    0.139    0.139
##     N11_Antagonism          0.029    0.003    9.800    0.000    0.073    0.073
##     N14_Disinhibtn          0.032    0.002   13.790    0.000    0.106    0.106
##   N10_Detachment ~~
##     N7_Thght_Dsrdr         -0.004    0.002   -2.211    0.027   -0.020   -0.020
##     N11_Antagonism          0.036    0.001   24.626    0.000    0.246    0.246
##     N14_Disinhibtn          0.006    0.001    5.636    0.000    0.051    0.051
##   N7_Thought_Disorder ~~
##     N11_Antagonism          0.059    0.002   35.653    0.000    0.339    0.339
##     N14_Disinhibtn          0.042    0.001   30.946    0.000    0.317    0.317
##   N11_Antagonism ~~
##     N14_Disinhibtn          0.047    0.001   38.415    0.000    0.472    0.472
##
## Variances:
##                    Estimate  Std.Err  z-value  P(>|z|)   Std.lv  Std.all
##    .PF_affect_comm    0.641    0.006  112.074    0.000    0.641    0.438
##    .PF_affect_diff    0.474    0.004  107.795    0.000    0.474    0.286
##    .PF_affct_tlrnc    0.518    0.005  109.728    0.000    0.518    0.339
##    .PF_frmng_rlshp    1.317    0.011  115.463    0.000    1.317    0.761
##    .PF_hlstc_prcpt    0.751    0.007  114.564    0.000    0.751    0.636
##    .PF_identity       0.248    0.002  103.588    0.000    0.248    0.214
##    .PF_impulse_reg    0.750    0.007  114.487    0.000    0.750    0.627
##    .PF_ntrn_rlshp_    0.691    0.006  113.021    0.000    0.691    0.497
##    .PF_self_rflctn    0.516    0.005  103.107    0.000    0.516    0.310
##    .PF_self_est_rg    0.586    0.005  113.217    0.000    0.586    0.511
##    .PF_anticipatin    0.908    0.008  113.935    0.000    0.908    0.571
##    .SAD_soc_inter     0.411    0.005   91.169    0.000    0.411    0.187
##    .SAD_avoidance     0.336    0.004   86.731    0.000    0.336    0.150
##    .SAD_performing    0.403    0.004  103.011    0.000    0.403    0.274
##    .SAD_devaluatin    0.417    0.004   95.686    0.000    0.417    0.161
##    .SAD_being_bsrv    0.605    0.006  107.169    0.000    0.605    0.323
##    .SAD_eating        0.565    0.005  108.977    0.000    0.565    0.352
##    .SAD_ttrctv_prs    1.269    0.011  113.072    0.000    1.269    0.571
##    .Spec_phob_avod    0.996    0.009  108.302    0.000    0.996    0.426
##    .Agoraphobia       0.722    0.006  113.790    0.000    0.722    0.543
##    .DEP_hopelessn     0.539    0.006   94.419    0.000    0.539    0.345
##    .DEP_depr_mood     0.312    0.003  102.534    0.000    0.312    0.258
##    .DEP_dmn_slfcnf    0.314    0.003   96.352    0.000    0.314    0.224
##    .DEP_worthlessn    0.560    0.006   94.826    0.000    0.560    0.302
##    .DEP_suicd_ntnt    0.800    0.007  113.097    0.000    0.800    0.603
##    .DEP_dimn_ntrst    0.301    0.003   90.868    0.000    0.301    0.206
##    .DEP_guilt         0.600    0.006  100.501    0.000    0.600    0.333
##    .DEP_crying        1.016    0.009  114.425    0.000    1.016    0.733
##    .DEP_decsn_prbl    0.655    0.006  112.041    0.000    0.655    0.377
##    .Panc_dctr_vsts    0.208    0.002   97.492    0.000    0.208    0.378
##    .Panic_attacks     0.408    0.004   92.408    0.000    0.408    0.260
##    .Panc_phys_symp    0.634    0.007   90.067    0.000    0.634    0.252
##    .Pnc_dysf_thght    0.359    0.004   92.655    0.000    0.359    0.236
##    .Panic_fer_rcrr    0.680    0.007  100.921    0.000    0.680    0.316
##    .Panic_avod_bhv    0.928    0.008  111.194    0.000    0.928    0.479
##    .Somatc_symptms    0.401    0.004  110.829    0.000    0.401    0.465
##    .Hypochondriass    1.103    0.010  115.136    0.000    1.103    0.718
##    .GAD_fre_flt_nx    0.248    0.004   67.694    0.000    0.248    0.121
##    .GAD_nervousnss    0.353    0.004   81.505    0.000    0.353    0.156
##    .GAD_irritablty    0.870    0.008  106.364    0.000    0.870    0.338
##    .GAD_excss_wrry    0.289    0.003   91.265    0.000    0.289    0.170
##    .DEP_rumination    0.872    0.008  112.901    0.000    0.872    0.489
##    .DEP_psycmot_gt    0.739    0.007  113.372    0.000    0.739    0.458
##    .OCD_sympt         0.598    0.005  113.500    0.000    0.598    0.582
##    .DI_perfectinsm    1.054    0.009  115.845    0.000    1.054    0.778
##    .EATA_restrictd    0.282    0.003   89.966    0.000    0.282    0.178
##    .EATA_omit_meal    0.069    0.001   57.563    0.000    0.069    0.080
##    .EATA_time_frms    0.185    0.002   98.657    0.000    0.185    0.223
##    .EATA_cont_clrs    0.254    0.002  104.424    0.000    0.254    0.309
##    .EATA_prgng_bhv    0.184    0.002  113.412    0.000    0.184    0.582
##    .EATA_excs_xrcs    0.402    0.004  113.552    0.000    0.402    0.587
##    .EATB_loss_ctrl    0.155    0.003   50.310    0.000    0.155    0.072
##    .EATB_cope_ngff    0.487    0.005   96.326    0.000    0.487    0.210
##    .EATB_overeatng    0.419    0.005   87.210    0.000    0.419    0.162
##    .EATB_ntrs_thgh    0.474    0.004  110.612    0.000    0.474    0.410
##    .DEP_incr_apptt    0.671    0.006  113.996    0.000    0.671    0.613
##    .EATB_soc_dsgsd    0.522    0.005  107.072    0.000    0.522    0.292
##    .PS_unusul_blfs    0.173    0.003   66.076    0.000    0.173    0.298
##    .PS_prcpt_dysrg    0.369    0.004  104.250    0.000    0.369    0.464
##    .Psychtc_symptm    0.234    0.003   83.205    0.000    0.234    0.357
##    .Dissociation      0.408    0.004  108.779    0.000    0.408    0.438
##    .PS_eccentricty    0.682    0.006  112.547    0.000    0.682    0.466
##    .Manic_sympt       0.248    0.002  110.158    0.000    0.248    0.700
##    .NA_sep_insecur    0.931    0.008  114.479    0.000    0.931    0.627
##    .DT_suspicisnss    1.141    0.010  114.283    0.000    1.141    0.605
##    .NA_perseveratn    0.646    0.006  113.958    0.000    0.646    0.573
##    .DEP_psycmt_nhb    0.550    0.005  107.876    0.000    0.550    0.575
##    .DEP_reduc_nrgy    0.179    0.003   68.546    0.000    0.179    0.159
##    .DEP_hypersomni    0.793    0.007  112.262    0.000    0.793    0.773
##    .DEP_crcd_md_fl    1.045    0.009  112.796    0.000    1.045    0.767
##    .DEP_concentrtn    0.576    0.005  107.868    0.000    0.576    0.370
##    .DEP_loss_apptt    0.541    0.005  114.563    0.000    0.541    0.746
##    .DEP_loss_libid    1.155    0.010  112.360    0.000    1.155    0.854
##    .DT_intimcy_avd    0.663    0.006  105.582    0.000    0.663    0.632
##    .DT_rstr_ffctvt    0.539    0.006   89.244    0.000    0.539    0.380
##    .DT_withdrawal     0.708    0.007  102.730    0.000    0.708    0.524
##    .DEP_emot_numbn    0.667    0.007   95.117    0.000    0.667    0.433
##    .DEP_crying_inb    0.965    0.009  110.554    0.000    0.965    0.604
##    .DT_anhedonia      0.550    0.006   94.425    0.000    0.550    0.364
##    .AN_callousness    0.227    0.002   93.006    0.000    0.227    0.565
##    .AN_hostility      0.535    0.005  105.361    0.000    0.535    0.585
##    .AN_deceitflnss    0.424    0.004   96.315    0.000    0.424    0.510
##    .AN_manipultvns    0.401    0.004   93.538    0.000    0.401    0.551
##    .AN_entitlement    0.415    0.004  105.454    0.000    0.415    0.746
##    .AN_grandiosity    0.748    0.007  106.237    0.000    0.748    0.641
##    .DI_impulsivity    0.739    0.006  114.086    0.000    0.739    0.585
##    .Spec_phob_fear    0.344    0.008   42.308    0.000    0.344    0.209
##    .Spec_phob_nmls    0.863    0.009   94.838    0.000    0.863    0.525
##    .Spec_phob_blod    0.616    0.006  102.410    0.000    0.616    0.623
##    .EATA_bdy_slfst    0.452    0.006   76.438    0.000    0.452    0.201
##    .EATB_worries      0.339    0.004   77.772    0.000    0.339    0.174
##    .EATA_fer_wghtg    0.487    0.006   82.976    0.000    0.487    0.206
##    .EATA_fer_wghtl    0.881    0.008  115.012    0.000    0.881    0.830
##    .Subst_use_drug    0.282    0.003  103.392    0.000    0.282    0.729
##    .DI_normviolatn    0.340    0.008   43.537    0.000    0.340    0.337
##    .Subst_use_alc     0.469    0.004  111.004    0.000    0.469    0.893
##    .DI_risk_taking    0.467    0.004  106.575    0.000    0.467    0.702
##    .DI_irrspnsblty    0.438    0.004  113.921    0.000    0.438    0.765
##     N8_PF_Problems    0.822    0.011   72.138    0.000    1.000    1.000
##     N9_Somatic_Dep    0.111    0.004   29.577    0.000    1.000    1.000
##     N2_Cognitiv_Dp    0.411    0.008   53.280    0.000    1.000    1.000
##     N4_GAD_OCD        0.850    0.010   84.798    0.000    1.000    1.000
##     N1_SAD_Agora      1.034    0.012   82.743    0.000    1.000    1.000
##     N3_Panic          0.310    0.004   72.753    0.000    1.000    1.000
##     N12_Spec_Phob     0.922    0.013   70.732    0.000    1.000    1.000
##     N13_ED_thoghts    1.530    0.054   28.578    0.000    1.000    1.000
##     N6_ED_Bulimic     1.868    0.018  106.288    0.000    1.000    1.000
##     N5_ED_Anorectc    1.203    0.013   94.702    0.000    1.000    1.000
##     N10_Detachment    0.159    0.005   31.773    0.000    1.000    1.000
##     N7_Thght_Dsrdr    0.232    0.004   60.001    0.000    1.000    1.000
##     N11_Antagonism    0.132    0.003   46.795    0.000    1.000    1.000
##     N14_Disinhibtn    0.075    0.002   33.139    0.000    1.000    1.000

library(semPlot)
semPaths(s1model,bifactor = "N8_PF_Problems",residuals = F, layout = "tree2", rotation=2,node.width = .8,cardinal = F,nCharNodes = 0, asize =  1, node.width = .5, reorder = T)

Include psychosocial impairment and affective well-being indicators as covariates in model to calculate latent regressions

addmood = "
mood_av ~~ N5_ED_Anorectic
mood_av ~~ N1_SAD_Agora
mood_av ~~ N2_Cognitive_Dep
mood_av ~~ N3_Panic
mood_av ~~ N4_GAD_OCD
mood_av ~~ N6_ED_Bulimic
mood_av ~~ N7_Thought_Disorder
mood_av ~~ N9_Somatic_Dep
mood_av ~~ N10_Detachment
mood_av ~~ N11_Antagonism
mood_av ~~ N12_Spec_Phob
mood_av ~~ N13_ED_thoughts
mood_av ~~ N14_Disinhibition
mood_av ~~ N8_PF_Problems


"

addpsf = "
psf_av ~~ N5_ED_Anorectic
psf_av ~~ N1_SAD_Agora
psf_av ~~ N2_Cognitive_Dep
psf_av ~~ N3_Panic
psf_av ~~ N4_GAD_OCD
psf_av ~~ N6_ED_Bulimic
psf_av ~~ N7_Thought_Disorder
psf_av ~~ N9_Somatic_Dep
psf_av ~~ N10_Detachment
psf_av ~~ N11_Antagonism
psf_av ~~ N12_Spec_Phob
psf_av ~~ N13_ED_thoughts
psf_av ~~ N14_Disinhibition
psf_av ~~ N8_PF_Problems
"



addmood_fu = "
mood_av_fu ~~ N5_ED_Anorectic
mood_av_fu ~~ N1_SAD_Agora
mood_av_fu ~~ N2_Cognitive_Dep
mood_av_fu ~~ N3_Panic
mood_av_fu ~~ N4_GAD_OCD
mood_av_fu ~~ N6_ED_Bulimic
mood_av_fu ~~ N7_Thought_Disorder
mood_av_fu ~~ N9_Somatic_Dep
mood_av_fu ~~ N10_Detachment
mood_av_fu ~~ N11_Antagonism
mood_av_fu ~~ N12_Spec_Phob
mood_av_fu ~~ N13_ED_thoughts
mood_av_fu ~~ N14_Disinhibition
mood_av_fu ~~ N8_PF_Problems


"

addpsf_fu = "
psf_av_fu ~~ N5_ED_Anorectic
psf_av_fu ~~ N1_SAD_Agora
psf_av_fu ~~ N2_Cognitive_Dep
psf_av_fu ~~ N3_Panic
psf_av_fu ~~ N4_GAD_OCD
psf_av_fu ~~ N6_ED_Bulimic
psf_av_fu ~~ N7_Thought_Disorder
psf_av_fu ~~ N9_Somatic_Dep
psf_av_fu ~~ N10_Detachment
psf_av_fu ~~ N11_Antagonism
psf_av_fu ~~ N12_Spec_Phob
psf_av_fu ~~ N13_ED_thoughts
psf_av_fu ~~ N14_Disinhibition
psf_av_fu ~~ N8_PF_Problems
"




#### estimate standard bifactor and correlated factors model
## ML
# basemodel.ml <- cfa(model = paste0(GPbase), sample.cov = scale_cov_pred, sample.nobs = 27173, std.lv = T)
s1model.ml.psf =  cfa(model = paste0(s1wn8asref, addpsf), sample.cov = scale_cov, sample.nobs = 27173)

s1model.ml.mood =  cfa(model = paste0(s1wn8asref, addmood), sample.cov = scale_cov, sample.nobs = 27173)

s1model.ml.psf_fu =  cfa(model = paste0(s1wn8asref, addpsf_fu), sample.cov = scale_cov, sample.nobs = 27173)

s1model.ml.mood_fu =  cfa(model = paste0(s1wn8asref, addmood_fu), sample.cov = scale_cov, sample.nobs = 27173)
#
# ##extract correlations
#
# scors.s1.psf = data.frame(inspect(s1model.ml.psf, what = "cor.all"))[c(vars,"psf_av"),c(vars,"psf_av")]
# scors.s1.mood = data.frame(inspect(s1model.ml.mood, what = "cor.all"))[c(vars,"mood_av"),c(vars,"mood_av")]


# cor.plot(scors.s1.mood, xlas = 2)
# cor.plot(scors.s1.psf, xlas = 2)

#
# predframe = data.frame(data_scale_pred[,c("id","psf_av","mood_av")])
# predframe = merge(predframe,fsnidextr.func(s1model.ml, data_scale_pred, vars), all.x = T )

#
# test = cor(predframe[c(vars,"mood_av")], use="pairwise.complete.obs")

Regression models and squared semipartial correlations (Delta Rsquare) for s-1 model based on latent correlation matrix: Average over one year after PF items were added

vars = c( "N8_PF_Problems",
 "N9_Somatic_Dep" ,      "N2_Cognitive_Dep"  ,
  "N4_GAD_OCD"  ,
 "N1_SAD_Agora"      ,   "N3_Panic"       ,     "N12_Spec_Phob"      ,
    "N13_ED_thoughts"  ,  "N6_ED_Bulimic"   ,
 "N5_ED_Anorectic"         ,   "N7_Thought_Disorder", "N10_Detachment",
      "N11_Antagonism"   ,      "N14_Disinhibition"       )


 predvarnames.s1 = vars


##correlations s1
cor_mat.m <- data.frame(inspect(s1model.ml.mood, what = "cor.all"))
cor_mat.m.s1 = cor_mat.m[c("mood_av",predvarnames.s1),
                            c("mood_av",predvarnames.s1)]

##correlations s1 psf
cor_mat.psf <- data.frame(inspect(s1model.ml.psf, what = "cor.all"))
cor_mat.psf.s1 = cor_mat.psf[c("psf_av",predvarnames.s1),
                            c("psf_av",predvarnames.s1)]


#### Regression model based on correlation matrix
library(metafor)
# # with cov
# regpsfs1 = matreg("psf_av", x=predvarnames.s1, cov_mat.psf.s1, cov = T, n=27173)

#with cor
regpsfs1 = matreg("psf_av", x=predvarnames.s1, cor_mat.psf.s1, n=27173, digits = 2,)
regmoods1 = matreg("mood_av", x=predvarnames.s1, cor_mat.m.s1, n=27173, digits = 2)



### calculate squared semipartial correlations for s factors and hierarchy
## average affective well-being up to 1y
# Extract the variances of the latent variables
cov_mat.m <- data.frame(inspect(s1model.ml.mood, what = "cov.all"))

cov_mat.m.s1 = cov_mat.m[c("mood_av",predvarnames.s1),
                            c("mood_av",predvarnames.s1)]

sqsemipcor.m.s1 <- round(data.frame(correlation::cor_to_spcor(cor = as.matrix(cor_mat.m.s1),cov = as.matrix(cov_mat.m.s1)))^2,digits = 4)

rownames(sqsemipcor.m.s1) = c("mood_av",predvarnames.s1)
colnames(sqsemipcor.m.s1) = c("mood_av",predvarnames.s1)


## average psychosocial functioning up to 1y
# Extract the variance covariance matrix of the latent and observed variables
cov_mat.psf <- data.frame(inspect(s1model.ml.psf, what = "cov.all"))

cov_mat.psf.s1 = cov_mat.psf[c("psf_av",predvarnames.s1),
                            c("psf_av",predvarnames.s1)]


sqsemipcor.psf.s1 <- round(data.frame(correlation::cor_to_spcor(cor = as.matrix(cor_mat.psf.s1) ,cov = as.matrix(cov_mat.psf.s1)))^2,digits = 4)

rownames(sqsemipcor.psf.s1) = c("psf_av",predvarnames.s1)
colnames(sqsemipcor.psf.s1) = c("psf_av",predvarnames.s1)


#
# #test with factor scores
# predframe = data.frame(data_scale_pred[,c("id","psf_av","mood_av")])
# predframe = merge(predframe,fsnidextr.func(s1model.ml, data_scale_pred, vars), all.x = T )
# temp = lm(paste0("mood_av ~ ", paste(vars[1:3], collapse = " + ")) , predframe)
# rockchalk::getDeltaRsquare(temp)




# correlation::cor_to_spcor(cortemp, covtemp)^2


## do tables
regpsf1tab = round(regpsfs1$tab, digits = 2)
regmoods1tab =round(regmoods1$tab, digits = 2)
regmoods1tab$df = NULL
regpsf1tab$df = NULL

regnsspc1 = cbind(regmoods1tab,SSPC = round(as.numeric(t(sqsemipcor.m.s1[1,2:15])), digits=3), regpsf1tab,SSPC = round(as.numeric(t(sqsemipcor.psf.s1[1,2:15])), digits=3))
rownames(regnsspc1) =  c("N8:PF Problems", "N9:Som. Depr. (resid.)", "N2:Cogn. Depr. (resid.)", "N4:GAD & OCD (resid.)", "N1:SAD & Agor. (resid.)", "N3:Panic & Som. (resid.)", "N12:Spec. Phob. (resid.)", "N13:ED Thoughts (resid.)", "N6:ED Bulimic (resid.)", "N5:ED Anorectic (resid.)", "N7:Thought Dis. (resid.)", "N10:Detachment (resid.)", "N11:Antagonism (resid.)", "N14:Disinhibition (resid.)")

kable(regnsspc1,caption = "Regression based on correlation and covariance matrix from bifactor-s-1 model. Prediction of average psychosocial functioning and affective well-being over one year 10 months after first assessment.") %>%
        kable_styling("striped") %>%
         add_header_above(c("   " = 1,"Prediction of average affective well-being over 1 year" = 7, "Prediction of average psychosocial functioning over 1 year" = 7)) %>% kable_paper %>% add_footnote(paste0("Table 5. ","N (affective well-being) = 25844",". N (psychosocial functioning) = 27173" ,". Total R2 affective well-being = ",round(regmoods1$R2,3),". Total R2 psychosocial functioning = ",round(regpsfs1$R2,3),". SSPC = squared semipartial correlation, i.e. the unique contribution of a specific predictor variable in a multiple regression analysis, while controlling for other predictors in the model. "))

Regression based on correlation and covariance matrix from bifactor-s-1 model. Prediction of average psychosocial functioning and affective well-being over one year 10 months after first assessment.

	Prediction of average affective well-being over 1 year							Prediction of average psychosocial functioning over 1 year
	beta	se	tval	pval	ci.lb	ci.ub	SSPC	beta	se	tval	pval	ci.lb	ci.ub	SSPC
N8:PF Problems	-0.55	0.00	-146.88	0.00	-0.56	-0.55	0.306	0.77	0	322.29	0.00	0.76	0.77	0.586
N9:Som. Depr. (resid.)	-0.12	0.00	-29.97	0.00	-0.13	-0.11	0.013	0.34	0	133.37	0.00	0.33	0.34	0.100
N2:Cogn. Depr. (resid.)	-0.52	0.00	-112.10	0.00	-0.53	-0.51	0.178	0.34	0	116.66	0.00	0.34	0.35	0.077
N4:GAD & OCD (resid.)	0.01	0.00	2.91	0.00	0.00	0.02	0.000	-0.02	0	-6.58	0.00	-0.02	-0.01	0.000
N1:SAD & Agor. (resid.)	0.01	0.00	2.10	0.04	0.00	0.02	0.000	0.08	0	30.66	0.00	0.08	0.09	0.005
N3:Panic & Som. (resid.)	0.00	0.00	-0.21	0.83	-0.01	0.01	0.000	0.01	0	5.43	0.00	0.01	0.02	0.000
N12:Spec. Phob. (resid.)	-0.02	0.00	-3.63	0.00	-0.03	-0.01	0.000	0.00	0	-0.72	0.47	-0.01	0.00	0.000
N13:ED Thoughts (resid.)	0.08	0.01	14.87	0.00	0.07	0.09	0.003	0.00	0	1.39	0.17	0.00	0.01	0.000
N6:ED Bulimic (resid.)	0.00	0.00	0.11	0.91	-0.01	0.01	0.000	0.01	0	5.57	0.00	0.01	0.02	0.000
N5:ED Anorectic (resid.)	-0.07	0.01	-13.24	0.00	-0.08	-0.06	0.002	-0.01	0	-2.94	0.00	-0.02	0.00	0.000
N7:Thought Dis. (resid.)	-0.10	0.00	-22.21	0.00	-0.11	-0.10	0.007	-0.01	0	-3.55	0.00	-0.02	0.00	0.000
N10:Detachment (resid.)	-0.14	0.00	-32.27	0.00	-0.15	-0.14	0.015	0.03	0	9.44	0.00	0.02	0.03	0.000
N11:Antagonism (resid.)	0.02	0.00	4.53	0.00	0.01	0.03	0.000	0.02	0	7.12	0.00	0.02	0.03	0.000
N14:Disinhibition (resid.)	0.03	0.00	7.53	0.00	0.03	0.04	0.001	0.01	0	4.98	0.00	0.01	0.02	0.000
a Table 5. N (affective well-being) = 25844. N (psychosocial functioning) = 27173. Total R2 affective well-being = 0.614. Total R2 psychosocial functioning = 0.847. SSPC = squared semipartial correlation, i.e. the unique contribution of a specific predictor variable in a multiple regression analysis, while controlling for other predictors in the model.

Regression models and squared semipartial correlations (Delta Rsquare) for s-1 model based on latent correlation matrix: Separated assessments

vars = c( "N8_PF_Problems",
 "N9_Somatic_Dep" ,      "N2_Cognitive_Dep"  ,
  "N4_GAD_OCD"  ,
 "N1_SAD_Agora"      ,   "N3_Panic"       ,     "N12_Spec_Phob"      ,
    "N13_ED_thoughts"  ,  "N6_ED_Bulimic"   ,
 "N5_ED_Anorectic"         ,   "N7_Thought_Disorder", "N10_Detachment",
      "N11_Antagonism"   ,      "N14_Disinhibition"       )


 predvarnames.s1 = vars

##correlations s1
cor_mat.m <- data.frame(inspect(s1model.ml.mood_fu, what = "cor.all"))
cor_mat.m.s1 = cor_mat.m[c("mood_av_fu",predvarnames.s1),
                            c("mood_av_fu",predvarnames.s1)]

##correlations s1 psf_fu
cor_mat.psf <- data.frame(inspect(s1model.ml.psf_fu, what = "cor.all"))
cor_mat.psf.s1 = cor_mat.psf[c("psf_av_fu",predvarnames.s1),
                            c("psf_av_fu",predvarnames.s1)]


#### Regression model based on correlatoin matrix
library(metafor)
# # with cov
# regpsfs1 = matreg("psf_av_fu", x=predvarnames.s1, cov_mat.psf.s1, cov = T, n=27173)

#with cor
regpsfs2 = matreg("psf_av_fu", x=predvarnames.s1, cor_mat.psf.s1, n=27173, digits = 2)
regmoods2 = matreg("mood_av_fu", x=predvarnames.s1, cor_mat.m.s1, n=27173, digits = 2)



### calculate squared semipartial correlations for s factors and hierarchy
## average affective well-being up to 1y
# Extract the variances of the latent variables
cov_mat.m <- data.frame(inspect(s1model.ml.mood_fu, what = "cov.all"))

cov_mat.m.s1 = cov_mat.m[c("mood_av_fu",predvarnames.s1),
                            c("mood_av_fu",predvarnames.s1)]

sqsemipcor.m.s1 <- round(data.frame(correlation::cor_to_spcor(cor = as.matrix(cor_mat.m.s1),cov = as.matrix(cov_mat.m.s1)))^2,digits = 4)

rownames(sqsemipcor.m.s1) = c("mood_av_fu",predvarnames.s1)
colnames(sqsemipcor.m.s1) = c("mood_av_fu",predvarnames.s1)


## average psychosocial functioning up to 1y
# Extract the variance covariance matrix of the latent and observed variables
cov_mat.psf <- data.frame(inspect(s1model.ml.psf_fu, what = "cov.all"))

cov_mat.psf.s1 = cov_mat.psf[c("psf_av_fu",predvarnames.s1),
                            c("psf_av_fu",predvarnames.s1)]


sqsemipcor.psf.s1 <- round(data.frame(correlation::cor_to_spcor(cor = as.matrix(cor_mat.psf.s1) ,cov = as.matrix(cov_mat.psf.s1)))^2,digits = 4)

rownames(sqsemipcor.psf.s1) = c("psf_av_fu",predvarnames.s1)
colnames(sqsemipcor.psf.s1) = c("psf_av_fu",predvarnames.s1)


## do tables
regpsf2tab = round(regpsfs2$tab, digits = 2)
regmoods2tab =round(regmoods2$tab, digits = 2)
regmoods2tab$df = NULL
regpsf2tab$df = NULL

regnsspc2 = cbind(regmoods2tab,SSPC = round(as.numeric(t(sqsemipcor.m.s1[1,2:15])), digits=3), regpsf2tab,SSPC = round(as.numeric(t(sqsemipcor.psf.s1[1,2:15])), digits=3))
rownames(regnsspc2) = c("N8:PF Problems", "N9:Som. Depr. (resid.)", "N2:Cogn. Depr. (resid.)", "N4:GAD & OCD (resid.)", "N1:SAD & Agor. (resid.)", "N3:Panic & Som. (resid.)", "N12:Spec. Phob. (resid.)", "N13:ED Thoughts (resid.)", "N6:ED Bulimic (resid.)", "N5:ED Anorectic (resid.)", "N7:Thought Dis. (resid.)", "N10:Detachment (resid.)", "N11:Antagonism (resid.)", "N14:Disinhibition (resid.)")

kable(regnsspc2,caption = "Regression based on correlation and covariance matrix from bifactor-s-1 model. Prediction of average psychosocial functioning and affective well-being in the first year of usage of the MindDoc App.") %>%
        kable_styling("striped") %>%
         add_header_above(c("   " = 1,"Prediction of average affective well-being over 1 year" = 7, "Prediction of average psychosocial functioning over 1 year" = 7)) %>% kable_paper %>% add_footnote(paste0("Table 4. ","N (affective well-being) = 10636",". N (psychosocial functioning) = 5342",". Total R2 affective well-being = ",round(regmoods2$R2,3),". Total R2 psychosocial functioning = ",round(regpsfs2$R2,3),". SSPC = squared semipartial correlation, i.e. the unique contribution of a specific predictor variable in a multiple regression analysis, while controlling for other predictors in the model. "))

Regression based on correlation and covariance matrix from bifactor-s-1 model. Prediction of average psychosocial functioning and affective well-being in the first year of usage of the MindDoc App.

	Prediction of average affective well-being over 1 year							Prediction of average psychosocial functioning over 1 year
	beta	se	tval	pval	ci.lb	ci.ub	SSPC	beta	se	tval	pval	ci.lb	ci.ub	SSPC
N8:PF Problems	-0.48	0.00	-104.57	0.00	-0.49	-0.47	0.232	0.58	0.00	132.36	0.00	0.57	0.58	0.331
N9:Som. Depr. (resid.)	-0.10	0.00	-20.93	0.00	-0.11	-0.09	0.009	0.29	0.00	62.55	0.00	0.28	0.30	0.074
N2:Cogn. Depr. (resid.)	-0.36	0.01	-63.79	0.00	-0.38	-0.35	0.086	0.21	0.01	38.46	0.00	0.20	0.22	0.028
N4:GAD & OCD (resid.)	0.06	0.01	11.33	0.00	0.05	0.07	0.003	-0.10	0.01	-20.14	0.00	-0.11	-0.09	0.008
N1:SAD & Agor. (resid.)	0.02	0.01	4.25	0.00	0.01	0.03	0.000	0.07	0.00	14.92	0.00	0.06	0.08	0.004
N3:Panic & Som. (resid.)	0.05	0.01	10.43	0.00	0.04	0.06	0.002	-0.02	0.00	-3.62	0.00	-0.03	-0.01	0.000
N12:Spec. Phob. (resid.)	-0.04	0.01	-6.43	0.00	-0.05	-0.02	0.001	0.03	0.01	5.13	0.00	0.02	0.04	0.000
N13:ED Thoughts (resid.)	0.06	0.01	8.51	0.00	0.04	0.07	0.002	0.02	0.01	2.56	0.01	0.00	0.03	0.000
N6:ED Bulimic (resid.)	-0.01	0.01	-2.16	0.03	-0.02	0.00	0.000	0.03	0.00	5.58	0.00	0.02	0.04	0.001
N5:ED Anorectic (resid.)	-0.07	0.01	-11.20	0.00	-0.08	-0.06	0.003	-0.02	0.01	-2.78	0.01	-0.03	0.00	0.000
N7:Thought Dis. (resid.)	-0.06	0.01	-10.95	0.00	-0.07	-0.05	0.002	-0.03	0.01	-4.83	0.00	-0.04	-0.02	0.000
N10:Detachment (resid.)	-0.12	0.01	-22.44	0.00	-0.13	-0.11	0.011	0.01	0.01	1.43	0.15	0.00	0.02	0.000
N11:Antagonism (resid.)	0.03	0.01	5.44	0.00	0.02	0.04	0.001	0.02	0.01	3.66	0.00	0.01	0.03	0.000
N14:Disinhibition (resid.)	0.08	0.01	14.38	0.00	0.07	0.09	0.004	-0.05	0.01	-8.63	0.00	-0.06	-0.04	0.001
a Table 4. N (affective well-being) = 10636. N (psychosocial functioning) = 5342. Total R2 affective well-being = 0.424. Total R2 psychosocial functioning = 0.487. SSPC = squared semipartial correlation, i.e. the unique contribution of a specific predictor variable in a multiple regression analysis, while controlling for other predictors in the model.

create table for paper

regnsspc.table = regnsspc1[,1:8]
names(regnsspc.table) = c("std. regr. weights", "SSPC","std. regr. weights", "SSPC","std. regr. weights", "SSPC","std. regr. weights", "SSPC" )
for (i in 1:nrow(regnsspc1)){
  regnsspc.table[i,3] = sub("0.", ".",paste0(regnsspc1[i,1], " [", sub("0.", ".",regnsspc1[i,5]), " , ", sub("0.", ".",regnsspc1[i,6]), "]"))
  regnsspc.table[i,1] = sub("0.", ".",paste0(regnsspc1[i,8], " [", sub("0.", ".",regnsspc1[i,12]), " , ", sub("0.", ".",regnsspc1[i,13]), "]"))
  regnsspc.table[i,4] = sub("0.", ".",paste0(regnsspc1[i,7]))
  regnsspc.table[i,2] = sub("0.", ".",paste0(regnsspc1[i,14]))
  regnsspc.table[i,7] = sub("0.", ".",paste0(regnsspc2[i,1], " [", sub("0.", ".",regnsspc2[i,5]), " , ", sub("0.", ".",regnsspc2[i,8]), "]"))
  regnsspc.table[i,5] = sub("0.", ".",paste0(regnsspc2[i,8], " [", sub("0.", ".",regnsspc2[i,12]), " , ", sub("0.", ".",regnsspc2[i,13]), "]"))
  regnsspc.table[i,8] = sub("0.", ".",paste0(regnsspc2[i,7]))
  regnsspc.table[i,6] = sub("0.", ".",paste0(regnsspc2[i,14]))


}
#
# regnsspc.table[i,1] = sub("0.", ".",paste0(regnsspc1[i,1], " [", sub("0.", ".",regnsspc1[i,5]), " , ", sub("0.", ".",regnsspc1[i,6]), "]"))
# regnsspc.table[i,3] = sub("0.", ".",paste0(regnsspc1[i,8], " [", sub("0.", ".",regnsspc1[i,12]), " , ", sub("0.", ".",regnsspc1[i,13]), "]"))
# regnsspc.table[i,2] = sub("0.", ".",paste0(regnsspc1[i,7]))
# regnsspc.table[i,4] = sub("0.", ".",paste0(regnsspc1[i,14]))
# regnsspc.table[i,5] = sub("0.", ".",paste0(regnsspc2[i,1], " [", sub("0.", ".",regnsspc2[i,5]), " , ", sub("0.", ".",regnsspc2[i,6]), "]"))
# regnsspc.table[i,7] = sub("0.", ".",paste0(regnsspc2[i,8], " [", sub("0.", ".",regnsspc2[i,12]), " , ", sub("0.", ".",regnsspc2[i,13]), "]"))
# regnsspc.table[i,6] = sub("0.", ".",paste0(regnsspc2[i,7]))
# regnsspc.table[i,8] = sub("0.", ".",paste0(regnsspc2[i,14]))

# # first year
# add_footnote(paste0("Table 4. ","N (affective well-being) = ",sum(!is.na(data_scale_pred$mood_av_fu)) ,". N (psychosocial functioning) = ", sum(!is.na(data_scale_pred$psf_av_fu)) ,". Total R2 affective well-being = ",round(regmoods1$R2,3),". Total R2 psychosocial functioning = ",round(regpsfs1$R2,3),". SSPC = squared semipartial correlation, i.e. the unique contribution of a specific predictor variable in a multiple regression analysis, while controlling for other predictors in the model. "))
#
# #10M
# add_footnote(paste0("Table 5. ","N (affective well-being) = ",sum(!is.na(data_scale_pred$mood_av)) ,". N (psychosocial functioning) = ", sum(!is.na(data_scale_pred$psf)) ,". Total R2 affective well-being = ",round(regmoods2$R2,3),". Total R2 psychosocial functioning = ",round(regpsfs2$R2,3),". SSPC = squared semipartial correlation, i.e. the unique contribution of a specific predictor variable in a multiple regression analysis, while controlling for other predictors in the model. "))


kable(regnsspc.table,caption = "Table 3. Multiple regression and semipartial correlations based on covariance matrix from bifactor-(S-1) model.") %>%
        kable_styling("striped") %>%
         add_header_above(c("   " = 1,"Psychosocial impairments" = 2, "Affective well-being" = 2,"Psychosocial impairments" = 2,"Affective well-being" = 2)) %>%
   add_header_above(c("   " = 1, "Prediction 1: Year one" = 4,"Prediction 2: Year two" = 4 ))%>%
  kable_paper %>% add_footnote(paste0("Note. N (affective well-being prediction 1) = 25844",". N (affective well-being prediction 2) = 10636",". N (psychosocial functioning prediction 1) = 27173", ". N (psychosocial functioning prediction 2) = 5342" ,". Total R2 affective well-being (prediction 1) = ",round(regmoods1$R2,3),". Total R2 psychosocial functioning (prediction 1) = ",round(regpsfs1$R2,3),". Total R2 affective well-being (prediction 2) = ",round(regmoods2$R2,3),". Total R2 psychosocial functioning (prediction 2) = ",round(regpsfs2$R2,3),". SSPC = squared semipartial correlation, i.e. the unique contribution of a specific predictor variable in a multiple regression analysis, while controlling for other predictors in the model. "))

Table 3. Multiple regression and semipartial correlations based on covariance matrix from bifactor-(S-1) model.

	Prediction 1: Year one				Prediction 2: Year two
	Psychosocial impairments		Affective well-being		Psychosocial impairments		Affective well-being
	std. regr. weights	SSPC	std. regr. weights	SSPC	std. regr. weights	SSPC	std. regr. weights	SSPC
N8:PF Problems	.77 [.76 , .77]	.586	-.55 [-.56 , -.55]	.306	.58 [.57 , .58]	.331	-.48 [-.49 , .58]	.232
N9:Som. Depr. (resid.)	.34 [.33 , .34]	.1	-.12 [-.13 , -.11]	.013	.29 [.28 , .3]	.074	-.1 [-.11 , .29]	.009
N2:Cogn. Depr. (resid.)	.34 [.34 , .35]	.077	-.52 [-.53 , -.51]	.178	.21 [.2 , .22]	.028	-.36 [-.38 , .21]	.086
N4:GAD & OCD (resid.)	-.02 [-.02 , -.01]	0	.01 [0 , .02]	0	-.1 [-.11 , -.09]	.008	.06 [.05 , -.1]	.003
N1:SAD & Agor. (resid.)	.08 [.08 , .09]	.005	.01 [0 , .02]	0	.07 [.06 , .08]	.004	.02 [.01 , .07]	0
N3:Panic & Som. (resid.)	.01 [.01 , .02]	0	.[-.01 , .01]	0	-.02 [-.03 , -.01]	0	.05 [.04 , -.02]	.002
N12:Spec. Phob. (resid.)	.[-.01 , 0]	0	-.02 [-.03 , -.01]	0	.03 [.02 , .04]	0	-.04 [-.05 , .03]	.001
N13:ED Thoughts (resid.)	.[0 , .01]	0	.08 [.07 , .09]	.003	.02 [0 , .03]	0	.06 [.04 , .02]	.002
N6:ED Bulimic (resid.)	.01 [.01 , .02]	0	.[-.01 , .01]	0	.03 [.02 , .04]	.001	-.01 [-.02 , .03]	0
N5:ED Anorectic (resid.)	-.01 [-.02 , 0]	0	-.07 [-.08 , -.06]	.002	-.02 [-.03 , 0]	0	-.07 [-.08 , -.02]	.003
N7:Thought Dis. (resid.)	-.01 [-.02 , 0]	0	-.1 [-.11 , -.1]	.007	-.03 [-.04 , -.02]	0	-.06 [-.07 , -.03]	.002
N10:Detachment (resid.)	.03 [.02 , .03]	0	-.14 [-.15 , -.14]	.015	.01 [0 , .02]	0	-.12 [-.13 , .01]	.011
N11:Antagonism (resid.)	.02 [.02 , .03]	0	.02 [.01 , .03]	0	.02 [.01 , .03]	0	.03 [.02 , .02]	.001
N14:Disinhibition (resid.)	.01 [.01 , .02]	0	.03 [.03 , .04]	.001	-.05 [-.06 , -.04]	.001	.08 [.07 , -.05]	.004
a Note. N (affective well-being prediction 1) = 25844. N (affective well-being prediction 2) = 10636. N (psychosocial functioning prediction 1) = 27173. N (psychosocial functioning prediction 2) = 5342. Total R2 affective well-being (prediction 1) = 0.614. Total R2 psychosocial functioning (prediction 1) = 0.847. Total R2 affective well-being (prediction 2) = 0.424. Total R2 psychosocial functioning (prediction 2) = 0.487. SSPC = squared semipartial correlation, i.e. the unique contribution of a specific predictor variable in a multiple regression analysis, while controlling for other predictors in the model.