如题，现在手上有一份区域的年度数据（数据样本量不好拓展），需要分析自变量对因变量的影响，打算采用的是多元线性回归。但是变量数k大于样本数n，使用逐步回归的方式最终的推荐结果是使用n-1个变量（有个常数项，加上就是n个变量）。相当于是n个方程联立求解，最终的结果表现特别好，r方等于1，p值为0。但是符号不符合预期。请问这种情况可以怎么处理呢？是否有其他更好的方法？