[R-br] Variáveis contínuas explicando Variáveis categóricas
Daniel Tiezzi
dtiezzi em usp.br
Quinta Setembro 18 16:47:38 BRT 2014
d <- dados[ ,-c( )]
cor(d)
Daniel Tiezzi, MD, PhD
Professor Associado
Departamento de Ginecologia e Obstetrícia
Setor de Mastologia e Oncologia Ginecológica
Faculdade de Medicina de Ribeirão Preto - USP
Tel.: 16 3602-2488
e-mail: dtiezzi em fmrp.usp.br
On Sep 18, 2014, at 4:41 PM, Jefferson Ferreira-Ferreira <jecogeo em gmail.com> wrote:
>
> Daniel,
>
> Tenho um erro:
>
> Error in cor(dados$agua_quente + dados$detergentes + dados$agua_trat + :
> forneça conjuntamente 'x' e 'y' ou algo semelhante a uma matriz 'x'
>
> Em 18 de setembro de 2014 16:31, Daniel Tiezzi <dtiezzi em usp.br> escreveu:
> Faça uma matrix de correlação das suas variáveis.
>
> > cor(dados$agua_quente + dados$detergentes + dados$agua_trat + dados$Filtr_leite + dados$Higie_resf + dados$pre_dipping + dados$toalha + dados$pos_dipping + dados$vacuo)
>
> Qual o retorno que você tem?
>
>
>
>
>
>
> Daniel Tiezzi, MD, PhD
> Professor Associado
> Departamento de Ginecologia e Obstetrícia
> Setor de Mastologia e Oncologia Ginecológica
> Faculdade de Medicina de Ribeirão Preto - USP
> Tel.: 16 3602-2488
> e-mail: dtiezzi em fmrp.usp.br
>
>
>
>
>
> On Sep 18, 2014, at 4:24 PM, Jefferson Ferreira-Ferreira <jecogeo em gmail.com> wrote:
>
>>
>> Pois então...
>>
>> Daniel ou mais algum colega da lista, teria alguma dica? Será que ANOVA ou lm seriam recomendadas para meus dados?
>>
>> Em 18 de setembro de 2014 16:04, <dtiezzi em usp.br> escreveu:
>> As suas variáveis não são colineares?
>>
>> Por isso está saindo o NA no modelo de regressão
>>
>> Daniel
>>
>>
>>
>> De: "Jefferson Ferreira-Ferreira" <jecogeo em gmail.com>
>> Para: r-br em listas.c3sl.ufpr.br
>> Enviadas: Quinta-feira, 18 de Setembro de 2014 15:48:04
>>
>> Assunto: Re: [R-br] Variáveis contínuas explicando Variáveis categóricas
>>
>>
>> Daniel e demais,
>>
>> Aparentemente não resolve, pois recebo NAs para diversos slopes. (ver summary do modelo abaixo). E acho que isso não pode estar certo.
>> Segue aqui um exemplo dos meus dados. O que eu gostaria é de ver o quanto as variáveis categóricas explicam minha media_CBT e minha media_CCS.
>>
>> Alguma ideia???
>>
>>
>> agua_quente detergentes agua_trat Filtr_leite Higie_resf pre_dipping toalha pos_dipping vacuo pulsador Aliment_pos media_CBT media_CCS
>> 1 0 0 1 0 1 1 1 1 1 1 931 623.3333333
>> 1 0 0 1 1 0 1 1 1 1 1 163.6666667 259
>> 1 0 0 1 1 1 1 1 0 0 0 690.3333333 306.3333333
>> 1 0 0 1 1 1 1 1 0 0 1 690.3333333 306.3333333
>>
>>
>>
>> summary(m)
>>
>> Call:
>> lm(formula = dados$media_CBT ~ dados$agua_quente + dados$detergentes +
>> dados$agua_trat + dados$Filtr_leite + dados$Higie_resf +
>> dados$pre_dipping + dados$toalha + dados$pos_dipping + dados$vacuo)
>>
>> Residuals:
>> 1 2 3 4
>> 0.000e+00 -3.155e-30 4.019e-14 -4.019e-14
>>
>> Coefficients: (7 not defined because of singularities)
>> Estimate Std. Error t value Pr(>|t|)
>> (Intercept) 4.043e+02 8.988e-14 4.499e+15 <2e-16 ***
>> dados$agua_quente NA NA NA NA
>> dados$detergentes NA NA NA NA
>> dados$agua_trat NA NA NA NA
>> dados$Filtr_leite NA NA NA NA
>> dados$Higie_resf -2.407e+02 6.962e-14 -3.457e+15 <2e-16 ***
>> dados$pre_dipping 5.267e+02 6.962e-14 7.565e+15 <2e-16 ***
>> dados$toalha NA NA NA NA
>> dados$pos_dipping NA NA NA NA
>> dados$vacuo NA NA NA NA
>> ---
>> Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
>>
>> Residual standard error: 5.684e-14 on 1 degrees of freedom
>> Multiple R-squared: 1, Adjusted R-squared: 1
>> F-statistic: 4.872e+31 on 2 and 1 DF, p-value: < 2.2e-16
>>
>>
>>
>>
>> Em 18 de setembro de 2014 15:36, <dtiezzi em usp.br> escreveu:
>> Jefferson
>>
>> Montando o modele assim
>>
>> model <- lm(test$x ~ test$c1 + test$c2 + test$c3 + test$c4)
>>
>> não resolve seu problema?
>>
>> Daniel
>>
>>
>> De: "Jefferson Ferreira-Ferreira" <jecogeo em gmail.com>
>> Para: r-br em listas.c3sl.ufpr.br
>> Enviadas: Quinta-feira, 18 de Setembro de 2014 13:33:57
>> Assunto: Re: [R-br] Variáveis contínuas explicando Variáveis categóricas
>>
>>
>>
>> Professor daniel,
>> Obrigado pela resposta. Ainda não testei essa alternativa. Mas a desvantagem dessa abordagem seria o grande número de modelos a serem ajustados, visto que tenho 11 variáveis explicativas categóricas para duas variáveis resposta numéricas contínuas. Isso me faria ajustar 22 modelos que teriam como premissa a independencia entre as variáveis explicativas.
>>
>> O que quero dizer é: o quanto a combinação das variáveis categóricas c1 c2 c3 c4 c5 c6... etc (biárias =0 ou 1) explicam meu x (numérico contínuo).
>>
>> Será que existe um modo de eu ajustar dois modelos? Tipo com a variável resposta X em função de todas as varíaveis explicativas categóricas e outro modelo com a variável resposta Y em função de todas as variáveis explicativas categóricas?
>>
>>
>>
>> c1 c2 c3 c4 x y
>> 0 1 1 1 931 623
>> 0 1 0 0 163 259
>> 1 0 1 0 690 306
>> 1 0 1 0 690 306
>>
>>
>> --
>>
>> Jefferson Ferreira-Ferreira
>> Geógrafo – GEOPROCESSAMENTO IDSM | Coordenadoria de TI
>>
>> Jefferson.ferreira em mamiraua.org.br
>> Instituto de Desenvolvimento Sustentável Mamirauá
>> Ministério da Ciência, Tecnologia e Inovação
>> Telefone: +55 97 3343-9710
>> Google Maps - Mapas deste e-mail:
>> Exibir mapa ampliado
>>
>>
>>
>> Contatos particulares:
>> (55) 9615-0100
>>
>>
>>
>>
>> Em 18 de setembro de 2014 12:56, Daniel Tiezzi <dtiezzi em usp.br> escreveu:
>> Você precisa fazer uma regressão.
>>
>> Segue um modelo
>>
>> # Regression analyses, standardized
>> model1.z <- lm(scale(PE$endurance) ~ scale(PE$age))
>> summary(model1.z)
>> confint(model1.z)
>>
>> model2.z <- lm(scale(PE$endurance) ~ scale(PE$activeyears))
>> summary(model2.z)
>> confint(model2.z)
>>
>> model3.z <- lm(scale(PE$endurance) ~ scale(PE$age) + scale(PE$activeyears))
>> summary(model3.z)
>> confint(model3.z)
>>
>> # Conduct a model comparison NHST to compare the fit of model2.z to the fit of model3.z
>> anova(model2.z, model3.z)
>>
>>
>> Acho que seria assim
>>
>>
>> Daniel
>>
>>
>>
>>
>> Daniel Tiezzi, MD, PhD
>> Professor Associado
>> Departamento de Ginecologia e Obstetrícia
>> Setor de Mastologia e Oncologia Ginecológica
>> Faculdade de Medicina de Ribeirão Preto - USP
>> Tel.: 16 3602-2488
>> e-mail: dtiezzi em fmrp.usp.br
>>
>>
>>
>>
>>
>> On Sep 18, 2014, at 12:47 PM, Jefferson Ferreira-Ferreira <jecogeo em gmail.com> wrote:
>>
>>
>> Prezados;
>>
>> Uma dúvida básica de um iniciante em análises estatísticas.
>> Tenho uma série de variáveis categóricas binárias e duas variáveis contínuas. Eu gostaria de saber o quando minhas variáveis categóricas explicam minhas duas variáveis contínuas. Por exemplo:
>>
>> c1 c2 c3 c4 x y
>> 0 1 1 1 931 623
>> 0 1 0 0 163 259
>> 1 0 1 0 690 306
>> 1 0 1 0 690 306
>>
>>
>> A pergunta é: o quanto a combinação das variáveis c explicam os valores de x e y? Ou, qual a correlação entre as variáveis c e as variáveis x e y?
>>
>> Podem me dar alguma ideia de análises possíveis?
>> Obrigado.
>>
>>
>> --
>>
>> Jefferson Ferreira-Ferreira
>> Geógrafo – GEOPROCESSAMENTO IDSM | Coordenadoria de TI
>>
>> Jefferson.ferreira em mamiraua.org.br
>> Instituto de Desenvolvimento Sustentável Mamirauá
>> Ministério da Ciência, Tecnologia e Inovação
>> Telefone: +55 97 3343-9710
>> Google Maps - Mapas deste e-mail:
>>
>> Exibir mapa ampliado
>>
>>
>>
>> Contatos particulares:
>> (55) 9615-0100
>> _______________________________________________
>> R-br mailing list
>> R-br em 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 em 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 em 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 em 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 em 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.
>>
>>
>>
>> --
>>
>> Jefferson Ferreira-Ferreira
>> Geógrafo – GEOPROCESSAMENTO IDSM | Coordenadoria de TI
>>
>> Jefferson.ferreira em mamiraua.org.br
>> Instituto de Desenvolvimento Sustentável Mamirauá
>> Ministério da Ciência, Tecnologia e Inovação
>> Telefone: +55 97 3343-9710
>> Google Maps - Mapas deste e-mail:
>> Exibir mapa ampliado
>>
>>
>>
>> Contatos particulares:
>> (55) 9615-0100
>> _______________________________________________
>> R-br mailing list
>> R-br em 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 em 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.
>
>
>
> --
>
> Jefferson Ferreira-Ferreira
> Geógrafo – GEOPROCESSAMENTO IDSM | Coordenadoria de TI
>
> Jefferson.ferreira em mamiraua.org.br
> Instituto de Desenvolvimento Sustentável Mamirauá
> Ministério da Ciência, Tecnologia e Inovação
> Telefone: +55 97 3343-9710
> Google Maps - Mapas deste e-mail:
> Exibir mapa ampliado
>
>
>
> Contatos particulares:
> (55) 9615-0100
> _______________________________________________
> R-br mailing list
> R-br em 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.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140918/2243bb7f/attachment.html>
Mais detalhes sobre a lista de discussão R-br