所以这是我的代码:
solfun<-function(par,data)
{
data1<-data[,1];
data2<-data[,2];
n<-length(data1)
LL<-1:n;
sigs<-1:(n+1);
z<-1:n;
xd<-1:n;
sigs[1]<-sqrt(mean(data1^2)) ;
for(i in 2:(n+1)) {
z[i-1]=data1[i-1]/sigs[i-1];
xd[i-1]=(data2[i-1]-0.43-log(sigs[i-1]))/0.29
sigs[i]=exp(par[1]/2+par[3]*log(sigs[i-1])+par[4]/2*xd[i-1]+par[2]/2*z[i-1])
} ;
for(i in 1:n) {
LL[i] <- (-0.5*log(sqrt(2*pi)*0.29)-(data2[i]-0.43-log(sigs[i]))^2/(2*0.29^2))
};
return(sum(-LL)) ;
}
solfun1<-function(par,data)
{
data1<-data[,1];
data2<-data[,2];
n<-length(data1)
LL<-1:n;
sigs<-1:(n+1);
z<-1:n;
xd<-1:n;
sigs[1]<-sqrt(mean(data1^2)) ;
for(i in 2:(n+1)) {
z[i-1]=data1[i-1]/sigs[i-1];
xd[i-1]=(data2[i-1]-0.43-log(sigs[i-1]))/0.29
sigs[i]=exp(par[1]/2+par[3]*log(sigs[i-1])+par[4]/2*xd[i-1]+par[2]/2*z[i-1])
} ;
for(i in 1:n) {
LL[i] <- (-0.5*log(sqrt(2*pi)*0.29)-(data2[i]-0.43-log(sigs[i]))^2/(2*0.29^2))
};
if (is.finite(sum(-LL))){
return(sum(-LL))
}
return(1000 + sum(par ^ 2)) ;
}
fit1<-optim(par<-c(0.19413280,0.02130018, 0.97327246, 0.13672007),solfun,data=data)
fit2<-optim(par<-c(0.19413280,0.02130018, 0.97327246, 0.13672007),solfun1,data=data,
method = "L-BFGS-B",lower=c(-10,-10, -1,-10),upper=c(10,10, 1,10))
根据我在互联网上收集的内容,这应该显示图像&#39; new.bmp&#39;在一个名为&#39; FRAME&#39;的窗口中。但所有这一切都是 this!。任何人都可以帮助我吗?
答案 0 :(得分:1)
这是一个愚蠢的安装问题。
我全心全意地道歉。