我想知道是否有人知道ASReml软件包中的pol()函数的R代码是什么(ASReml是由VSNi开发的付费非开源软件包)。
当拟合随机回归模型时,可以调用pol()函数从数值矢量生成正交多项式。
我从VSNi(https://forum.vsni.co.uk/viewtopic.php?t=1075)在线论坛上找到了pol()函数的以下代码。
尝试代码时,我注意到如果重新排序输入向量,则输出多项式将发生变化。但是,如果我重新排序用于拟合随机回归模型的数据框,则模型输出(包括方差分量和随机效应系数)根本不会改变。因此,我想知道pol()函数的实际代码是否在进行转换之前将输入向量从低到高排序?