SSI has enjoyed great success over the years in the development and publishing of statistical software and is proud to announce the release of LISREL 10.1.
In an effort to meet the growing demands of our LISREL 8 and 9 user community, SSI has developed LISREL 10, which is on the cutting edge of current technology. The program has been tested extensively on the Microsoft Windows platform with Windows 7 and Windows 10 operating systems.
Background
Modern structural equation modeling is based on raw data. With LISREL 10, if raw data is available in a LISREL data system file or in a text file, one can read the data into LISREL and formulate the model using either SIMPLIS syntax or LISREL syntax.
LISREL 10 contains fixes to all bugs reported by users of LISREL 9. The new LISREL features are summarized next.
Multiple group analyses using a single data file
In practice, many multivariate data sets are observations from several groups. Examples of these groups are genders, languages, political parties, countries, faculties, colleges, schools, etc.
LISREL may be used to fit multiple group structural equation models to multiple group data. Traditional statistical methods such as Maximum Likelihood (ML), Robust Maximum Likelihood (RML), Weighted Least Squares (WLS), Diagonally Weighted Least Squares (DWLS), Generalized Least Squares (GLS) and Un-weighted Least Squares (ULS) are available for complete multiple group data while the Full Information Maximum Likelihood (FIML) method is available for incomplete multiple group data.
In previous versions of LISREL, the user was required to create separate data files for each group. Suppose that the groups to be analyzed consisted of data collected in eight countries, the implication is that eight datasets must to be created in order to fit a multiple group structural equation model.
A new feature implemented in LISREL 10 allows researchers to use a single dataset that contains a group variable that can be defined by
- Using the Data menu when a LISREL system file (.lsf) is opened
- By inserting the line $GROUPS=<group variable name> anywhere in the syntax file.
Consider the dataset efficacy_4countries.lsf shown below. There are 4 countries and portion of the data from countries 2 and 3 are shown below.
To use this dataset in a multiple group analysis, use the Data menu from the main menu bar and select the Group Variable… option (see below)
Select COUNTRY from the list of variables and when done click the OK button.
The LISREL Examples folder contains a sub-folder named MGROUPS that contains examples for the following statistical procedures:
For a detailed example, see the Assessment of Invariance, (Section 2 in the “Additional Topics Guide.pdf”) that can be accessed via the Help option on the main menu bar:
Models for grouped- and discrete-tie survival data
For data that are clustered and/or repeated, models including random effects provide a convenient way of accounting for association in correlated survival data.
Several authors have noted the relationship between ordinal regression models (using complementary log-log and logistic link functions) and survival analysis models for grouped and discrete time. In LISREL 10 a generalization of an ordinal random-effects regression model to handle correlated grouped-time survival data is implemented. This model accommodates multivariate normally-distributed random effects, and additionally, allows for a general form for model covariates.
Assuming a proportional or partial proportional, hazards or odds model, a maximum marginal likelihood solution is implemented using multi-dimensional quadrature to numerically integrate over the distribution of random-effects. The reference guide “Survival Models for grouped data.pdf” contains examples and references and is accessible via the online Help menu.
Models for ordinal outcomes and the proportional odds versus non-proportional odds assumption
Extensive work on the development of methods for the analysis of ordinal response data has been undertaken by numerous researchers. These developments have focused on the extension of methods for dichotomous variables to ordinal response data, and have been mainly in terms of logistic and probit regression models. The proportional odds model is a common choice for the analysis of ordinal data. In LISREL 10, it is possible to fit both proportional and non-proportional odds models to verify the proportional odds assumption using a chi-square difference test. The reference guide “Models for proportional and non-proportional odds.pdf” contains examples and references and is accessible via the online Help menu.
P-value for C1 statistic under normality
LISREL and PRELIS functionality
LISREL 10 estimates the d eigenvalues and the p-value of the linear combination. Essentially, this makes the other chi-squares C2, C3, C4 and C5 less important, since they are based on approximations of the distribution using only the mean and variance of the eigenvalues.
Data conversion using stat/transfer
Selection of new features available Stat/Transfer Version 14
Version 14 has added support for the following formats:
- Stata 15/MP
- BayesiaLab (Write Only)
- JSON-Stat (Read Only)Version 14 has larger limits for:
- Excel Files > 4 GB
- SAS files > 32K variables
- Stata Files > 32K variables
- dBASE > 2GB
FIML for ordinal and continuous variables
The LISREL implementation allows for the use of design weights to fit SEM models to a mixture of continuous and ordinal manifest variables with or without missing values with optional specification of stratum and/or cluster variables. It also deals with the issue of robust standard error estimation and the adjustment of the chi-square goodness of fit statistic.
This method is based on adaptive quadrature and a user can specify any one of the following four link functions:
- Logit
- Probit
- Complementary Log-log
- Log-Log
Examples to illustrate this feature are given in the folders \orfimlex and \ls9ex.
Three-level Multilevel Generalized Linear Models
Cluster or multi-stage samples designs are frequently used for populations with an inherent hierarchical structure. Ignoring the hierarchical structure of data has serious implications. The use of alternatives such as aggregation and disaggregation of information to another level can induce an increase in co-linearity among predictors and large or biased standard errors for the estimates.
The collection of models called Generalized Linear Models (GLIMs) have become important, and practical, statistical tools. The basic idea of GLIMs is an adaption of standard regression to quite different kinds of data. The variables may be dichotomous, ordinal (as with a 5-point Likert scale), counts (number of arrest records), or nominal. The motivation is to tailor the regression relationship connecting the outcome to relevant independent variables so that it is appropriate to the properties of the dependent variable. The statistical theory and methods for fitting Generalized Linear Models (GLIMs) to survey data was implemented in LISREL 8.8.
Researchers from the social and economic sciences are often applying these methods to multilevel data and consequently, inappropriate results are obtained. The LISREL statistical module for the analysis of multilevel data allows for design weights. Two estimation methods, MAP (maximization of the posterior distribution) and QUAD (adaptive quadrature) for fitting generalized linear models to multilevel data are available. The LISREL Multilevel Generalized Linear models module (MGLIM) allows for a wide variety of sampling distributions and link functions.
The LISREL 10 MGLIM module also include zero-inflated Poisson and zero-inflated Negative-Binomial models and prints results for unit-specific and population-average estimates of the fixed effects.
Examples in the folder \mglimex illustrate these features.
Four and Five-level Multilevel Models for continuous outcome variables
Social science research often entails the analysis of data with a hierarchical structure. A frequently cited example of multilevel data is a dataset containing measurements on children nested within schools, with schools nested within education departments.
The need for statistical models that take account of the sampling scheme is well recognized and it has been shown that the analysis of survey data under the assumption of a simple random sampling scheme may give rise to misleading results.
Multilevel models are particularly useful in the modeling of data from complex surveys. Cluster or multi-stage samples designs are frequently used for populations with an inherent hierarchical structure. Ignoring the hierarchical structure of data has serious implications. The use of alternatives such as aggregation and disaggregation of information to another level can induce an increase in co-linearity among predictors and large or biased standard errors for the estimates. In order to address concerns regarding the appropriate analyses of survey data, the LISREL multilevel module for continuous data now also handles up to five levels and features an option for users to include design weights on levels 1, 2 , 3, 4 or 5 of the hierarchy.
Examples are given in the \mlevelex folder.
Filename extensions
To ensure backwards compatibility, users can still run previously created syntax files using a .psf file, but to open an existing .psf file using the graphical user’s interface, the user has to rename it to .lsf.
MVABOOK examples
This book can be used by Master and PhD students and researchers in the economic, social, behavioral, and many other sciences that need to have a basic understanding of multivariate statistical theory and methods for their analysis of multivariate data. It can also be used as a text book for courses on multivariate statistical analysis. All examples are listed in the Table of Contents. All the syntax and data files for these examples are distributed with LISREL 10.1 and are located in
- LISREL Examples\MVABOOK\CHAPTER1
- LISREL Examples\MVABOOK\CHAPTER2
- LISREL Examples\MVABOOK\CHAPTER3
- LISREL Examples\MVABOOK\CHAPTER4
- LISREL Examples\MVABOOK\CHAPTER5
- LISREL Examples\MVABOOK\CHAPTER6
- LISREL Examples\MVABOOK\CHAPTER7
- LISREL Examples\MVABOOK\CHAPTER8
- LISREL Examples\MVABOOK\CHAPTER9
- LISREL Examples\MVABOOK\CHAPTER10
Running LISREL in batch mode
"c:\program files (x86)\LISREL10\MLISREL64_10"
where
Program name = LISREL, PRELIS, MULTILEV, MAPGLIM or SURVEYGLIM
Example:
Syntax File = "c:\LISREL Examples\ls9ex\npv1a.spl"
Output File = "c:\LISREL Examples\ls9ex\npv1a.out"
Examples of batch files (RunLISREL.bat and RunSIMPLIS.bat) are given in the \ls9ex folder. These batch files will run all the LISREL and SIMPLIS syntax files in this folder.
Documentation
A list of PDF guides, accessible via the online Help menu is given below.
- New features in LISREL
- Graphical User’s Interface
- PRELIS Examples Guide
- LISREL Examples Guide
- Multilevel Modeling Guide
- Complex Survey Sampling
- Generalized Linear Modeling Guide
- Multilevel Generalized Linear Modeling Guide
- Models for Proportional and Non-proportional Odds
- Survival Models for Grouped Data
- LISREL Syntax Guide
- SIMPLIS Syntax Guide
- PRELIS Syntax Guide
- Additional Topics Guide
Documentation of the LISREL graphical user’s interface is also available as an online Help file. The Help file has features that simplify navigation across topics:
The Complex Survey Sampling Guide includes structural equation modeling (SEM) for continuous variables and SEM for a mixture of ordinal and continuous variables. LISREL uses full information maximum likelihood under complex survey data with data missing at random.
The Additional Topics Guide includes sections on assessment of invariance, multiple imputation, multilevel structural equation modeling and multilevel non-linear regression.
Examples
A selection of examples, illustrating some of the new features is given below.
Analysis of ordinal data using quadrature (\ls9ex)
Path Diagram Representation
Path Diagram (Standardized Solution)
Portion of output file
The last part of the output file is shown below.
For the moment we note the value of the deviance statistic −2 ln L = 19858.06. Since there is no value of −2 ln L for a saturated model, it is impossible to say whether this is large or small in some absolute sense. The deviance statistic can therefore only be used to compare different models for the same data.
To illustrate, the difference between the deviance statistic for this model and the deviance statistic for a model with one latent variable (Efficacy2a.spl) is 19934.57-19858.06 =76.51, which suggests that the two-dimensional model fits the data much better than the uni-dimensional model.
The output also gives estimates of the thresholds, their standard errors and z-values. The thresholds are parameters of the model but are seldom useful in analysis of a single sample.
Analysis of ordinal data using imputation and ACM (\ls9ex)
Path Diagram Representation
Descriptive statistics
Parameter Estimates
Goodness of Fit Statistics
The last portion of the output file is a summary of fit statistics and confidence intervals. These statistics are discussed in the Appendix of the New Features in LISREL 9 guide, available in PDF format via the LISREL online Help menu.
Fit Statistics
Three-level Generalized Linear Model (\mglimex)
Selected portions of the output file are displayed below.
Parameter Estimates and Odds Ratios
Estimated variance components
Estimates of the variance components on levels 2 and 3 and the associated p-values indicate that the PreTHKS coefficient does not vary significantly over classes. Note however, that the covariance term is almost significant. The level-3 intercept effect is also not significant. These results seem to indicate a level-2 model random intercept model as being more appropriate.
A level-4 Model with Continuous Outcome Variable (\mlevelex)
Data for the first 10 participants on most of the variables are shown below in the form of a LISREL spreadsheet file, named therapist_L4.lsf.
The variables of interest are:
- site is the level-4 identification variable (49 units in total).
- therapis is the level-3 identification variable (187 units in total).
- particip is the level-2 identification variable (1192 units in total).
- assesmt is a score assigned by a therapist to a particular participant on occasion 0, 1 or 2.
- gender is a gender indicator, with a value of 0 indicating a male participant and 1 a female participant.
- occasion is a predictor variable coded 0, 1 and 2.
- thera1 - thera4 are dummy coded variables indicating four types of therapy.Only selected parts of the output are shown. The output describing the estimated fixed effects after convergence is shown first. From the z-values and associated exceedance probabilities, we see that except for the coefficient associated with gender, the remaining coefficients are all highly significant.
A study of the random part of the model shows that all the intercept effects are highly significant, except for the level-3 (therapists) intercept. From this, we conclude that intercept estimates vary significantly over sites, but not over therapists.
Observed Residuals (\obsresex)
The LSFfile command is to create a new lsf file containing latent variable scores. The Estimate Residuals command is to add estimated residuals to the new lsf file. The name of the new file is POLIDEMstdnew.LSF. The first ten cases of this file are shown below.