projects
/
morpheus.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
Fix accuracy.R + add postTreatment draft
[morpheus.git]
/
reports
/
accuracy.R
diff --git
a/reports/accuracy.R
b/reports/accuracy.R
index
96f3a4a
..
63cd5aa
100644
(file)
--- a/
reports/accuracy.R
+++ b/
reports/accuracy.R
@@
-10,7
+10,7
@@
optimBeta <- function(N, n, p, beta, b, link, ncores)
K <- ncol(fargs$beta)
M <- computeMoments(fargs$X, fargs$Y)
mu <- computeMu(fargs$X, fargs$Y, list(K=K, M=M))
K <- ncol(fargs$beta)
M <- computeMoments(fargs$X, fargs$Y)
mu <- computeMu(fargs$X, fargs$Y, list(K=K, M=M))
- op <- optimParams(fargs$X, fargs$Y, K, fargs$link, M
)
+ op <- optimParams(fargs$X, fargs$Y, K, fargs$link, M
, 1) #only 1 OpenMP core
x_init <- list(p=rep(1/K,K-1), beta=mu, b=rep(0,K))
res2 <- NULL
tryCatch({
x_init <- list(p=rep(1/K,K-1), beta=mu, b=rep(0,K))
res2 <- NULL
tryCatch({
@@
-27,14
+27,14
@@
optimBeta <- function(N, n, p, beta, b, link, ncores)
dat <- as.data.frame( cbind(fargs$Y,fargs$X) )
res2 <- NULL
tryCatch({
dat <- as.data.frame( cbind(fargs$Y,fargs$X) )
res2 <- NULL
tryCatch({
- fm <- flexmix( cbind(V1, 1-V1) ~ .
-V1
, data=dat, k=K,
+ fm <- flexmix( cbind(V1, 1-V1) ~ ., data=dat, k=K,
model = FLXMRglm(family = binomial(link = link)) )
model = FLXMRglm(family = binomial(link = link)) )
- p
<- mean(fm@posterior[["scaled"]][,1
])
+ p
f <- colMeans(fm@posterior[["scaled"]
])
out <- refit(fm)
beta_b <- sapply( seq_len(K), function(i) {
as.double( out@components[[1]][[i]][,1] )
} )
out <- refit(fm)
beta_b <- sapply( seq_len(K), function(i) {
as.double( out@components[[1]][[i]][,1] )
} )
- res2 <- rbind(p, beta_b[2:nrow(beta_b),], beta_b[1,])
+ res2 <- rbind(p
f
, beta_b[2:nrow(beta_b),], beta_b[1,])
}, error = function(e) {
res2 <- NA
})
}, error = function(e) {
res2 <- NA
})