คุณจะสร้างเมทริกซ์สุ่มใน R ด้วยค่า colSum คงที่ได้อย่างไร

คุณสามารถทำเช่นนี้ นี่คือฟังก์ชันที่สร้างเมทริกซ์ที่เต็มไปด้วยตัวเลขสุ่มจากการแจกแจงแบบปกติ จากนั้นจะจัดคอลัมน์ทั้งหมดให้อยู่ตรงกลางเป็น 0 และเพิ่มจำนวนที่ต้องการ

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)

หากคุณต้องการให้แน่ใจว่ามีทั้งค่าบวกและลบในเมทริกซ์ คุณสามารถเพิ่มค่าเบี่ยงเบนมาตรฐานของ rnorm ได้โดยการเพิ่มอาร์กิวเมนต์ sd=... ในการเรียกใช้ฟังก์ชัน หรือบังคับให้ sd เท่ากับค่าเฉลี่ยด้วย