Я работаю над генетическим проектом в 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? (К какому выводу я буду обращаться)