Bioestatística

 

​Coeficiente de correlação

​

Como exemplo, clicar aqui para acessar o banco de dados em Excel

​

Podemos abrir o banco de dados no programa R usando, por exemplo, a função read_excel() do pacote readxl

​

# Usando o pacote readxl

​

install.packages("readxl")

 

library(readxl)

​

dados <- read_excel("E:\\Aulas\\R\\Dados aula 3 ex 2.xlsx")

​

# A função summary() descreve as variáveis do banco de dados

​

summary(dados) 

​

# A função attach() permite que as variáveis do banco de dados sejam acessadas simplesmente digitando seus nomes

       
attach(dados)

​

# A função cor(x,y) retorna a correlação entre as variáveis x e y

​        

cor(Peso,Altura)                     # Retorna a correlação (Pearson) entre Peso e Altura, -0.07115096

​

cor(Peso,Altura,method="pearson")    # Correlação de Pearson,  -0.07115096

​

cor(Peso,Altura,method="spearman")   # Correlação de Spearman, -0.02004432

​

cor(Peso,Altura,method="kendall")    # Correlação de Kendall,  -0.01317936

​

​

A função cor.test() também calcula o coeficiente de correlação, mas ela traz mais detalhes.

​

cor.test(Peso,Altura)

​

        Pearson's product-moment correlation

 

data:  Peso and Altura

t = -0.43972, df = 38, p-value = 0.6626

alternative hypothesis: true correlation is not equal to 0

95 percent confidence interval:

 -0.3743629  0.2458065

sample estimates:

        cor

-0.07115096

​

Neste exemplo, perceba que a função cor.test() também retorna um valor p de um teste de hipóteses que tem por hipótese nula que o coeficiente de correlação é igual a zero na população correspondente, e um intervalo de confiança 95% para o coeficiente de correlação.

​

# Vamos armazenar no objeto "a" o resultado da função cor.test(Peso,Altura)

​

a <- cor.test(Peso,Altura)

​

names(a)

​

[1] "statistic"   "parameter"   "p.value"     "estimate"    "null.value"

[6] "alternative" "method"      "data.name"   "conf.int" 

​

a$p.value    # retorna o valor p do teste de significância do coeficiente de correlação

​

a$estimate   # retorna a estimativa do coeficiente de correlação

​

a$conf.int   # retorna o intervalo de confiança para o coeficiente de  correlação

​

​

# Coeficiente de correlação de Spearman

​

cor.test(Peso,Altura,method="spearman")   

​

        Spearman's rank correlation rho

​

data:  Peso and Altura

S = 10874, p-value = 0.9023

alternative hypothesis: true rho is not equal to 0

sample estimates:

        rho

-0.02004432

​

​

​