===================================================================
RCS file: /home/cvs/OpenXM/src/hgm/mh/src/mh-r.c,v
retrieving revision 1.4
retrieving revision 1.5
diff -u -p -r1.4 -r1.5
--- OpenXM/src/hgm/mh/src/mh-r.c	2013/03/05 06:00:52	1.4
+++ OpenXM/src/hgm/mh/src/mh-r.c	2013/03/05 07:03:37	1.5
@@ -1,4 +1,4 @@
-/* $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.4 2013/03/05 06:00:52 takayama Exp $
  R interface module
 */
 
@@ -17,4 +17,20 @@ 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;
+    }
+  }
+
 }