Wie generieren Sie eine zufällige Matrix in R mit einem konstanten colSum-Wert?

Question

Ich suche nach einer Möglichkeit, eine Reihe von Matrizen zu generieren, bei denen alle Spaltensummen einem bestimmten Wert entsprechen. Ich dachte, ich könnte Zufallsmatrizen generieren und überprüfen, ob colSums = Wert benötigt, aber das scheint nicht zu funktionieren. Wie generieren wir zufällige Matrizen in R mit positiven und negativen Werten?

Answer 1

Sie können so tun. Dies ist eine Funktion, die eine Matrix mit Zufallszahlen aus einer Normalverteilung erstellt. Dann zentriert es alle Spalten auf 0 und summiert die gewünschte Zahl.

RandomMatrixColsum = function(mycolsum=1, dim=c(10, 10), ...){ 
    m <- matrix(rnorm(prod(dim), ...), ncol=dim[2]) 
    ms <- scale(m, center=T, scale=F) + mycolsum/dim[1] 
    return(ms) 
} 
RandomMatrixColsum(8) 

Wenn Sie sicherstellen möchten, sowohl positive als auch negative Werte in Ihrer Matrix haben, können Sie die Standardabweichung der rnorm erhöhen, indem ein Argument sd=... zum Funktionsaufruf hinzufügen, oder auch die sd Kraft zu sein gleich dem Mittelwert.