1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
# last modified 2014-10-27 by J. Fox
read.moments <- function(...){
.Deprecated("readMoments", package="sem")
readMoments(...)
}
readMoments <- function(file="", text, diag=TRUE,
names=as.character(paste("X", 1:n, sep=""))){
elements <- if (!missing(text)) scan(text=text) else scan(file=file)
m <- length(elements)
d <- if (diag) 1 else -1
n <- floor((sqrt(1 + 8*m) - d)/2)
if (m != n*(n + d)/2)
stop("wrong number of elements (cannot make square matrix)")
if (length(names) != n) stop("wrong number of variable names")
X <- diag(n)
X[upper.tri(X, diag=diag)] <- elements
rownames(X) <- colnames(X) <- names
t(X)
}
|
