Novembro 2016 - R-br

Heterocedasticidade e teste de média
by Wagner Bonat 04 Nov '16

04 Nov '16

Caros, Alguém postou esse conjunto de dados com problema de pressupostos, principalmente heterocedasticidade. Agora a pouco veio outro e-mail com um problema similar. Fiz um exemplo um pouco mais detalhado de como isso pode ser facilmente resolvido e mostrando o efeito disso no modelo. # Example 2 ------------------------------------------------------------ Fenois = c(337.311, 344.874, 342.353, 325.546, 333.950, 330.588, 328.067, 328.067, 318.824, 331.429, 333.950, 334.790, 336.471, 338.151, 342.353, 259.160, 252.437, 268.403, 265.882, 266.723, 287.731, 88.571, 88.571, 90.252, 41.513, 52.437, 49.076, 88.571, 88.571, 90.252, 64.202, 60.000, 61.681) Cor = factor(c(rep("ambar",6), rep("ambar_claro",3), rep("ambar",6), rep("ambar_claro",6),rep("branco",6), rep("extra_ambar_claro",3),rep("branco",3))) # Exploratory analysis boxplot(Fenois ~ Cor) tapply(Fenois, Cor, sd) dados <- data.frame(Fenois, Cor) dados$id <- 1 # Fitting fit1 <- mcglm(c(Fenois ~ Cor), list(mc_id(dados)), data = dados) fit2 <- mcglm(c(Fenois ~ Cor), list(mc_dglm(~ Cor, id = "id", data = dados)), covariance = "expm", data = dados) # Goodness-of-fit gof(fit1) gof(fit2) # Comparing estimates and standard errors coef(fit1, type = "beta", std.error = TRUE) coef(fit2, type = "beta", std.error = TRUE) O interessante é que a estimativa pontual é exatamente a mesma, porém olha a enorme diferença nos erros padrões dos betas. -- Wagner Hugo Bonat ---------------------------------------------------------------------------------------------- Department of Mathematics and Computer Science (IMADA) University of Southern Denmark (SDU) and Laboratório de Estatística e Geoinformação (LEG) Universidade Federal do Paraná (UFPR)

3 2

heterocedasticidade
by Luiz Leal 04 Nov '16

04 Nov '16

PrezadosNos casos em que a transformação box-cox não é capaz remover a heterocedasticidade há alguma outra solução? Eu utilizei a sugestão do Leg/ UFPR nos meus dados, mas não funcionou.(exemplo: <http://www.leg.ufpr.br/Rpira/Rpira/node17.html>) Segue em anexo o script e os dados. Desde já agradeço Luiz

2 2

Sorteios aleatórios de um cadastro (AAS ou amostra estratificada)
by Andre Oliveira 02 Nov '16

02 Nov '16

Bom dia, estou buscando informações de como realizar sorteios aleatórios de um cadastro (AAS ou amostra estratificada) com objetivo de inferir para média ou para a proporção. Alguma dica de pacote para tal trabalho? obg André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

4 4

mutate vs within
by Leonard Mendonça de Assis 02 Nov '16

02 Nov '16

Boa tarde a todos Surgiu uma dúvida conceitual aqui: Eu tenho um data.frame que precisa unir com vários outros, além de fazer várias transformações. Estou utilizando a seguinte sequencia para fazer tudo bhz.data <- read_csv2( "Belo Horizonte 20161027.csv", col_types=colunas) %>% filter(Filtro == 1 & Filtro.Votar==1) %>% select(-Filtro, -Filtro.Votar) %>% left_join(label.Sexo,by=c("Sexo"="value")) %>% select(-Sexo) %>% left_join(label.Faixa.Idade, by=c("Faixa.Idade"="value")) %>% select(-Faixa.Idade) %>% left_join(label.Renda, by=c("Renda"="value")) %>% select(-Renda) %>% left_join(label.classe, by=c("Renda.Classe"="value")) %>% select(-Renda.Classe) %>% left_join(label.Escolaridade, by=c("Escolaridade"="value")) %>% select(-Escolaridade) %>% left_join(label.Etinia, by=c('Etinia'="value")) %>% select(-Etinia) %>% left_join(label.Religião, by=c("Religião"="value")) %>% select(-Religião) %>% left_join(label.Trabalha, by=c("Trabalha"="value")) %>% select(-Trabalha) %>% mutate(Area.1 = ifelse(is.na(Area.1),0,Area.1), Area.2 = ifelse(is.na(Area.2),0,Area.2), Area.3 = ifelse(is.na(Area.3),0,Area.3)) %>% mutate(area.Educação = ifelse(Area.1 == 1,1,0) + ifelse(Area.2 == 1,1,0) + ifelse(Area.3 == 1,1,0), area.Habitação = ifelse(Area.1 == 2,1,0) + ifelse(Area.2 == 2,1,0) + ifelse(Area.3 == 2,1,0), area.Limpeza = ifelse(Area.1 == 3,1,0) + ifelse(Area.2 == 3,1,0) + ifelse(Area.3 == 3,1,0), area.Obras.Mobilidade = ifelse(Area.1 == 4,1,0) + ifelse(Area.2 == 4,1,0) + ifelse(Area.3 == 4,1,0), area.Obras.Pavimentação = ifelse(Area.1 == 5,1,0) + ifelse(Area.2 == 5,1,0) + ifelse(Area.3 == 5,1,0), area.Saneamento = ifelse(Area.1 == 6,1,0) + ifelse(Area.2 == 6,1,0) + ifelse(Area.3 == 6,1,0), area.Saúde = ifelse(Area.1 == 7,1,0) + ifelse(Area.2 == 7,1,0) + ifelse(Area.3 == 7,1,0), area.Transpoorte = ifelse(Area.1 == 8,1,0) + ifelse(Area.2 == 8,1,0) + ifelse(Area.3 == 8,1,0), area.Conservação = ifelse(Area.1 == 9,1,0) + ifelse(Area.2 == 9,1,0) + ifelse(Area.3 == 9,1,0), area.Acessibilidade = ifelse(Area.1 == 10,1,0) + ifelse(Area.2 == 10,1,0) + ifelse(Area.3 == 10,1,0), area.Combate.Mosquito = ifelse(Area.1 == 11,1,0) + ifelse(Area.2 == 11,1,0) + ifelse(Area.3 == 11,1,0), area.MeioAmbiente = ifelse(Area.1 == 12,1,0) + ifelse(Area.2 == 12,1,0) + ifelse(Area.3 == 12,1,0), area.Social = ifelse(Area.1 == 13,1,0) + ifelse(Area.2 == 13,1,0) + ifelse(Area.3 == 13,1,0), area.Esporte.Lazer = ifelse(Area.1 == 14,1,0) + ifelse(Area.2 == 14,1,0) + ifelse(Area.3 == 14,1,0), area.Geração.Empregos = ifelse(Area.1 == 15,1,0) + ifelse(Area.2 == 15,1,0) + ifelse(Area.3 == 15,1,0)) %>% select(-Area.1, -Area.2, -Area.3) A questão é a seguinte: No caso da variável área, que é recodificada de 3 para 15 variáveis, quando eu substituo esta sequencia enorme de mutate por um within, as variáveis ficam em uma ordem diferente da que eu preciso. A questão é: como deixar o código mais limpo nesta situação, visto que esta alternativa funciona, mas parece ser meio gambiarra? Leonard

3 3

Teste de independência dos resíduos
by Maurício Sangiogo 02 Nov '16

02 Nov '16

Bom dia listeiros! Estou iniciando no R, e gostaria de pedir ajuda a vocês sobre uma dúvida... Como eu faço no R para testar o pressuposto da anova de "independência dos resíduos"? Eu li algumas citações ao uso do Durbin-watson test para análise de resíduos em regressões, porém não sei se o mesmo se encaixaria (ou se ao menos está correto) para análise de experimentos com fatores qualitativos (em DBC ou DIC). Também vi algo em plotar um gráfico dos resíduos, mas não entendi o como. .. Se alguma boa alma estatística puder me ajudar, muito obrigado, se não puderem agradeço igualmente. Att .................. Eng. Agr. Maurício Sangiogo TEL:(55) 96822030

4 3

comparação de médias - problema com os pressupostos
by Maurício Lordêlo 02 Nov '16

02 Nov '16

Caros, Preciso fazer um teste de comparação de médias mas estou tendo problemas com os pressupostos. Tentei uma transformação Box Cox mas não consegui êxito. Alguma sugestão? Fenois = c(337.311, 344.874, 342.353, 325.546, 333.950, 330.588, 328.067, 328.067, 318.824, 331.429, 333.950, 334.790, 336.471, 338.151, 342.353, 259.160, 252.437, 268.403, 265.882, 266.723, 287.731, 88.571, 88.571, 90.252, 41.513, 52.437, 49.076, 88.571, 88.571, 90.252, 64.202, 60.000, 61.681) Cor = factor(c(rep("ambar",6),rep("ambar_claro",3),rep("ambar",6),rep("ambar_claro",6),rep("branco",6), rep("extra_ambar_claro",3),rep("branco",3))) Fenois; Cor boxplot(Fenois~Cor,ylab=expression(paste("Fenois"," ","(",mg," ",kg^-1,")")), xlab="cor",names=c("Âmbar", "Âmbar Claro", "Branco", "Extra Âmbar Claro")) tapply(Fenois, Cor, mean) tapply(Fenois, Cor, sd) tapply(Fenois, Cor, var) #Teste de homogeneidade de variâncias bartlett.test(Fenois ~ Cor) #Gráficos para avaliar a variável resposta qqnorm(Fenois) qqline(Fenois) hist(Fenois) boxplot(Fenois) #Transformação Box-Cox require(MASS) boxcox(Fenois~Cor, lambda = seq(-2, 2, 1/10)) boxcox(Fenois~Cor, lambda = seq(0.5, 1.5, 1/100)) #Optando por lambda = 0.6 Fenois2 = ((Fenois^0.6) - 1)/0.6 par(mfrow=c(1,2)) boxplot(Fenois~Cor,ylab=expression(paste("Fenois"," ","(",mg," ",kg^-1,")")), xlab="cor",names=c("Âmbar", "Âmbar Claro", "Branco", "Extra Âmbar Claro")) boxplot(Fenois2~Cor,ylab=expression(paste("Fenois"," ","(",mg," ",kg^-1,")")), xlab="cor",names=c("Âmbar", "Âmbar Claro", "Branco", "Extra Âmbar Claro")) bartlett.test(Fenois2 ~ Cor)

3 2

Sorteios aleatórios de um cadastro (AAS ou amostra estratificada)
by Andre Oliveira 02 Nov '16

02 Nov '16

Bom dia, estou buscando informações de como realizar sorteios aleatórios de um cadastro (AAS ou amostra estratificada) com objetivo de inferir para média ou para a proporção. Alguma dica de pacote para tal trabalho? obg André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

1 0

Barplot
by Emerson Cotta Bodevan 01 Nov '16

01 Nov '16

Prezados, bom dia. Estou tentando construir um barplot horizontal tal que as barras representem a frequencia de uso de dez medicamentos. Gostaria de colocar os nomes de todos os medicamentos a esquerda das barras. Não estou conseguindo colocar todos os nomes. Somente alguns deles aparecem. Alguma sugestão? Agradeço qualquer ajuda. Att., *Emerson*

4 4