search for: insem

...- p[3]*log(p[3])-log(gamma(p[3]))+sum(lnhaz)+ log(gamma(di+p[3]))-(p[3]+di)*log(p[3]+cumhaz) -2*lik } di <-sum(stat) initial<-c(500000,2.8,.5,1.8) # initial 1/exp(constant) t<-nlm(likelihood.weibul,initial,print.level=1,hessian=T) Best Regards, Denis Lalountas University of Patras insem<-read.table("Melanoma.txt",header=T) id <- insem$id timeto<-insem$time stat<-insem$status x1<-insem$ulcer x2 <-insem$thickness x2 <- log(x2) id time status sex age year thickness ulcer 1 1 10 0 1 76 1972 6.76 1 2 2 30 0 1 56 1...