Skip to contents

Cornejo2025_analysis

Source: vignettes/Cornejo2025_analysis.Rmd
Cornejo2025_analysis.Rmd
processedTongue = parafac4microbiome::processDataCube(Cornejo2025$Tongue, sparsityThreshold=0.5, considerGroups=TRUE, groupVariable="GenderID", centerMode=1, scaleMode=2)

processedSaliva = parafac4microbiome::processDataCube(Cornejo2025$Saliva, sparsityThreshold=0.5, considerGroups=TRUE, groupVariable="GenderID", centerMode=1, scaleMode=2)

processedCytokine = Cornejo2025$Cytokines
modelTongue = parafac4microbiome::parafac(processedTongue$data, nfac=2, nstart=100)
modelSaliva = parafac4microbiome::parafac(processedSaliva$data, nfac=2, nstart=100)
modelCytokine = parafac4microbiome::parafac(processedCytokine$data, nfac=2, nstart=100)
Y = as.numeric(as.factor(Cornejo2025$Tongue$mode1$GenderID))
Ycnt = Y - mean(Y)
tongue_NPLS = triPLS1(processedTongue$data, Ycnt, 1)

Y = as.numeric(as.factor(Cornejo2025$Saliva$mode1$GenderID))
Ycnt = Y - mean(Y)
saliva_NPLS = triPLS1(processedSaliva$data, Ycnt, 1)

Y = as.numeric(as.factor(Cornejo2025$Cytokines$mode1$GenderID))
Ycnt = Y - mean(Y)
cytokine_NPLS = triPLS1(processedCytokine$data, Ycnt, 1)