Displaying 20 results from an estimated 600 matches similar to: "MplusAutomation"
2012 Nov 01
2
SEM validation: Cross-Validation vs. Bootstrapping
Hello All,
Recently, I was asked to help out with an SEM cross-validation analysis. Initially, the project was based on "sample-splitting" where half of cases were randomly assigned to a training sample and half to a testing sample. Attempts to replicate a model developed in the training sample using the testing sample were not entirely successful. A number of parameter estimates were
2010 Dec 11
2
remove quotes from the paste output
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
2010 Mar 14
3
CRAN (and crantastic) updates this week
CRAN (and crantastic) updates this week
New packages
------------
* apcluster (1.0.1)
Ulrich Bodenhofer
http://crantastic.org/packages/apcluster
The apcluster package implements Frey's and Dueck's Affinity
Propagation clustering in R. The algorithms are analogous to the
Matlab code published by Frey and Dueck.
* BioPhysConnectoR (1.6-1)
Franziska Hoffgaard
2012 Aug 10
1
Lavaan: Immediate non-positive definite matrix
Hi,
I recently tried to estimate a linear unconditional latent growth curve on
7 repeated measures using lavaan (most recent version):
modspec='
alpha =~ 1*read_g0 + 1*read_g1 + 1*read_g2 + 1*read_g3 + 1*read_g4 +
1*read_g5 + 1*read_g6
beta =~ 0*read_g0 + 1*read_g1 + 2*read_g2 + 3*read_g3 + 4*read_g4 +
5*read_g5 + 6*read_g6
'
gmod=lavaan(modspec, data=math, meanstructure=T,
2013 Mar 16
3
Running other programs from R
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
2007 May 23
2
saving datafreame object problem
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 hint on this.
2008 Jul 27
2
Link functions in SEM
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,
2011 Mar 17
2
Incorrect degrees of freedom in SEM model using lavaan
I have been trying to use lavaan (version 0.4-7) for a simple path model,
but the program seems to be computing far less degrees of freedom for my
model then it should have. I have 7 variables, which should give (7)(8)/2 =
28 covariances, and hence 28 DF. The model seems to only think I have 13
DF. The code to reproduce the problem is below. Have I done something
wrong, or is this something I
2011 Jul 27
1
Inserting weights in ltm package
Afternoon R help,
I want to run Rasch/IRT analyses using the ltm package, however, I am
using large scale survey data which requires weighting for accurate
results. I attempted to create a weighted object to insert into the
formulae of the ltm packages, however, the survey data only includes
30 replicate weights and a sampling weight. The svrepdesign requires
additional information such as
2012 Dec 07
1
Polychor() - why does it take that long?
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
2011 Dec 15
2
fundamental guide to use of numerical optimizers?
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 believe that's right, and if some
2007 May 21
1
can I get same results using lme and gls?
Hi All
I was wondering how to get the same results with gls and lme. In my lme, the
design matrix for the random effects is (should be) a identity matrix and
therefore G should add up with R to produce the R matrix that gls would report
(V=ZGZ'+R). 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
2012 Jan 23
1
R not giving significance tests for coefficients/estimates?
> 3x4
Error: unexpected symbol in "3x4"
R has no idea that you equate "x" as multiplication.. use an astrix
> 3*4
[1] 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
2012 Apr 06
2
Multivariate Multilevel Model: is R the right software for this problem
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.
2002 Feb 21
2
Re: Factor analysis of categorical or mixed categorical/continuousdata in
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-----
2007 Nov 01
0
Mplus and R
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. Miller, PhD
HYPERLINK
2008 Nov 05
0
latent class analysis of nominal and continuous indicators
Dear R-help listers,
I am a new convert to R. I am trying to use a r package to conduct latent
class analysis as a triangulation check of my cluster analysis 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
2011 Mar 05
0
extractModelParameters HELP!!!
Hello Gurus,
I'm using one of the function, i.e. extractModelParameters, from the
MplusAutomation package to read parameters fromMplus out files into R. I have
360 cells and each cell has several out files (max. 100). I'm extracting the
model parameters from each cell and doing subsequent calculations, primarily, on
'est.' column. The following code demonstrates how I'm
2008 Dec 11
0
Equivalent to Full Information Maximum Likelihood (FIML) in R?
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 technology is discussed on page 25. I have asked the
developer of the
2024 May 22
1
confint Attempts to Use All Server CPUs by Default
Following up on this -- on my system, I have 69 packages installed
that appear to provide something like a confint() method:
h <- help.search("confint", agrep = FALSE)
p <- sort(unique(h$matches$Package))
length(p)
## [1] 69
p
[1] "bamlss" "bbmle" "binom" "brglm2"
[5] "broom"