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: [email protected] twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== 2011/9/10 Fernando Antonio de souza <[email protected]> > 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 >
_______________________________________________ R-br mailing list [email protected] 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.
