A strange problem with sem:
I downloaded the sem library and then, I specified my simple measurement model
(below). I highlighted it and ran it. It ran, but it did NOT tell me "22
lines read". And nothing works after that - it looks like it runs, but it
does not produce anything...
Did I make a mistake somewhere in the model? (notice, TIME has only 1 indicator
- t1, and I fixed t1's error variance at 0.414.)
Thank you!
model1 <- specify.model()
NECESSITY -> n1, NA, 1
NECESSITY -> n2, lam_nec_2, NA
NECESSITY -> n3, lam_nec_3, NA
NECESSITY -> n4, lam_nec_4, NA
FRUGALITY -> f1, NA, 1
FRUGALITY -> f2, lam_frug_2, NA
FRUGALITY -> f3, lam_frug_3, NA
TIME -> t1, NA, 1
n1 <-> n1, theta_n1, NA
n2 <-> n2, theta_n2, NA
n3 <-> n3, theta_n3, NA
n4 <-> n4, theta_n4, NA
f1 <-> f1, theta_f1, NA
f2 <-> f2, theta_f2, NA
f3 <-> f3, theta_f3, NA
t1 <-> t1, NA, 0.414
NECESSITY <-> NECESSITY, phi_NN, NA
FRUGALITY <-> FRUGALITY, phi_FF, NA
TIME <-> TIME, phi_TT, NA
NECESSITY <-> TIME, phi_NT, NA
NECESSITY <-> FRUGALITY, phi_NF, NA
FRUGALITY <-> TIME, phi_FT, NA
____________________________________________________________________________________
Finding fabulous fares is fun.
Dear Dimitri, This works find for me. It's not possible to know what you did wrong without more information, but one guess (probably wrong) is that you didn't enter a blank line to exit from specify.model(), which simply uses scan() to read the model:> library(sem) > model1 <- specify.model()1: NECESSITY -> n1, NA, 1 2: NECESSITY -> n2, lam_nec_2, NA 3: NECESSITY -> n3, lam_nec_3, NA 4: NECESSITY -> n4, lam_nec_4, NA 5: FRUGALITY -> f1, NA, 1 6: FRUGALITY -> f2, lam_frug_2, NA 7: FRUGALITY -> f3, lam_frug_3, NA 8: TIME -> t1, NA, 1 9: n1 <-> n1, theta_n1, NA 10: n2 <-> n2, theta_n2, NA 11: n3 <-> n3, theta_n3, NA 12: n4 <-> n4, theta_n4, NA 13: f1 <-> f1, theta_f1, NA 14: f2 <-> f2, theta_f2, NA 15: f3 <-> f3, theta_f3, NA 16: t1 <-> t1, NA, 0.414 17: NECESSITY <-> NECESSITY, phi_NN, NA 18: FRUGALITY <-> FRUGALITY, phi_FF, NA 19: TIME <-> TIME, phi_TT, NA 20: NECESSITY <-> TIME, phi_NT, NA 21: NECESSITY <-> FRUGALITY, phi_NF, NA 22: FRUGALITY <-> TIME, phi_FT, NA 23: Read 22 records>Regards, John On Wed, 11 Apr 2007 11:41:25 -0700 (PDT) John Smith <ld7631 at yahoo.com> wrote:> A strange problem with sem: > > I downloaded the sem library and then, I specified my simple > measurement model (below). I highlighted it and ran it. It ran, but > it did NOT tell me "22 lines read". And nothing works after that - it > looks like it runs, but it does not produce anything... > Did I make a mistake somewhere in the model? (notice, TIME has only 1 > indicator - t1, and I fixed t1's error variance at 0.414.) > Thank you! > > model1 <- specify.model() > NECESSITY -> n1, NA, 1 > NECESSITY -> n2, lam_nec_2, NA > NECESSITY -> n3, lam_nec_3, NA > NECESSITY -> n4, lam_nec_4, NA > FRUGALITY -> f1, NA, 1 > FRUGALITY -> f2, lam_frug_2, NA > FRUGALITY -> f3, lam_frug_3, NA > TIME -> t1, NA, 1 > n1 <-> n1, theta_n1, NA > n2 <-> n2, theta_n2, NA > n3 <-> n3, theta_n3, NA > n4 <-> n4, theta_n4, NA > f1 <-> f1, theta_f1, NA > f2 <-> f2, theta_f2, NA > f3 <-> f3, theta_f3, NA > t1 <-> t1, NA, 0.414 > NECESSITY <-> NECESSITY, phi_NN, NA > FRUGALITY <-> FRUGALITY, phi_FF, NA > TIME <-> TIME, phi_TT, NA > NECESSITY <-> TIME, phi_NT, NA > NECESSITY <-> FRUGALITY, phi_NF, NA > FRUGALITY <-> TIME, phi_FT, NA > > > >____________________________________________________________________________________> Finding fabulous fares is fun. > > ______________________________________________ > R-help at stat.math.ethz.ch 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.-------------------------------- John Fox Department of Sociology McMaster University Hamilton, Ontario, Canada http://socserv.mcmaster.ca/jfox/
Thanks a lot - I did not realize one had to hit "return" after running it (or leave an extra empty line below it). It works! Dimitri ----- Original Message ---- From: John Fox <jfox at mcmaster.ca> To: John Smith <ld7631 at yahoo.com> Cc: r-help at stat.math.ethz.ch Sent: Wednesday, April 11, 2007 5:20:11 PM Subject: Re: [R] sem is not "taking" the model Dear Dimitri, This works find for me. It's not possible to know what you did wrong without more information, but one guess (probably wrong) is that you didn't enter a blank line to exit from specify.model(), which simply uses scan() to read the model:> library(sem) > model1 <- specify.model()1: NECESSITY -> n1, NA, 1 2: NECESSITY -> n2, lam_nec_2, NA 3: NECESSITY -> n3, lam_nec_3, NA 4: NECESSITY -> n4, lam_nec_4, NA 5: FRUGALITY -> f1, NA, 1 6: FRUGALITY -> f2, lam_frug_2, NA 7: FRUGALITY -> f3, lam_frug_3, NA 8: TIME -> t1, NA, 1 9: n1 <-> n1, theta_n1, NA 10: n2 <-> n2, theta_n2, NA 11: n3 <-> n3, theta_n3, NA 12: n4 <-> n4, theta_n4, NA 13: f1 <-> f1, theta_f1, NA 14: f2 <-> f2, theta_f2, NA 15: f3 <-> f3, theta_f3, NA 16: t1 <-> t1, NA, 0.414 17: NECESSITY <-> NECESSITY, phi_NN, NA 18: FRUGALITY <-> FRUGALITY, phi_FF, NA 19: TIME <-> TIME, phi_TT, NA 20: NECESSITY <-> TIME, phi_NT, NA 21: NECESSITY <-> FRUGALITY, phi_NF, NA 22: FRUGALITY <-> TIME, phi_FT, NA 23: Read 22 records>Regards, John On Wed, 11 Apr 2007 11:41:25 -0700 (PDT) John Smith <ld7631 at yahoo.com> wrote:> A strange problem with sem: > > I downloaded the sem library and then, I specified my simple > measurement model (below). I highlighted it and ran it. It ran, but > it did NOT tell me "22 lines read". And nothing works after that - it > looks like it runs, but it does not produce anything... > Did I make a mistake somewhere in the model? (notice, TIME has only 1 > indicator - t1, and I fixed t1's error variance at 0.414.) > Thank you! > > model1 <- specify.model() > NECESSITY -> n1, NA, 1 > NECESSITY -> n2, lam_nec_2, NA > NECESSITY -> n3, lam_nec_3, NA > NECESSITY -> n4, lam_nec_4, NA > FRUGALITY -> f1, NA, 1 > FRUGALITY -> f2, lam_frug_2, NA > FRUGALITY -> f3, lam_frug_3, NA > TIME -> t1, NA, 1 > n1 <-> n1, theta_n1, NA > n2 <-> n2, theta_n2, NA > n3 <-> n3, theta_n3, NA > n4 <-> n4, theta_n4, NA > f1 <-> f1, theta_f1, NA > f2 <-> f2, theta_f2, NA > f3 <-> f3, theta_f3, NA > t1 <-> t1, NA, 0.414 > NECESSITY <-> NECESSITY, phi_NN, NA > FRUGALITY <-> FRUGALITY, phi_FF, NA > TIME <-> TIME, phi_TT, NA > NECESSITY <-> TIME, phi_NT, NA > NECESSITY <-> FRUGALITY, phi_NF, NA > FRUGALITY <-> TIME, phi_FT, NA > > > >____________________________________________________________________________________> Finding fabulous fares is fun. > > ______________________________________________ > R-help at stat.math.ethz.ch 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.-------------------------------- John Fox Department of Sociology McMaster University Hamilton, Ontario, Canada http://socserv.mcmaster.ca/jfox/ ____________________________________________________________________________________ The fish are biting.