OpenXM/src/hgm/mh/src/mh-r.c - diff

Return to mh-r.c CVS log

Up to [local] / OpenXM / src / hgm / mh / src

-version 1.4, 2013/03/05 06:00:52
+version 1.7, 2013/03/08 04:54:01
 Line 1
 Line 1
 Line 1
- /* $OpenXM: OpenXM/src/hgm/mh/src/mh-r.c,v 1.3 2013/03/01 07:20:41 takayama Exp $
+ /* $OpenXM: OpenXM/src/hgm/mh/src/mh-r.c,v 1.6 2013/03/07 05:23:31 takayama Exp $
   R interface module
  */
 Line 17  int Rmh_cwishart_gen(int *mp,int *np,double *beta,doub
 Line 17  int Rmh_cwishart_gen(int *mp,int *np,double *beta,doub
 Line 17  int Rmh_cwishart_gen(int *mp,int *np,double *beta,doub
    cw = mh_cwishart_gen(*mp,*np,beta,*x0p,*approxDegp,*hp,*dpp,*xp,modep);
    xy[0] = cw->x;
    for (i=1; i<=rank; i++) xy[i] = (cw->f)[i-1];
+   if ((modep[2] > 0) && cw->aux) {
+     struct SFILE *sfp3;
+     char *s3;
+     char str[1024];
+     double x;
+     s3 = (char *)cw->aux;
+     sfp3 = mh_fopen(s3,"r",0);
+     for (i=cw->rank+1; i<modep[2]+(cw->rank)+1; i++) xy[i] = 0.0;
+     for (i=cw->rank+1; i<modep[2]+(cw->rank)+1; i++) {
+       if (!mh_fgets(str,1024,sfp3)) break;
+       sscanf(str,"%lg",&x);
+       xy[i] = x;
+     }
+   }
+   return(0);
+ }
+ int Rmh_set_strategy(int *m,double *err,double *ans) {
+   mh_set_strategy(*m,err);
+   ans[0] = 0.0;
+   return(0);
  }