Last updated: 2017-08-24
Code version: 6089d1d
par( mfrow = c(2,2), mar =c(0.1,0.1,0.66,0.1))
layout(matrix(c(1,1,1,1,1,0,2,2,2,2,2,0,0,0,3,3,3,3,3,0,0,0), 2, 11, byrow = TRUE))
In this case, we simulate data set with 100 train samples and 100 test samples.
#dsc-query omega -o methodsP300.ipynb --title "compare methods in P 300" --target toeplitz.P toeplitz.n hub_graph.P hub_graph.n dense_edors_renyi_graph.P dense_edors_renyi_graph.n sparse_edors_renyi_graph.P sparse_edors_renyi_graph.n cluster_graph.P cluster_graph.n scale_free_graph.P scale_free_graph.n band_graph.P band_graph.n diagonal.P diagonal.n loglik.score flash SFAmix tiger glasso clime hybrid_SFAmix hybrid_flash --condition "hub_graph.P = 300" "hub_graph.n = 200" "band_graph.P = 300" "band_graph.n = 200" "cluster_graph.P = 300" "cluster_graph.n = 200" "scale_free_graph.P = 300" "scale_free_graph.n = 200" "diagonal.P = 300" "diagonal.n = 200" "toeplitz.P = 300" "toeplitz.n = 200" "dense_edors_renyi_graph.P=300" "dense_edors_renyi_graph.n=200" "sparse_edors_renyi_graph.P=300" "sparse_edors_renyi_graph.n=200" --language R
source("../testdata/plotfunction.R")
data = readRDS("../testdata/omega-0811/output/methodsP300.rds")
methods = c("flash", "SFAmix", "tiger", "glasso", "clime", "hybrid_SFAmix", "hybrid_flash")
score = "loglik_score"
title = "P=300 N = 100 loglik"
plot_per_method(data, methods, score,title,T)
#dsc-query omega -o methodsP300mse.ipynb --title "compare methods in P 300 in MSE" --target toeplitz.P toeplitz.n hub_graph.P hub_graph.n dense_edors_renyi_graph.P dense_edors_renyi_graph.n sparse_edors_renyi_graph.P sparse_edors_renyi_graph.n cluster_graph.P cluster_graph.n scale_free_graph.P scale_free_graph.n band_graph.P band_graph.n diagonal.P diagonal.n mse.score flash SFAmix tiger glasso clime hybrid_SFAmix hybrid_flash --condition "hub_graph.P = 300" "hub_graph.n = 200" "band_graph.P = 300" "band_graph.n = 200" "cluster_graph.P = 300" "cluster_graph.n = 200" "scale_free_graph.P = 300" "scale_free_graph.n = 200" "diagonal.P = 300" "diagonal.n = 200" "toeplitz.P = 300" "toeplitz.n = 200" "dense_edors_renyi_graph.P=300" "dense_edors_renyi_graph.n=200" "sparse_edors_renyi_graph.P=300" "sparse_edors_renyi_graph.n=200" --language R
source("../testdata/plotfunction.R")
data = readRDS("../testdata/omega-0811/output/methodsP300mse.rds")
methods = c("flash", "SFAmix", "tiger", "glasso", "clime", "hybrid_SFAmix", "hybrid_flash")
score = "mse_score"
title = "P=300 N = 100 MSE"
plot_per_method(data, methods, score,title,F)
#dsc-query omega -o methodsP300hyvarinen.ipynb --title "compare methods in P 300 in hyvarinen" --target toeplitz.P toeplitz.n hub_graph.P hub_graph.n dense_edors_renyi_graph.P dense_edors_renyi_graph.n sparse_edors_renyi_graph.P sparse_edors_renyi_graph.n cluster_graph.P cluster_graph.n scale_free_graph.P scale_free_graph.n band_graph.P band_graph.n diagonal.P diagonal.n hyvarinen.score flash SFAmix tiger glasso clime hybrid_SFAmix hybrid_flash --condition "hub_graph.P = 300" "hub_graph.n = 200" "band_graph.P = 300" "band_graph.n = 200" "cluster_graph.P = 300" "cluster_graph.n = 200" "scale_free_graph.P = 300" "scale_free_graph.n = 200" "diagonal.P = 300" "diagonal.n = 200" "toeplitz.P = 300" "toeplitz.n = 200" "dense_edors_renyi_graph.P=300" "dense_edors_renyi_graph.n=200" "sparse_edors_renyi_graph.P=300" "sparse_edors_renyi_graph.n=200" --language R
source("../testdata/plotfunction.R")
data = readRDS("../testdata/omega-0811/output/methodsP300hyvarinen.rds")
methods = c("flash", "SFAmix", "tiger", "glasso", "clime", "hybrid_SFAmix", "hybrid_flash")
score = "hyvarinen_score"
title = "P=300 N = 100 hyvarinen"
plot_per_method(data, methods, score,title,F)
#dsc-query omega -o methodsP1000loglik.ipynb --title "compare methods in P 1000 in loglik" --target toeplitz.P toeplitz.n hub_graph.P hub_graph.n dense_edors_renyi_graph.P dense_edors_renyi_graph.n sparse_edors_renyi_graph.P sparse_edors_renyi_graph.n cluster_graph.P cluster_graph.n scale_free_graph.P scale_free_graph.n band_graph.P band_graph.n diagonal.P diagonal.n loglik.score flash SFAmix tiger glasso clime hybrid_SFAmix hybrid_flash --condition "hub_graph.P = 1000" "hub_graph.n = 200" "band_graph.P = 1000" "band_graph.n = 200" "cluster_graph.P = 1000" "cluster_graph.n = 200" "scale_free_graph.P = 1000" "scale_free_graph.n = 200" "diagonal.P = 1000" "diagonal.n = 200" "toeplitz.P = 1000" "toeplitz.n = 200" "dense_edors_renyi_graph.P=1000" "dense_edors_renyi_graph.n=200" "sparse_edors_renyi_graph.P=1000" "sparse_edors_renyi_graph.n=200" --language R
source("../testdata/plotfunction.R")
data = readRDS("../testdata/omega-0811/output/methodsP1000loglik.rds")
methods = c("flash", "SFAmix", "tiger", "glasso", "clime", "hybrid_SFAmix", "hybrid_flash")
score = "loglik_score"
title = "P=1000 N = 100 loglik"
# plot_per_method(data, methods, score,title,T,T)
plot_per_method(data, methods, score,title,T,F)
#dsc-query omega -o methodsP1000mse.ipynb --title "compare methods in P 1000 in MSE" --target toeplitz.P toeplitz.n hub_graph.P hub_graph.n dense_edors_renyi_graph.P dense_edors_renyi_graph.n sparse_edors_renyi_graph.P sparse_edors_renyi_graph.n cluster_graph.P cluster_graph.n scale_free_graph.P scale_free_graph.n band_graph.P band_graph.n diagonal.P diagonal.n mse.score flash SFAmix tiger glasso clime hybrid_SFAmix hybrid_flash --condition "hub_graph.P = 1000" "hub_graph.n = 200" "band_graph.P = 1000" "band_graph.n = 200" "cluster_graph.P = 1000" "cluster_graph.n = 200" "scale_free_graph.P = 1000" "scale_free_graph.n = 200" "diagonal.P = 1000" "diagonal.n = 200" "toeplitz.P = 1000" "toeplitz.n = 200" "dense_edors_renyi_graph.P=1000" "dense_edors_renyi_graph.n=200" "sparse_edors_renyi_graph.P=1000" "sparse_edors_renyi_graph.n=200" --language R
source("../testdata/plotfunction.R")
data = readRDS("../testdata/omega-0811/output/methodsP1000mse.rds")
methods = c("flash", "SFAmix", "tiger", "glasso", "clime", "hybrid_SFAmix", "hybrid_flash")
score = "mse_score"
title = "P=1000 N = 100 MSE"
# plot_per_method(data, methods, score,title,F,T)
plot_per_method(data, methods, score,title,F,F)
#dsc-query omega -o methodsP1000hyvarinen.ipynb --title "compare methods in P 1000 in hyvarinen" --target toeplitz.P toeplitz.n hub_graph.P hub_graph.n dense_edors_renyi_graph.P dense_edors_renyi_graph.n sparse_edors_renyi_graph.P sparse_edors_renyi_graph.n cluster_graph.P cluster_graph.n scale_free_graph.P scale_free_graph.n band_graph.P band_graph.n diagonal.P diagonal.n hyvarinen.score flash SFAmix tiger glasso clime hybrid_SFAmix hybrid_flash --condition "hub_graph.P = 1000" "hub_graph.n = 200" "band_graph.P = 1000" "band_graph.n = 200" "cluster_graph.P = 1000" "cluster_graph.n = 200" "scale_free_graph.P = 1000" "scale_free_graph.n = 200" "diagonal.P = 1000" "diagonal.n = 200" "toeplitz.P = 1000" "toeplitz.n = 200" "dense_edors_renyi_graph.P=1000" "dense_edors_renyi_graph.n=200" "sparse_edors_renyi_graph.P=1000" "sparse_edors_renyi_graph.n=200" --language R
source("../testdata/plotfunction.R")
data = readRDS("../testdata/omega-0811/output/methodsP1000hyvarinen.rds")
methods = c("flash", "SFAmix", "tiger", "glasso", "clime", "hybrid_SFAmix", "hybrid_flash")
score = "hyvarinen_score"
title = "P=1000 N = 100 hyvarinen"
plot_per_method(data, methods, score,title,F)
#dsc-query omega -o methodsP100loglik.ipynb --title "compare methods in P 100 in loglik" --target toeplitz.P toeplitz.n hub_graph.P hub_graph.n dense_edors_renyi_graph.P dense_edors_renyi_graph.n sparse_edors_renyi_graph.P sparse_edors_renyi_graph.n cluster_graph.P cluster_graph.n scale_free_graph.P scale_free_graph.n band_graph.P band_graph.n diagonal.P diagonal.n loglik.score flash SFAmix tiger glasso clime hybrid_SFAmix hybrid_flash --condition "hub_graph.P = 100" "hub_graph.n = 200" "band_graph.P = 100" "band_graph.n = 200" "cluster_graph.P = 100" "cluster_graph.n = 200" "scale_free_graph.P = 100" "scale_free_graph.n = 200" "diagonal.P = 100" "diagonal.n = 200" "toeplitz.P = 100" "toeplitz.n = 200" "dense_edors_renyi_graph.P=100" "dense_edors_renyi_graph.n=200" "sparse_edors_renyi_graph.P=100" "sparse_edors_renyi_graph.n=200" --language R
source("../testdata/plotfunction.R")
data = readRDS("../testdata/omega-0811/output/methodsP100loglik.rds")
methods = c("flash", "SFAmix", "tiger", "glasso", "clime", "hybrid_SFAmix", "hybrid_flash")
score = "loglik_score"
title = "P=100 N = 100 loglik"
plot_per_method(data, methods, score,title,T)
#dsc-query omega -o methodsP100mse.ipynb --title "compare methods in P 100 in MSE" --target toeplitz.P toeplitz.n hub_graph.P hub_graph.n dense_edors_renyi_graph.P dense_edors_renyi_graph.n sparse_edors_renyi_graph.P sparse_edors_renyi_graph.n cluster_graph.P cluster_graph.n scale_free_graph.P scale_free_graph.n band_graph.P band_graph.n diagonal.P diagonal.n mse.score flash SFAmix tiger glasso clime hybrid_SFAmix hybrid_flash --condition "hub_graph.P = 100" "hub_graph.n = 200" "band_graph.P = 100" "band_graph.n = 200" "cluster_graph.P = 100" "cluster_graph.n = 200" "scale_free_graph.P = 100" "scale_free_graph.n = 200" "diagonal.P = 100" "diagonal.n = 200" "toeplitz.P = 100" "toeplitz.n = 200" "dense_edors_renyi_graph.P=100" "dense_edors_renyi_graph.n=200" "sparse_edors_renyi_graph.P=100" "sparse_edors_renyi_graph.n=200" --language R
source("../testdata/plotfunction.R")
data = readRDS("../testdata/omega-0811/output/methodsP100mse.rds")
methods = c("flash", "SFAmix", "tiger", "glasso", "clime", "hybrid_SFAmix", "hybrid_flash")
score = "mse_score"
title = "P=100 N = 100 MSE"
plot_per_method(data, methods, score,title,F)
#dsc-query omega -o methodsP100hyvarinen.ipynb --title "compare methods in P 100 in hyvarinen" --target toeplitz.P toeplitz.n hub_graph.P hub_graph.n dense_edors_renyi_graph.P dense_edors_renyi_graph.n sparse_edors_renyi_graph.P sparse_edors_renyi_graph.n cluster_graph.P cluster_graph.n scale_free_graph.P scale_free_graph.n band_graph.P band_graph.n diagonal.P diagonal.n hyvarinen.score flash SFAmix tiger glasso clime hybrid_SFAmix hybrid_flash --condition "hub_graph.P = 100" "hub_graph.n = 200" "band_graph.P = 100" "band_graph.n = 200" "cluster_graph.P = 100" "cluster_graph.n = 200" "scale_free_graph.P = 100" "scale_free_graph.n = 200" "diagonal.P = 100" "diagonal.n = 200" "toeplitz.P = 100" "toeplitz.n = 200" "dense_edors_renyi_graph.P=100" "dense_edors_renyi_graph.n=200" "sparse_edors_renyi_graph.P=100" "sparse_edors_renyi_graph.n=200" --language R
source("../testdata/plotfunction.R")
data = readRDS("../testdata/omega-0811/output/methodsP100hyvarinen.rds")
methods = c("flash", "SFAmix", "tiger", "glasso", "clime", "hybrid_SFAmix", "hybrid_flash")
score = "hyvarinen_score"
title = "P=100 N = 100 hyvarinen"
plot_per_method(data, methods, score,title,F)
#dsc-query omega -o methodsP10loglik.ipynb --title "compare methods in P 10 in loglik" --target toeplitz.P toeplitz.n hub_graph.P hub_graph.n dense_edors_renyi_graph.P dense_edors_renyi_graph.n sparse_edors_renyi_graph.P sparse_edors_renyi_graph.n cluster_graph.P cluster_graph.n scale_free_graph.P scale_free_graph.n band_graph.P band_graph.n diagonal.P diagonal.n loglik.score flash SFAmix tiger glasso inv_cov clime hybrid_SFAmix hybrid_flash --condition "hub_graph.P = 10" "hub_graph.n = 200" "band_graph.P = 10" "band_graph.n = 200" "cluster_graph.P = 10" "cluster_graph.n = 200" "scale_free_graph.P = 10" "scale_free_graph.n = 200" "diagonal.P = 10" "diagonal.n = 200" "toeplitz.P = 10" "toeplitz.n = 200" "dense_edors_renyi_graph.P=10" "dense_edors_renyi_graph.n=200" "sparse_edors_renyi_graph.P=10" "sparse_edors_renyi_graph.n=200" --language R
source("../testdata/plotfunction.R")
data = readRDS("../testdata/omega-0811/output/methodsP10loglik.rds")
methods = c("flash", "SFAmix", "tiger", "glasso", "clime", "hybrid_SFAmix", "hybrid_flash","inv_cov")
score = "loglik_score"
title = "P=10 N = 100 loglik"
plot_per_method(data, methods, score,title,T)
#dsc-query omega -o methodsP10mse.ipynb --title "compare methods in P 10 in MSE" --target toeplitz.P toeplitz.n hub_graph.P hub_graph.n dense_edors_renyi_graph.P dense_edors_renyi_graph.n sparse_edors_renyi_graph.P sparse_edors_renyi_graph.n cluster_graph.P cluster_graph.n scale_free_graph.P scale_free_graph.n band_graph.P band_graph.n diagonal.P diagonal.n mse.score flash SFAmix tiger glasso inv_cov clime hybrid_SFAmix hybrid_flash --condition "hub_graph.P = 10" "hub_graph.n = 200" "band_graph.P = 10" "band_graph.n = 200" "cluster_graph.P = 10" "cluster_graph.n = 200" "scale_free_graph.P = 10" "scale_free_graph.n = 200" "diagonal.P = 10" "diagonal.n = 200" "toeplitz.P = 10" "toeplitz.n = 200" "dense_edors_renyi_graph.P=10" "dense_edors_renyi_graph.n=200" "sparse_edors_renyi_graph.P=10" "sparse_edors_renyi_graph.n=200" --language R
source("../testdata/plotfunction.R")
data = readRDS("../testdata/omega-0811/output/methodsP10mse.rds")
methods = c("flash", "SFAmix", "tiger", "glasso", "clime", "hybrid_SFAmix", "hybrid_flash","inv_cov")
score = "mse_score"
title = "P=10 N = 100 MSE"
plot_per_method(data, methods, score,title,F)
# plot_per_method(data, methods, score,title,F,T)
#dsc-query omega -o methodsP10hyvarinen.ipynb --title "compare methods in P 10 in hyvarinen" --target toeplitz.P toeplitz.n hub_graph.P hub_graph.n dense_edors_renyi_graph.P dense_edors_renyi_graph.n sparse_edors_renyi_graph.P sparse_edors_renyi_graph.n cluster_graph.P cluster_graph.n scale_free_graph.P scale_free_graph.n band_graph.P band_graph.n diagonal.P diagonal.n hyvarinen.score flash SFAmix tiger glasso clime inv_cov hybrid_SFAmix hybrid_flash --condition "hub_graph.P = 10" "hub_graph.n = 200" "band_graph.P = 10" "band_graph.n = 200" "cluster_graph.P = 10" "cluster_graph.n = 200" "scale_free_graph.P = 10" "scale_free_graph.n = 200" "diagonal.P = 10" "diagonal.n = 200" "toeplitz.P = 10" "toeplitz.n = 200" "dense_edors_renyi_graph.P=10" "dense_edors_renyi_graph.n=200" "sparse_edors_renyi_graph.P=10" "sparse_edors_renyi_graph.n=200" --language R
source("../testdata/plotfunction.R")
data = readRDS("../testdata/omega-0811/output/methodsP10hyvarinen.rds")
methods = c("flash", "SFAmix", "tiger", "glasso", "clime", "hybrid_SFAmix", "hybrid_flash","inv_cov")
score = "hyvarinen_score"
title = "P=10 N = 100 hyvarinen"
plot_per_method(data, methods, score,title,F)
for the same scale, there are some problems. we can not tell the difference from some relavant small change. one can just set same_scale=T
to see the result.
sessionInfo()
R version 3.3.0 (2016-05-03)
Platform: x86_64-apple-darwin13.4.0 (64-bit)
Running under: OS X 10.12.6 (unknown)
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] fmsb_0.6.1 bindrcpp_0.2 dplyr_0.7.2 reshape2_1.4.2
[5] cowplot_0.8.0 ggplot2_2.2.1 workflowr_0.4.0 rmarkdown_1.6
loaded via a namespace (and not attached):
[1] Rcpp_0.12.12 bindr_0.1 knitr_1.17 magrittr_1.5
[5] munsell_0.4.3 colorspace_1.3-2 R6_2.2.2 rlang_0.1.2
[9] stringr_1.2.0 plyr_1.8.4 tools_3.3.0 grid_3.3.0
[13] gtable_0.2.0 git2r_0.19.0 htmltools_0.3.6 assertthat_0.2.0
[17] yaml_2.1.14 lazyeval_0.2.0 rprojroot_1.2 digest_0.6.12
[21] tibble_1.3.3 glue_1.1.1 evaluate_0.10.1 labeling_0.3
[25] stringi_1.1.5 scales_0.4.1 backports_1.1.0 pkgconfig_2.0.1
This R Markdown site was created with workflowr