powered by
partialcor(formula, data, level = 0.95)
m <- lvm(c(y1,y2,y3)~x1+x2) covariance(m) <- c(y1,y2,y3)~y1+y2+y3 d <- sim(m,500) partialcor(~x1+x2,d)
Run the code above in your browser using DataLab