Revision 1.1, Sun Mar 23 10:17:40 2014 UTC (10 years, 3 months ago) by tkoyama
Branch: MAIN
CVS Tags: RELEASE_1_3_1_13b, HEAD

hgm_ko_orthant.c can work without gsl,
and hgm_ko_orthant.so generated by the make command
can be called from R. For example,
> source("ko_orthant.R")
> x <- matrix(c(1,0.5,0.5,1),nrow=2)
> y <- c(1.0,0.5)
> HgmKoNcorthant(x,y)
[1] 0.6301772
We also use it from command line:
  ./hgm_ko_orthant sample_data1.txt

However, hgm_ko_orthant.c may have some new bugs,
and I wil test it more.