powered by
The average weighted estimator and the unknown weighted estimator of the PQL in Poisson-GLMS through damped Gauss-Newton
pqlLogist(data,G,nk)
estimation value of betaW
estimation value of betaA
estimation of MSEW
estimation of MSEA
data is a highly correlated data set
G is the number of nodes
n1 is the length of each data subset
p<- 5;G<- 20;n<- 1000;nk=200 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) pqlLogist(data,G,nk)
Run the code above in your browser using DataLab