Как интерпретировать результат коэффициента из функции MFAmix в пакете PCAmixdata - PullRequest
0 голосов
/ 08 июля 2019

Я работаю над генетическим проектом в R. В этом наборе данных 497 строк и 11 226 столбцов. Каждая строка представляет собой серию генетических маркеров для определенного гибрида, а каждый столбец представляет собой генетический маркер («G1», «G2» и т. Д.) Со значениями 1, 0, -1 и NA. Я пытаюсь уменьшить размерность, используя MFAmix, чтобы продолжить свои прогнозы.

Вот пример из набора данных.

Hybrid  G1  G5  G8  G9  G10
P1000:P2030 0   -1  0   1   0
P1006:P1384 0   0   0   0   1
P1006:P1401 0   NA  NA  NA  1
P1006:P1412 0   0   0   0   1
P1006:P1594 0   0   0   0   1
P1013:P1517 0   0   0   0   1
res.mfamix<-MFAmix(data=geno3,groups=index,
                  name.groups=names(geno3),ndim=50,rename.level=TRUE,graph=FALSE)
coeff <- res.mfamix$coef

Выходные данные res.mfamix содержат координаты отдельных лиц и количественные переменные. Поскольку теперь у меня есть только 50 измерений, как я могу продолжить свои прогнозы с выводом из res.mfamix? (К какому выводу я буду обращаться)

...