Ich habe eine n x p Matrix und möchte die als definierte n x n Matrix B berechnen

B[i, j] = f(A[i,], A[j,])

Dabei ist f eine Funktion, die Argumente mit der entsprechenden Dimensionalität akzeptiert. Gibt es einen tollen Trick, um dies in R zu berechnen? f ist symmetrisch und positiv-definit (wenn dies bei der Berechnung helfen kann).

EDIT: Praneet gebeten, f anzugeben. Das ist ein guter Punkt. Obwohl ich denke, dass es interessant wäre, eine effiziente Lösung für jede Funktion zu haben, würde ich in dem wichtigen Fall, in dem f (x, y) base :: norm ist (xy, type = 'F' ).