我正在将一些LAB代码移植到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)))
它会产生相同的数字，直到列符号(如LAB文档所预期):
| 1 | 0 | 0 | 0 |
|---|---|---|---|
| 1 | 1 | 0 | 0 |
| 1 | 2 | 1 | 0 |
| 1 | 3 | 3 | 1 |
可以产生具有正确列符号的预期结果的等效R代码是什么？
我正在将一些LAB代码移植到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)))
它会产生相同的数字，直到列符号(如LAB文档所预期):
| 1 | 0 | 0 | 0 |
|---|---|---|---|
| 1 | 1 | 0 | 0 |
| 1 | 2 | 1 | 0 |
| 1 | 3 | 3 | 1 |
可以产生具有正确列符号的预期结果的等效R代码是什么？
您可以try pracma套餐中的pascal个
> 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