Re: [R-br] Heterocedasticidade e teste de média

Tentei rodar, e ele não encontra as funções, tem alguma bibliotéca para isto?      fit1 <- mcglm(c(Fenois ~ Cor), list(mc_id(dados)), data = ""> Erro: não foi possível encontrar a função "mcglm" > fit2 <- mcglm(c(Fenois ~ Cor), list(mc_dglm(~ Cor, id = "id", data = ""> +               covariance =

Re: [R-br] Heterocedasticidade e teste de média

Olá Wagner, Obrigado pela ilustração do exemplo. Não consegui reproduzir a linha do fit2, especificamente o comando mc_dglm. Você sugeriria algum texto para maior aprofundamento? *Paulo Dick* Estatístico / Epidemiologia em Saúde Pública Tel.: (55 21) 99591-2716 Em 4 de novembro de 2016 16:46,

[R-br] Heterocedasticidade e teste de média

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

Re: [R-br] heterocedasticidade

Complementando minha resposta. dados <- read.table("dados.csv", header = TRUE, sep = ";", dec = ",") plot(dados$y ~ dados$x) # Veja que o A33 nao tem variabilidade vou remover dados <- dados[which(dados$x != "A33"),] dados$x <- droplevels(dados$x) tapply(dados$y, dados$x, sd) # Aparentemente o

Re: [R-br] heterocedasticidade

Pq simplesmente não tenta modelar a variancia? Coloca um modelo de regressão pro parâmetro de variancia. Ou tenta um distribuição que seja naturalmente heterocedástica tipo Gama ou mais geral Tweedie. Em 4 de novembro de 2016 18:06, Luiz Leal via R-br escreveu: >

[R-br] heterocedasticidade

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:  )  Segue em anexo o script e os dados. Desde já