1

Windows 7 で 32 ビット R 3.0.2 と Java jdk1.7.0_45 を使用して、NetBeans 7.4 で実行する JRI をインストールしました。

次の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使用バージョンに再度ダウンロードします

1 つのオプションは、このライブラリを以前のバージョンの R にダウンロードしたことです。

于 2014-04-06T13:06:36.617 に答える