O melhor, Giselle, é vc tornar uma amostra do seu arquivo disponível (mesmo que seja com valores trocados, etc). Eu usei como exemplo o arquivo que vc passou como referência no sei primeiro email. On Nov 20, 2013 1:44 PM, "Giselle Davi" <[email protected]> wrote:
> Isso foi a primeira coisa que conferi !!!! Olha só a saída: > > Sys.setenv(HTTP_PROXY="http://giselle:selecaogenomic@2013@ > proxy.cnpms.embrapa.br:3128") > > getwd() > [1] "C:/GS_tutorial" > > dir(getwd()) > [1] "DatasetS2_29SAWSN_gbs.csv" > [2] "DatasetS2_29SAWSN_gbs.txt" > [3] "Geno_RRBLUP_All.txt" > [4] "RData" > > GBS = read.csv('DatasetS2_29SAWSN_gbs.csv', header=TRUE, > stringsAsFactors=FALSE) > > parse.GBS <- function(x) { > + unique.x <- unique(x) > + alleles <- setdiff(unique.x,union("H","N")) > + y <- rep(0,length(x)) > + y[which(x==alleles[1])] <- -1 > + y[which(x==alleles[2])] <- 1 > + y[which(x=="N")] <- NA > + return(y) > + } > > X <- apply(GBS[,-c(1:3)],1,parse.GBS) > > dim(X) > NULL > > E agora ??? > > > > Em Quarta-feira, 20 de Novembro de 2013 13:31, Leandro Marino < > [email protected]> escreveu: > O problema está aqui: > > cannot open file 'tpg12-06-0006-dataset-s2': No such file or directory > o arquivo não existe no local que está. > > veja o diretorio que vc está com getwd() e liste os arquivos do diretorio > com dir(getwd()), você verá que não há nenhum arquivo com esse nome. > > > <http://www.showmetech.com.br/> > *LEANDRO MARINO* | Showmetech > Estatístico | Fotógrafo > Mobile: 55 21 9845-7707 > Email: [email protected] > Profiles SMT: facebook <http://www.facebook.com/showmetech> | > twitter<http://www.twitter.com/showmetech> > | Google+ <https://plus.google.com/115440786884851522659> > > <http://www.leandromarino.com.br/> > > > >
_______________________________________________ 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.
