帮我看看这个R程序的错误

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

在书上看到的程序，结果运行且有好多错误，望高手帮我解答一下，问题出在什么地方。谢谢。

exam=function(Y,ep=1e-5){

  phi=function(x)1/sqrt(2*pi)*exp(-x^2/2)

  p=0.6;mu=3.5;n=length(Y)

  k=1

  repeat{

     k=k+1

p[k]=1/n*p[k-1]*sum(phi(Y-mu[k-1])/p[k-1]*phi(Y-mu[k-1])+(1-p[k-1])*phi(Y)))

mu[k]=sum(phi(Y-mu[k-1])*Y/(p[k-1]*phi(Y-mu[k-1])

       +(1-p[k-1])*phi(Y)))/sum(phi(Y-mu[k-1])/(p[k-1]*phi(Y-mu[k-1])

       +(1-p[k-1])*phi(Y)))

if(abs(p[k]-p[k-1])<=ep&abs(mu[k]-mu[k-1])<=ep)

    break

         }

list(p=p,mu=mu,iter=k-1)

}

错误提示如下：

p[k]=1/n*p[k-1]*sum(phi(Y-mu[k-1])/p[k-1]*phi(Y-mu[k-1])+(1-p[k-1])*phi(Y)))"

> mu[k]=sum(phi(Y-mu[k-1])*Y/(p[k-1]*phi(Y-mu[k-1])

+        +(1-p[k-1])*phi(Y)))/sum(phi(Y-mu[k-1])/(p[k-1]*phi(Y-mu[k-1])

+        +(1-p[k-1])*phi(Y)))

错误: 没有"phi"这个函数

> if(abs(p[k]-p[k-1])<=ep&abs(mu[k]-mu[k-1])<=ep)

+     break

错误: 找不到对象'ep'

>          }

错误: 意外的'}'在"         }"里

> list(p=p,mu=mu,iter=k-1)

$p

[1] 0.6

$mu

[1] 3.5

$iter

[1] 0

> }

错误: 意外的'}'在"}"里

>

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：R程序 function length repeat leng 程序

没有人能帮你的。除非你指出是哪一本书，第几页，可以帮你看看。

错误很简单，你的p[k]那一行(第七行)最后多了一个括号

qoiqpwqr 发表于 2011-12-2 22:16
错误很简单，你的p[k]那一行(第七行)最后多了一个括号

谢谢谢谢，现在成功了

R程序的结果在哪显示呢