在python中求拟合模型的协方差矩阵(相当于python中的vcov()(R function))
lmfit <- lm(formula = Y ~ X, data=Data_df)
lmpred <- predict(lmfit, newdata=Data_df, se.fit=TRUE, interval = "prediction")
std_er <- sqrt(((X0) %*% vcov(lmfit)) %*% t(X0))
尝试用python转换上面的代码。为此,我需要找到拟合模型的协方差矩阵,即vcov。 我不能使用np.cov公司()试图找到模型的协方差矩阵。在
我已经用过了statsmodels.回归.线性模型.OLSResults.cov_参数(),但我得到的值与R中的值不同
scipy ODR代码可以独立计算参数协方差矩阵,这里是从my的源代码中提取的一个示例尊尊网在线曲线拟合员:
请提供您正在使用的具体细节。在
假设您对数据使用numpy数组,这里有^{} 估计器
相关问题 更多 >
编程相关推荐