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 <giselle_davi@yahoo.com.br>:
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 <dmarcelino@live.com> 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 <giselle_davi@yahoo.com.br> 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 <dmarcelino@live.com> 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 <nandodesouza@gmail.com>:
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 <giselle_davi@yahoo.com.br> 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 R-br@listas.c3sl.ufpr.br 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 R-br@listas.c3sl.ufpr.br 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 R-br@listas.c3sl.ufpr.br 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.