search for: mplus

Has anyone successfully linked Mplus and R and have code on how to do it? It would be great to use R to simulate datasets and then ship to Mplus for SEM analysis. I know one researcher who does this quite regularly but I can’t get ahold of him right now. Sincerely, Jeff Miller **************************************** Jeffrey M....

Dear list, I want to run a statistical program (using its .exe file) from R by writing a script. I know there are some packages that call WinBUGS, Mplus etc. form R. I just want to call the .exe extension of this program and run several times writing a code in R. Thus, I want to have the output inside R. I just don't know where to start. Does anyone have any idea about that? Is there a universal package to call application files of other stat...

...tely returned the following error message: Error in chol.default(S) : the leading minor of order 5 is not positive definite Error in lavSampleStatsFromData(Data = lavaanData, rescale = (lavaanOptions$estimator == : lavaan ERROR: sample covariance can not be inverted I ran the same model in Mplus and found that it has low covariance coverage for the 7 repeated measures, but all coverage is greater than 0. The model ran fine in Mplus once I set the covariance coverage limit to .001. That provided me starting values to try in lavaan, but again it immediately returned the same error message....

Is it possible to fit a structural equation model with link functions in R? I am trying to build a logistic-regression-like model in sem, because incorporating the dichotomous variables linearly seems inappropriate. Mplus can do something similar by specifying a 'link' parameter, but I would like to be able to do it in R, ofcourse. I have explored the 'sem' package from John Fox, but it does not seem to be able to fit non-linear relations. Is there some R-package or way to get this done? I have als...

Do I miss here something? dtaa = read.table("http://www.ats.ucla.edu/stat/mplus/examples/ma_snijders/mlbook1.dat", sep=",") head(dtaa) # shows the data as it should be save(dtaa,"dtaa",file="c:/dtaa") d = load("c:/dtaa") head(d) # all data is lost, it only shows [1] "dtaa" "dtaa" Thanks for your hi...

...;, "HHincome", "Race", "Age", "Gender", "Stress1"))) model = ' Internet ~ HHincome + Race + Age + Gender + Stress1 Stress3 ~ Internet + HHincome + Race + Age + Gender + Stress1 ' fit=sem(model, sample.nobs=161, sample.cov=myCov, mimic.Mplus=FALSE) #check the number of parameters being estimated inspect(fit, what="free") #Note the DF for the Chi-square is 0, when it should be 28-13 = 15 summary(fit) Andrew Miles [[alternative HTML version deleted]]

...n pursuing this approach? 2. Does anyone know of some substantial published applications of this approach using SEM? 3. Would any of the available R packages for SEM (e.g., lavaan, sem, OpenMx) be particulary straightforward to use in doing the bootstrapping? Thus far, the SEM has been done using MPLUS. I've not tried SEM in R yet, but would be interested in giving it a shot.?The SEM itself is relatively straightforward. Four latent variables, one with 7 indicators and the others with 4 indicators each. A couple of indirect paths involving mediation. Some pretty non-normal data though.??Lots...

...is using the cluster package in R. I have about 30 cases and 6 indicators, some of which are binary indicators and others are ratio-level variables (percentages). I looked around for information in flexmix, lca, and poLCA, and couldn't find anything relevant or probably I missed them. I checked MPlus, and it seems to me MPlus can handle latent class analysis with indicators of different measurement levels (continous or nominal). And my intuition is that R ought to be able to it. Any ideas about which package I should use and any examples that I can model after. As a long-time user of Stata, I...

R list- i have begun using the "MplusAutomation" while piloting a large-scale simulation (~200,000 replications). since the package takes advantage of the DOS batch mode available in Mplus, each replication starts and activates a new instance of a command prompt window. this effectively locks me out of my computer for the duration...

Is there an equivalent to MPlus's Full Information Maximum Likelihood (FIML) missing data estimator for R? If so, is there a way to take covariance structures produced by such a package and perform multiple regression with these? If you are unfamiliar with Mplus' FIML below is a link to their manual. Their estimation tec...

...y data only includes 30 replicate weights and a sampling weight. The svrepdesign requires additional information such as jacknife weights in an external object (for scale, rscale), however, I do not have these. I typically just use the sampling weight when doing analyses in other platforms such as MPlus, however there does not appear to a way of referring to this directly in the ltm formulae. My question is: am I able to insert sampling weights directly into an ltm analysis, or is there a way of creating a weighted data object with the information available which I could then insert into the ltm...

Hello. Using the polychor function > polychor(data[c(s1,s2)] ) for polychoric correlations of two ordinal variables in R takes a long time for N=7000 (20 minutes+) and significantly slows down my computer. Now, I have a pretty old computer, but it takes about 20 seconds for MPLUS to print out the complete polychoric correlation matrix for all 16 variables, while I am running the R function only for 2 of the 16 variables. Why is that? Can that process be speeded up? What makes a polychoric correlation so much more computationally intensive than a pearson correlation? Thank...

I was in a presentation of optimizations fitted with both MPlus and SAS yesterday. In a batch of 1000 bootstrap samples, between 300 and 400 of the estimations did not converge. The authors spoke as if this were the ordinary cost of doing business, and pointed to some publications in which the nonconvergence rate was as high or higher. I just don't belie...

...Added complexity is that I have 3 levels, so I have R, G and say H (V=WHW'+ZGZ'+R). The lme is giving me the correct results, I am having trouble finding the right corresponding specification for the gls. Thanks for your help. Toby dtaa = read.table("http://www.ats.ucla.edu/stat/mplus/examples/ma_snijders/mlbook1.dat", sep=",") dta1 = reshape(dtaa, list(c("V10","V12")), "score", direction="long", drop=c("V2","V3","V4","V5","V6","V7","V8","V9",&...

...12 dominic wrote > > This is basically my code: > > library(MASS) > lmsreg(formula = b0 ~ b1 + b3 + b1xb2, data=mydata) > > b1xb2 is an interaction but it was the centered value for a continuous > variable times a categorical variable. > > I am used to using Mplus statistical software and SAS for robust > regression, and they'll usually give you a t-test, p-value as well as the > SE and/or CI. But I am not getting it here even though I have sent that is > the correct code to get it on some websites that have posted their output. > I > &gt...

Hello, I've been trying to answer a problem I have had for some months now and came across multivariate multilevel modeling. I know MPLUS and SPSS quite well but these programs could not solve this specific difficulty. My problem: 9 correlated dependent variables (medical symptoms; categorical, 0-3), 5 measurement points, 10 time-varying covariates (life events; dichotomous, 0-1), N ~ 900. Up to 35% missing values on some variables,...

Hi, I'm generating the name of the variable with paste function and then using that variable name further to get the specific position value from the data.frame, here is the snippet from my code: modelResults <- extractModelParameters("C:/PilotStudy/Mplus_Input/Test", recursive=TRUE) #extractModelParameters reads all the output files from the Test folder and create the following variables in R for each file read: #C..PilotStudy.Mplus_Input.Test.rep1.out.unstandardized.est #C..PilotStudy.Mplus_Input.Test.rep2.out.unstandardized.est #C..PilotS...

I am looking to fit one or more latent categorical variables to data that is a mixture of categorical and continuous variables. Factor analysis would work for continuous data, latent class analysis for categorical data. I understand that in a package such as MPlus I could perform a single analysis of both data types. Are there similar routines available in R? Stuart -----Original Message----- From: Prof Brian Ripley <ripley at stats.ox.ac.uk> To: Dr Stuart Leask <stuart.leask at nottingham.ac.uk> Cc: r-help at stat.math.ethz.ch <r-help at st...

I have been working through the book "Applied longitudinal data analysis: modeling change and event occurrence" by Judith D. Singer and John B. Willett. I have been working examples using SAS and also using it as an opportunity for learning to use R for statistical analysis. I ran into some difficulties in chapter 8 which deals with using structural equation modeling. I have tried to

Dears, I'm a college student and In doing my statistics homework. I use R with SEM package as my tool for sem analysis, but my teacher told me AMOS is more suitable for such analysis. Could someone help tell me whether it is true that some commercial software is better accepted in academic fields? Sorry if I should not post such topics here. -- Best Regards, Reeyarn T. Lee Accounting