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
|
# semPaths.loadings <- function(object,...)
# {
# invisible(semPaths(semPlotModel(object),...))
# }
#
### SINGLE GROUP MODEL ###
semPlotModel.loadings <- function(object, ...)
{
# Check if object is of class "sem":
if (!"loadings"%in%class(object)) stop("Input must be a 'factanal' object")
manNames <- rownames(object)
latNames <- colnames(object)
# Define Pars:
Pars <- data.frame(
label = "",
lhs = rep(latNames,each=length(manNames)),
edge = "--",
rhs = rep(manNames,times=length(latNames)),
est = c(object),
std = c(object),
group = "",
fixed = FALSE,
par = 1:length(object),
stringsAsFactors=FALSE)
# Variable dataframe:
Vars <- data.frame(
name = c(manNames[order(apply(abs(object),1,which.max))],latNames),
manifest = c(rep(TRUE,nrow(object)),rep(FALSE,ncol(object))),
exogenous = NA,
stringsAsFactors=FALSE)
semModel <- new("semPlotModel")
semModel@Pars <- Pars
semModel@Vars <- Vars
semModel@Computed <- FALSE
semModel@Original <- list(object)
semModel@ObsCovs <- list()
semModel@ImpCovs <- list()
return(semModel)
}
|
