GOLL_W<-function(tempo,alpha,theta,a,b){ gt5=(a*tempo^(a-1)/b^a)*(exp(-(tempo/b)^a)) Gt5=1-exp(-(tempo/b)^a) ft5=(alpha*theta*gt5*(Gt5^(alpha*theta-1))*((1-Gt5^(theta))^(alpha-1)))/((Gt5^(alpha*theta)+(1-Gt5^theta)^alpha)^2) Ft5=(Gt5^(alpha*theta))/(Gt5^(alpha*theta)+ (1-Gt5^(theta))^alpha) St5=1-Ft5 ht5=ft5/St5 return(ht5) } tempo=c(seq(0,10,0.01)) v_lty <-c('longdash ', 'dotdash', 'dotted', 'dashed', 'solid')