[R-br] Transformação BoxCox no R

[Helder Gramacho]

Olá Walmes,

Obrigado pela resposta, minha dúvida ainda é bem básica olha só como estou
fazendo:

> dados<-read.table(file="DADOS_N_NORMAIS.csv",sep=",",header=T,dec=".")
> dados
                     REG     VTN1     VTN2
1               Andradina  7713.50  9090.91
2              Araraquara 10743.80 10927.46
3               Araçatuba  8677.69  9917.36
4                   Assis     0.00  5268.60
5                   Avaré  8677.69  8790.38
6                Barretos 10330.58 10330.58
7                   Bauru     0.00     0.00
8                Botucatu  9445.10 10909.09
9       Bragança Paulista 11386.59 12267.56
10              Campinas   9228.65 13739.67
11              Catanduva  8795.75 10743.80
12                Dracena  4700.41  5920.80
13          Fernandópolis  7988.98  8677.69
14                Franca   8801.65  9691.96
15        General Salgado  8719.01  9356.55
16          Guaratinguetá  2561.98  3246.75
17           Itapetininga  9447.78 10800.15
18               Itapeva   5406.34  6060.61
19            Jaboticabal  8264.46 11983.47
20                 Jales      0.00 10123.97
21                    Jaú  5475.21  9194.21
22                Limeira 12534.44 13016.53
23                  Lins   8884.30 10468.32
24               Marília   5454.55  6311.04
25             Mogi Mirim 13085.40 15289.26
26       Mogi das Cruzes   8942.52  9513.09
27               Orlândia     0.00 13774.10
28              Ourinhos   9958.68 10440.77
29       Pindamonhangaba   4958.68  7914.81
30             Piracicaba  7300.28 10212.51
31   Presidente Prudente   3751.15  4309.33
32   Presidente Venceslau     0.00     0.00
33              Registro   2654.96  1977.11
34         Ribeirao Preto 11404.96 11511.22
35               Sorocaba 14462.81 13688.02
36 Sao José do Rio Preto  13429.75 14003.67
37  Sao Joao da Boa Vista  9435.26  9729.53
38              Sao Paulo  6466.94 10789.26
39                   Tupa  2768.60  3140.50
40           Votuporanga   8884.30  9366.39
> attach(dados)
> shapiro.test(VTN1)

        Shapiro-Wilk normality test

data:  VTN1
W = 0.9369, p-value = 0.02713

> library(MASS)
> boxcox(VTN1)
Erro: $ operator is invalid for atomic vectors

Chega neste ponto estou fazendo algo muito errado só não sei o que é...

Hélder


2013/6/25 walmes . <walmeszeviani em gmail.com>

> No primeiro caso indica que uma transformação log pode ser considerada
> (aceita a hipótese de que lambda==0).
>
> m0 <- lm(Volume~log(Height)+log(Girth), data=trees)
> par(mfrow=c(2,2))
> plot(m0)
> layout(1)
>
> require(MASS)
> boxcox(m0, lambda=seq(-0.25, 0.25, length=10))
> abline(v=0, col=2)
>
> m1 <- lm(log(Volume)~log(Height)+log(Girth), data=trees)
> par(mfrow=c(2,2))
> plot(m1)
> layout(1)
>
> À disposição.
> Walmes.
>
> ==========================================================================
> Walmes Marques Zeviani
> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
> Departamento de Estatística - Universidade Federal do Paraná
> fone: (+55) 41 3361 3573
> VoIP: (3361 3600) 1053 1173
> e-mail: walmes em ufpr.br
> skype: walmeszeviani
> twitter: @walmeszeviani
> homepage: http://www.leg.ufpr.br/~walmes
> linux user number: 531218
> ==========================================================================
>
>
> _______________________________________________
> 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.
>



-- 
Hélder Gramacho dos Santos
Engenheiro Agrônomo
Mestrando Ciênc. Geod. Tec. Geoinformação-UFPE
*
agrohelder em gmail.com <agrohelder em hotmail.com>
*
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130625/33f43356/attachment.html>