推 mangogogo:可以用用看Tinn-R 03/12 20:20
推 KnightX:你有檢查過到第三行為止的 k 值嗎？建議你先搞懂指令吧... 03/13 01:36
→ wulingking:啊~@@ 謝謝~ 03/13 01:49
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)