%%%%%%%%% Τα δεδομένα αναφέρονται σε ώρες επιβίωσης που μετρούν την αντοχή συγκεκριμένων %%%%%%%%%%%%%%πλοίων σε συνθήκες πίεσης


lifetimes<-c(1051,1337,1389,1921,1942,2322,3629,4006,4012,4063, 4921,5445,5620,5817,5905,5956,6068,6121,6473,7501,7886,8108,8546,8666,8831,9106,9711,9806,10205,10396, 10861,11026,11214,11362,11604,11608,11745,11762,11895,12044,13520,13670,14110,14496,15395,16179,17092,17568,17568)




x=seq(0,10, length=100)
plot(x, dweibull(x, shape=2, scale=1), type="l", ylab="Density")
lines(x, dweibull(x, shape=2, scale=2), lty=2)
title(main="Weibull Distribution, shape=2")
leg.names<-c("scale=1","scale=2")



loglikelihood <- function(data, theta, lambda=2)
 { logl<-sum((lambda-1)*log(data)+log(lambda)-lambda*log(theta)-(data/theta)^{lambda})
   return(logl)
 }


theta1 <- seq(7000, 15000, by=100)
 loglik=rep(NA, length(theta1))
 for (i in 1:length(theta1))
 {
 loglik[i]=loglikelihood(lifetimes, theta1[i])
 }