search for: loglik1

Hi All, Does anyone know how to compare the means of several circular-data samplings? Any related websites, references and softwares? Thanks Xiaohua -- Xiaohua Dai, Dr. [[alternative HTML version deleted]]

...the other part doesn't. I don't see anything wrong in my code. Any hint would be greatly appreciated. Thank you. Z<-rbinom(20,1,0.5) #covariate time<-rgamma(20,shape=2,scale=0.5) #failure time del<-c(0,0,1,1,0,1,1,1,0,1,0,0,1,1,0,1,1,1,0,1) #event indicator ll<-function(p){ loglik1<-del*log(exp(Z*p[3])) #log of the first part of the likelihood function haz<-(time*exp(Z*p[3]))*p[1] ll2<-vector("numeric",length(1:20)) for (s in 1:20){ integrand<-function(x) x^(del[s]+1/p[2]-1)*exp(-x/p[2])*exp(-x*haz[s])/(gamma(1+1/p[2])*p[2]^(1/p[2])) ll2[s]<-log(in...

Hello, I have found a problem with using the glm function with a gamma family. I have a vector of data, assumed to be generated by a gamma distribution. The parameters of this gamma distribution are estimated in two ways (i) using the glm() function, (ii) "by hand", using the optim() function. I find that the -2*likelihood at the maximum found by (i) is substantially larger than that

...PCONVERGE(0.000001) SINGULAR(0.00000001) /MODEL /STEPWISE=PIN(.05) POUT(0.1) MINEFFECT(0) RULE(SINGLE) ENTRYMETHOD(LR) REMOVALMETHOD(LR) /INTERCEPT=INCLUDE /PRINT=FIT PARAMETER SUMMARY LRT CPS STEP MFI IC. ', file=syntaxfile, sep="", append=F) -> Loglik0: 135.02 -> Loglik1: 129.80 Thanks, S?ren