Klinedata <- read.moments(diag=FALSE, names=c('LowMorale', 'IllnessSymptoms', 'NeurologicalDysfunction', 'PoorRelationships', 'DiminishedSES'))
     .53
     .15  .18
     .52  .29  -.05
     .30  .34  .23   .09

pathmodel1= specify.model()
	NeurologicalDysfunction -> LowMorale, p1, NA
	DiminishedSES -> LowMorale, p2, NA
	LowMorale -> IllnessSymptoms, p3, NA
	LowMorale -> PoorRelationships, p4, NA
	LowMorale <-> LowMorale, d1, NA
	IllnessSymptoms <-> IllnessSymptoms, d2, NA
	PoorRelationships <-> PoorRelationships, d3, NA

psychosomaticmodel = sem(pathmodel1, Klinedata, 469, fixed.x=c('DiminishedSES', 'NeurologicalDysfunction'))

summary(psychosomaticmodel)


#Conventional Medical Model

pathmodel2= specify.model()
	IllnessSymptoms -> DiminishedSES , p1a, NA
	IllnessSymptoms -> LowMorale, p1b, NA
	NeurologicalDysfunction -> DiminishedSES, p2a, NA
	NeurologicalDysfunction -> PoorRelationships, p2b, NA
	DiminishedSES -> LowMorale, p3, NA
	LowMorale -> PoorRelationships, p4, NA
	DiminishedSES  <-> DiminishedSES , d1, NA
	LowMorale <-> LowMorale, d2, NA
	PoorRelationships <-> PoorRelationships, d3, NA

conventionalmedicalmodel=sem(pathmodel2, Klinedata, 469, fixed.x=c('IllnessSymptoms', 'NeurologicalDysfunction'))

summary(conventionalmedicalmodel)