Gisele, boa noite! Rode o código abaixo e havendo erros, copie e cole a mensagem mostrada no console em seu próximo email. ### <code r> install.packages(c("synbreed", "synbreedData"), dep=TRUE) ### instala pacotes library(synbreed) ### requisita pacote library(synbreedData) ### requisita pacote (.packages()) ### verifica pacotes carregados! data(maize) newDHpheno <- data.frame(Trait=1000,row.names="newDH"); newDHpheno #simulating genotypic data newDHgeno <- matrix(sample(c(0,1), ncol(maize$geno), replace=TRUE), nrow=1) rownames(newDHgeno) <- "newDH" #new pedigree newDHpedigree <- data.frame(ID="newDH", Par1=0, Par2=0, gener=0) #new covar information newDHcovar <- data.frame(family=NA, DH=1, tbv=1000, row.names="newDH") # add individual maize2 <- add.individuals(maize,newDHpheno,newDHgeno,newDHpedigree,newDHcovar) summary(maize2) ### </code> Aqui roda sem problemas... Veja a saída de summary(maize2)logo abaixo:
summary(maize2) object of class 'gpData' covar No. of individuals 1611 phenotyped 1251 genotyped 1251 pheno No. of traits: 1
Trait Min. : 120.7 1st Qu.: 142.8 Median : 148.9 Mean : 149.6 3rd Qu.: 154.9 Max. :1000.0 geno No. of markers 1117 genotypes 0 1 frequencies 0.340134 0.659866 NA's 0.000 % map No. of mapped markers 1117 No. of chromosomes 10 markers per chromosome 1 2 3 4 5 6 7 8 9 10 76 96 99 122 85 106 154 130 121 128 pedigree Number of individuals 1611 Par 1 219 Par 2 221 generations 15
Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]