绘制威布尔的可能性

由 发布于

我想用一系列形状参数theta绘制大小为1000的Weibull样本的似然函数。我使用了标准化的weibull，因此比例lambda为1。但是输出是水平直线。

n<-1000
lik <- function(theta, x){
  K<- length(theta)
  n<- length(x)
  out<- rep(0,K)
  for(k in 1:K){
    out[k] <- prod(dweibull(x, shape= theta[k], scale=1))   
  }
  return(out)
}
theta<-seq(0.01, 10, by = 0.01)
x <- rweibull(n, shape= 0.5, scale= 1)
plot(theta, lik(theta, x), type="l", lwd=2)