我想得到分位数回归系数的95%置信区间。您可以在R中使用rq包的quantreg函数计算分位数回归(与OLS模型相比):LM<-lm(mpg~disp, data = mtcars)我可以用限制函数得到线性模型的95%置信区间:当我使用分位数回归时,我理解以下代码会产生引导标准错误:
summary.rq(QR,se="boot"
本文尝试用QR迭代法求解特征值和特征向量,代码非常简单。但是,与linalg.eigs的答案相比,QR迭代的答案总是有一些相反的或不正确的值。import numpy as np for i in range(100):Q, R = linalg.qr(A) return np.diag(R), Q
c