projects
/
morpheus.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
Fix multistart script
[morpheus.git]
/
reports
/
accuracy.R
diff --git
a/reports/accuracy.R
b/reports/accuracy.R
index
96f3a4a
..
24191d8
100644
(file)
--- a/
reports/accuracy.R
+++ b/
reports/accuracy.R
@@
-29,12
+29,12
@@
optimBeta <- function(N, n, p, beta, b, link, ncores)
tryCatch({
fm <- flexmix( cbind(V1, 1-V1) ~ .-V1, data=dat, k=K,
model = FLXMRglm(family = binomial(link = link)) )
tryCatch({
fm <- flexmix( cbind(V1, 1-V1) ~ .-V1, data=dat, k=K,
model = FLXMRglm(family = binomial(link = link)) )
- p <- mean(fm@posterior[["scaled"]][,1])
+ p
f
<- mean(fm@posterior[["scaled"]][,1])
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
})