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