程序报错但是找不到错误，求大神帮忙！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

generate_data=function(T,theta.true){
epsilon=rnorm(T,0,theta.true[1])
y=epsilon
y[1]=epsilon[1]
for(t in 2:T)
   {y[t]=theta.true[2]*epsilon[t-1]+epsilon[t]}
return(y)}

L=function(data,x1,x2){
dt=0
for(t in 1:T)
{dt=dt+log(x1*(1-x2^(2*t+2))/(1-x2^(2*t)))}
Y=rep(0,T)
Y[1]=data[1]
for(t in 2:T)
{Y[t]=data[t]-x2*(1-x2^(2*t-2))*Y[t-1]/(1-x2^(2*t))}
Sum=data[1]^2/((1+x2^2)*x1)
for(t in 2:T)
{Sum=Sum+Y[t]^2*(1-x2^(2*t))/(x1*(1-x2^(2*t+2)))}
l=-0.5*dt-0.5*Sum
return(l)
}

likeli=function(data,f,j){
deta=1e-5
g1=(L(f+deta,j)-L(f,j))/deta
g2=(L(f,j+deta)-L(f,j))/deta
g=c(g1,g2)
return(g)}

NR=function(data,epsilon){
A=diag(2)
theta.temp=c(0.2,0.6)
Gtemp=likeli(data,theta.temp[1],theta.temp[2])  
for(iter in 1:T){
theta=theta.temp+A%*%Gtemp
G=likeli(data,theta[1],theta[2])
A=A-(A%*%(G-Gtemp)%*%t(G-Gtemp)%*%A)/(t(G-Gtemp)%*%A%*%(G-Gtemp))-((theta-theta.temp)%*%t(theta-theta.temp))/((t(G-Gtemp)%*%(theta-theta.temp))
dif=max(abs(theta-theta.temp))
if(dif<=epsilon)break
theta.temp=theta
Gtemp=likeli(data,theta.temp[1],theta.temp[2])}
return(theta)
}
近段时间在学时间序列分析，老师要求采取改进的牛顿迭代方法来进行极大似然估计
是对Ma（1）过程，但是问题出现在最后一个NR函数，
A=A-(A%*%(G-Gtemp)%*%t(G-Gtemp)%*%A)/(t(G-Gtemp)%*%A%*%(G-Gtemp))-((theta-theta.temp)%*%t(theta-theta.temp))/((t(G-Gtemp)%*%(theta-theta.temp))
dif=max(abs(theta-theta.temp))
之后会报错显示错误: unexpected symbol in:
"A=A-(A%*%(G-Gtemp)%*%t(G-Gtemp)%*%A)/(t(G-Gtemp)%*%A%*%(G-Gtemp))-((theta-theta.temp)%*%t(theta-theta.temp))/((t(G-Gtemp)%*%(theta-theta.temp))
dif"

A的表达式是按照书上的来的，检查了好多遍应该没有输错，但一直找不到错误原因，求各位大神帮帮忙！！感激不尽！！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：function generate RETURN Theta Temp 程序

应该是这句" if(dif<=epsilon)break " 的问题，如果条件true，没有返回的值，会报错。
例如：
fun =function(x){
    if(max(x)<2) break
   return(max(x))
}
> a = c(0,1, 0.2, 0.3, 0.5)
> fun(a)
Error in fun(a) : no loop for break/next, jumping to top level
>
> a = c(0, 1, 3, 0.5)
> fun(a)
[1] 3
建议改成：
fun =function(x){
    if(max(x)<2){return(max(x));  break}
   return(max(x))
}