best.rsa2 {RSAtools} | R Documentation |
Compare two polynomial models against the data
Description
Compare two polynomial models, for example to test parametric constraints in STEP2 of the 3-step identification strategy (see RSAmodel
).
Usage
best.rsa2(RSA_object, m1, m2, order = c("wAIC"), robust = TRUE)
Arguments
RSA_object |
x an object of class "RSA_object" generated by RSAmodel() |
m1 |
First model to be compared (contained in RSA_object$models) |
m2 |
Second model to be compared (contained in RSA_object$models) |
order |
Fit index used to determine best-fitting model. |
robust |
A boolean stating whether robust fit indices should be extracted (default= TRUE) |
Value
A table containing fit indices for each model
Examples
##### Test a variant within a family (e.g., FM26_PARALLELASYMWEAK)
##Define variant as constraints
list_variant <- list()
list_variant[["variant1"]] <- c('
####First-order polynomials: variant-specific
b1 == -1/4*b2
####Second-order polynomials: variant-specific
b5 == b3/2
b4 == 0
####Third-order polynomials: FM26_PARALLELASYMWEAK
b6 == 0
b7 == 0
b9 == b8/-3
')
RSA_NSfit <- RSAmodel(formula= engagement ~ needs*supplies,
data= sim_NSfit, model= c("FM26_PARALLELASYMWEAK","USER"),
user_model= list_variant)
##Compare variant to best-fitting family (e.g., LRT_pvalue p > .05)
best.rsa2(RSA_NSfit,m1="variant1",m2="FM26_PARALLELASYMWEAK")[2,1:3]
[Package RSAtools version 0.1.1 Index]