Prezados(as)

Alguém poderia  identificar o que está errado  na seguinte função:

diag_gee_binomial <- function(modelo, dados, umaJanela=T, selGraficos = c(1,1,1,1),

                              identifica=c(0,0,0,0))

{

  X <- model.matrix(as.formula(paste("~ ", modelo$call$formula[3])), dados)

  y <- modelo$y

  beta <- coef(modelo)

  R <- modelo$work

  mi <- fitted(modelo)

  individuo <- modelo$id

   repet <- dim(modelo$work)[1]

  ue <- modelo$nobs/repet

   N <- nrow(X)

  p <- ncol(X)

  # Matriz C <- A * Delta                #Ligação canonica -> Delta=Identidade

  A <- diag(mi*(1-mi),N)

Etc...........

Na execução,

diag_gee_binomial(model.ar1, d1, umaJanela=T, selGraficos = c(1,1,1,1),

                  identifica=c(0,0,0,0))

Ocorre a seguinte mensagem de erro:

Error in diag(mi * (1 - mi), N) :

  'nrow' or 'ncol' cannot be specified when 'x' is a matrix

O link para a o CMR da função completa é:

https://dl.dropboxusercontent.com/u/73337918/diag_gee_binomial.txt

O banco de dados d1 pode ser baixado do link

https://dl.dropboxusercontent.com/u/73337918/d1.RData

Desde já agradeço

Att.