我有一张如下表:Type Actual Predicted
A 4 3
A 10 18
A 13 11
B 3 10
B 4 2
B 8 33
C 20 17
C 40 33
C 87 80
C 32 30
我想计算每种类型的R^2和RMSE。执行此操作的代码如下:
^{pr2}$
示例输出表(值为假设值):Type R^2 RMSE
A 0.66 4
B 1.00 6
C 0.03 1
上面的代码有效,并给出了我想要的输出。但是现在我想在类型级别的表中添加置信/预测间隔。我真的在网上搜遍了如何在没有运气的情况下做到这一点。在
概念性问题:如果我想要以95%的置信度捕获实际值的值范围,我是对实际列还是对预测列运行置信区间?在
下面是我想要的样本表:Type Conf_Int_90% Conf_Int_80%
A (21, 100) (5, 55)
B (10, 46) (3, 14)
C (1, 19) (12, 19)
我觉得置信区间代码是这样的:st.t.interval(0.95, len(a)-1, loc=np.mean(a), scale=st.sem(a)) BUT ...
我应该将哪些特定代码合并到现有代码中(如上所示),以便获得所需的表输出?在