Eu acho que os dados que você passou não reproduz o procedimento que você está tentando implementar. Contudo, repare que uma das bases de dados tem NA como row.names, enquanto as outras não. Isso pode ter acontecido na geração dos dados a partir de dput() sem controlar a opção "control", mas pode ser que os dados que você tem já tenham herdado isso também.
names(gen); names(gen1); names(gen2) > table(gen %in% colnames(gen1)) FALSE 10 > table(colnames(gen) %in% colnames(gen1)) FALSE 10 2015-01-08 11:31 GMT-02:00 Giselle Davi <[email protected]>: > Então, segue: > dput(sample(gen, 10)): >> dput(sample(gen, 10)) > structure(list(S1_5719699 = c(0L, 1L, 0L, 1L, 0L, 0L, 0L, 1L, > 0L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, > 0L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, > 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 1L, 2L, > 1L, 2L, 1L, 1L, 1L, 2L, 1L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, > 0L, 1L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, > 1L, 0L), S1_228254954 = c(1L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, > 1L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, > 0L, 0L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 1L, > 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, > 0L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, > 1L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, > 0L), S1_35929087 = c(0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 1L, 1L, > 1L, 1L, 2L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, > 2L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, > 1L, 1L, 2L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, > 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, > 0L, 1L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 1L > ), S1_271311840 = c(0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, > 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, > 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, > 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, > 1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, > 0L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L), > S1_110780932 = c(1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, > 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, > 0L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 1L, 0L, > 0L, 0L, 0L, 0L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, > 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, > 0L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 1L, 0L, > 0L, 0L, 0L, 0L, 1L), S1_29391767 = c(0L, 0L, 1L, 0L, 0L, > 0L, 1L, 0L, 0L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 0L, 0L, > 1L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, > 1L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 0L, 0L, 1L, 0L, 0L, > 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, > 1L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, > 1L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L), S1_55455793 = c(1L, > 2L, 2L, 2L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, > 2L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 1L, 0L, 1L, 1L, 1L, > 1L, 1L, 0L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, > 1L, 1L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, > 1L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 0L, 1L, 2L, 2L, 2L, > 2L, 2L, 1L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 1L), > S1_35135950 = c(1L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 0L, 0L, > 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, > 0L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 0L, 0L, 1L, 0L, 0L, > 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 1L, > 2L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, > 1L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, > 1L, 2L, 1L, 1L), S1_5516961 = c(0L, 0L, 0L, 0L, 0L, 1L, 0L, > 1L, 0L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 1L, 0L, 0L, 0L, 0L, > 0L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, > 0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, > 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, > 0L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, > 0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L), S1_30262521 = c(0L, 0L, > 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, > 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, > 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, > 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, > 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, > 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L)), .Names = > c("S1_5719699", > "S1_228254954", "S1_35929087", "S1_271311840", "S1_110780932", > "S1_29391767", "S1_55455793", "S1_35135950", "S1_5516961", "S1_30262521" > ), row.names = c(NA, -90L), class = "data.frame") >> >> dput(sample(gen1, 10)) >> dput(sample(gen1, 10)) > structure(list(S1_55500762 = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, > 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, > 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, > 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, > 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, > 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, > 0L, 0L), S1_93331343 = c(0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, > 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, > 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, > 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, > 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, > 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, > 0L), S1_76339723 = c(2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, > 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, > 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, > 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, > 0L, 0L, 0L, 0L, 0L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, > 0L, 0L, 0L, 0L, 0L, 0L, 0L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L > ), S1_217586127 = c(0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 1L, > 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, > 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, > 1L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, > 1L, 1L, 1L, 2L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, > 0L, 1L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 1L), > S1_289338857 = c(1L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 2L, > 1L, 1L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 0L, 0L, 0L, 1L, 0L, 1L, > 0L, 0L, 2L, 1L, 1L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 0L, 0L, 0L, > 1L, 0L, 1L, 0L, 0L, 2L, 1L, 1L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, > 1L, 1L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 2L, > 1L, 1L, 1L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 2L, 1L, 1L, 1L, > 2L, 1L, 2L, 1L, 1L), S1_189190087 = c(1L, 1L, 1L, 2L, 2L, > 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, > 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, > 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 0L, 0L, 0L, 1L, 1L, > 1L, 1L, 0L, 0L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 0L, 0L, > 0L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 0L, > 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 0L, 0L), S1_82633392 = c(0L, > 0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, > 0L, 0L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 1L, 1L, > 0L, 0L, 0L, 0L, 0L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, > 1L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, > 1L, 1L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, > 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L), > S1_211264709 = c(0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, > 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, > 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, > 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, > 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, > 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, > 1L, 1L, 2L, 1L, 1L), S1_281166059 = c(2L, 2L, 2L, 2L, 1L, > 1L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, > 1L, 1L, 0L, 0L, 1L, 1L, 0L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, > 1L, 1L, 1L, 1L, 1L, 0L, 0L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 0L, > 0L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 0L, 0L, 1L, 1L, 0L, 2L, 2L, > 2L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, > 1L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 1L), S1_250893506 = c(0L, > 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, > 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, > 0L, 0L, 0L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, > 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, > 2L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, > 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L)), .Names = > c("S1_55500762", > "S1_93331343", "S1_76339723", "S1_217586127", "S1_289338857", > "S1_189190087", "S1_82633392", "S1_211264709", "S1_281166059", > "S1_250893506"), row.names = c("D 1xF 1", "D 1xF 2", "D 1xF 3", > "D 1xF 4", "D 1xF 5", "D 1xF 6", "D 1xF 7", "D 1xF 8", "D 1xF 9", > "D 2xF 1", "D 2xF 2", "D 2xF 3", "D 2xF 4", "D 2xF 5", "D 2xF 6", > "D 2xF 7", "D 2xF 8", "D 2xF 9", "D 3xF 1", "D 3xF 2", "D 3xF 3", > "D 3xF 4", "D 3xF 5", "D 3xF 6", "D 3xF 7", "D 3xF 8", "D 3xF 9", > "D 4xF 1", "D 4xF 2", "D 4xF 3", "D 4xF 4", "D 4xF 5", "D 4xF 6", > "D 4xF 7", "D 4xF 8", "D 4xF 9", "D 5xF 1", "D 5xF 2", "D 5xF 3", > "D 5xF 4", "D 5xF 5", "D 5xF 6", "D 5xF 7", "D 5xF 8", "D 5xF 9", > "D 6xF 1", "D 6xF 2", "D 6xF 3", "D 6xF 4", "D 6xF 5", "D 6xF 6", > "D 6xF 7", "D 6xF 8", "D 6xF 9", "D 7xF 1", "D 7xF 2", "D 7xF 3", > "D 7xF 4", "D 7xF 5", "D 7xF 6", "D 7xF 7", "D 7xF 8", "D 7xF 9", > "D 8xF 1", "D 8xF 2", "D 8xF 3", "D 8xF 4", "D 8xF 5", "D 8xF 6", > "D 8xF 7", "D 8xF 8", "D 8xF 9", "D 9xF 1", "D 9xF 2", "D 9xF 3", > "D 9xF 4", "D 9xF 5", "D 9xF 6", "D 9xF 7", "D 9xF 8", "D 9xF 9", > "D 10xF 1", "D 10xF 2", "D 10xF 3", "D 10xF 4", "D 10xF 5", "D 10xF 6", > "D 10xF 7", "D 10xF 8", "D 10xF 9"), class = "data.frame") >> > dput(sample(gen2, 10)) >> dput(sample(gen2, 10)) > structure(list(S1_95703374 = c(0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, > 0L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, > 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 1L, 2L, 1L, > 1L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 0L, 0L, > 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, > 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, > 1L, 0L), S1_19948955 = c(1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, > 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, > 0L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 0L, 0L, 0L, 0L, > 1L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 1L, 0L, 0L, > 0L, 0L, 1L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 2L, > 1L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, > 2L), S1_64972675 = c(1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, > 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, > 2L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, > 0L, 0L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 0L, 0L, 0L, > 0L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, > 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L > ), S1_212259765 = c(1L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 2L, 1L, > 1L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, > 1L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 1L, > 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, > 0L, 1L, 0L, 0L, 2L, 1L, 1L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 0L, 0L, > 0L, 1L, 0L, 1L, 0L, 0L, 2L, 1L, 1L, 1L, 2L, 1L, 2L, 1L, 1L), > S1_55454770 = c(0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, > 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, > 0L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, > 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, > 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, > 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, > 0L, 1L, 0L, 0L), S1_49131053 = c(1L, 1L, 1L, 2L, 2L, 2L, > 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 0L, 0L, 0L, > 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, > 0L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, > 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, > 1L, 1L, 1L, 0L, 0L, 0L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, > 1L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 1L), S1_195416162 = c(0L, > 1L, 0L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 1L, 1L, 0L, > 1L, 0L, 0L, 1L, 0L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, > 1L, 1L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, > 1L, 0L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 1L, 1L, 0L, > 1L, 0L, 1L, 2L, 1L, 2L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, > 2L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 2L, 2L, 1L, 2L, 1L), > S1_293795584 = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, > 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, > 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, > 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, > 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, > 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, > 1L, 1L, 1L, 1L, 1L), S1_43059180 = c(0L, 0L, 0L, 0L, 0L, > 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, > 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, > 0L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, > 2L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, > 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, > 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L), S1_298690372 = c(1L, > 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, > 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, > 0L, 0L, 0L, 0L, 0L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, > 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, > 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, > 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L)), .Names = > c("S1_95703374", > "S1_19948955", "S1_64972675", "S1_212259765", "S1_55454770", > "S1_49131053", "S1_195416162", "S1_293795584", "S1_43059180", > "S1_298690372"), row.names = c("D 1xF 1", "D 1xF 2", "D 1xF 3", > "D 1xF 4", "D 1xF 5", "D 1xF 6", "D 1xF 7", "D 1xF 8", "D 1xF 9", > "D 2xF 1", "D 2xF 2", "D 2xF 3", "D 2xF 4", "D 2xF 5", "D 2xF 6", > "D 2xF 7", "D 2xF 8", "D 2xF 9", "D 3xF 1", "D 3xF 2", "D 3xF 3", > "D 3xF 4", "D 3xF 5", "D 3xF 6", "D 3xF 7", "D 3xF 8", "D 3xF 9", > "D 4xF 1", "D 4xF 2", "D 4xF 3", "D 4xF 4", "D 4xF 5", "D 4xF 6", > "D 4xF 7", "D 4xF 8", "D 4xF 9", "D 5xF 1", "D 5xF 2", "D 5xF 3", > "D 5xF 4", "D 5xF 5", "D 5xF 6", "D 5xF 7", "D 5xF 8", "D 5xF 9", > "D 6xF 1", "D 6xF 2", "D 6xF 3", "D 6xF 4", "D 6xF 5", "D 6xF 6", > "D 6xF 7", "D 6xF 8", "D 6xF 9", "D 7xF 1", "D 7xF 2", "D 7xF 3", > "D 7xF 4", "D 7xF 5", "D 7xF 6", "D 7xF 7", "D 7xF 8", "D 7xF 9", > "D 8xF 1", "D 8xF 2", "D 8xF 3", "D 8xF 4", "D 8xF 5", "D 8xF 6", > "D 8xF 7", "D 8xF 8", "D 8xF 9", "D 9xF 1", "D 9xF 2", "D 9xF 3", > "D 9xF 4", "D 9xF 5", "D 9xF 6", "D 9xF 7", "D 9xF 8", "D 9xF 9", > "D 10xF 1", "D 10xF 2", "D 10xF 3", "D 10xF 4", "D 10xF 5", "D 10xF 6", > "D 10xF 7", "D 10xF 8", "D 10xF 9"), class = "data.frame") >> > > > Desde já agradeço a atenção. > > Att, > > Giselle > > > > Em Quinta-feira, 8 de Janeiro de 2015 11:15, Daniel Marcelino > <[email protected]> escreveu: > > > Não resolve muito. O que voce mandou sao apenas os nomes das colunas, o que > diz absolutamente nada sobre como seus dados estão salvos internamente. > O comando é o dput mesmo. O que voce pode fazer é selecionar uma pequena > amostra dos dados tipo: dput(sample(dados, 10)) > > Daniel > > Please, excuse me for typos and brevity. Sent from my iPhone > > On Jan 8, 2015, at 9:01 AM, Giselle Davi <[email protected]> wrote: > > Prezados, > > > Utilizei o comando (dput()) sugerido, porém a saída vira tudo texto imagino > que deve ser por causa do tamanho do arquivo.Com o objetivo de tornar o > código reproduzível achei melhor salvar os arquivos com os head de cada > objeto e anexá-los ao e-mail. Imagino que assim tenha-se uma visão melhor > dos dados. Desde já agradeço a atenção. > > Att, > > Giselle > > > > > Em Quarta-feira, 7 de Janeiro de 2015 22:28, Daniel Marcelino > <[email protected]> escreveu: > > > Sem uma amostra dos dados ou informações mais detalhadas a estrutura > do banco é difícil oferecer uma solução. Mas eu acredito que possa ser > problema envolvendo os níveis dos fatores. Veja ?droplevels . > > Daniel > > 2015-01-07 17:24 GMT-02:00 Fernando Antonio de souza > <[email protected]>: >> Sua questão é interessante mas seu código não é reproduzível, o que >> difulcuta a ajuda. faça um dput() de parte de seu banco de dados copie e >> cole no e-mail para que possamos ter facil acesso e assim poder lhe >> ajudar. >> >> Seja mais clara sobre o seu problema para obter ajuda mais rapidamente. >> >> Em 7 de janeiro de 2015 15:13, Giselle Davi <[email protected]> >> escreveu: >>> >>> Prezados, >>> >>> Alguém poderia me ajudar a selecionar um dataframe, para este caso o nome >>> é list que corresponde à Lista de SNP's selecionados, a partir de de >>> outro >>> dataframe, que para este caso corresponde ao gen1 ( sem precisar >>> indicar a >>> posição das colunas )? >>> >>> Utilizei a rotina especificada abaixo, mas o comando não me disponibiliza >>> a seleção de todas as 22894 colunas, ou seja, de todos os SNP's. Existem >>> 7 >>> deles que não conseguem ser selecionados por serem declarados como >>> falsos. >>> Alguém poderia me ajudar ? Desde já agradeço. >>> >>> Att, >>> >>> Giselle >>> >>> >>> ROTINA: >>> #1)Lendo arquivo genotipos >>> gen=read.table("matriz_genotipos.txt",h=T)#lendo arquivo genotipos >>> dim(gen)#90 x 24165. Tem que colocar o nome dos cruzamentos na primerira >>> coluna >>> head(gen) >>> >>> #2)Gerando o arquivo com o nome dos SNP's selecionados >>> nomes<- read.table("D_F_nomes.txt",h=F, sep=",",row.names=1,quote="") >>> dim(nomes)#90 x 0 >>> >>> #3)Arquivo de genótipos com os indivíduos nas linhas e com todos SNP's >>> nas >>> colunas >>> gen1<-cbind(nomes,gen) >>> head(gen1)#juntou ! >>> dim(gen1)#90 x 24165 >>> >>> #4)Lendo o arquivo de frequência alélica >>> freq=read.table("chr1_alefreq.txt",h=T) #lendo arquivo freq alel >>> dim(freq)#22894 x 4 >>> head(freq) >>> >>> #5) Gerando o arquivo com os SNP's selecionados pela MAF >>> lista_SNP_sel<- freq$nun_snp_freq >>> head(lista_SNP_sel)#Está como fator >>> list<-as.character(lista_SNP_sel) >>> head(list) >>> length(list) >>> dim(gen1) >>> >>> #6)Verificando se os SNP's selecionados são realmente os 22894 >>> selecionados pela MAF >>> table(list%in%colnames(gen1)) >>> >>> #SURPRESA:Existem 7 SNP's lidos como falsos e por isso não completam os >>> 22894 ! >>> logico<-colnames(gen1)%in%list >>> head(logico) >>> head(gen1) >>> head(list) >>> gen2<- gen1[,logico] >>> head(gen2) >>> dim(gen2)#90 x 22887 ( Existem 7 snps que estão sendo declarados como >>> falsos, por que ?) >>> >>> >>> _______________________________________________ >>> R-br mailing list >>> [email protected] >>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >>> código mínimo reproduzível. >> >> >> >> _______________________________________________ >> R-br mailing list >> [email protected] >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >> código >> mínimo reproduzível. > > <head(gen).txt> > > <head(gen1).txt> > > <head(freq).txt> > > <head(list).txt> > > <head(gen2).txt> > > > > > _______________________________________________ > R-br mailing list > [email protected] > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código > mínimo reproduzível. _______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
