Matrixmultiplikation mit Rcpp - verschiedene Werte bei der Zuweisung der Ausgabe

8

Ich schrieb die folgende Funktion, um eine Matrix zu berechnen (Informationsmatrix des Weibull-Modells)

%Vor%

Nun möchte ich die Matrix für einige Werte in R berechnen:

%Vor%

Diese Matrix ist in Ordnung. Nun, wenn ich zuerst die Ausgabe zuweisen und dann drucken, werde ich eine andere Matrix haben (einige der Elemente werden in einige große Werte geändert werden!).

%Vor%

Wie zu sehen ist, erhalten die Elemente (1, 2), (1, 4), (2, 2) und (3, 2) andere große Werte. Ich würde jede Hilfe schätzen.

Update1 als Antwort auf @Ronald:

  • version.string R Version 3.1.2 (2014-10-31)
  • Windows 7 x64
  • Plattform x86_64-w64-mingw32
  • Rcpp_0.11.3, RcppEigen_0.3.2.3.0, tools_3.1.2
  • Rtools Version 3.2.0.1948
Ehsan Masoudi 22.04.2015, 11:30
quelle

1 Antwort

9

Sie haben vergessen, Fisher_mat mit Nullen vor der Summierung zu initialisieren, so dass es manchmal zufälligen Müll enthält. Hinzufügen von z.B. Fisher_mat.setZero(); vor einer for Schleife gibt Ihnen eine konsistente Ausgabe.

Überprüfen Sie diese Referenz: RcppEigen Einführung .

    
tonytonov 22.04.2015, 12:56
quelle

Tags und Links