library(sem)

holzinger = read.table('holzall.dat',header=T)
names(holzinger)

S =  var(holzinger[,c(7:15,20,21,23,25)])

model.holz = specify.model()
    spatial   -> visual,   lam011, NA
    spatial   -> cubes,    lam021, NA
    spatial   -> paper,    lam031, NA
    spatial   -> lozenges, lam041, NA
    spatial   -> numberr,  lam061, NA
    spatial   -> object,   lam071, NA
    spatial   -> figurew,  lam081, NA
    spatial   -> general,  lam091, NA
    spatial   -> paragrap, lam101, NA
    spatial   -> wordc,    lam121, NA
    spatial   -> wordm,    lam131, NA
    memory    -> visual,   lam012, NA
    memory    -> cubes,    lam022, NA
    memory    -> paper,    lam032, NA
    memory    -> wordr,    lam052, NA
    memory    -> numberr,  lam062, NA
    memory    -> object,   lam072, NA
    memory    -> figurew,  lam082, NA
    memory    -> general,  lam092, NA
    memory    -> paragrap, lam102, NA
    memory    -> wordc,    lam122, NA
    memory    -> wordm,    lam132, NA
    verbal    -> visual,   lam013, NA
    verbal    -> cubes,    lam023, NA
    verbal    -> paper,    lam033, NA
    verbal    -> numberr,  lam063, NA
    verbal    -> object,   lam073, NA
    verbal    -> figurew,  lam083, NA
    verbal    -> general,  lam093, NA
    verbal    -> paragrap, lam103, NA
    verbal    -> sentence, lam113, NA
    verbal    -> wordc,    lam123, NA
    verbal    -> wordm,    lam133, NA
    visual   <-> visual,   th01,   NA
    cubes    <-> cubes,    th02,   NA
    paper    <-> paper,    th03,   NA
    lozenges <-> lozenges, th04,   NA
    wordr    <-> wordr,    th05,   NA
    numberr  <-> numberr,  th06,   NA
    object   <-> object,   th07,   NA
    figurew  <-> figurew,  th08,   NA
    general  <-> general,  th09,   NA
    paragrap <-> paragrap, th10,   NA
    sentence <-> sentence, th11,   NA
    wordc    <-> wordc,    th12,   NA
    wordm    <-> wordm,    th13,   NA
    spatial  <-> spatial,    NA,    1
    memory   <-> memory,     NA,    1
    verbal   <-> verbal,     NA,    1
    spatial  <-> memory,   psi12,  NA
    spatial  <-> verbal,   psi13,  NA
    memory   <-> verbal,   psi23,  NA

obs.vars.holz = row.names(S)

sem.holz = sem(model.holz, S, 145, obs.vars.holz,maxiter=1000)

summary(sem.holz,digits=3)
std.coef(sem.holz,digits=3)
mod.indices(sem.holz)