Displaying 1 result from an estimated 1 matches for "pcor4".
Did you mean:
pcor
2007 May 18
1
partial correlation significance
...1, resid.sd, "*"), 2, resid.sd, "*")
#colnames(pcc) <- rownames(pcc) <- colnames(x)
if (test) {
t.df <- ndata - nvar
t <- pcc/sqrt((1 - pcc^2)/t.df)
print(t);
pcc <- list(coefs = pcc, sig = t > qt(1 - (p/2), df = t.df))
}
return(pcc)
}
pcor4 <- function(x, y, z) {
return(cor.test(lm(x~z)$resid,lm(y~z)$resid));
}
[[alternative HTML version deleted]]