projects
/
morpheus.git
/ commitdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
| commitdiff |
tree
raw
|
patch
|
inline
| side by side (parent:
0e0fa9f
)
Fix mistake in multistart.R
author
Benjamin Auder
<benjamin.auder@somewhere>
Fri, 20 Dec 2019 22:30:22 +0000
(23:30 +0100)
committer
Benjamin Auder
<benjamin.auder@somewhere>
Fri, 20 Dec 2019 22:30:22 +0000
(23:30 +0100)
reports/multistart.R
patch
|
blob
|
blame
|
history
diff --git
a/reports/multistart.R
b/reports/multistart.R
index
28d471f
..
e8a9300
100644
(file)
--- a/
reports/multistart.R
+++ b/
reports/multistart.R
@@
-9,8
+9,9
@@
testMultistart <- function(N, n, p, beta, b, link, nstart, ncores)
# 1 start
library(morpheus)
K <- ncol(fargs$beta)
# 1 start
library(morpheus)
K <- ncol(fargs$beta)
+ mu <- computeMu(io$X, io$Y, list(K=K, M=fargs$M))
op <- optimParams(fargs$X, fargs$Y, K, fargs$link, fargs$M)
op <- optimParams(fargs$X, fargs$Y, K, fargs$link, fargs$M)
- x_init <- list(p=rep(1/K,K-1), beta=
fargs$
mu, b=rep(0,K))
+ x_init <- list(p=rep(1/K,K-1), beta=mu, b=rep(0,K))
res2 <- NULL
tryCatch({
res2 <- do.call(rbind, op$run(x_init))
res2 <- NULL
tryCatch({
res2 <- do.call(rbind, op$run(x_init))
@@
-52,8
+53,6
@@
testMultistart <- function(N, n, p, beta, b, link, nstart, ncores)
library(morpheus)
io = generateSampleIO(fargs$n, fargs$p, fargs$beta, fargs$b, fargs$link)
fargs$M <- computeMoments(io$X, io$Y)
library(morpheus)
io = generateSampleIO(fargs$n, fargs$p, fargs$beta, fargs$b, fargs$link)
fargs$M <- computeMoments(io$X, io$Y)
- mu <- computeMu(io$X, io$Y, list(M=fargs$M))
- fargs$mu <- mu
fargs$X <- io$X
fargs$Y <- io$Y
fargs
fargs$X <- io$X
fargs$Y <- io$Y
fargs