作者luckybp5 (低赔率也可以杀~很~大~喔)
看板Statistics
标题[问题] 可以帮我看哪里有错误吗??
时间Sat Jun 20 22:54:44 2009
E<-function(m11,m12,m21,m22,y11,y12,y21,y22,y1a,y1b,y2a,y2b)
{
i<-1
w<-y11+(m11/(m11+m12))*y1a+(m11/(m11+m21))*y1b
x<-y12+(m12/(m11+m12))*y1a+(m12/(m12+m22))*y2b
y<-y21+(m21/(m21+m22))*y2a+(m21/(m11+m21))*y1b
z<-y22+(m22/(m21+m22))*y2a+(m22/(m12+m22))*y2b
w_hat<-y11+(w/(w+x))*y1a+(w/(w+y))*y1b
x_hat<-y12+(x/(w+x))*y1a+(x/(x+z))*y2b
y_hat<-y21+(y/(y+z))*y2a+(y/(w+y))*y1b
z_hat<-y22+(z/(y+z))*y2a+(z/(x+z))*y2b
while(abs(w_hat-w)>=10^(-6))
while(abs(x_hat-x)>=10^(-6))
while(abs(y_hat-y)>=10^(-6))
while(abs(z_hat-z)>=10^(-6))
{
i<i+1
w=w_hat
x=x_hat
y=y_hat
z=z_hat
w<-y11+(m11/(m11+m12))*y1a+(m11/(m11+m21))*y1b
x<-y12+(m12/(m11+m12))*y1a+(m12/(m12+m22))*y2b
y<-y21+(m21/(m21+m22))*y2a+(m21/(m11+m21))*y1b
z<-y22+(m22/(m21+m22))*y2a+(m22/(m12+m22))*y2b
w_hat<-y11+(w/(w+x))*y1a+(w/(w+y))*y1b
x_hat<-y12+(x/(w+x))*y1a+(x/(x+z))*y2b
y_hat<-y21+(y/(y+z))*y2a+(y/(w+y))*y1b
z_hat<-y22+(z/(y+z))*y2a+(z/(x+z))*y2b
}
return(w_hat)
return(x_hat)
return(y_hat)
return(z_hat)
}
QQ<-E(100,50,75,75,100,50,75,75,30,28,60,60)
请问一下
我return那部份这样子写可以吗??????????
还是可以省略不用写阿????
请高手帮我解答一下~
谢啦
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 122.125.75.47
※ 编辑: luckybp5 来自: 122.125.75.47 (06/20 22:56)
1F：推 Jordan23:根本没宣告阿 06/20 23:07
※ 编辑: luckybp5 来自: 122.125.75.47 (06/20 23:40)
2F：→ clickhere:宣告?! 06/21 02:57
3F：→ clickhere:他有吧. 只是整个从while开始就不知道要做什了. 06/21 03:01
4F：推 Jordan23:文章有改过了 = = 06/21 03:03
5F：→ clickhere:return(c(w_hat, x_hat, y_hat, z_hat))? 06/21 03:13
6F：推 KnightX:return 其实可以省略不写, 合起来组成向量就好... 06/21 10:37
7F：→ luckybp5:感谢楼上大大们解释 谢谢 06/21 11:03