R 相当于 MATLAB 的 pascal(n, 1) 不同的列符号

Question

我正在将一些 MATLAB 代码移植到 R，并且遇到了这个函数：

pascal(4, 1)

根据文档，“返回 Pascal 矩阵的下三角 Cholesky 因子（直到列的符号）”并产生

1	0	0	0
1	-1	0	0
1	-2	1	0
1	-3	3	-1

我认为等效的 R 代码是

t(chol(pascal(4)))

它会产生相同的数字，直到列符号为止，如 MATLAB 文档所预期的那样：

1	0	0	0
1	1	0	0
1	2	1	0
1	3	3	1

产生具有正确列符号的预期结果的等效 R 代码是什么？

Answer 1

您可以尝试包装中的

pascal

pracma

> library(pracma)

> pascal(4, 1)
     [,1] [,2] [,3] [,4]
[1,]    1    0    0    0
[2,]    1   -1    0    0
[3,]    1   -2    1    0
[4,]    1   -3    3   -1