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
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.