Dear Sarah,
It's generally a good idea to include a reproducible example if you want to get
help with a problem, but in this case it's a safe bet that the problem is that
the model you specified has no variance or covariance parameters for the
variables x1 and x2, which, I assume, you mean to be exogenous. The easiest way
to include these variances and covariance in the model is to specify the
argument fixed.x=c("x1", "x2") in the call to sem().
In addition:
(1) Your model is fully recursive (guessing that all the x's and y's are
observed variables), and so it amounts to four OLS regressions. You could just
use lm() to fit the model.
(2) It's generally easier in the sem package to use specifyEquations() than
specifyModel() for model specification.
(3) If you have the original data set, as you do, it's generally preferable to
use the data argument to sem() than to pass it the covariance matrix for the
observed variables.
I hope that this helps,
John
------------------------------------------------
John Fox
McMaster University
Hamilton, Ontario, Canada
http://socserv.mcmaster.ca/jfox/
On Sat, 2 Nov 2013 11:02:31 +0100
Sarah Rogers <rogerssara...@gmail.com> wrote:
> Hello,
>
> I have just started to work on a path analysis (see attached image for the
> diagram), but I have encountered an error message.
>
>
>
>
> This is the code I have used:
>
> cov_matrix<-var(xdata)
>
> library(sem)
> model.xdata<-specifyModel()
> x1 -> y2, xy12, NA
> x2 -> y1, xy21, NA
> y1 -> y2, yy12, NA
> y2 -> y3, yy23, NA
> y2 -> y4, yy24, NA
> y3 -> y4, yy34, NA
> y2 <-> y2, y2error, NA
> y1 <-> y1, y1error, NA
> y3 <-> y3, y3error, NA
> y4 <-> y4, y4error, NA
>
> model.xdata.sem <- sem(model.xdata, cov_matrix, nrow(xdata))
>
> and the error message is:
> Error in csem(model = model.description, start, opt.flag = 1, typsize =
> typsize, :
> The matrix is non-invertable.
>
> I fear to have a problem in the data.
> I would be very grateful if you could help me to solve this problem and
> proceed with my analyses.
>
> thank you in advance for your help!
> Sarah
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
