powered by
The average weighted estimator and the unknown weighted estimator of the PQL in Poisson-GLMS through damped Gauss-Newton
pqlBpoisson2(data,G,nk)
betaBA, betaBW, MSEA, MSEW
is a design matrix with uniform distribution and the response vector
is the number of subsets.
is the size of subsets.
p<- 5;G<- 20;n<- 1000;nk=50 X<- matrix(runif(n * p, 0, 0.5), ncol = p) beta =matrix(runif(p, 0, 1),nrow=p) L=X%*%beta y<- rpois(n, exp(L)) data=cbind(y,X) pqlBpoisson2(data,G,nk)
Run the code above in your browser using DataLab