| Title: | R/EQS Interface |
|---|---|
| Description: | Contains the function run.eqs() which calls an EQS script file, executes the EQS estimation, and, finally, imports the results as R objects. These two steps can be performed separately: call.eqs() calls and executes EQS, whereas read.eqs() imports existing EQS outputs as objects into R. It requires EQS 6.2 (build 98 or higher). |
| Authors: | Patrick Mair [aut, cre], Eric Wu [aut] |
| Maintainer: | Patrick Mair <[email protected]> |
| License: | GPL-3 |
| Version: | 0.8-13 |
| Built: | 2024-10-25 02:43:38 UTC |
| Source: | https://github.com/cran/REQS |
This function calls an EQS script file (.eqs) and executes it.
call.eqs(EQSpgm, EQSmodel, serial, Rmatrix = NA, datname = NA, LEN = 2000000)call.eqs(EQSpgm, EQSmodel, serial, Rmatrix = NA, datname = NA, LEN = 2000000)
EQSpgm |
String containing path including program name where EQS is located (see details) |
EQSmodel |
String containing path where .eqs script file is located (see details) |
serial |
EQS serial number as character |
Rmatrix |
Optional matrix argument if data or covariances are stored in R |
datname |
If |
LEN |
Integer containing number of working array units. By default, it is 2000000 8 bytes units |
If the path in EQSpgm and EQSmodel contains a blank, single quotes and double quotes
are required in argument. See EQSpgm argument in examples. The last statement in the EQSpgm argument refers
to the name of the executable program file. Under Windows it is ".../WINEQS" (referring to WINEQS.exe), under Mac ".../MACEQS" and
under Linux ".../EQS". When specifying the path, use slash instead of backslash.
The .ETS, .CBK and .ETP files are written in the directory where the .eqs file is located. Unless another path is provided within in the .eqs script file.
The argument datname must match with the input data specified in the corresponding .eqs file.
Returns TRUE is the estimation was succesfully and FALSE otherwise.
Patrick Mair, Eric Wu
Bentler, P. M. (1995). EQS Program Manual. Encino, CA: Multivariate Software Inc.
## Not run: ##not executable, valid serial number has to be provided res <- call.eqs(EQSpgm = "C:/Program Files/EQS61/WINEQS.EXE", EQSmodel = "c:/eqs61/examples/manul7.eqs", serial = "1234") ## End(Not run)## Not run: ##not executable, valid serial number has to be provided res <- call.eqs(EQSpgm = "C:/Program Files/EQS61/WINEQS.EXE", EQSmodel = "c:/eqs61/examples/manul7.eqs", serial = "1234") ## End(Not run)
This function reads EQS output files (.ets, .CBK and .ETP) into R and stores the results as objects.
read.eqs(file)read.eqs(file)
file |
The name (string) of the .ets file or the full path which the data are to be read from. If it does not contain an absolute path, the file name is relative to the current working directory, 'getwd()'. A .CBK and .ETP file have to be of the same name and in the same directory. |
The value list below provides objects for the full EQS output. If in EQS some objects are not computed, the corresponding values in R are NA.
Returns a list with the following objects:
model.info |
General model information |
pval |
p-values for various test statistics |
fit.indices |
Variuos fit indices |
model.desc |
Descriptive measures |
Phi |
Phi matrix |
Gamma |
Gamma matrix |
Beta |
Beta matrix |
par.table |
Parameter table (with standard errors) |
sample.cov |
Sample covariance matrix |
sigma.hat |
Model covariance matrix |
inv.infmat |
Inverse information matrix |
rinv.infmat |
Robust inverse information matrix |
cinv.infmat |
Corrected inverse information matrix |
derivatives |
First derivatives |
moment4 |
Matrix with 4th moments |
ssolution |
Standardized elements |
Rsquared |
R-squared measures |
fac.means |
Factor means |
var.desc |
Descriptive measures for the variables (univariate statistics) |
indstd |
Independent variable standardization vector |
depstd |
Dependent variable standardization vector |
Patrick Mair, Eric Wu
Bentler, P. M. (2008). EQS Program Manual. Encino, CA: Multivariate Software Inc.
## Not run: ##not executable eqsout <- read.eqs("c:/home/user/eqs/eqsfile.ets") ## End(Not run)## Not run: ##not executable eqsout <- read.eqs("c:/home/user/eqs/eqsfile.ets") ## End(Not run)
Calls an EQS script file from R, executes EQS, and imports the results into R. Basically
it is a wrapper function of call.eqs and the subsequent read.eqs.
run.eqs(EQSpgm, EQSmodel, serial, Rmatrix = NA, datname = NA, LEN = 2000000)run.eqs(EQSpgm, EQSmodel, serial, Rmatrix = NA, datname = NA, LEN = 2000000)
EQSpgm |
String containing path where EQS is located (see details) |
EQSmodel |
String containing path where .eqs script file is located (see details) |
serial |
EQS serial number as integer value |
Rmatrix |
Optional matrix argument if data or covariances are stored in R |
datname |
If |
LEN |
Integer containing number of working array units. By default, it is 2000000 8 bytes units |
If the path in EQSpgm and EQSmodel contains a blank, single quotes and double quotes
are required in argument. See EQSpgm argument in examples. The last statement in the EQSpgm argument refers
to the name of the executable program file. Under Windows it is ".../WINEQS" (referring to WINEQS.exe), under Mac ".../MACEQS" and
under Linux ".../EQS". When specifying the path, use slash instead of backslash.
The .ETS, .CBK and .ETP files are written in the directory where the .eqs file is located. Note that these 3 files must be in the same directory than the .eqs file.
The argument datname must match with the input data specified in the corresponding .eqs file.
This option can be used for simulations: Generate data in R, run.eqs() on with the corresponding
data argument, pick out the relevant return values.
The value list below provides objects for the full EQS output. If in EQS some objects are not computed, the corresponding values in R are NA.
Returns a list with the following objects:
success |
|
model.info |
General model information |
pval |
p-values for various test statistics |
fit.indices |
Variuos fit indices |
model.desc |
Descriptive measures |
Phi |
Phi matrix |
Gamma |
Gamma matrix |
Beta |
Beta matrix |
par.table |
Parameter table (with standard errors) |
sample.cov |
Sample covariance matrix |
sigma.hat |
Model covariance matrix |
inv.infmat |
Inverse information matrix |
rinv.infmat |
Robust inverse information matrix |
cinv.infmat |
Corrected inverse information matrix |
derivatives |
First derivatives |
moment4 |
Matrix with 4th moments |
ssolution |
Standardized elements |
Rsquared |
R-squared measures |
fac.means |
Factor means |
var.desc |
Descriptive measures for the variables (univariate statistics) |
indstd |
Independent variable standardization vector |
depstd |
Dependent variable standardization vector |
Patrick Mair, Eric Wu
Bentler, P. M. (1995). EQS Program Manual. Encino, CA: Multivariate Software Inc.
## Not run: ##not executable, valid serial number has to be provided res <- run.eqs(EQSpgm = "C:/Program Files/EQS61/WINEQS.EXE", EQSmodel = "c:/eqs61/examples/manul7.eqs", serial = "1234") ##For instance, to extract the parameter table you can do res$par.table ##simulation example: not executable, provide serial number and proper eqs script file ##simulated 100 replications, extract CFI cfivec <- NULL for (i in 1:100) { X <- matrix(rnorm(1000), ncol = 10, nrow = 100) res <- run.eqs(EQSpgm = "C:/Program Files/EQS61/WINEQS.EXE", EQSmodel = "c:/eqs61/examples/manul7.eqs", data = X, datname = "manul7.dat", serial = "1234") cfivec <- c(cfivec, res.run$fit.indices[9,]) } ## End(Not run)## Not run: ##not executable, valid serial number has to be provided res <- run.eqs(EQSpgm = "C:/Program Files/EQS61/WINEQS.EXE", EQSmodel = "c:/eqs61/examples/manul7.eqs", serial = "1234") ##For instance, to extract the parameter table you can do res$par.table ##simulation example: not executable, provide serial number and proper eqs script file ##simulated 100 replications, extract CFI cfivec <- NULL for (i in 1:100) { X <- matrix(rnorm(1000), ncol = 10, nrow = 100) res <- run.eqs(EQSpgm = "C:/Program Files/EQS61/WINEQS.EXE", EQSmodel = "c:/eqs61/examples/manul7.eqs", data = X, datname = "manul7.dat", serial = "1234") cfivec <- c(cfivec, res.run$fit.indices[9,]) } ## End(Not run)