#' ---
#' title: "LCE0602 - Estatística Experimental - Aula01"
#' author: "Renata Alcarde Sermarini"
#' date: "08/2023"
#' ---

#' # Funções básicas e objetos do R 
#' 
#' ## R como uma calculadora

# comentário

#' # Exemplo
#'
#' Os dados a seguir correspondem ao peso de 30 colmos de cana-de-açúcar (em g). 
#'

#' ## Conjunto de dados
#' 
x <- c(11.5,	18.7,	15,	
       17.6,	12.4,	21,	
       14.4,	15.2,	16.9,	14.5,
       16.2,	14.1,	15.7,	16,	
       15.7,	13.1,	15.2,	18,	
       15.1,	19.6,
       18.1,	14.7,	16.4,	13.8,	
       18.8,	15.5,	17.6,	16.8,	
       15.2,	13.7)

x[19]
sort(x)
order(x)

#' ## Calculando medidas de posição e dispersão
#' 
#' 
#' - Média
mean(x) 
mean(c(11.5,	18.7,	15,	
       17.6,	12.4,	21,	
       14.4,	15.2,	16.9,	14.5,
       16.2,	14.1,	15.7,	16,	
       15.7,	13.1,	15.2,	18,	
       15.1,	19.6,
       18.1,	14.7,	16.4,	13.8,	
       18.8,	15.5,	17.6,	16.8,	
       15.2,	13.7))

#' - Variância
#' 
var(x)

#' - Desvio-padrão
#' 
sd(x)
#' ### Exercício: calcule o coeficiente de variação
100*sd(x)/mean(x)


#' 
#' - Quantis (ou percentis)
#' 
quantile(x)
quantile(x, type = 5)
quantile(x,0.975)
quantile(x,0.025)

#' ### Exercício: Calcule a amplitude interquartílica
# q3 - q1

quantile(x,0.75) - quantile(x,0.25)


#' ### Exercício: Observe e interprete o resultado dos seguintes comandos
#' 
min(x)
max(x)
summary(x)
range(x)



#' ## Representação gráfica 
#' 
#' ### Histograma
#' 
hist(x)

#' Verificando argumentos da função hist()
#' 
?hist
#' Algumas modificações
#' 
hist(x, 
     col="red",
     border="black")
hist(x, 
     col="red", 
     border="blue",
     freq=FALSE)

#' ### Diagrama de ramo-e-folhas
stem(x)
#' ### Boxplot
boxplot(x)

#' texto.....
#' 
#' 
#' ## Intervalo de confiança (variância desconhecida)
#' ### 1-$\alpha$ = 0,95
t.test(x,
       conf.level=0.95)
#' ### 1-$\alpha$ = 0,99
t.test(x,
       conf.level=0.99)