| version 1.8, 2014/03/16 03:11:07 |
version 1.9, 2016/02/13 05:55:09 |
|
|
| /* $OpenXM: OpenXM/src/hgm/mh/src/mh-r.c,v 1.7 2013/03/08 04:54:01 takayama Exp $ |
/* $OpenXM: OpenXM/src/hgm/mh/src/mh-r.c,v 1.8 2014/03/16 03:11:07 takayama Exp $ |
| R interface module |
R interface module |
| */ |
*/ |
| |
|
| Line 41 int Rmh_cwishart_gen(int *mp,int *np,double *beta,doub |
|
| Line 41 int Rmh_cwishart_gen(int *mp,int *np,double *beta,doub |
|
| int Rmh_set_strategy(int *m,double *err,double *ans) { |
int Rmh_set_strategy(int *m,double *err,double *ans) { |
| mh_set_strategy(*m,err); |
mh_set_strategy(*m,err); |
| ans[0] = 0.0; |
ans[0] = 0.0; |
| |
return(0); |
| |
} |
| |
|
| |
int Rmh_pFq_gen(int *mp, |
| |
int *pp, double *a, |
| |
int *qp, double *b, |
| |
int *ef_typep, |
| |
double *beta,double *x0p, |
| |
int *approxDegp,double *hp, int *dpp, double *xp, |
| |
int *modep,int *rankp, |
| |
int *automaticp,double *assigned_series_errorp,int *verbosep, |
| |
double *xy) { |
| |
struct cWishart *cw; |
| |
int rank; |
| |
int i; |
| |
rank = *rankp; |
| |
cw = mh_pFq_gen(*mp, |
| |
*pp,a, |
| |
*qp,b, |
| |
*ef_typep, |
| |
beta,*x0p,*approxDegp,*hp,*dpp,*xp,modep, |
| |
*automaticp,*assigned_series_errorp,*verbosep); |
| |
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); |
return(0); |
| } |
} |
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to mh-r.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM / src / hgm / mh / src