RSA (version 0.5.1)

RSA: Performs several RSA model tests on a data set with two predictors

Description

Performs several RSA model tests on a data set with two predictors

Main functions are: RSA, demoRSA, plotRSA, RSA.ST, compare.

Usage

RSA(formula, data = NULL, sample.cor = NULL,
    sample.nobs = NULL, center = FALSE, scale = FALSE,
    na.rm = FALSE, out.rm = TRUE, breakline = FALSE,
    models = "default", verbose = TRUE, ...)

Arguments

formula

A formula in the form z ~ x*y, specifying the variable names used from the data frame, where z is the name of the response variable, and x and y are the names of the predictor variables.

data

A data frame with the variables

sample.cor

A matrix with the sample correlation matrix - you have to provide either raw data in the data argument, or a sample.cor and a sample.nobs

sample.nobs

Number of observations in the sample (goes along with sample.cor)

center

Should predictor variables be centered on the sample mean before analyses?

scale

Should predictor variables be scaled to SD = 1 before analyses?

na.rm

Remove missings before proceeding?

out.rm

Should outliers according to Bollen & Jackman (1980) criteria be excluded from analyses?

breakline

Should the breakline in the unconstrained absolute difference model be allowed (the breakline is possible from the model formulation, but empirically rather unrealistic ...)

verbose

Should additional information during the computation process be printed?

models

A vector with names of all models that should be computed. Should be any from c("absdiff", "absunc", "diff", "additive", "IA", "sqdiff", "SSD", "SRSD", "full"). For models="all", all models are computed, for models="default"

...

Additional parameters passed to the lavaan sem function. For example: se="boot"

Details

No details so far.

References

...

See Also

demoRSA, plotRSA, RSA.ST

Examples

# Compute response surface from a fake data set
set.seed(0xBEEF)
n <- 300
err <- 15
x <- rnorm(n, 0, 5)
y <- rnorm(n, 0, 5)
df <- data.frame(x, y)
df <- within(df, {
	diff <- x-y
	absdiff <- abs(x-y)
	sqdiff <- (x-y)^2
	z.diff <- diff + rnorm(n, 0, err)
	z.abs <- absdiff + rnorm(n, 0, err)
	z.sq <- sqdiff + rnorm(n, 0, err)
	z.add <- diff + 0.4*x + rnorm(n, 0, err)
	z.complex <- 0.4*x + - 0.2*x*y + + 0.1*x^2 - 0.03*y^2 + rnorm(n, 0, err)
})

r1 <- RSA(z.sq~x*y, df)
print(r1)
compare(r1)
plot(r1)
plot(r1, model="SRSD")
plot(r1, model="full", type="c")
getPar(r1, "coef")	# print model parameters including SE and CI
RSA.ST(r1)	# get surface parameters

# Motive congruency example
data(motcon)
r.m <- RSA(negAct~EM*IM, motcon)