-1

単純な線形モデルを仮定します: Stay ~ Diet + Smoke + Diet:Smoke 問題は、係数 (境界線) の行列で次元名を渡すことです。

Rのskrypt: (syntaks.txt)

BaseFrame=read.csv("DietSmoke.txt")

modd<-glm(Stay ~ Diet + Smoke + Diet:Smoke , BaseFrame, family = gaussian) 

coeff=summary(modd)$coef

lista<-list(r1=coeff,r2=dimnames(coeff))

RGui コンソールでの結果:

リスト
$r1

              Estimate Std. Error       t value  Pr(>|t|)
(Intercept)   3.666667e+00   2.441159  1.502019e+00 0.1538496
DietB         1.833333e+00   2.989797  6.131965e-01 0.5489321
DietC         2.666667e+00   3.452321  7.724273e-01 0.4518611
SmokeY        5.333333e+00   3.452321  1.544855e+00 0.1432147
DietB:SmokeY -5.166667e+00   4.566991 -1.131306e+00 0.2756896
DietC:SmokeY  2.756121e-15   4.882319  5.645107e-16 1.0000000

$r2
$r2[[1]]

[1] "(Intercept)" "DietB"   "DietC"   "SmokeY" "DietB:SmokeY" "DietC:SmokeY"

$r2[[2]]

[1] "Estimate"   "Std. Error" "t value"    "Pr(>|t|)"

RGui コンソール オブジェクト (r1、r2) に表示されます。オブジェクト r2 は、行名の文字列配列です。R 構文を Rcaller に配置すると、オブジェクト r2 は Java に渡されません。リストの代わりに、オブジェクト dimnames(coeff) を単独で配置すると、同じです。

Java コード:

RCaller caller = new RCaller();
caller.setRscriptExecutable("C:\\R\\R-3.1.2\\bin\\Rscript.exe");
RCode code = new RCode();
code.addRCode(syntaks);
caller.setRCode(code);
caller.runAndReturnResult("lista");
System.out.println(caller.getParser().getNames());

Eclipse コンソールに「[]」が表示される

オブジェクトのリストを使用すると、一部のオブジェクトのみが表示される別の例を追加できます。

リクエストに感謝します。

4

1 に答える 1

0

あなたの例の問題は、その部分が

lista <- list(r1 = coeff, r2 = dimnames(coeff))

list オブジェクトを list に含めます。つまり、

dimnames(coeff)

はリスト自体であり、RCaller で返されるオブジェクトはプリミティブ (行列とベクトル) のリストでなければなりません。lista の定義を変更する

lista <- list(r1=coeff, r2=dimnames(coeff)[[1]], r3=dimnames(coeff)[[2]])

lista を行列と 2 つの文字列ベクトルのリストにします。私はあなたが使用するデータを持っていないので、いくつかのランダムなデータであなたの例を複製します:

    RCaller caller = RCaller.create();
    RCode code = RCode.create();

    code.addDoubleArray("Stay", getRandomDouble(30));
    code.addDoubleArray("Diet", getRandomDouble(30));
    code.addDoubleArray("Smoke", getRandomDouble(30));
    code.addRCode("BaseFrame <- as.data.frame(cbind(Stay,Diet,Smoke))");
    code.addRCode("modd <- glm(Stay ~ Diet + Smoke + Diet:Smoke, BaseFrame, family = gaussian)");
    code.addRCode("coeff <- summary(modd)$coef");
    code.addRCode("lista<-list(r1 = coeff, r2 = dimnames(coeff)[[1]], r3 = dimnames(coeff)[[2]])");
    caller.setRCode(code);
    caller.runAndReturnResult("lista");

出力は次の方法で処理できます。

    System.out.println("r1:");
    double[][] data = caller.getParser().getAsDoubleMatrix("r1");
    for (int i = 0; i < data.length; i++){
        for (int j = 0; j < data[0].length; j++){
            System.out.print(data[i][j] + " ");
        }
        System.out.println();
    }

    System.out.println("r2:");
    String[] names = caller.getParser().getAsStringArray("r2");
    for (int i = 0; i < names.length; i++){
        System.out.println(names[i]);
    }

    System.out.println("r3:");
    names = caller.getParser().getAsStringArray("r3");
    for (int i = 0; i < names.length; i++){
        System.out.println(names[i]);
    }

出力は次のとおりです。

r1:
59.2591703775371 0.12534103762444 0.0438737218457358 -0.0042978710538064 
18.7562448759218 0.3291208843117 0.31313594046636 0.00547569333255683 
3.1594368046245 0.38083586791087 0.140110783132699 -0.784899882587023 
0.00398304744186597 0.706419085801913 0.889652249263241 0.439606023811926 
r2:
(Intercept)
Diet
Smoke
Diet:Smoke
r3:
Estimate
Std. Error
t value
Pr(>|t|)
于 2016-09-05T10:54:51.670 に答える