Вызов rchisq с параметрами вектора

Question

Что происходит, когда я звоню rchisq(100,1:100,1:100):

• генерирует ли R 100 чисел с df = 1, ncp = 1?
• или он генерирует 100 чисел с df = k, ncp = k для k = 1 ... 100?

Я хочу знать, могут ли df и ncp быть векторами или нет. Это не ясно в документации (по сравнению с rnorm). Я подозреваю, что они также могут быть векторами, и переработка происходит, если длины различаются (?)

Answer 1

Просто протестируйте его, используя set.seed():

R> set.seed(42); rchisq(2, 1:2, 1:2)              # base case
[1] 8.676 1.653
R> set.seed(42); rchisq(1, 1, 1); rchisq(1, 2, 2) # matches
[1] 8.676
[1] 1.653
R> set.seed(42); rchisq(1, 1, 1); rchisq(1, 1, 1) # does not match
[1] 8.676
[1] 0.5874

, чтобы было похоже, что вы получаете N рисунков с использованием степеней свободы и нецентральности из соответствующего значения в предоставленном векторе.

С другой стороны, рециркуляция, как на самом деле правила, так как скалярные значения для df и ncp рециркулируются до длины вектора, что имеет некоторый смысл.