[問題] 請問有關R的程式

看板Statistics作者 (等的好辛苦)時間18年前 (2008/03/12 20:10), 編輯推噓2(201)
留言3則, 3人參與, 最新討論串3/3 (看更多)
k=array(99) k[1]=0.01 for(i in 2:99)k[i]=k[i]+0.01 a=read.table("1.txt") b=read.table("2.txt") a=as.matrix(a) b=as.matrix(b) c=array(101) d=array(101) for(i in 2:101){c[i]=(a[i]-a[i-1])/a[i-1];d[i]=(b[i]-b[i-1])/b[i-1];} c[1]=0 d[1]=0 e=array(101) for(i in 1:101)e[i]=(1/3)*c[i]+(2/3)*d[i] f=array(99) f[1]=100000 ccc=function(k,x,y,i) (1-k)*sum(x*y*k^(i-1)) for(i in 1:99){for(j in 1:101){f[i]=sum(e^2-ccc(k[i],a[j],b[j],j))}} 我是希望找出f的最小值在哪 可是寫到最後一行跑出來的f很怪@@" 想請問是哪裡出問題了 另外R不能按Enter..不能分行好痛苦@@ -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 59.116.44.140 ※ 編輯: wulingking 來自: 59.116.44.140 (03/12 20:12) ※ 編輯: wulingking 來自: 59.116.44.140 (03/12 20:20)

03/12 20:20, , 1F
可以用用看Tinn-R
03/12 20:20, 1F

03/13 01:36, , 2F
你有檢查過到第三行為止的 k 值嗎?建議你先搞懂指令吧...
03/13 01:36, 2F

03/13 01:49, , 3F
啊~@@ 謝謝~
03/13 01:49, 3F
文章代碼(AID): #17ryWq00 (Statistics)
文章代碼(AID): #17ryWq00 (Statistics)