# MODEL II
# Install the package 'sem'
# carry data library
library (sem)
# Change R's directory to that in which the data base
.CSV file is
# Load database
teste<- read.csv("BaseSEM.csv",header=T)
teste
# sem package model specification
mod.wh.1 <- specify.model()# Type these values that
specify the model's relations (just use de Ctrl+R over each
relation).
Support->SP,NA,1
Support->SM,sm1,NA
Rejection->RP,rp1,NA
Rejection->RM,NA,1
Overprot->OP,op1,NA
Overprot->OM,NA,1
MotI->MIAO,NA,1
MotI->MIEE,moti1,NA
MotI->MIC,moti2,NA
MotE->MERE,NA,1
MotE->MEIn,mote1,NA
MotE->MEId,mote2,NA
Support->MotI,NA,1
Rejection->MotI,rej1,NA
Overprot->MotI,over1,NA
Support->MotE,sup,NA
Rejection->MotE,rej2,NA
Overprot->MotE,NA,1
MotI->SCI,misci,NA
MotE->SCI,mesci,NA
Support<->Support,ss1,NA
Rejection<->Rejection,rr1,NA
Overprot<->Overprot,oo1,NA
OP<->OP,opcov,NA
OM<->OM,omcov,NA
SP<->SP,spcov,NA
SM<->SM,smcov,NA
RP<->RP,rpcov,NA
RM<->RM,rmcov,NA
MotI<->MotI,micov,NA
MotE<->MotE,mecov,NA
MIAO<->MIAO,miaocov,NA
MIEE<->MIEE,mieecov,NA
MIC<->MIC,miccov,NA
MEIn<->MEIn,meincov,NA
MEId<->MEId,meidcov,NA
MERE<->MERE,merecov,NA
SCI<->SCI,scicov,NA
Support<->Rejection,suprej,NA
Support<->Overprot,supoover,NA
Rejection<->Overprot,rejover,NA
MotI<->MotE,mots1,NA
# Insert de covariance matrix
cov.teste <- cov(teste, y = NULL, use = "everything",
method = c("pearson", "kendall", "spearman"))
cov.teste
# Estimate de model (Here is where we have been finding
difficulties)
sem.wh.1 <- sem(mod.wh.1, cov.teste, N=126)
summary(sem.wh.1)
Prof. Ms. Faculdade Ingá
Tel. 4499298078/4399368588