search for: llne

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...