Re: [R-br] Teste Paralelismo no r

Geovane Assumindo que as regressões propostas nos gráficos (p8) fazem sentido e satisfazem os pressupostos do modelo de regressão linear poderia fazer o seguinte: library(car) escol='Superior/Pós graduação' model2=lm(Risco~0+(Idade*Quintis),subset(dados,Escolaridade3cat==escol)) summary(model2) # o teste para ‘Idade:Quintis5‘ na verdade está testando se os coeficientes angulares são iguais ( ou a diferença deles é zero) #para testar a diferença dos interceptos linearHypothesis(model2, "Quintis1 - Quintis5 =0") linearHypothesis(model2, "Idade:Quintis5 =0") Aqui mostrei os testes para uma escolaridade (escol), na verdade vc tem mais informação, teria que rodar um só modelo com todas as informações e fazer os testes nele: model1=lm(Risco~0+Idade*Escolaridade3cat*Quintis,data=dados) summary(model1) Boa sorte Robert From: Benilton Carvalho Sent: Thursday, March 06, 2014 2:15 PM To: r-br ; geovane barbosa Subject: Re: [R-br] Teste Paralelismo no r Geovane, Para realizar uma analise de dados, o estatistico precisa entender a (hi/e)storia dos dados. Por exemplo, sem saber do que se tratam os dados, e' dificil acreditar que a regressao linear simples que voce apresenta no seu exemplo seja adequada. A outra coisa: como/onde entra o fator "Quintis" existente no seu conjunto de dados? Enfim, supondo adequacidade em todos os sentidos (cabera' a vc verificar os pressupostos), se o interesse e' averiguar a identidade das regressoes nos tres grupos escolares (vou ignorar Quintis aqui, pois realmente nao sei o que e', nem como entra na sua analise), eu ajustaria os modelos com interacao (vis-a-vis Escolaridade) e compararia com o modelo simplificado (que contem apenas Idade): fit1 = lm(Risco ~ Idade*Escolaridade3cat, data=dados) fit0 = lm(Risco ~ Idade, data=dados) anova(fit0, fit1) benilton Em 6 de março de 2014 10:43, geovane barbosa <geovanecb@yahoo.com.br> escreveu: Olá pessoal , segue em anexo um exemplo reprodutível. Quero verificar se os interceptos proporcionados pelas retas de regressão bem como a inclinação são estatisticamente diferente e significativas, como faço isso? segue abaixo os codigos e banco num link. O teste do paralelismo consiste em verificar se as inclinações são estatisticamente significantes? dados<-read.csv("bancoteste.csv",sep=";",dec=",",header=T) dados$Escolaridade3cat<-factor(dados$Escolaridade3cat,label=c("Fundamental","Médio","Superior/Pós graduação"),levels=1:3) dados$Quintis<-factor(dados$Quintis) c <- ggplot(dados, aes(y=Risco, x=Idade, linetype=Quintis)) p8 <-c + stat_smooth(method=lm, se=FALSE,color="black",size=1.5)+ xlab("Idade")+ylab("Risco")+scale_colour_hue(name="")+theme_bw()+facet_grid(.~ Escolaridade3cat,space="free",scales="free_x") https://www.dropbox.com/s/ymyzh7p0xhmumda/bancoteste.csv Em Quinta-feira, 6 de Março de 2014 10:12, geovane barbosa <geovanecb@yahoo.com.br> escreveu: Olá Prof, segue em anexo um exemplo reprodutível. Como eu lhe disse anteriormente, queria verificar se os interceptos proporcionados pelas retas de regressão bem como a inclinação são estatisticamente diferente e significativas, como faço isso? segue abaixo os codigos e o banco de dados ficticio. O teste do paralelismo consiste em verificar se as inclinações são estatisticamente significantes? dados<-read.csv("bancoteste.csv",sep=";",dec=",",header=T) dados$Escolaridade3cat<-factor(dados$Escolaridade3cat,label=c("Fundamental","Médio","Superior/Pós graduação"),levels=1:3) dados$Quintis<-factor(dados$Quintis) c <- ggplot(dados, aes(y=Risco, x=Idade, linetype=Quintis)) p8 <-c + stat_smooth(method=lm, se=FALSE,color="black",size=1.5)+ xlab("Idade")+ylab("Risco")+scale_colour_hue(name="")+theme_bw()+facet_grid(.~ Escolaridade3cat,space="free",scales="free_x") Em Quarta-feira, 5 de Março de 2014 16:54, Benilton Carvalho <beniltoncarvalho@gmail.com> escreveu: A gente pede que exemplos reprodutíveis sejam dados para que possamos dar uma resposta definitiva sem a necessidade de idas e vindas. Como eu não entendi quais são os três cenários, fica a tarefa p vc: o caso a que vc se refere consiste em testar se as interações das variáveis de interesse com os referidos grupos são conjuntamente iguais a zero. b On Mar 5, 2014 11:00 AM, "geovane barbosa" <geovanecb@yahoo.com.br> wrote: Olá pessoal bom dia. Estou com uma questão que não consigo soluciona-la. Realize um regressão linear múltipla para comprar um grau de risco por idade e por escolaridade. Quero verificar se a inclinação proporcionada por essas retas é estatisticamente significativa, ou seja comparando a inclinação nos três cenários. bem como tambem seu intercepto. Como eu faço isso? Seria o test de Brant? Não consigo encontrar tal ferramenta no r. abraços _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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. _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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. _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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. -------------------------------------------------------------------------------- _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.
participantes (1)
-
Robert Iquiapaza