setwd("d:\\Estatística")
data = read.table("wine.txt", header = TRUE)
sample1 = read.table("sample1.txt", header = TRUE)
sample2 = read.table("sample2.txt", header = TRUE)

print("Tamanho da amostra:")
print(length(data$pH))
print()

print("Dados de quality:")
print("Media:")
print(mean(data$quality))
print("Mediana:")
print(median(data$quality))
print("Desvio:")
print(sd(data$quality))
print("Variancia:")
print(var(data$quality))
print("Quartis:")
print(quantile(data$quality))
hist(data$quality)
print()

print("Dados de pH:")
print("Media:")
print(mean(data$pH))
print("Mediana:")
print(median(data$pH))
print("Desvio:")
print(sd(data$pH))
print("Variancia:")
print(var(data$pH))
print("Quartis:")
print(quantile(data$pH))
hist(data$pH)
print()

print("Dados de citric_acid:")
print("Media:")
print(mean(data$citric_acid))
print("Mediana:")
print(median(data$citric_acid))
print("Desvio:")
print(sd(data$citric_acid))
print("Variancia:")
print(var(data$citric_acid))
hist(data$citric_acid)
print("Quartis:")
print(quantile(data$citric_acid))
print()

print("Dados de residual_sugar:")
print("Media:")
print(mean(data$residual_sugar))
print("Mediana:")
print(median(data$residual_sugar))
print("Desvio:")
print(sd(data$residual_sugar))
print("Variancia:")
print(var(data$residual_sugar))
print("Quartis:")
print(quantile(data$residual_sugar))
hist(data$residual_sugar)
print()

print("Dados de alcohol:")
print("Media:")
print(mean(data$alcohol))
print("Mediana:")
print(median(data$alcohol))
print("Desvio:")
print(sd(data$alcohol))
print("Variancia:")
print(var(data$alcohol))
print("Quartis:")
print(quantile(data$alcohol))
hist(data$alcohol)
print()

print("Dados de sample1:")
print("Media:")
print(mean(sample1$alcohol))
print("Mediana:")
print(median(sample1$alcohol))
print("Desvio:")
print(sd(sample1$alcohol))
print("Variancia:")
print(var(sample1$alcohol))
print("Quartis:")
print(quantile(sample1$alcohol))
hist(sample1$alcohol)
print()

print("Dados de sample2:")
print("Media:")
print(mean(sample2$alcohol))
print("Mediana:")
print(median(sample2$alcohol))
print("Desvio:")
print(sd(sample2$alcohol))
print("Variancia:")
print(var(sample2$alcohol))
hist(sample2$alcohol)
print("Quartis:")
print(quantile(sample2$alcohol))
print()