data(hidden_fac.dat)
fit <- svpls(10,10,hidden_fac.dat,pmax = 5)
fit$genes
Y.corrected <- fit$Y.corr
data(hidden_fac.dat)
gen <- c(1:15,50:65)
sub <- c(1:5,11:17)
hfp(fit,gen,sub,hidden_fac.dat)Run the code above in your browser using DataLab