Я пытаюсь подогнать модель PARAFA C к 4-мерному набору FTIR-спектров, организованных в матрицу:
dim(matrix)
6844, 34, 147
Когда я запускаю
library(multiway)
parafac<-parafac(matrix, nfac=3, nstart=1)
cor<-corcondia(matrix, parafac)
cor
[1] NA
То же самое происходит с любым значением nfac аргумент
Если я правильно понял теорию, лежащую в основе parafa c, этого не должно произойти с тестом согласованности ядра.