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
|
\name{Random.Start}
\alias{Random.Start}
\title{Generate a Random Orthogonal Rotation}
\usage{
Random.Start(k)
}
\arguments{
\item{k}{An integer indicating the dimension of the square matrix.}
}
\description{
Random orthogonal rotation to use as Tmat matrix to start GPForth or GPFoblq.
}
\value{An orthogonal matrix.}
\details{
The random start function produces an orthogonal matrix with columns
of length one based on the QR decompostion.
}
\seealso{
\code{\link{GPForth}},
\code{\link{GPFoblq}},
\code{\link{oblimin}}
}
\examples{
Global.min <- function(A,method,B=10){
fv <- rep(0,B)
seeds <- sample(1e+7, B)
for(i in 1:B){
cat(i," ")
set.seed(seeds[i])
gpout <- GPFoblq(A=A, Random.Start(ncol(A)), method=method)
dtab <- dim(gpout$Table)
fv[i] <- gpout$Table[dtab[1],2]
cat(fv[i], "\n")
}
cat("Min is ",min(fv),"\n")
set.seed(seeds[order(fv)[1]])
ans <- GPFoblq(A=A, Random.Start(ncol(A)), method=method)
ans
}
data("Thurstone", package="GPArotation")
Global.min(box26,"simplimax",10)
}
\author{Coen A. Bernaards and Robert I. Jennrich
with some R modifications by Paul Gilbert
}
\concept{rotation}
\keyword{multivariate}
|
