[R-br] ajuda soft R

[Walmes Zeviani]

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 em ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================


2011/9/10 Fernando Antonio de souza <nandodesouza em 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
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110910/8b3297fa/attachment.html>