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@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@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list R-br@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@gmail.com <agrohelder@hotmail.com> *