1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
|
semPlotModel_Onyx <- function(object)
{
# Parse Onyx model:
doc <- xmlParse(object)
# Get Nodes and Edges:
Nodes <- getNodeSet(doc, "/model/graph/node")
Edges <- getNodeSet(doc, "/model/graph/edge")
Const <- as.logical(sapply(Nodes, function(n) xmlGetAttr(n, "constant")))
# Get NodeNames:
NodeNames <- sapply(Nodes, function(n) xmlGetAttr(n, "caption"))
NodeNames[Const] <- ""
# Get edgelist:
Edgelist <- data.frame(
From = as.numeric(as.character(sapply(Edges, function(n) xmlGetAttr(n, "sourceNodeId")))),
To = as.numeric(as.character(sapply(Edges, function(n) xmlGetAttr(n, "targetNodeId")))),
stringsAsFactors=FALSE) + 1
# Define Pars:
Pars <- data.frame(
label = sapply(Edges, function(n) xmlGetAttr(n, "parameterName")),
lhs = NodeNames[Edgelist$From],
edge = ifelse(as.logical(sapply(Edges, function(n) xmlGetAttr(n, "doubleHeaded"))),"<->","->"),
rhs = NodeNames[Edgelist$To],
est = as.numeric(as.character(sapply(Edges, function(n) xmlGetAttr(n, "value")))),
std = NA,
group = "",
fixed = as.logical(sapply(Edges, function(n) xmlGetAttr(n, "fixed"))),
par = 0,
stringsAsFactors=FALSE)
Pars$edge[Pars$lhs==""] <- "int"
Pars$par <- 1:nrow(Pars)
# Vars:
Vars <- data.frame(
name = NodeNames,
manifest = !as.logical(sapply(Nodes, function(n) xmlGetAttr(n, "latent"))),
exogenous = NA,
stringsAsFactors=FALSE)
Vars <- Vars[c(which(Vars$manifest),which(!Vars$manifest)),]
# Return:
semModel <- new("semPlotModel")
semModel@Pars <- Pars
semModel@Vars <- Vars
semModel@Computed <- FALSE
semModel@Original <- list(doc)
semModel@ObsCovs <- list()
semModel@ImpCovs <- list()
# semModel@Thresholds <- Thresh
return(semModel)
}
|
