## Leitura e manipulação de dados
### Ana Carolina Luchetta


##Distância entre cidades
cidades <- c("Atenas", "Madri", "Paris", "Estocolmo")     #concatenando os rótulos das cidades envolvidas
matriz.cidades <- matrix (,4,4)                           #matriz 4 por quarto. A primeira parte está em branco pois eu não tenho valores ainda para preencher.
rownames(matriz.cidades) <- cidades                       #inserindo rótulo de linha na matriz
colnames(matriz.cidades) <- cidades                       #inserindo rótulo de coluna na matriz
matriz.cidades                                            #verifico a matriz com os rótulos e nos valores vejo que consta NA (não inserimos ainda os valores)
dist = c(3949, 3000, 3927, 1273, 3188, 1827)              #distâncias concatenadas
matriz.cidades[lower.tri(matriz.cidades)] <- dist         #colocar os valores na matriz
diag(matriz.cidades) <- 0                                 #atribuindo 0 a diagonal
matriz.cidades                                            #verifico
require(datasets)
eurodist                                                  #retorna as distâncias listadas pelo R em uma matriz semelhante a produzida por mmim anteriormente (onde os valores de Atenas a Estocolmo não são repetidos quando falamos de Estocolmo para Atenas)

##Criação de um data frame
cor <- rep(c("claro", "escuro"), each=9)                  #repetição dos fenótipos
dieta <- rep(c("A","B","C"), each=3,2)                    #repetição das dietas
peso <- c(0.1 , 1.1 , 3.7,  5.7, -1.2, -1.5,  3.0, -0.4, 0.6, 1.5, -0.1, 2.0,  0.6, -3.0, -0.3,  -0.2, 0.3, 1.5)
experimento <- data.frame(cor, dieta, peso)               #tabela com cores, dieta e variação de peso de cada um dos ratos
média.cor <- tapply(peso, cor, mean)                      #na função tapply, colocar(quem tem o mesmo tamanho e é valor numérico, a média de quem, função que quero executar
média.cor                                                 #conferindo a média para cada cor
média.dieta <- tapply(peso, dieta, mean)
média.dieta

##Criando uma matriz
valores=rnorm(15, 10, sqrt(3.6))                          #dostribuição normal, de 15 elementos, com média 10 e desvio padrão (igual raiz quadrada da variância)
valores
matriz <- matrix (valores, ncol=5)                        #monta a tabela com os resultados
matriz
colnames(matriz) <- paste("C",1:5)
colnames(matriz)                                          #modificou o nome das colunas
rownames(matriz) <-paste("L", 1:3)
rownames (matriz)                                         #modificou o nome das linhas
média.coluna <- apply(matriz,2, mean)
média.linha <- apply (matriz, 1, mean)
média.linha                                               #verifico a média da coluna
média.coluna                                              #verifico a média da linha
var.coluna <- apply(matriz,2, var)
var.linha <- apply (matriz, 1, var)
var.linha                                                 #verifico a variância da coluna
var.coluna                                                #verifico a variância da linha
framecol <- data.frame(média.coluna, var.coluna)
framelin <- data.frame(média.linha, var.linha)
framecol                                                  #média e variância da coluna
framelin                                                  #média e variância da linha
junto <- list(framecol, framelin)
junto                                                     #tudo junto em uma tabela

##