<font face="trebuchet ms,sans-serif">Fernando,<br><br>Ao declarar a fórmula do modelo na função glm(), a matriz de delineamento (X) é construída de acordo com uma restrição paramétrica (por padrão zerar o efeito associado ao primeiro nível do fator) que torna não singular a matriz X'X. Você não precisa se preocupar com isso. É só declarar a fórmula do modelo na glm(), conhecer a restrição paramétrica usada que fica fácil interpretar a solução pedindo um summary(). Veja esse exemplo mediocre mas que ilustra o que falei,<br>
<br><span style="font-family: courier new,monospace;">tr <- gl(3,4)</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">y <- rbinom(length(x), size=50, prob=0.5)</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">m0 <- glm(cbind(y, 50-y)~tr, family=binomial)</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">summary(m0)      # solução</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">model.matrix(m0) # matriz de delinamento X</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">contrasts(tr)    # tipo de restrição</span><br style="font-family: courier new,monospace;">
<br>Aconselho que envie suas mensagens para a nossa lista de discussão sobre R. Mais pessoas podem contribuir com você, com respostas objetivas e rápidas. Informações de cadastro, envio e etiqueta disponíveis em<br><br><a href="http://www.leg.ufpr.br/doku.php/software:rbr">http://www.leg.ufpr.br/doku.php/software:rbr</a><br>
<br>À disposição.<br>Walmes.<br><br clear="all"></font><span style="font-family: trebuchet ms,sans-serif;">==========================================================================</span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">Walmes Marques Zeviani</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">Departamento de Estatística - Universidade Federal do Paraná</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">fone: (+55) 41 3361 3573</span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">VoIP: (3361 3600) 1053 1173</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">e-mail: <a href="mailto:walmes@ufpr.br" target="_blank">walmes@ufpr.br</a></span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">twitter: @walmeszeviani</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">homepage: <a href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">linux user number: 531218</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">==========================================================================</span><br>

<br><br><div class="gmail_quote">2011/9/10 Fernando Antonio de souza <span dir="ltr"><<a href="mailto:nandodesouza@gmail.com">nandodesouza@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Caro Prof. Walmer,<br><br>Eu preciso fazer uma análise de regressão linear com efeitos fixos (Y~X+Z+X:Z), porém a matriz  X é singular e eu necessito fazer a inversa generalizada para que eu possa encontrar uma solução. Estou utilizando a função glm do software R. A pouco tempo descobri que no pacote (MASS) há uma função chamada GINV o qual gera a inversa generalizada porém não sei como colocá-la dentro do modelo para obter a solução. Você pode me orientar, como fazer isto?<br>

<br>Grato pela atenção<br><br>Fernando<br>Zootecnista, Mestre em Nutrição de Ruminantes<br>
</blockquote></div><br>