library(readxl)
novilho <- read_excel("G:/Meu Drive/Dropbox/Regressão minha/novilho.xlsx")
attach(novilho)
novilho
# Análise de correlação
pairs(novilho)
# matriz de covariancia dos dados
as.matrix(sig<-round(cov(novilho),2))
sig
# matriz de correlação linear simples dos dados
as.matrix(cs<-round(cor(novilho),2))
cs
# ajuste dos modelos
m1<-lm(PT~Peso, data=novilho)
r1<-residuals(m1)
m2<-lm(AC~Peso, data=novilho)
r2<- residuals(m2)
# correlação parcial entre PT e AC, ajustadas
# para o Peso
cor.test(r1,r2)

# correlação parcial entre PT e AC, ajustadas
# para o Peso, matriz covariância
sig11<-sig[1:2,1:2];as.matrix(sig11)
sig12<-sig[1:2,4];as.matrix(sig12)
sig22<-sig[4,4];as.matrix(sig22)
sig21<-sig[4,1:2];as.matrix(sig21)

sig_ajs<-sig11-(sig12%*%solve(sig22)%*%sig21);
sig_ajs

cparcial<-sig_ajs[1,2]/sqrt(sig_ajs[1,1]*sig_ajs[2,2])
cparcial

# Correlação parcial por meio da matriz de correlação
cs
inv_cs<-solve(cs)
inv_cs
roij<-inv_cs[1,2]/sqrt(inv_cs[1,1]*inv_cs[2,2])
roij

sig_t<-sig[-3,-3];as.matrix(sig_t)
inv_cst<-solve(sig_t)
roijt<-inv_cst[1,2]/sqrt(inv_cst[1,1]*inv_cst[2,2])
roijt