求助修改R脚本--小问题

1. e3=numeric(length(x0))
   
2. for(h in 1:(length(x0))){
   
3. e3[h]<- abs(e[h]-mean(e))
   
4. } #求偏差（残差离差的绝对值）
   
5. cat("偏差（残差离差的绝对值）:",'\n',e3,'\n')
   
6. fazhi<-0.6745*sx
   
7. cat("小误差阀值:",fazhi,'\n')

复制代码

这是灰色预测脚本的一小部分，执行的结果是：

我的目标很简单：
在屏幕显示：
p=
p是概率，是偏差小于阀值的概率。比如，上面的6个偏差全部小于阀值0.18，此时p=1。
谢谢！

p<-length(e3[e3<fazhi])/length(e3)
cat("p=",p,'\n')

x<-e3-fazhi
p<-length(x[x<0])/length(e3)
cat("p=",p,'\n')

syslon 发表于 2014-7-29 00:35
p

非常感谢！

URSIMON 发表于 2014-7-29 01:24
x

非常感谢您帮助！