Pvalue_ith_coeff <-
function(coef_vec,R,normr,dfSJ,i) {
#VC <- (solve(R)) %*% t(solve(R)) * (normr^2)/dfSJ
VC <- (pinv(R)) %*% t(pinv(R)) * (normr^2)/dfSJ
VCi <- sqrt(VC[i,i])
statistic <- (coef_vec[i])/VCi
p <- 2 * (1-pt(abs(statistic),dfSJ))
return(p)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.