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:
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 .