<div dir="ltr"><p class="MsoNormal">Prezados(as)</p>

<p class="MsoNormal">Alguém poderia  identificar
o que está errado  na seguinte função:</p>

<p class="">diag_gee_binomial <- function(modelo, dados,
umaJanela=T, selGraficos = c(1,1,1,1), </p>

<p class="">                             
identifica=c(0,0,0,0))</p>

<p class="">{</p>

<p class="">  </p>

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

<p class="">  y <- modelo$y</p>

<p class="">  beta <-
coef(modelo)</p>

<p class="">  <span lang="EN-US">R <- modelo$work</span></p>

<p class=""><span lang="EN-US">  mi <- fitted(modelo)</span></p>

<p class=""><span lang="EN-US">  </span>individuo <- modelo$id</p>

<p class="">   repet <-
dim(modelo$work)[1]</p><p class=""></p>

<p class="">  ue <-
modelo$nobs/repet</p>

<p class="">   N <- nrow(X)</p><p class=""></p>

<p class="">  p <- ncol(X)</p>

<p class="">  </p>

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

<p class="">  A <-
diag(mi*(1-mi),N)</p>

<p class="MsoNormal">Etc...........</p>

<p class="MsoNormal">Na execução,</p>

<p class="">diag_gee_binomial(model.ar1, d1, umaJanela=T, selGraficos
= c(1,1,1,1), </p>

<p class="">                 
identifica=c(0,0,0,0))</p>

<p class=""> </p>

<p class="">Ocorre a seguinte mensagem de erro:</p>

<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-US" style="font-size:10pt;font-family:'Lucida Console';color:rgb(197,6,11);background-color:rgb(225,226,229)">Error in diag(mi *
(1 - mi), N) : </span></p>

<p class="MsoNormal" style="margin-bottom:0.0001pt"><span lang="EN-US" style="font-size:10pt;font-family:'Lucida Console';color:rgb(197,6,11);background-color:rgb(225,226,229)">  'nrow' or 'ncol' cannot be specified when 'x'
is a matrix</span><span lang="EN-US" style="font-size:10pt;font-family:'Lucida Console';color:black;background-color:rgb(225,226,229)"></span></p>

<p class=""><span lang="EN-US"> </span></p>

<p class="">O link para a o CMR da função completa é:</p>

<p class=""><a href="https://dl.dropboxusercontent.com/u/73337918/diag_gee_binomial.txt">https://dl.dropboxusercontent.com/u/73337918/diag_gee_binomial.txt</a></p>

<p class=""> </p>

<p class="">O banco de dados d1 pode ser baixado do link</p>

<p class=""> </p>

<p class=""><a href="https://dl.dropboxusercontent.com/u/73337918/d1.RData">https://dl.dropboxusercontent.com/u/73337918/d1.RData</a></p>

<p class=""> </p>

<p class="">Desde já agradeço</p>

<p class="">Att.</p>

<p class=""> </p><div><br></div>-- <br><div dir="ltr">Gilenio Borges Fernandes<div>Professor Associado IV<br>Universidade Federal da Bahia<br>Instituto de Matemática<br>Departamento de Estatística<br>Av. Adhemar de Barros, s/n – Ondina.<br>
40.170-110 - Salvador - BA, Brasil<br>Tel.: (071)3283-6280/6336  Fax:  (071)3283-6276<br>URL: <a href="http://lattes.cnpq.br/6764860618464860" target="_blank">http://lattes.cnpq.br/6764860618464860</a><br></div></div>
</div>