# Redirecionar o diretório
setwd("C:\\Documents and Settings\\kcc.MOB-PESQUISA\\Meus documentos\\Mestrado\\Dissertação\\Resultados e discussao\\R\\Estudo1")

#leitura dos dados (sem denominá-los)
read.csv2("aaps.csv")

# Leitura do arquivo de dados (dando o nome de "aaps")
aaps<-read.csv2("aaps.csv")

# Exibicao do conteudo dos dados
aaps

# Verificar os nomes das variaveis
names(aaps)

# Resumo das informacoes nos dados
summary(aaps)

# Media do peso das sementes 
mean(aaps$pesosementes)

#media do peso dos cones no ano de 2005
mean(pesocones[ano==2005])

#sempre que for colocar uma condicao lógica usar dois sinais de =!

# Anexando os nomes das variaveis
attach(aaps)

# Desvio padrao do peso do cone em 2006
sd(pesocones)

# Boxplot do peso dos cones em funcao do ano, com nomeacao das coordenadas e coloracao
boxplot(pesocones~ano,xlab="ano", ylab="Peso(g)",main="Peso dos cones em função do ano",col=colors()[637])

# Diagrama de dispersao
plot(pesocones,pesosementes)

#Histograma colorido
hist(pesocones, xlab="Peso dos cones", ylab="Frequencia", main="Histograma da frequencia de peso dos cones",col="steelblue")

#sequencia para salvar o gráfico (sem ve-lo)
x11()
jpeg("aaps.jpg")
boxplot(pesocones~ano)
dev.off()

#para fechar todos os graficos
graphics.off()

# comando teste "t" Considerando amostras pareadas (já considera o teste de homogeneidade das variancias falso e o nivel de confiaca 0,95)
t.test(pesocones[ano==2005],pesocones[ano==2006],paired=T)

# comando teste "t" Considerando amostras nao pareadas
t.test(pesocones[ano==2005],pesocones[ano==2006])