Fernando, 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, tr <- gl(3,4) y <- rbinom(length(x), size=50, prob=0.5) m0 <- glm(cbind(y, 50-y)~tr, family=binomial) summary(m0) # solução model.matrix(m0) # matriz de delinamento X contrasts(tr) # tipo de restrição 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 http://www.leg.ufpr.br/doku.php/software:rbr À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== 2011/9/10 Fernando Antonio de souza <nandodesouza@gmail.com>

Caro Prof. Walmer,

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?

Grato pela atenção

Fernando Zootecnista, Mestre em Nutrição de Ruminantes