ManonMartin / limpca

https://manonmartin.github.io/limpca/
1 stars 2 forks source link

Trout Vignette #14

Closed bgovaerts closed 1 year ago

bgovaerts commented 1 year ago

In the Trout Vignette on the GitHub web site, the last 3 graphics do not appear after the last chance below.

Effects <- c("Day","Treatment","Day:Treatment") for(i in 1:3){ Pval_log <- -log10(mat_pval[,Effects[i]]) resA=resASCA[[Effects[i]]] PC12Load <- as.vector(sqrt(resA$loadings[,1:2]^2%*%resA$singvar[1:2]^2)) matres <- cbind(PC12Load,Pval_log) A[[i]] <- plotScatter(Y = matres, xy = c(1,2), points_labs = rownames(matres),xlab = "PC1&2 ASCA loadings", ylab = "-log10(p-values)", title = paste(Effects[i],"effect")) } A[[1]] A[[2]] A[[3]]