Learn R Programming

pql (version 0.1.0)

pqlBpoisson1: The weight Gauss-Newton estimators of the PQL in Poisson-GLMS

Description

The average weighted estimator and the unknown weighted estimator of the PQL in Poisson-GLMS through damped Gauss-Newton

Usage

pqlBpoisson1(data,G,nk)

Value

betaBA, betaBW, MSEA, MSEW

Arguments

data

is a design matrix with uniform distribution and the response vector

G

is the number of subsets.

nk

is the size of subsets

Examples

Run this code
G <- 20;n=1000;p=5; nk=50
X<- matrix(runif(1000* 5, 0, 0.5), ncol = 5)
beta =matrix(runif(p, 0, 1),nrow=p)
L=X%*%beta
y<- rpois(1000, exp(L))
data=cbind(y,X)
pqlBpoisson1(data,G,nk)

Run the code above in your browser using DataLab