作者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