similar to: (Statistics question) - Nonlinear regression and simultaneous equation

I have a student that I'm encouraging to use R rather than SAS or Stata and within just 2 weeks he has come up with a question that stumps me. What does a person do about endogeneity in generalized linear models? Suppose Y1 and Y2 are 5 category ordinal dependent variables. I see that MASS has polr for estimation of models like that, as long as they are independent. But what if the

Does anyone know if someone is developing full-information maximum likelihood (FIML) estimation algorithms for basic regression functions, like glm()? I think that having FIML options for workhorse functions that typically use ML would give R an edge over other statistical software, given how well FIML performs in missing data situations compared to ML. While my current level of programming

Dear R users/developers I requested help to solve the problem of formulating Multivariate Sample selection model by using Full Information Maximum Likelihood (FIML)estimation method. I could not get any response. I formulated the following code of FIML to analyse univariate sample selection problem. Would you please advise me where is my problem library (sem) library(nrmlepln) Selection

Hi there, I am really new to statistics in R and statistics itself as well. My situation: I ran a lot of OLS regressions with different independent variables. (using the lm() function). After having done that, I know there is endogeneity due to omitted variables. (or perhaps due to any other reasons). And here comes the Hausman test. I know this test is used to identify endogeneity. But what I

Is there a way to use full information maximum likelihood (FIML) to estimate missing data in the sem package? For example, suppose I have a dataset with complete information on X1-X3, but missing data (MAR) on X4. Is there a way to use FIML in this case? I know lavaan and openmx allow you to do it, but I couldn't find anything in the documentation for the sem package. Thanks! -- Dustin Fife

Dear all, I am trying to estimate simultaneous equation model concerning growth in russian regions. I run the analysis by means of FIML in R sem package. I am not familiar with SEM yet, but I've just got several suitable estimated specifications. Nevertheless, sometimes R gives the following warning message: Warning message: Negative parameter variances. Model is probably underidentified.

Hi, can anybody tell me how the Hausman test for endogenty works? I have a simulated model with three correlated predictors (X1-X3). I also have an instrument W for X1 Now I want to test for endogeneity of X1 (i.e., when I omit X2 and X3 from the equation). My current approach: library(systemfit) fit2sls <- systemfit(Y~X1,data=data,method="2SLS",inst=~W) fitOLS <-

Hi I was asked if lme can use FIML (Full Information Maximum Likelihood) instead of REML or ML but I don't know the answer. Does anybody know if this is implemented in R? Thanks Francisco

Hi, I am wondering if anyone knows of an easy way to fit a saturated model using the sem package on raw data? Say the data were: mtcars[, c("mpg", "hp", "wt")] The model would estimate the three means (intercepts) of c("mpg", "hp", "wt"). The variances of c("mpg", "hp", "wt"). The covariance of mpg with

Hello! I am trying to reproduce (for a publication) analyses that I ran several months ago using lavaan, I'm not sure which version, probably 0.4-12. A sample model is given below: pathmod='mh30days.log.w2 ~ mh30days.log + joingroup + leavegroup + alwaysgroup + grp.partic.w2 + black + age + bivoc + moved.conf + local.noretired + retired + ds + ministrytime + hrswork + nomoralescore.c +

Hello, I have a set of 100 variables with 1560 observations. I did an O.L.S regression of three of these variables on a fourth. But there are problems of endogeneity... So I look in my dataset for instruments to do an IV. I can't find a good instrument because their correlation with my endogeneous variables are too low. But I see that when I create a combined variable composed of 12 variables

Hi there, as a student sociology, I'm starting to learn about SEM. The course I follow is based on LISREL, but I want to use the SEM-package on R parallel to it. Using LISREL, I found it to be very usable to be able to see the total direct and total indirect effects (standardized and unstandardized) in the output. Can I create these effects using R? I know how to calculate them

Hello- I am using R to build my initial models, but need to use AMOS to compare the models of two groups (adults vs. kids). The problem is I am getting different results with R and AMOS for the initial models of the separate groups (and the R results make more sense). The parameter estimates (path coefficients and variances) from both programs are nearly identical, but the model chi-squares

Hi - I am confronting a situation where I have a set of structural equation and one or two of my responses are multinomial. I understand that sem would not deal with the unordered response. So I am thinking of the following two ways: 1. Expanding my response to a new set of binary variables corresponding to each label of my multinomial response. Then use each of these as a separate response in my

Hi, I am using tsls function from sem package to estimate a model which includes large number of data. Among its predictors, it has a nominal data which has about 10 possible values. So I expand this parameter into 9-binary-value predictors with the coefficient of base value equals 0. I also have another continuous predictor. The problem is that, whenever I run the tsls, I will get 'System

Hi Frances, I have not touched the system.fit package for quite some time, but to solve your problem the following two pointers might be helpful: 1) Recast your model in the revised form, i.e., include your identity directly into your reaction functions, if possible. 2) For solving your model, you can employ the Gau?-Seidel method (see https://en.wikipedia.org/wiki/Gauss%E2%80%93Seidel_method).

Can R (which as I am quickly discovering is a awsome package) easily perform Two-Stage, Three-Stage, and Seemingly Unrelated Regression? Jeff D. Hamann -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the

Who was speaking about non-linear models in the first place??? The Klein-Model(s) and pretty much all simultaneous equation models encountered in macro-econometrics are linear and/or can contain linear approximations to non-linear relationships, e.g., production functions of the Cobb-Douglas type. Best, Bernhard -----Urspr?ngliche Nachricht----- Von: Berend Hasselman [mailto:bhh at xs4all.nl]

Frances, I would not advise Gauss-Seidel for non linear models. Can be quite tricky, slow and diverge. You can write your model as a non linear system of equations and use one of the nonlinear solvers. See the section "Root Finding" in the task view NumericalMathematics suggesting three packages (BB, nleqslv and ktsolve). These package are certainly able to handle medium sized models.

Dear R-users, A new package called `lavaan' (for latent variable analysis) has been uploaded to CRAN. The current version of lavaan (0.3-1) can be used for path analysis, confirmatory factor analysis, structural equation modeling, and growth curve modeling. More information can be found on the website: http://lavaan.org Some notable features of lavaan: - the 'lavaan model