Desdobrar modelos de regresão em fatorial duplo (um fator quantitativo)
Estou com dificuldades em desdobrar a interação em uma regressão em fatorial. Segue um exemplo abaixo. Como obter um modelo de regressão de r em função de b para a1 e a2?????? a<-as.factor(c(rep(1,9),rep(2,9))) b<-c(rep(0,3),rep(5,3),rep(10,3), rep(0,3),rep(5,3),rep(10,3)) r<-c(1,1.2,1.3,2,2.3,2.25,4,3.9,3.8,5,4.8,4.7,3.5,3,3.4,1.7,1.5,1.2) b<-as.factor(b) m<-lm(r~a*b) anova(m) b<-as.numeric(b) regressão<-lm(r~b) regressão anova(regressão) Como desdobrar b dentro de a1 e a2????? Obrigado. Marcelo.
Marcelo, É só modificar a fórmula, m <- lm(r~a/b) anova(m) Para obter as estimativas da forma mais interpretável, ou seja, separado por nível de a, use m <- lm(r~-1+a/b) summary(m) À 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 ==========================================================================
Ok Walmes. Resolvido. Obrigado. --- Em qui, 5/5/11, Walmes Zeviani <walmeszeviani@gmail.com> escreveu: De: Walmes Zeviani <walmeszeviani@gmail.com> Assunto: Re: [R-br] Desdobrar modelos de regresão em fatorial duplo (um fator quantitativo) Para: r-br@listas.c3sl.ufpr.br Data: Quinta-feira, 5 de Maio de 2011, 20:14 Marcelo, É só modificar a fórmula, m <- lm(r~a/b) anova(m) Para obter as estimativas da forma mais interpretável, ou seja, separado por nível de a, use m <- lm(r~-1+a/b) summary(m) À 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 ========================================================================== -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Walmes, mais uma dúvida. É possível desdobrar a anova??? Semelhante ao resultado do pacote ExpDes?? Tentei usando contrast mais dá erro por usar um fator quantitativo. É possível? --- Em qui, 5/5/11, Walmes Zeviani <walmeszeviani@gmail.com> escreveu: De: Walmes Zeviani <walmeszeviani@gmail.com> Assunto: Re: [R-br] Desdobrar modelos de regresão em fatorial duplo (um fator quantitativo) Para: r-br@listas.c3sl.ufpr.br Data: Quinta-feira, 5 de Maio de 2011, 20:14 Marcelo, É só modificar a fórmula, m <- lm(r~a/b) anova(m) Para obter as estimativas da forma mais interpretável, ou seja, separado por nível de a, use m <- lm(r~-1+a/b) summary(m) À 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 ========================================================================== -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Bom dia pessoal, Fui rodar uma rotina que eu já havia realizado antes, mas dessa vez começou a aparecer o seguinte erro: model.eme<-aov(parasitismo~tratamento+tempo,data=liberacao) Erro em storage.mode(y) <- "double" : invalid to change the storage mode of a factor Além disso: Warning message: In model.response(mf, "numeric") : using type="numeric" with a factor response will be ignored Alguem faz alguma idéis o que seja isso, Obrigado, Alexandre dos Santos Ingenieur forestier, Msc. INRA- Biostatistique et Processus Spatiaux (BioSP) Domaine Saint-Paul Site Agroparc 84914 - Avignon - France Tél. : +33 (0)6 87 95 16 29
sim... sua variavel resposta 'parasitismo' eh um fator, o que certamente nao era da primeira vez. b 2011/6/12 Alexandre dos Santos <alexandresantosbr@yahoo.com.br>:
Bom dia pessoal,
Fui rodar uma rotina que eu já havia realizado antes, mas dessa vez começou a aparecer o seguinte erro:
model.eme<-aov(parasitismo~tratamento+tempo,data=liberacao)
Erro em storage.mode(y) <- "double" :
invalid to change the storage mode of a factor
Além disso: Warning message:
In model.response(mf, "numeric") :
using type="numeric" with a factor response will be ignored
Alguem faz alguma idéis o que seja isso,
Obrigado,
Alexandre dos Santos
Ingenieur forestier, Msc.
INRA- Biostatistique et Processus Spatiaux (BioSP)
Domaine Saint-Paul Site Agroparc 84914 - Avignon - France Tél. : +33 (0)6 87 95 16 29
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- Successful people ask better questions, and as a result, they get better answers. (Tony Robbins)
Obrigado Benilton, Mas o problema é que se eu tentava fazer as.numeric(), a função transformava os zeros do vetor em 1. Abraço, -----Original Message----- From: r-br-bounces@listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] On Behalf Of Benilton Carvalho Sent: domingo, 12 de junho de 2011 11:17 To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Erro storage(y) sim... sua variavel resposta 'parasitismo' eh um fator, o que certamente nao era da primeira vez. b 2011/6/12 Alexandre dos Santos <alexandresantosbr@yahoo.com.br>:
Bom dia pessoal,
Fui rodar uma rotina que eu já havia realizado antes, mas dessa vez começou a aparecer o seguinte erro:
model.eme<-aov(parasitismo~tratamento+tempo,data=liberacao)
Erro em storage.mode(y) <- "double" :
invalid to change the storage mode of a factor
Além disso: Warning message:
In model.response(mf, "numeric") :
using type="numeric" with a factor response will be ignored
Alguem faz alguma idéis o que seja isso,
Obrigado,
Alexandre dos Santos
Ingenieur forestier, Msc.
INRA- Biostatistique et Processus Spatiaux (BioSP)
Domaine Saint-Paul Site Agroparc 84914 - Avignon - France Tél. : +33 (0)6 87 95 16 29
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- Successful people ask better questions, and as a result, they get better answers. (Tony Robbins) _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
vc precisa ler com cuidado a ajuda do 'factor', pois ela explica como transformar (aproximadamente) um fator para numerico.... b 2011/6/12 Alexandre dos Santos <alexandresantosbr@yahoo.com.br>:
Obrigado Benilton,
Mas o problema é que se eu tentava fazer as.numeric(), a função transformava os zeros do vetor em 1. Abraço,
-----Original Message----- From: r-br-bounces@listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] On Behalf Of Benilton Carvalho Sent: domingo, 12 de junho de 2011 11:17 To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Erro storage(y)
sim... sua variavel resposta 'parasitismo' eh um fator, o que certamente nao era da primeira vez.
b
2011/6/12 Alexandre dos Santos <alexandresantosbr@yahoo.com.br>:
Bom dia pessoal,
Fui rodar uma rotina que eu já havia realizado antes, mas dessa vez começou a aparecer o seguinte erro:
model.eme<-aov(parasitismo~tratamento+tempo,data=liberacao)
Erro em storage.mode(y) <- "double" :
invalid to change the storage mode of a factor
Além disso: Warning message:
In model.response(mf, "numeric") :
using type="numeric" with a factor response will be ignored
Alguem faz alguma idéis o que seja isso,
Obrigado,
Alexandre dos Santos
Ingenieur forestier, Msc.
INRA- Biostatistique et Processus Spatiaux (BioSP)
Domaine Saint-Paul Site Agroparc 84914 - Avignon - France Tél. : +33 (0)6 87 95 16 29
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- Successful people ask better questions, and as a result, they get better answers. (Tony Robbins) _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- Successful people ask better questions, and as a result, they get better answers. (Tony Robbins)
Obrigado Benilton, Funcionou com as.numeric(as.character(parasitismo)), -----Original Message----- From: r-br-bounces@listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] On Behalf Of Benilton Carvalho Sent: domingo, 12 de junho de 2011 14:20 To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Erro storage(y) vc precisa ler com cuidado a ajuda do 'factor', pois ela explica como transformar (aproximadamente) um fator para numerico.... b 2011/6/12 Alexandre dos Santos <alexandresantosbr@yahoo.com.br>:
Obrigado Benilton,
Mas o problema é que se eu tentava fazer as.numeric(), a função transformava os zeros do vetor em 1. Abraço,
-----Original Message----- From: r-br-bounces@listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] On Behalf Of Benilton Carvalho Sent: domingo, 12 de junho de 2011 11:17 To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Erro storage(y)
sim... sua variavel resposta 'parasitismo' eh um fator, o que certamente nao era da primeira vez.
b
2011/6/12 Alexandre dos Santos <alexandresantosbr@yahoo.com.br>:
Bom dia pessoal,
Fui rodar uma rotina que eu já havia realizado antes, mas dessa vez começou a aparecer o seguinte erro:
model.eme<-aov(parasitismo~tratamento+tempo,data=liberacao)
Erro em storage.mode(y) <- "double" :
invalid to change the storage mode of a factor
Além disso: Warning message:
In model.response(mf, "numeric") :
using type="numeric" with a factor response will be ignored
Alguem faz alguma idéis o que seja isso,
Obrigado,
Alexandre dos Santos
Ingenieur forestier, Msc.
INRA- Biostatistique et Processus Spatiaux (BioSP)
Domaine Saint-Paul Site Agroparc 84914 - Avignon - France Tél. : +33 (0)6 87 95 16 29
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- Successful people ask better questions, and as a result, they get better answers. (Tony Robbins) _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- Successful people ask better questions, and as a result, they get better answers. (Tony Robbins) _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
;-) qdo vc for ler conjuntos de dados relativamente pequenos (nesse caso e' comum o R converter para fator), vc pode se beneficiar com o uso de stringsAsFactors=FALSE no read.table() e variantes. Olhe tambem os argumentos 'as.is' e 'colClasses'. b 2011/6/12 Alexandre dos Santos <alexandresantosbr@yahoo.com.br>:
Obrigado Benilton,
Funcionou com as.numeric(as.character(parasitismo)),
-----Original Message----- From: r-br-bounces@listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] On Behalf Of Benilton Carvalho Sent: domingo, 12 de junho de 2011 14:20 To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Erro storage(y)
vc precisa ler com cuidado a ajuda do 'factor', pois ela explica como transformar (aproximadamente) um fator para numerico....
b
2011/6/12 Alexandre dos Santos <alexandresantosbr@yahoo.com.br>:
Obrigado Benilton,
Mas o problema é que se eu tentava fazer as.numeric(), a função transformava os zeros do vetor em 1. Abraço,
-----Original Message----- From: r-br-bounces@listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] On Behalf Of Benilton Carvalho Sent: domingo, 12 de junho de 2011 11:17 To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Erro storage(y)
sim... sua variavel resposta 'parasitismo' eh um fator, o que certamente nao era da primeira vez.
b
2011/6/12 Alexandre dos Santos <alexandresantosbr@yahoo.com.br>:
Bom dia pessoal,
Fui rodar uma rotina que eu já havia realizado antes, mas dessa vez começou a aparecer o seguinte erro:
model.eme<-aov(parasitismo~tratamento+tempo,data=liberacao)
Erro em storage.mode(y) <- "double" :
invalid to change the storage mode of a factor
Além disso: Warning message:
In model.response(mf, "numeric") :
using type="numeric" with a factor response will be ignored
Alguem faz alguma idéis o que seja isso,
Obrigado,
Alexandre dos Santos
Ingenieur forestier, Msc.
INRA- Biostatistique et Processus Spatiaux (BioSP)
Domaine Saint-Paul Site Agroparc 84914 - Avignon - France Tél. : +33 (0)6 87 95 16 29
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- Successful people ask better questions, and as a result, they get better answers. (Tony Robbins) _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- Successful people ask better questions, and as a result, they get better answers. (Tony Robbins) _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- Successful people ask better questions, and as a result, they get better answers. (Tony Robbins)
participantes (4)
-
Alexandre dos Santos -
Benilton Carvalho -
Marcelo Cardoso mello -
Walmes Zeviani