#------Raykov's Covariance Matrix--------
Variable Y1 Y2 Y3 Y4
Y1 1.8870
Y2 1.5744 4.3173
Y3 2.1767 3.1596 7.5345
Y4 0.1015 0.2136 0.2690 2.9428

#-------------------------------------




#----------------My Code-----------------

library(OpenMx)
manifests <- c(names(data))
latents <- c("F1","Theta1", "Theta2", "Theta3", "Theta4")

myRaykovModel <- mxModel("Raykov Model", type = "RAM", 
mxData(observed = dataCov, type = "cov", numObs = 400),

manifestVars = manifests,
latentVars = latents,
mxPath(
from = manifests,
arrows = 2,
free = c(TRUE,TRUE,TRUE,TRUE),
labels = c("e1","e2","e3","e4")),

mxPath(from = "F1",
arrows = 2,
free = FALSE, 
values = 1.0,
labels = "varF1"
),

mxPath(
from = "F1",
to = manifests,
arrows = 1,
free= c(TRUE,TRUE,TRUE,TRUE),
values = 1,
labels = c("L1","L2","L3","L4")),


mxAlgebra((e1/(L1+L2+L3+L4)),
name = "Eta1Constraint"),

mxPath(from = "Theta1",
free = FALSE,
labels = "Eta1Constraint[1,1]"
)


)


myRaykovFit <- mxRun(myRaykovModel)