1

我已安装 JRI 以在 Windows 7 上使用 32 位 R 3.0.2 和 Java jdk1.7.0_45 与 NetBeans 7.4 一起运行。

我正在使用以下 Java 代码

     REXP load=re.eval("source('C:\\\\SearchPath\\\\gammaDistAnova.r')");
     String errStr=load.asString();
     REXP stats=re.eval("gammaDistAnova.getStats(ref, target)");

调用以下 R 脚本。

gammaDistAnova.getStats<-function(ref, target){
  library("MASS") # Library containing fitdistr
  library("mgcv")
  library("stats")

  # Get mean and SD of reference
  gr=fitdistr(ref+0.00001,"gamma")  # Ref is a vector of (reference) values.  The     0.00001 is to prevent errors due to zero values
  meanRef<-gr[1]$estimate['shape']/gr[1]$estimate['rate']  # Mean of ref vector  
  SDref<-sqrt(gr[1]$estimate['shape']/(gr[1]$estimate['rate']^2)) # SD of ref vector

  # Get mean and SD of target
  gt=fitdistr(target+0.00001,"gamma")  # target is a vector of (target) values.  The 0.00001 is to prevent errors due to zero values
  meanTarget<-gt[1]$estimate['shape']/gt[1]$estimate['rate']  # Mean of target vector  
  SDTarget<-sqrt(gt[1]$estimate['shape']/(gt[1]$estimate['rate']^2)) # SD of target vector

  # Analysis of variance between the distributions
  n=300
  x=rgamma(n, shape=gr[1]$estimate['shape'], scale=1/gr[1]$estimate['rate'])
  y=rgamma(n, shape=gt[1]$estimate['shape'], scale=1/gt[1]$estimate['rate'])
  random1 <- sample(c("level1","level2","level3"), n, replace=TRUE)
  debug=list(random1 = ~1)
#  glmm1 <- gamm(y ~ x, random=list(random1 = ~1))
#  anova(glmm1$gam)

  out=list("refMean"=meanRef, "refSD"=SDref, "targetMean"=meanTarget, "targetSD"=SDTarget, "lx"=length(x), "ly"=length(y))
  out
}

一切运行良好（列表中返回的所有内容似乎都是有效的，并且是我所期望的），直到我取消注释

glmm1 <- gamm(y ~ x, random=list(random1 = ~1))

在这种情况下，函数返回 null，表示失败。

虽然该函数在 JRI 上失败，但它在 RStudio 版本 0.98.501 上运行没有任何问题。

编辑：

我试过了

  glmm1 <- gamm(y ~ x, random=list(random1 = ~1))
  errStr=geterrmessage()
  errStr

但

 re.eval("gammaDistAnova.getStats(ref, target)");

仍然返回 null

4

1 回答 1

1

你应该检查你的“lib”目录....\R\win-library\3.0.2

并检查您是否可以查看该mgcv软件包，如果没有，请再次将其下载到您使用的 R 版本中

一种选择是您将此库下载到早期版本的 R 中。

于 2014-04-06T13:06:36.617 回答