[R-br] Duvida Regressão

[Paulo Justiniano]

Se a sua resposta é continua, ao inves de cair neste buuaco negro de 
testes de aderencia
sugoiri que voce examene alguma tranformacaod e variavel
como 0or exmeplo MASS:boxcox()

a funcao vai achar o valor do parametro (verossmihanca) que faz com que os 
dados mais se aproximem da normal.
A ventagf34m é que voce pode declarar o modelo
sem necessidadse de fazer em 2 passos (estimar residuos depois testar 
nomalidade)

algo como

boxcox(Y ~bloco + tratamento, ...


de acordo com seu modelo.

Melhor que isto ainda ;e encotnrar um modelo adequado para sua variavel
(Gamma, beta, etc) e usar a analise com esta distribuicao





On Sun, 29 Apr 2012, Sérgio Henrique almeida da silva ju wrote:

> Olá
>  
> Estou trabalhando com a taxa de uma doença e gostaria de rodar um modelo de regressão.
>  
> Meus dados não apresentaram distribuição normal e por teste de aderência deu uma distribuição lognomal
>  
> -> fitdistr(tx_esqui, "lognormal")
>      meanlog       sdlog
>   6.25725753   1.90754231
>  (0.11630491) (0.08223999)
>  
> ->ks.test(tx_esqui, "plnorm", sdlog=1.90754231, meanlog=6.25725753,alternative=c("two.sided"))
>         One-sample Kolmogorov-Smirnov test
> data:  tx_esqui
> D = 0.0701, p-value = 0.1419
> alternative hypothesis: two-sided
>  
> Tem como eu rodar um modelo de regressão com essa variável? Como fazê-lo? Qual a melhor forma de trabalhar com essa
> variável?
>  
> Abraços
>  
>  
> Segue a minha variável dependente:
>  
> structure(list(tx_esqui = c(3L, 3L, 4L, 4L, 4L, 4L, 8L, 8L, 8L,
> 9L, 10L, 10L, 19L, 22L, 22L, 23L, 23L, 25L, 25L, 26L, 27L, 28L,
> 29L, 30L, 31L, 32L, 34L, 36L, 39L, 41L, 42L, 45L, 46L, 49L, 49L,
> 56L, 57L, 58L, 58L, 59L, 61L, 62L, 62L, 68L, 70L, 74L, 76L, 80L,
> 82L, 82L, 83L, 83L, 86L, 101L, 102L, 107L, 107L, 112L, 113L,
> 118L, 120L, 121L, 123L, 126L, 135L, 138L, 140L, 153L, 156L, 158L,
> 179L, 191L, 200L, 201L, 204L, 207L, 209L, 212L, 215L, 223L, 227L,
> 228L, 233L, 247L, 248L, 249L, 255L, 257L, 258L, 258L, 259L, 264L,
> 274L, 277L, 279L, 287L, 297L, 309L, 310L, 313L, 314L, 326L, 328L,
> 329L, 331L, 332L, 332L, 340L, 354L, 358L, 365L, 372L, 379L, 383L,
> 395L, 396L, 418L, 423L, 447L, 449L, 451L, 453L, 493L, 532L, 535L,
> 553L, 554L, 565L, 573L, 581L, 584L, 599L, 603L, 604L, 611L, 672L,
> 703L, 708L, 727L, 756L, 799L, 809L, 812L, 815L, 816L, 862L, 864L,
> 866L, 932L, 946L, 951L, 964L, 969L, 983L, 985L, 989L, 1032L,
> 1102L, 1102L, 1112L, 1112L, 1199L, 1218L, 1245L, 1252L, 1283L,
> 1302L, 1327L, 1327L, 1328L, 1340L, 1367L, 1382L, 1413L, 1441L,
> 1468L, 1487L, 1493L, 1512L, 1523L, 1567L, 1576L, 1585L, 1598L,
> 1620L, 1623L, 1634L, 1664L, 1704L, 1813L, 1816L, 1824L, 1842L,
> 1892L, 1939L, 1977L, 2010L, 2035L, 2060L, 2100L, 2152L, 2184L,
> 2210L, 2248L, 2282L, 2288L, 2358L, 2462L, 2473L, 2527L, 2541L,
> 2555L, 2576L, 2603L, 2705L, 2711L, 2895L, 2916L, 2958L, 2963L,
> 3088L, 3110L, 3304L, 3427L, 3489L, 3535L, 3629L, 3638L, 3682L,
> 3704L, 3796L, 3932L, 4160L, 4293L, 4294L, 4366L, 4516L, 4715L,
> 4746L, 5152L, 5263L, 5307L, 5823L, 6025L, 6306L, 6313L, 6330L,
> 6649L, 7114L, 7158L, 7340L, 7488L, 7665L, 8055L, 8394L, 8773L,
> 8850L, 8918L, 9273L, 9356L, 10110L, 10281L, 10310L, 10519L, 11111L,
> 13525L, 14465L, 14955L, 16629L)), .Names = "tx_esqui", class = "data.frame", row.names = c(NA,
> -269L))
> 
> 
> --
> Sérgio Henrique Almeida da Silva Junior
> Doutorando em Epidemiologia em Saúde Pública
> Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ
> http://lattes.cnpq.br/1611345552843383
> Tel: (21) 94429486/78101651 id: 123*20942
> 
> 
> 
>