Displaying 1 result from an estimated 1 matches for "llne".
Did you mean:
line
2008 Apr 08
1
Weibull maximum likelihood estimates for censored data
...g((a^b)*b*(D[2:(NT-1),1]^(b-1))*exp(-((a*D[2:(NT-1),1])^b))))
fd1 = (a^b)*b*(D[1,1]^(b-1))*exp(-(a*D[1,1])^b)
fdn = (a^b)*b*(D[NT,1]^(b-1))*exp(-(a*D[NT,1])^b)
S1 = exp(-(a*D[1,1])^b)
SN = exp(-(a*D[NT,1])^b)
-(D[1,2]*log(S1)+(1-D[1,2])*log(fd1)+ f + D[NT,2]*log(SN)+
(1-D[NT,2])*log(fdn))
}
lLne <- function(A){
D = CDMatrix
NT = nrow(D)
b=1
f = sum(log(A*b*(D[2:(NT-1),1]^(b-1))*exp(-(A^(1/b)*D[2:(NT-1),1])^b)))
fd1 = A*b*(D[1,1]^(b-1))*exp(-(A^(1/b)*D[1,1])^b)
fdn = A*b*(D[NT,1]^(b-1))*exp(-(A^(1/b)*D[NT,1])^b)
S1 = exp(-(A^(1/b)*D[1,1])^b)
SN = exp(-(A^(1/b)*D[NT,1])^b)
lLw...