[R-br] Ajuda numa função para análise de resíduos em GEE

[Gilenio Borges Fernandes]

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.



-- 
Gilenio Borges Fernandes
Professor Associado IV
Universidade Federal da Bahia
Instituto de Matemática
Departamento de Estatística
Av. Adhemar de Barros, s/n – Ondina.
40.170-110 - Salvador - BA, Brasil
Tel.: (071)3283-6280/6336  Fax:  (071)3283-6276
URL: http://lattes.cnpq.br/6764860618464860
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131201/3f67cdce/attachment.html>