[R] Lavaan

David Barron dnbarron at gmail.com
Wed Jun 24 10:39:14 CEST 2015


Does the package semPlot not do what you want?  I notice that you got
an error when you used library(semPlot) because you don't have all the
dependencies installed. When you have semPlot working, it should be
able to produce a graphical output of the results, including
'covariance arrows'.

David

On 23 June 2015 at 17:47, deva d <devazresearch at gmail.com> wrote:
> i am attaching a .csv file, and the associated code worked out in R Studio.
>
> i used the lavaan and sem packages, and conducted it.
>
> now, i wish to draw the SEM model, as is available in AMOS other packages
> and how does one draw the covariance arrows in R.
>
> ONE STATISTICS oriented question - how can one provide interpretation for
> negative coefficients.
>
> kindly guide.
>
> thanks and regds,
>
>
>
>
>
> taxliability <- read.csv("~/R WORK SPACE/taxliability.csv")
>>   View(taxliability)
>> model <-'tax~ inc + exp + svg + inv'
>> fit <- sem(model, data = taxliability)
> Error: could not find function "sem"
>> library("lavaan", lib.loc="~/R/win-library/3.2")
> This is lavaan 0.5-18
> lavaan is BETA software! Please report any bugs.
>> fit <- sem(model, data = taxliability)
> Warning message:
> In lav_model_vcov(lavmodel = lavmodel, lavsamplestats = lavsamplestats,  :
>   lavaan WARNING: could not compute standard errors!
>   lavaan NOTE: this may be a symptom that the model is not identified.
>
>> library("sem", lib.loc="~/R/win-library/3.2")
> Error in loadNamespace(i, c(lib.loc, .libPaths()), versionCheck = vI[[i]])
> :
>   there is no package called ‘htmlwidgets’
> Error: package or namespace load failed for ‘sem’
>> fit <- sem(model, data = taxliability)
> Warning message:
> In lav_model_vcov(lavmodel = lavmodel, lavsamplestats = lavsamplestats,  :
>   lavaan WARNING: could not compute standard errors!
>   lavaan NOTE: this may be a symptom that the model is not identified.
>
>> summary(fit,rsq=T, fit.measures=TRUE)
> lavaan (0.5-18) converged normally after   1 iterations
>
>   Number of observations                            66
>
>   Estimator                                         ML
>   Minimum Function Test Statistic                0.000
>   Degrees of freedom                                 0
>
> Model test baseline model:
>
>   Minimum Function Test Statistic              160.444
>   Degrees of freedom                                 4
>   P-value                                        0.000
>
> User model versus baseline model:
>
>   Comparative Fit Index (CFI)                    1.000
>   Tucker-Lewis Index (TLI)                       1.000
>
> Loglikelihood and Information Criteria:
>
>   Loglikelihood user model (H0)              -3441.453
>   Loglikelihood unrestricted model (H1)      -3441.453
>
>   Number of free parameters                          5
>   Akaike (AIC)                                6892.905
>   Bayesian (BIC)                              6903.854
>   Sample-size adjusted Bayesian (BIC)         6888.113
>
> Root Mean Square Error of Approximation:
>
>   RMSEA                                          0.000
>   90 Percent Confidence Interval          0.000  0.000
>   P-value RMSEA <= 0.05                          1.000
>
> Standardized Root Mean Square Residual:
>
>   SRMR                                           0.000
>
> Parameter estimates:
>
>   Information                                 Expected
>   Standard Errors                             Standard
>
>                    Estimate  Std.err  Z-value  P(>|z|)
> Regressions:
>   tax ~
>     inc               0.103
>     exp              -0.023
>     svg              -0.073
>     inv               0.222
>
> Variances:
>     tax           4662558.169
>
> R-Square:
>
>     tax               0.912
>
>  semPlot
> Error: object 'semPlot' not found
>> library("semPlot", lib.loc="~/R/win-library/3.2")
> Error in loadNamespace(i, c(lib.loc, .libPaths()), versionCheck = vI[[i]])
> :
>   there is no package called ‘htmlwidgets’
> Error: package or namespace load failed for ‘semPlot’
>> semPlot::
> Error: unexpected end of line in "semPlot::"
>>
>
>
>
>
> *....*
>
> *Deva*
>
>
> ...............
>
>
>
> *in search of knowledge, everyday something is added ....*
>
> *in search of wisdom, everyday something is dropped  ... an old Chinese
> Proverb*
> :::::::::::::::::::::::::
>
> On Tue, Jun 23, 2015 at 9:03 PM, Sarah Goslee <sarah.goslee at gmail.com>
> wrote:
>
>> Hi,
>>
>> There are various tutorials for lavaan online, and even an entire book
>> on the R package. Have you worked through those examples and tutorials
>> successfully? If so, a clearer description of what you've tried and
>> what failed is required to be able to help you.
>>
>> Without a reproducible example that includes some sample data (fake is
>> fine), the code you used, and some clear idea of what output you
>> expect, it's impossible to figure out how to help you. Here are some
>> suggestions for creating a good reproducible example:
>>
>> http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example
>>
>> Sarah
>>
>> On Tue, Jun 23, 2015 at 10:48 AM, DzR <devazresearch at gmail.com> wrote:
>> > Dear Senior users of R/R Studio,
>> >
>> > I am very new to this environment hence am unable to plot the SEM models
>> including use of graphic package ggplot.
>> >
>> > Request for some help in getting the plots please.
>> >
>> > Thanks,
>> >
>> > -----
>> > Deva
>>
>> --
>> Sarah Goslee
>> http://www.functionaldiversity.org
>>
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> 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.



More information about the R-help mailing list