[問題] R 可以幫我看哪裡程式有錯誤嗎??

看板Statistics作者luckybp5 (低賠率也可以殺~很~大~喔)時間16年前 (2009/06/20 22:54)推噓3(3推 0噓 4→)

留言7則, 4人參與討論串1/1

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)

推

Jordan23

06/20 23:07, , 1^F

06/20 23:07, 1^F

※ 編輯: luckybp5 來自: 122.125.75.47 (06/20 23:40)

→

clickhere

06/21 02:57, , 2^F

06/21 02:57, 2^F

→

clickhere

06/21 03:01, , 3^F