作者wulingking (等的好辛苦)
看板Statistics
标题[问题] 请问有关R的程式
时间Wed Mar 12 20:10:27 2008
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)
1F：推 mangogogo:可以用用看Tinn-R 03/12 20:20
2F：推 KnightX:你有检查过到第三行为止的 k 值吗？建议你先搞懂指令吧... 03/13 01:36
3F：→ wulingking:啊~@@ 谢谢~ 03/13 01:49