Проблемы построения / редактирования родословной с родством2 и optiSel - PullRequest
0 голосов
/ 28 июня 2019

У меня есть данные, состоящие из примерно 750 собак, которых я хочу собрать в родословную.Данные включают ID, DADID, MOMID, SEX, коэффициенты инбридинга и цвет шерсти.Я следовал рекомендациям по правильному формату родословных данных и рекомендациям по коду, представленным в пакетах kinship2 и optiSel, но, похоже, не могу понять это.

Я уже создал суб-родословную, котораяотображает несколько человек и их предков / предшественников.Я могу построить его с помощью функции pedPlot из пакета optiSel, но имена отдельных лиц перекрываются и, хотя пакет говорит, что значения родства 2 можно использовать для изменения эстетики, я не смог использовать эти значения для изменения длины ветви,размер текста или даже ширина графика.Я пытаюсь использовать функцию kinship2 plot.pedigree, но у меня возникают еще большие проблемы с этим

library(dplyr)
> library(tidyr)
> library(pedigree)
> library(kinship2)
> library(optiSel)
> library(gridExtra)
> library(ggplot2)
> library(ggthemes)
> View(pedigree)

> ajaxped<-subPed(pedigree, keep = c("Customs Ajax AAD"), prevGen = 
1,succGen = 2)
> ajaxped<-subPed(pedigree, keep = c("Customs Ajax AAD"), prevGen = 1, 
succGen = 2)
> pedplot(ajaxped)
> pedplot(ajaxped, cex = 0.5, width = 8)
> pedplot(ajaxped, cex = 0.5)
#these 2 are the same; the width doesn't change anything
> plot.pedigree(ajaxped,id= ajaxped$ID, Sire = ajaxped$Sire, Dame = 
 ajaxped$Dam, sex = ajaxped$Sex)
Error in matrix(0, nrow = n) : data is too long
> plot.pedigree(ajaxped)       
 Error in matrix(0, nrow = n) : data is too long
> plot.pedigree(pedigree, id = pedigree$ID)       
Error in plot.pedigree(pedigree, id = pedigree$ID) : Wrong length for id

Сюжет, созданный из pedPlot;ajaxped

1 Ответ

0 голосов
/ 04 июля 2019

Если вы посмотрите на пакет Igraph , вы можете легко построить свою родословную

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...