例子的话——
a=-0.618034
b=1.618034
l=function(x){
s=2*log(1+x^2)
t=2*log(1+(x-1)^2)
s-t
}
de=function(x){
4*x/(1+x^2)-4*(x-1)/(1+(x-1)^2)
}
x=seq(from=-4,to=4,by=0.01)
lstar=ifelse(x>a,ifelse(x<b,l(x)*2,de(b)*(x-b)+l(b)),de(a)*(x-a)+l(a))
y=lstar
plot(x,y,type='l',ylab= "log-likelihood ratio")
扫码加好友，拉您进群
全部版块
我的主页
收藏
