Load Libraries

library(foreign)
library(tidyLPA)
library(dplyr)

Load Data

data <- read.spss("data.sav", to.data.frame = T)
data.sav: Very long string record(s) found (record type 7, subtype 14), each will be imported in consecutive separate variables

Find Best Fitting Model

It appears as though model 2 with 2 profiles fits best (as indicated by low BIC). Model 3 with 3 profiles has the lowest BIC but I’m suspicious of it due to large variance.

Model 2 with 2 Profiles

model2 <- estimate_profiles(data, 
                        RTA_factor_1, RTA_factor_2, RTA_factor_3, RTA_factor_4,
                        model = 2,
                        n_profiles = 2) 
Fit varying means, equal variances and covariances (Model 2) model with 2 profiles.
LogLik is 5493.057
BIC is 11094.924
Entropy is 0.925

Entropy is above .8 suggesting distinct profiles.

Estimated Means and Variance for Model 2

y-axis = mean

y-axis = mean centered

y-axis = standard deviation

Model 3 with 3 Profiles

model3 <- estimate_profiles(data, 
                        RTA_factor_1, RTA_factor_2, RTA_factor_3, RTA_factor_4,
                        model = 3,
                        n_profiles = 3) 
Fit varying means and variances, covariances fixed to 0 (Model 3) model with 3 profiles.
LogLik is 5403.342
BIC is 10955.582
Entropy is 0.928

y-axis = mean

y-axis = mean centered

y-axis = standard deviation

LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OgogIGh0bWxfbm90ZWJvb2s6IGRlZmF1bHQKICBodG1sX2RvY3VtZW50OiBkZWZhdWx0CiAgcGRmX2RvY3VtZW50OiBkZWZhdWx0CiAgd29yZF9kb2N1bWVudDogZGVmYXVsdAotLS0KCiMgTG9hZCBMaWJyYXJpZXMKYGBge3J9CmxpYnJhcnkoZm9yZWlnbikKbGlicmFyeSh0aWR5TFBBKQpsaWJyYXJ5KGRwbHlyKQpgYGAKCiMgTG9hZCBEYXRhCmBgYHtyfQpkYXRhIDwtIHJlYWQuc3BzcygiZGF0YS5zYXYiLCB0by5kYXRhLmZyYW1lID0gVCkKCmRhdGEgPC0gc2VsZWN0KGRhdGEsIEFidXNlLCBSVEFfZmFjdG9yXzEsIFJUQV9mYWN0b3JfMiwKICAgICAgICAgICAgICAgUlRBX2ZhY3Rvcl8zLCBSVEFfZmFjdG9yXzQsIFNleFZpY3RpbWl6YXRpb25fbm8xLAogICAgICAgICAgICAgICBQaHlzaWNhbF92aWN0aW1pemF0aW9uLCBTZXhWaWN0aW1pemF0aW9uX2FnZ19ubzEsCiAgICAgICAgICAgICAgIFBoeXNpY2FsX0FnZ3Jlc3Npb24pCgpkYXRhIDwtIGZpbHRlcihkYXRhLCBBYnVzZSA9PSAicmVwb3J0ZWQgYW4gYWJ1c2UiKQoKZGF0YSRpZCA8LSAxOm5yb3coZGF0YSkKYGBgCgojIEZpbmQgQmVzdCBGaXR0aW5nIE1vZGVsCmBgYHtyfQpjb21wYXJlX3NvbHV0aW9ucyhkYXRhLCBSVEFfZmFjdG9yXzEsIFJUQV9mYWN0b3JfMiwgUlRBX2ZhY3Rvcl8zLCBSVEFfZmFjdG9yXzQpCmBgYAoKSXQgYXBwZWFycyBhcyB0aG91Z2ggbW9kZWwgMiB3aXRoIDIgcHJvZmlsZXMgZml0cyBiZXN0IChhcyBpbmRpY2F0ZWQgYnkgbG93IEJJQykuIE1vZGVsIDMgd2l0aCAzIHByb2ZpbGVzIGhhcyB0aGUgbG93ZXN0IEJJQyBidXQgSSdtIHN1c3BpY2lvdXMgb2YgaXQgZHVlIHRvIGxhcmdlIHZhcmlhbmNlLgoKIyBNb2RlbCAyIHdpdGggMiBQcm9maWxlcwpgYGB7cn0KbW9kZWwyIDwtIGVzdGltYXRlX3Byb2ZpbGVzKGRhdGEsIAogICAgICAgICAgICAgICAgICAgICAgICBSVEFfZmFjdG9yXzEsIFJUQV9mYWN0b3JfMiwgUlRBX2ZhY3Rvcl8zLCBSVEFfZmFjdG9yXzQsCiAgICAgICAgICAgICAgICAgICAgICAgIG1vZGVsID0gMiwKICAgICAgICAgICAgICAgICAgICAgICAgbl9wcm9maWxlcyA9IDIsCiAgICAgICAgICAgICAgICAgICAgICAgIHJldHVybl9vcmlnX2RmID0gVFJVRSkKYGBgCgpFbnRyb3B5IGlzIGFib3ZlIC44IHN1Z2dlc3RpbmcgZGlzdGluY3QgcHJvZmlsZXMuIAoKIyBFc3RpbWF0ZWQgTWVhbnMgYW5kIFZhcmlhbmNlIGZvciBNb2RlbCAyCgojIyB5LWF4aXMgPSBtZWFuCmBgYHtyfQpwbG90X3Byb2ZpbGVzKG1vZGVsMikgI3RvX2NlbnRlciA9IFRSVUUsIHRvX3NjYWxlID0gVFJVRQpgYGAKCiMjIHktYXhpcyA9IG1lYW4gY2VudGVyZWQKYGBge3J9CnBsb3RfcHJvZmlsZXMobW9kZWwyLCB0b19jZW50ZXIgPSBUKSAjdG9fY2VudGVyID0gVFJVRSwgdG9fc2NhbGUgPSBUUlVFCmBgYAoKIyMgeS1heGlzID0gc3RhbmRhcmQgZGV2aWF0aW9uCmBgYHtyfQpwbG90X3Byb2ZpbGVzKG1vZGVsMiwgdG9fY2VudGVyID0gVCwgdG9fc2NhbGUgPSBUKSAjdG9fY2VudGVyID0gVFJVRSwgdG9fc2NhbGUgPSBUUlVFCmBgYAoKIyBNb2RlbCAzIHdpdGggMyBQcm9maWxlcwpgYGB7cn0KbW9kZWwzIDwtIGVzdGltYXRlX3Byb2ZpbGVzKGRhdGEsIAogICAgICAgICAgICAgICAgICAgICAgICBSVEFfZmFjdG9yXzEsIFJUQV9mYWN0b3JfMiwgUlRBX2ZhY3Rvcl8zLCBSVEFfZmFjdG9yXzQsCiAgICAgICAgICAgICAgICAgICAgICAgIG1vZGVsID0gMywKICAgICAgICAgICAgICAgICAgICAgICAgbl9wcm9maWxlcyA9IDMpIApgYGAKCiMjIHktYXhpcyA9IG1lYW4KYGBge3J9CnBsb3RfcHJvZmlsZXMobW9kZWwzKSAjdG9fY2VudGVyID0gVFJVRSwgdG9fc2NhbGUgPSBUUlVFCmBgYAoKIyMgeS1heGlzID0gbWVhbiBjZW50ZXJlZApgYGB7cn0KcGxvdF9wcm9maWxlcyhtb2RlbDMsIHRvX2NlbnRlciA9IFQpICN0b19jZW50ZXIgPSBUUlVFLCB0b19zY2FsZSA9IFRSVUUKYGBgCgojIyB5LWF4aXMgPSBzdGFuZGFyZCBkZXZpYXRpb24KYGBge3J9CnBsb3RfcHJvZmlsZXMobW9kZWwzLCB0b19jZW50ZXIgPSBULCB0b19zY2FsZSA9IFQpICN0b19jZW50ZXIgPSBUUlVFLCB0b19zY2FsZSA9IFRSVUUKYGBgCgo=