- R-br - DINF UFPR

Como construir gráficos para analisar o tempo de execução.
by Jackeline Bonetti Campos 02 Abr '13

02 Abr '13

Boa noite! Tenho 5 algoritmos quaisquer e seus respectivos tempos de execução dependentes da entrada n: T(n)= (n/1000)^2; T(n)= (n/1000)^3 ; T(n)= (n/1000)^0,5; T(n)= (n/1000)^0,2; e T(n)= log2 * (n/1000). Como posso esboçar um gráfico para comparar os tempos de execução? Agradeço desde já! Att,Jackeline.

2 2

[OFF-TOPIC] Vagas de doutorado no programa "Scaling problems in statistics" (Alemanha)
by Tatiane Beduschi 02 Abr '13

02 Abr '13

Caros colegas, O programa de doutorado do qual faco parte está com vagas abertas para doutorado nas áreas de ecologia, matemática, economia, estatística e genética. Peco desculpas por usar esse espaco para propaganda, mas, como o programa é muito voltado para estatística, achei que poderia ser interessante para muitos participantes desta lista. As inscricoes vao até 15 de Maio. Mais informacoes podem ser encontradas neste link: http://www.uni-goettingen.de/en/156593.html Atenciosamente, Tatiane -- MSc. Tatiane Beduschi RTG 1644 "Scaling Problems in Statistics" University of Goettingen DNPW, Agroecology Grisebachstr. 6 D-37077 Goettingen Germany tatibed(a)gmail.com

1 0

Grafico de trajetórias - zap plot
by Pedro Emmanuel Alvarenga Americano do Brasil 01 Abr '13

01 Abr '13

Amigos de R, Estou fazendo umas análises iniciais de trajetórias de biomarcadores ao longo de um tempo em que pacientes receberam um tratamento. Eu gostaria de ter uma percepção do conjunto de trajetórias sem ajuste e procurei por ai formas de faze-lo. Encontrei stats::interaction.plot e epicalc::followup.plot. Mas por algum motivo desconhecido, essas duas funções estavam truncando as trajetórias em 120 dias, mesmo solicitando para mostrar as trajetórias em 180 dias que eu sei o que é o maximo de tempo observado contido nos dados que tenho. Assim, depois de bater cabeça um tanto investi um dia inteiro para fazer a função abaixo e compartilho com voces. Essa função é bem semelhante a followup.plot, mas, quando se deseja não plotar todos os sujeitos para não poluir o gráfico, ao inves de selecionar sujeitos aleatórios, eu coloquei um argumento chamado zap que seleciona n sujeitos como valores iniciais, minimo, maximo e mediano. # data = dataset in long format # id = character indicating the name of an unique id variable for each patient in data (tested with numeric variable) # time = a character indicating the name of a numeric variable representing the time of the measure # outcome = a character indicating the name of a numeric variable representing the outcome # plot.arg = a list of arguments to pass to plot.default # line.arg = one or more lists of arguments to pass to lines (one list to each strata of by) # zap = if auto then the n subjects with the lower, median, higher initial outcome values will be ploted, anything different from "auto" will plot all subjects # n = number of subject to use in zap # by = a character indicating the name of a factor/character variable representing strata for ploting # think about making conditions for zap such as only the median, only the lowest or only the highest, the final highest, final lowest, final median etc zapplot <- function(data,id,time,outcome, plot.arg=list(xlab='time',ylab=outcome), line.arg=list(list(lty=2,col='blue'),list(lty=3,col='red')), leg.args=list(x='top',inset=c(0,-.15), horiz=TRUE,xjust=.5,col=c('blue','red'),lty=2:3, bty='n',xpd=T), zap='auto',n=1,by=NULL){ # if(!is.null(by)){ # if(length(line.arg)!=nlevels(by)){ # stop('If by is porvided, number of lists in line.arg and number of levels in by should be the same! Try improving the lists in line.arg.')} # } data <- data[which(!is.na(data[,time]) & !is.na(data[,outcome])),] plot.arg <- c(list(x = data[,time],y = data[,outcome],type = 'n'),plot.arg) do.call(plot,plot.arg) if(!is.null(by)){ if(!is.factor(data[,by])){ data[,by] <- as.factor(data[,by]) } tmp2 <- data for(j in 1:nlevels(data[,by])){ # i=2 data <- tmp2 data <- data[which(data[,by]==levels(data[,by])[j]),] tmp <- data[which(data[,time]==min(data[,time])),] tmp <- tmp[order(tmp[,outcome]),] # tmp[,c(id,time,outcome)] id.max <- tmp[(nrow(tmp)-(n-1)):nrow(tmp),id] # id.min <- tmp[(1:n),id] # data[which(data$id==id.min),outcome] tmp$med <- abs(tmp[,outcome]-median(tmp[,outcome],TRUE)) tmp <- tmp[order(tmp$med),] id.med <- tmp[1:n,id] if(zap=='auto'){ for(i in 1:length(id.min)){ line.arg.min <- c(list(x = data[which(data[,id]==id.min[i]),time],y = data[which(data[,id]==id.min[i]),outcome]),line.arg[[j]]) do.call(lines,line.arg.min) } for(i in 1:length(id.med)){ line.arg.med <- c(list(x = data[which(data[,id]==id.med[i]),time],y = data[which(data[,id]==id.med[i]),outcome]),line.arg[[j]]) do.call(lines,line.arg.med) } for(i in 1:length(id.max)){ line.arg.max <- c(list(x = data[which(data[,id]==id.max[i]),time],y = data[which(data[,id]==id.max[i]),outcome]),line.arg[[j]]) do.call(lines,line.arg.max) } } else{ uni <- unique(data[,id]) for(i in 1:length(uni)){ line.arg.all <- c(list(x=data[which(data[,id]==uni[i]),time],y=data[which(data[,id]==uni[i]),outcome]),line.arg[[j]]) do.call(lines,line.arg.all) } } } # COLOCAR AQUI OS COMANDOS DE LEGENDA para identificar as diferentes categorias do by leg.args <- c(legend=list(levels(data[,by])),leg.args) do.call(legend,leg.args) } else{ tmp <- data[which(data[,time]==min(data[,time])),] tmp <- tmp[order(tmp[,outcome]),] # tmp[,c(id,time,outcome)] id.max <- tmp[(nrow(tmp)-(n-1)):nrow(tmp),id] # id.min <- tmp[(1:n),id] # data[which(data$id==id.min),outcome] tmp$med <- abs(tmp[,outcome]-median(tmp[,outcome],TRUE)) tmp <- tmp[order(tmp$med),] id.med <- tmp[1:n,id] if(zap=='auto'){ for(i in 1:length(id.min)){ line.arg.min <- c(list(x = data[which(data[,id]==id.min[i]),time],y = data[which(data[,id]==id.min[i]),outcome]),line.arg[[1]]) do.call(lines,line.arg.min) } for(i in 1:length(id.med)){ line.arg.med <- c(list(x = data[which(data[,id]==id.med[i]),time],y = data[which(data[,id]==id.med[i]),outcome]),line.arg[[1]]) do.call(lines,line.arg.med) } for(i in 1:length(id.max)){ line.arg.max <- c(list(x = data[which(data[,id]==id.max[i]),time],y = data[which(data[,id]==id.max[i]),outcome]),line.arg[[1]]) do.call(lines,line.arg.max) } } else{ uni <- unique(data[,id]) for(i in 1:length(uni)){ line.arg.all <- c(list(x=data[which(data[,id]==uni[i]),time],y=data[which(data[,id]==uni[i]),outcome]),line.arg[[1]]) do.call(lines,line.arg.all) } } } } # zapplot(b,'pront','dia','HLA',zap='all',n=2,by='CD4cod') # Abraço forte e que a força esteja com voces, Dr. Pedro Emmanuel A. A. do Brasil Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365, CEP 21040-360, Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

1 0

Deixar as barras e linhas mais próximas dos eixos y no lattice.
by Alisson Lucrecio 01 Abr '13

01 Abr '13

Caros colegas r-br, Boa noite. Em um gráfico xyplot da lattice com sequencia de números contínuos no eixo x posso inserir xlim = (0, 100), isso aproxima o conteúdo do gráfico para as linhas do eixo y. Entretanto, em um gráfico onde tenho no eixo x data não tenho como usar xlim, assim queria colocar as barras e as linhas mais próximas dos eixos y e que não fique espeço vazio entre o inicio e final dos dados. Obrigado. Alisson Lucrécio da Costa

1 0

Repartir matriz....
by andrebvs＠bol.com.br 01 Abr '13

01 Abr '13

3 5

salvar grafico do R em 600dpi
by Ana Carolina Martins 01 Abr '13

01 Abr '13

Oi colegas, preciso enviar para publicação figuras geradas no R com formato TIFF e resolução de 600 dpi.porém quando salvo em tiff ele automaticamente salva com 300dpi.alguém sabe como mudar essa resolução? agradeço desde já!Carol

2 1

Cruzamento entre variáveis
by Wellington em Cristo Jesus Voltará 01 Abr '13

01 Abr '13

Caros, Gostaria de saber se dar para o R fazer cruzamentos entre variáveis, e se houver essapossibilidade como poderia ser feito?E ainda estou com outro problema, os valores ausentes (missing), como poderia ser feitoesses cruzamentos sem perder as demais informações? Ex: (Cruzar a variável Sexo com a variável X1? Semestre com a variável X3?) Sexo Semestre X1 X2 X31 3 2 3 21 3 1 3 12 1 2 31 3 1 22 2 4 1 12 2 1 3 1 Att. Wellington ReinaldoEstatístico

3 3

Re: [R-br] Distribuição normal truncada
by Hélio Gallo Rocha 01 Abr '13

01 Abr '13

Luiz, achei runif não achei o pacote msn se puder dar uma opinião sobre meu problema, agradeceria... muito Abraço Hélio Em 1 de abril de 2013 11:57, Luiz Roberto Martins Pinto < luizroberto.uesc(a)gmail.com> escreveu: > Hélio, > > 1) Agradeço por responder e compartilhar seu problema. > 2) Para encontrar o pacote que contém tnorm faça o seguinte: > library.packages('sos') > library(sos) > ???tnorm > Vc verá que o tnorm está no pacote msn. Então use > library.packages('msn') > library(msn) > pronto para usar tnorm. > > > Luiz Roberto Martins Pinto > Prof. Pleno/DCET/UESC > Laboratório de Estatística Computacional > Universidade Estadual de Santa Cruz > Ilhéus-Bahia > > luizroberto.uesc(a)gmail.com > skype: lrmpinto > http://lattes.cnpq.br/2732314327604831 > > > > > Em 1 de abril de 2013 11:29, Hélio Gallo Rocha <heliogallorocha(a)gmail.com>escreveu: > >> Luiz Roberto. >> >> Primeiro, não sou expert nem em estatistica e nem em R, vou tentar >> resumir em itens pra ficar mais claro >> >> 1. estou fazendo previsão de safra em café; >> 2. testando 4 modelos matemáticos, que levam em consideração: altura, >> diametro e contagem de frutos; >> 3. agora estou na época de amostragem, a partir de maio colho e verei >> qual é o melhor método.; >> 4. sendo o café uma planta "perene", a variação de dados é alta, não >> tanto altura e diametro, Cv +/- 10%, mas na contagem de frutos o CV >> costumar dar acima de 40%; >> 5. então é natural a minha preocupação com amostragem seja grande, não >> quero que minha conclusão do experimento seja "inconclusiva" por falta de >> amostragem; >> 6. Já amostrei 86 pontos, variando amostradores, metodologia de >> amostragem dos frutos ( aleatória e não); >> 7. então tenho sd e mean; >> 9. tem um post da semana passada com um script do Walmes usando rnorm, >> que justamente é sobre o n a ser amostrado; >> 10. com relação a n coloco a quantidade total de plantas do talhão, não >> sei se poso fazer isso, ainda estou estudando; >> 11. rodei e verifiquei que ela gera dados com máximo e minino maiores do >> que já tenho, então to querendo que estes dados sejam mais próximos da >> realidade; >> >> tentei rodar o rtnorm mas ele deve pertencer a alguma pacote, pois o R dá >> o erro: >> Erro: não foi possível encontrar a função "rtnorm" >> >> estou pesquisando >> >> abraço >> >> Hélio >> >> >> Em 1 de abril de 2013 09:57, Luiz Roberto Martins Pinto < >> luizroberto.uesc(a)gmail.com> escreveu: >> >> Ola Hélio, >>> >>> Interessante seu interesse pela geração de dados com distribuição normal >>> truncada. >>> Você pode compartilhar o motivo/aplicação? >>> Luiz Roberto >>> >>> Luiz Roberto Martins Pinto >>> Prof. Pleno/DCET/UESC >>> Laboratório de Estatística Computacional >>> Universidade Estadual de Santa Cruz >>> Ilhéus-Bahia >>> >>> luizroberto.uesc(a)gmail.com >>> skype: lrmpinto >>> http://lattes.cnpq.br/2732314327604831 >>> >>> >>> >> >> >> -- >> Hélio Gallo Rocha >> IFSULDEMINAS - Câmpus Muzambinho >> >> _______________________________________________ >> R-br mailing list >> R-br(a)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. >> > > -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

1 0

Re: [R-br] Distribuição normal truncada
by Hélio Gallo Rocha 01 Abr '13

01 Abr '13

Luiz Roberto. Primeiro, não sou expert nem em estatistica e nem em R, vou tentar resumir em itens pra ficar mais claro 1. estou fazendo previsão de safra em café; 2. testando 4 modelos matemáticos, que levam em consideração: altura, diametro e contagem de frutos; 3. agora estou na época de amostragem, a partir de maio colho e verei qual é o melhor método.; 4. sendo o café uma planta "perene", a variação de dados é alta, não tanto altura e diametro, Cv +/- 10%, mas na contagem de frutos o CV costumar dar acima de 40%; 5. então é natural a minha preocupação com amostragem seja grande, não quero que minha conclusão do experimento seja "inconclusiva" por falta de amostragem; 6. Já amostrei 86 pontos, variando amostradores, metodologia de amostragem dos frutos ( aleatória e não); 7. então tenho sd e mean; 9. tem um post da semana passada com um script do Walmes usando rnorm, que justamente é sobre o n a ser amostrado; 10. com relação a n coloco a quantidade total de plantas do talhão, não sei se poso fazer isso, ainda estou estudando; 11. rodei e verifiquei que ela gera dados com máximo e minino maiores do que já tenho, então to querendo que estes dados sejam mais próximos da realidade; tentei rodar o rtnorm mas ele deve pertencer a alguma pacote, pois o R dá o erro: Erro: não foi possível encontrar a função "rtnorm" estou pesquisando abraço Hélio Em 1 de abril de 2013 09:57, Luiz Roberto Martins Pinto < luizroberto.uesc(a)gmail.com> escreveu: > Ola Hélio, > > Interessante seu interesse pela geração de dados com distribuição normal > truncada. > Você pode compartilhar o motivo/aplicação? > Luiz Roberto > > Luiz Roberto Martins Pinto > Prof. Pleno/DCET/UESC > Laboratório de Estatística Computacional > Universidade Estadual de Santa Cruz > Ilhéus-Bahia > > luizroberto.uesc(a)gmail.com > skype: lrmpinto > http://lattes.cnpq.br/2732314327604831 > > > -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

1 0

duvidas na função factanal
by ynha.thays＠gmail.com 01 Abr '13

01 Abr '13

Eu estou com dúvidas na função factanal que está dando o erro de matriz singular. Por favor, alguém poderia me ajudar? Preciso dessa função para analisar os dados do meu Doutorado. Thays

1 0

limitando rnorm
by Hélio Gallo Rocha 01 Abr '13

01 Abr '13

Caros Teria como impor limites, máx e min para os valores gerados no uso de rnorm? Agradeço desde já -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 2

Re: [R-br] Digest R-br, volume 25, assunto 41
by Pamela Chiroque 01 Abr '13

01 Abr '13

Disculpem alguem sabe se existe a funcao pFq (Hypergeometrica?), ajuda por favor. Obrigada 2013/3/31 <r-br-request(a)listas.c3sl.ufpr.br> > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > > Tópicos de Hoje: > > 1. duvidas sobre o r (Juraci Mendes Moreira) > 2. Re: duvidas sobre o r (Daniel Marcelino) > 3. Re: duvidas sobre o r (Tiago Souza Marçal) > 4. Problemas com grid.arrange (Bernardo Rangel Tura) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sat, 30 Mar 2013 11:06:32 -0700 (PDT) > From: Juraci Mendes Moreira <juracimendesmoreira(a)yahoo.com.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] duvidas sobre o r > Message-ID: > <1364666792.97559.YahooMailClassic(a)web161005.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > como fazer analise de variancia em quadrado latino. qual o comando de anana > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130330/db11dc93/atta… > > > > ------------------------------ > > Message: 2 > Date: Sat, 30 Mar 2013 15:15:06 -0400 > From: Daniel Marcelino <dmarcelino(a)live.com> > To: R-br Lista <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] duvidas sobre o r > Message-ID: > <CAAjbdiK= > KN2rZWdymDiXUpHnavEfeJpDKQp0RJZ4S8i64ge88w(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Juraci, dê uma olhada nas funções ?aov, ?lme e ?lmer. > Além disso a sua pergunta não esclarece o que você precisa. O bom seria > você enviar um exemplo dos seus dados. O experimento que você está > realizando tem medidas repetidas ou não? > Por exemple, para analisar a variância em um quadrado latino simples, você > pode usar assim: > > out <- aov(time~participante+modo) > summary(out) > > > > 2013/3/30 Juraci Mendes Moreira <juracimendesmoreira(a)yahoo.com.br> > > > como fazer analise de variancia em quadrado latino. qual o comando de > anana > > _______________________________________________ > > R-br mailing list > > R-br(a)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. > > > > > > -- > \begin{signature} > Daniel Marcelino > ? dm.silva(a)umontreal.ca > ? (514) 343 6111 #3799 > Skype: d.marcelino > ? 3200 Jean Brillant, Office C5071 > Montreal, QC; H3T 1N8 > Canada > \end{signature} > > "Small steps toward a much better world" > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130330/1ae3897e/atta… > > > > ------------------------------ > > Message: 3 > Date: Sat, 30 Mar 2013 23:55:24 +0300 > From: Tiago Souza Marçal <tiagosouzamarcal(a)hotmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] duvidas sobre o r > Message-ID: <BLU170-W124F2B19F369AEAFA04AE7ED9DC0(a)phx.gbl> > Content-Type: text/plain; charset="utf-8" > > Para fazer uma anova em dql é só usar os seguntes comantos: > trat<-factor(trat) linha<-factor(linha) coluna<-factor(coluna) > av<-aov(resp~linha+coluna+trat) # Como pode ver a linha e a coluna são > somadas aos tratamentos. Caso você queira fazer uma analise completa de um > experimento dql eu te indico a usar a função dql()::ExpDes Att. Tiago. > Date: Sat, 30 Mar 2013 11:06:32 -0700 > From: juracimendesmoreira(a)yahoo.com.br > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] duvidas sobre o r > > como fazer analise de variancia em quadrado latino. qual o comando de anana > _______________________________________________ > R-br mailing list > R-br(a)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. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130330/41e0f417/atta… > > > > ------------------------------ > > Message: 4 > Date: Sun, 31 Mar 2013 09:02:37 -0300 > From: Bernardo Rangel Tura <tura(a)centroin.com.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Problemas com grid.arrange > Message-ID: <515825DD.70306(a)centroin.com.br> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > > Boa tarde, > > estou precisando fazer um gráfico chamado forest plot porem não posso > usar os já prontos no R. Já esta quase pronto utilizando este script. > > require(gridExtra) > nomes <- sample(letters,5) > logOR <- abs(rnorm(5,0,.3)) > sd_logOR <-rbeta(5,10,10) > OR <- round(exp(logOR),2) > lci_OR <- round(exp(logOR-1.96*sd_logOR),2) > uci_OR <- round(exp(logOR+1.96*sd_logOR),2) > OR <- paste(OR,' (',lci_OR,' - ',uci_OR,')',sep='') > df <- data.frame(nomes,OR,lci_OR,uci_OR) > p1 <- tableGrob(df[,1:2],cols=c('Estudos','OR (CI)'),row.just > ='left',col.just ='center', core.just ='right') > require(ggplot2) > OR <- round(exp(logOR),2) > df <- data.frame(nomes,OR,lci_OR,uci_OR) > p2 <- p2 <- ggplot(df, aes(x=nomes, y=OR, ymin=lci_OR, > ymax=uci_OR))+geom_hline(yintercept=1,lty=2)+expand_limits(y = > c(0,max(uci_OR)*1.1))+geom_pointrange()+coord_flip()#+ > geom_point(aes(size = sd_logOR)) > grid.arrange(p1,p2,ncol=2) > > Porém aqui no fim há uma desproporção entre os dados da tabela com os > dados do gráfico. > > Porventura alguém sabe como resolver isto ? > > []s > Tura > > > ------------------------------ > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > Fim da Digest R-br, volume 25, assunto 41 > ***************************************** >

3 2

Problemas com grid.arrange
by Bernardo Rangel Tura 31 Mar '13

31 Mar '13

Boa tarde, estou precisando fazer um gráfico chamado forest plot porem não posso usar os já prontos no R. Já esta quase pronto utilizando este script. require(gridExtra) nomes <- sample(letters,5) logOR <- abs(rnorm(5,0,.3)) sd_logOR <-rbeta(5,10,10) OR <- round(exp(logOR),2) lci_OR <- round(exp(logOR-1.96*sd_logOR),2) uci_OR <- round(exp(logOR+1.96*sd_logOR),2) OR <- paste(OR,' (',lci_OR,' - ',uci_OR,')',sep='') df <- data.frame(nomes,OR,lci_OR,uci_OR) p1 <- tableGrob(df[,1:2],cols=c('Estudos','OR (CI)'),row.just ='left',col.just ='center', core.just ='right') require(ggplot2) OR <- round(exp(logOR),2) df <- data.frame(nomes,OR,lci_OR,uci_OR) p2 <- p2 <- ggplot(df, aes(x=nomes, y=OR, ymin=lci_OR, ymax=uci_OR))+geom_hline(yintercept=1,lty=2)+expand_limits(y = c(0,max(uci_OR)*1.1))+geom_pointrange()+coord_flip()#+ geom_point(aes(size = sd_logOR)) grid.arrange(p1,p2,ncol=2) Porém aqui no fim há uma desproporção entre os dados da tabela com os dados do gráfico. Porventura alguém sabe como resolver isto ? []s Tura

1 0

duvidas sobre o r
by Juraci Mendes Moreira 30 Mar '13

30 Mar '13

como fazer analise de variancia em quadrado latino. qual o comando de anana

3 2

Fwd: Redes neurais
by Bernardo Silva 29 Mar '13

29 Mar '13

Boa tarde, estou a trabalhar em redes neurais utilizando o pacote nnet. Gostaria de saber qual o comando para encontrar o melhor modelo, por exemplo qd trabalhamos em glms utilizamos o stepwise, mas em redes neurais nao sei, podem ajudar? li que tb pode ser o stepwise mas depois nao consigo implementar obrigado

1 0

nomeando pontos
by Hélio Gallo Rocha 28 Mar '13

28 Mar '13

Caros, Busquei em msg anteriores sobre plotar nomes de atributos num mapa, na msg de Alexandro Yahoo parece interessante, até consegui baixar os shapes, mas no script faz referencia a arquivos do excel que ele não anexou, assim fica dificil reproduzir. Estou usando o pacote geoR, Alguém poderia dar uma força neste sentido de plotar o nome de um atributo, tratamento por exemplo, no mapa? grato -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

6 11

Considerar ângulo da área para distribuir coordenadas com expand.grid ()
by ASANTOS 28 Mar '13

28 Mar '13

Boa noite pessoal, Gostaria de saber se existe uma maneira de controlar a inclinação da distribuição dos pontos em expand.grid() em função de um limite de área, como por exemplo no meu CRM abaixo, onde no exemplo 1 esta tudo Ok, mas no exemplo 2 as coordenadas não consideram a inclinação da área. Seria possível usar algum argumento para controlar isto? # Exemplo 1 x<-c(0,0,100,100,0) y<-c(0,100,100,0,0) plot(x,y, type="l", col="red") # p0<-c(10,10) # # c0 <- seq(p0[1],p0[1]+95,by=20) d0 <- seq(p0[2],p0[2]+95,by=20) d1 <- expand.grid(x=c0, y=d0) points(d1[,1],d1[,2]) # #Exemplo 2 x<-c(0,0,50,50,0) y<-c(0,100,150,50,0) plot(x,y, type="l", col="red") # p0<-c(10,10) # # c0 <- seq(p0[1],p0[1]+95,by=20) d0 <- seq(p0[2],p0[2]+95,by=20) d1 <- expand.grid(x=c0, y=d0) points(d1[,1],d1[,2]) # Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br ======================================================================

2 2

Unique com menos variáveis
by Fátima Lima Paula 28 Mar '13

28 Mar '13

Prezados, tenho um banco de dados com muitas variáveis e com linhas idênticas em relação a algumas variáveis. Quero retirar as linhas idênticas, mas para que isso fique melhor é necessário fazer o unique apenas com algumas variáveis. Quando faço isso são eliminadas mais linhas. Mas, por outro lado o unique me retorna um banco com menos variáveis. Como posso fazer para ficar com o banco com todas as colunas e sem as linhas idênticas. Exemplo: nome=c("A","B","C","D","E","F","F","A","C","K","L","M") idade=c("1","2","3","4","5","6","6","1","3","11","12","13") cor=c("azul","preto","verde","branco","amarelo","cinza","cinza", "marrom", "amarelo","bege","azul","azul") material=c("lápis","ignorado","caneta","papel","borracha","lápis","lápis", "lápis","ignorado","caneta","apontador","pilha") banco1=cbind(nome,idade,cor,material) # banco que tenho dim(banco1) banco2=cbind(nome,idade) # banco para fazer o unique dim(banco2) b1unique=unique(banco1) # se faço o unique com o banco todo retornam 11 linhas e 4 colunas dim(b1unique) b2unique=unique(banco2) # fazendo o unique com o banco com as 2 variáveis retornam 9 linhas e 2 colunas dim(b2unique) Para fazer o unique, de forma a eliminar o maior número de repetidas, quero utilizar apenas as variáveis nome e idade. O b1unique (com todas as variáveis) retorna 11 linhas e o b2unique retorna 9. Quero o banco com as nove linhas, mas com todas as colunas do banco1. Verifiquei que existe um argumento incomparables no unique, mas vi que não é aplicável para bancos grandes. Alguém poderia me ajudar, por favor? Obrigada Fátima

2 6

Re: [R-br] Teste para decisão do número ótimo de amostra [RESOLVIDO]
by Hélio Gallo Rocha 27 Mar '13

27 Mar '13

Caros Alexandre e Walmes. Estou usando este teste de decisão de número de amostras para ver se minha amostragem está ok. Mas fiquei na dúvida sobre o valor 10000 da linha pop.total <- rnorm(10000, mean=20, sd=15) Vou exemplificar com dados de previsão de safra em café, no caso, são relativos a altura da planta de café. Estas medidas referem a coleta em 3 pontos da lavoura, cada ponto é composto por 5 plantas, então sd e mean são de 15 plantas. Estes 3 pontos amostrais estão dentro de uma área de 1.42 ha com 1775 plantas. sd=0.2587516 mean= 3.013333 assim acho que o valor 10000 seria 1775, ok? Posso determinar o numero de amostras para para uma área de 20 ha com 37085 plantas? Grato Hélio Em 27 de março de 2013 12:30, Hélio Gallo Rocha <heliogallorocha(a)gmail.com>escreveu: > Caros Alexandre e Walmes. > > Estou usando este teste de decisão de número de amostras para ver se minha > amostragem está ok. > > Mas fiquei na dúvida sobre o valor 10000 da linha > > pop.total <- rnorm(10000, mean=20, sd=15) > > > > Vou exemplificar com dados de previsão de safra em café, no caso, são > relativos a altura da planta de café. > Estas medidas referem a coleta em 3 pontos da lavoura, cada ponto é > composto por 5 plantas, então sd e mean são de 15 plantas. > Estes 3 pontos amostrais estão dentro de uma área de 1.42 ha com 1775 > plantas. > > sd=0.2587516 > mean= 3.013333 > > assim acho que o valor 10000 seria 1775, ok? > > Posso determinar o numero de amostras para para uma área de 20 ha > com 37085 plantas? > > Grato > > Hélio > > > > Em 20 de março de 2013 18:49, ASANTOS [via R-br] < > ml-node+s2285057n4658689h27(a)n4.nabble.com> escreveu: > > Perfeito Walmes, >> >> Era isto mesmo que eu buscava, >> >> Abraço, >> >> Alexandre >> >> >> Em 20/03/2013 17:31, Walmes Zeviani escreveu: >> >> Alexandre, >> >> Segue um CMR simples: >> >> # tamanho de amostra para medida de interesse que >> # é o erro absoluto ao estimar a média da população >> # erro = abs(média_estimada-média_verdadeira) >> >> pop.total <- rnorm(10000, mean=20, sd=15) # população >> m <- mean(20) # média verdadeira >> i <- unique(floor(2^seq(1, 12, by=0.5))); i # tamanhos de amostra >> nr <- 100 # número de execuções de >> reamostragem >> fun <- function(size, pop, m){ >> # função que calcula o erro a partir de uma amostra >> x <- sample(pop, size, replace=TRUE) >> # retorna a diferença entre estimado e real >> abs(mean(x)-m) >> } >> >> results1 <- sapply(i, >> function(j){ >> replicate(nr, fun(size=j, pop=pop.total, m=m)) >> }) >> str(results1) # matriz: colunas correspondem à i >> results2 <- colMeans(results1) # tira a média para cada i >> results2 >> >> plot(results2~i) >> plot(results2~log2(i)) >> >> # qual o tamanho de amostra para ter um erro no ponto >> # médio entre o máximo (quando n=1) e o mínimo (n->infinito, erro=0) >> >> n0 <- nls(results2~A*B^(-log2(i)+1), start=list(A=1, B=1)) >> coef(n0) >> >> perc <- 0.3 # erro tolerado de 30% do máximo >> (usuário define) >> xi <- -log(perc)/log(coef(n0)["B"])+1 # tamanho de amostra correspondente >> >> plot(results2~log2(i), >> xlim=c(1, max(log2(i))), ylim=c(0,max(results2, coef(n0)["A"])), >> xlab="Tamanho da amostra (log2 unidades)", >> ylab="Erro absoluto médio ao estimar a média") >> abline(v=1, lty=2) >> with(as.list(coef(n0)), curve(A*B^(-x+1), add=TRUE)) >> with(as.list(coef(n0)), abline(h=A*c(1,perc), v=xi, lty=2, col=2)) >> >> ceiling(2^xi) # tamanho de amostra à ser considerado >> >> À disposição. >> Walmes. >> >> >> ========================================================================== >> Walmes Marques Zeviani >> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) >> Departamento de Estatística - Universidade Federal do Paraná >> fone: (+55) 41 3361 3573 >> VoIP: (3361 3600) 1053 1173 >> e-mail: [hidden email]<http://user/SendEmail.jtp?type=node&node=4658689&i=0> >> skype: walmeszeviani >> twitter: @walmeszeviani >> homepage: http://www.leg.ufpr.br/~walmes >> linux user number: 531218 >> ========================================================================== >> >> >> >> _______________________________________________ >> R-br mailing list[hidden email] <http://user/SendEmail.jtp?type=node&node=4658689&i=1>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. >> >> >> -- >> ====================================================================== >> Alexandre dos Santos >> Proteção Florestal >> Coordenador do curso Técnico em Florestas >> Vice Coordenador do curso de Engenharia Florestal >> IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso >> Campus Cáceres >> Caixa Postal 244 >> Avenida dos Ramires, s/n >> Bairro: Distrito Industrial >> Cáceres - MT CEP: 78.200-000 >> Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)[hidden email] <http://user/SendEmail.jtp?type=node&node=4658689&i=2> >> [hidden email] <http://user/SendEmail.jtp?type=node&node=4658689&i=3> >> ====================================================================== >> >> >> _______________________________________________ >> R-br mailing list >> [hidden email] <http://user/SendEmail.jtp?type=node&node=4658689&i=4> >> 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. >> >> ------------------------------ >> If you reply to this email, your message will be added to the >> discussion below: >> >> http://r-br.2285057.n4.nabble.com/R-br-Teste-para-decisao-do-numero-otimo-d… >> To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…> >> . >> NAML<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_view…> >> > > > > -- > Hélio Gallo Rocha > IFSULDEMINAS - Câmpus Muzambinho > -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 2

Re: [R-br] text mining
by Vinicius Brito Rocha 27 Mar '13

27 Mar '13

tura , estou usando esse. parece bom. Obrigado. aliais se tiver material sobre o assunto poderia me mandar? abs Vinicius Enviado por Samsung Mobile -------- Mensagem original -------- De : Bernardo Rangel Tura <tura(a)centroin.com.br> Data: Para: r-br(a)listas.c3sl.ufpr.br Assunto: Re: [R-br] text mining On 03/17/2013 03:49 PM, Vinicius Brito Rocha wrote: > Senhores, > > estou começando a estudar text mining. Quais pacotes, os senhores > recomendam? > > -- > /Vinicius Brito Rocha./ Vinicius, No R tem um pacote chamado tm que uso com bons resultados. Existem vários outros que você pode ver em http://cran.r-project.org/web/views/NaturalLanguageProcessing.html []s Tura _______________________________________________ R-br mailing list R-br(a)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.

2 1

Re: [R-br] PCA ou RDA?
by Marcos Vital 27 Mar '13

27 Mar '13

Olá, Cézar A sugestão do Luis Paulo, de usar uma CCA pode ser uma boa saída. Mas não sei se seria uma saída muito segura: dê uma olhada nas críticas feitas à CCA nos livros 'Numerical Ecology' ou no 'Numerical Ecology with R'. Os autores vão chamar a atenção para os pressupostos da análise, que raramente são atendidos (especialmente o de multinormalidade), o que pode comprometer o resultado. Por outro lado, eu sou da opinião de que você não necessariamente deveria ir pelo caminho das ordenações. Se você tem hipóteses claras sobre o efeito das variáveis explicativas (as características das ilhas) sobre a riqueza de espécie, então seria muito mais interessante partir para análises de significância! Você pode fazer regressões múltiplas para cada grupo funcional, por exemplo. Considero que esta saída seria mais interessante pelo fato de você poder ter em mãos não só a significância do efeito (que a CCA te daria também) mas uma medida do tamanho do efeito. E em um contexto ecológico como o seu, a medida do efeito das variáveis explicativas sobre a riqueza é infinitamente mais importante do que uma simples medida de sua significância. É claro que a significância estatística é fundamental, mas sozinha ela dirá muito pouco sobre o fenômeno que você está estudando. Atenciosamente Marcos V C Vital Message: 8 > Date: Wed, 27 Mar 2013 05:10:12 +0300 > From: Cézar Borges <cezar-fonseca(a)hotmail.com> > To: Lista_R <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] PCA ou RDA? > Message-ID: <COL126-W40F5F8DECC86B6190A749AFED10(a)phx.gbl> > Content-Type: text/plain; charset="iso-8859-1" > > > Olá lista,estou com um banco de dados com 9 variáveis, divididos em > atributos espacias de ilhas (área, forma e distância) e número de espécies > de grupos funcionais vegetais (árvores, arbustos, lianas...). Neste > sentindo, gostaria de realizar uma análise de ordenação, mas estou em > dúvida de qual seria mais apropriada, PCA ou RDA, uma vez que, quero > verificar o efeito das características espacias das ilhas sobre a > vegetação. Qual pacote e comando posso estar usando para essas > análises???Abraços! > > > Att, > > > > Cézar Fonseca Borges > > > > > > > > Mestrando em Ciências Ambientais > > UFPA/MPEG/EMBRAPA - > Amazônia Oriental > > Geógrafo > > Email: cezar-fonseca(a)hotmail.com/ caborges(a)museu-goeldi.br > > Lattes: http://lattes.cnpq.br/9805721240137485 > > > > -- Marcos Vinícius Carneiro Vital Universidade Federal de Alagoas Instituto de Ciências Biológicas e da Saúde Setor de Biodiversidade e Ecologia

1 0

PCA ou RDA?
by Cézar Borges 26 Mar '13

26 Mar '13

Olá lista,estou com um banco de dados com 9 variáveis, divididos em atributos espacias de ilhas (área, forma e distância) e número de espécies de grupos funcionais vegetais (árvores, arbustos, lianas...). Neste sentindo, gostaria de realizar uma análise de ordenação, mas estou em dúvida de qual seria mais apropriada, PCA ou RDA, uma vez que, quero verificar o efeito das características espacias das ilhas sobre a vegetação. Qual pacote e comando posso estar usando para essas análises???Abraços! Att, Cézar Fonseca Borges Mestrando em Ciências Ambientais UFPA/MPEG/EMBRAPA - Amazônia Oriental Geógrafo Email: cezar-fonseca(a)hotmail.com/ caborges(a)museu-goeldi.br Lattes: http://lattes.cnpq.br/9805721240137485

2 1

Listar variáveis relevantes
by Aline Lipsky 26 Mar '13

26 Mar '13

Boa tarde. Eu tenho um data.frame com 6000 observações e 55 variáveis. Preciso saber quais variáveis diferentes de zero ocorrem para cada observação. Se alguém puder me ajudar. Agradeço desde já a atenção. Aline F. Lipsky

2 4

Interpretação RDA-Parcial
by Luciano F. Sgarbi 26 Mar '13

26 Mar '13

Boa tarde, Estou com um problema na interpretação (ou é um erro) dos resultados da RDA-Parcial usando a função "varpart" {vegan}. É possível ter valores de R² Ajustado menores que 0 (negativos) ou maiores que 1 (neste caso 1.986 para a porção não explicada)? Qualquer ajuda será muito bem vinda. Alguém conhece outra função/pacote que faça o mesmo O CMR segue abaixo: (y<-read.csv("http://ge.tt/api/1/files/51HX3Dc/0/blob?download",h=T)) (y.phys<-y[,1:5]) (y.chem<-y[,6:10]) (y.comu<-y[,11:ncol(y)]) require(vegan) varpart(y.comu,y.phys,y.chem,transfo="hell") Abraços, Luciano -- *Luciano F. Sgarbi* Mestrando em Ecologia e Evolução Laboratório de Ecologia de Insetos (sl. 222) Departamento de Ecologia Instituto de Ciências Biológicas - ICB 1 Universidade Federal de Goiás, campus II Goiânia-GO 74001-970 Brazil

1 0

TwinSpan
by Robson 26 Mar '13

26 Mar '13

1 0

Ajuste em multcompBoxplot
by Marcelo Laia 26 Mar '13

26 Mar '13

Estou a usar o codigo abaixo > xzx <-multcompBoxplot(DESC_COMP ~ CLONETRAT,data = data.cerato, + sortFn = "median", decreasing=TRUE, + horizontal=FALSE, compFn = "TukeyHSD", + plotList=list( + boxplot=list(fig=c(0, 0.75, 0, 1), + las=3, cex.axis=1.0), + multcompLetters=list( + fig=c(0.87, 0.97, 0.03, 0.98), + type='Letters') )) > para gerar um grafico com multiplas comparacoes. O grafico é gerado, mas com as letras fora de alinhamento com os nomes no eixo x. Voce saberia me dizer onde eu tenho que ajustar para alinha-los? Veja o "a" está desalinhado com o primeiro item do eixo x. O mesmo acontece com o ultimo. Já tentei algumas configuracoes par, mas, esse boxplot parece ser especial e nao aceita alguns parametros normalmente passados para a funcao plot(). Eis o grafico gerado: https://dl.dropbox.com/u/34009642/boxplot_DESC_S_interacao.jpg -- Marcelo Brasil (Brazil, for English Speakers) Linux user number 487797

2 1

Medidas repetidas no tempo: aov, Anova, lme ou multcomp?
by Marcelo Laia 26 Mar '13

26 Mar '13

Colega, A dias que venho lendo a respeito de como realizar uma análise de medidas repetidas no tempo. Encontrei uma farta bibliografia e diversas sugestões de como fazer. Eis um dado mínimo: datafilename="http://dl.dropbox.com/u/34009642/Dados_Omissao.csv" data.min <- read.table(datafilename, header=T, sep="\t", dec=",") head(data.min) data.min <- within(data.min, { Clone <- factor(Clone) Tempo <- factor(Tempo) Trat <- factor(Trat) Planta <- factor(Planta) }) summary(data.min) A partir desses dados e do que eu li, tentei: dados_omissao_h.aov <- aov(h ~ Clone * Trat * Tempo + Error(Planta), data = data_min) # http://tinyurl.com/bjlas86 Esse http://pealco.net/2009/01/30/repeated-measures-anova-in-r.html dá outra abordagem. Eu também rodei o código abaixo. library(nlme) lme(h ~ Clone+Trat+Tempo,random=~1|Planta,data=data.min) # Ainda li os textos seguintes: http://egret.psychol.cam.ac.uk/statistics/R/anova.html#8.4 http://egret.psychol.cam.ac.uk/statistics/R/enteringdata.html#reshape http://www.psych.yorku.ca/cribbie/6130/ex14_rscript.R No entanto, sinceramente, não sei qual abordagem utilizar nos nossos dados. Desejamos verificar se há diferenças entre clones, entre tratamentos, a interação clone tratamento e a evolução ao longo do tempo das variáveis. Ainda, há plantas que morreram ao longo do experimento devido ao próprio tratamento. E ainda tem a dica do Wlames no ridiculas para dados binomial! Alguma sugestão? -- Marcelo Brasil (Brazil, for English Speakers) Linux user number 487797

4 5

Roy e Bose
by Ze Henrique 24 Mar '13

24 Mar '13

Boa noite a todos. Gostaria de saber se alguém conhece alguma função e o pacote que faz a análise de Roy e Bose (estatística multivariada) Desde já obrigado. -- *José Henrique Soler Guilhen* Graduando em Agronomia Universidade Federal do Espírito Santo - CCAUFES

1 0

erro
by Pablo . 24 Mar '13

24 Mar '13

o que significa esse erro: Erro em next[1] = -1 : lado esquerdo de atribuição inválido (NULL).. já tentei várias coisas e não muda esse erro...

2 1

Gráfico Biplot via psych para Analise de fatores!
by Leandro Peixouto 24 Mar '13

24 Mar '13

Bom dia Colegas R's Estou com dificuldade em separar os gráficos oriundos de uma analise de fatores com mais de 3 eixos via função "biplot". Ao realizar a análise gráfica o R apresenta a janela com todos os gráficos dos fatores combinados dois a dois e na diagonal outros gráficos que não me interessam. O que realmente preciso é plotar os gráficos PC1 vsPC2, PC1 vsPC3 e PC2 vsPC3. Como sempre vejo aqui, segue abaixo um CMR. Library(psych) data(bfi) fa3 <- fa(bfi[1:200,1:15], 3, fm="ml", scores=TRUE, rotate="varimax") biplot(fa3) Na maioria dos exemplo encontrado sempre o fazem com dois eixos. Só que minha análise demonstrou ser necessario 5 fatores. Desde já agradeço a ajuda! Leandro Santos Peixouto Prof. IFBaiano Campus Guanambi-BA Doutorando em Genética e Mlhoramento de Plantas -UFLA

1 0

Exercicios de Regressão Linear em R.
by Bruno Monteiro Dantas 24 Mar '13

24 Mar '13

Pessoal, Estou fazendo uma cadeira de estatistica na faculdade e gostaria de entender mais da linguagem R com questões relacionadas a regressão... Alguem poderia me recomendar algum livro com exercicios resolvidos para que eu possa resolve-los e logo em seguida verificar se estão corretos? Acho pessimo resolver uma questão e não saber se está certa ou não. Peço isso, pois terei que fazer uma prova em laboratorio com questões de regressão. Diante disto gostaria de me familiarizar logo com a resolução de questões deste tipo. Desde Já, muito Obrigado. Bruno Monteiro Dantas, Estudante de Ciencias da Computação - UFRN.

4 3

Re: [R-br] nomeando pontos
by Hélio Gallo Rocha 23 Mar '13

23 Mar '13

Valeu Alexandro Verei mais tarde, depois te falo grato Hélio Em 23 de março de 2013 12:57, Alexandro (Yahoo) <vl.alexandro(a)yahoo.com.br>escreveu: > > OI Hélio, segue anexo tudo bonitinho > Instale as bibliotecas do arquivo > plotar_mapa_R > e execute a rotina, depende de cada shape file > a ordem correta do arquivo do excel para as respectivas cidades! > > Espero ter ajudado! > Alexandro > > > ------------------------------ > *De:* Hélio Gallo Rocha <heliogallorocha(a)gmail.com> > *Para:* r-br <r-br(a)listas.c3sl.ufpr.br> > *Enviadas:* Sábado, 23 de Março de 2013 12:16 > *Assunto:* [R-br] nomeando pontos > > Caros, > > Busquei em msg anteriores sobre plotar nomes de atributos num mapa, na msg > de > Alexandro Yahoo parece interessante, até consegui baixar os shapes, mas no > script faz referencia a arquivos do excel que ele não anexou, assim fica > dificil reproduzir. > > Estou usando o pacote geoR, > Alguém poderia dar uma força neste sentido de plotar o nome de um > atributo, tratamento por exemplo, no mapa? > > grato > -- > Hélio Gallo Rocha > IFSULDEMINAS - Câmpus Muzambinho > > _______________________________________________ > R-br mailing list > R-br(a)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. > > -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

1 0

Falha na instalação de pacotes
by Jobenil Júnior 23 Mar '13

23 Mar '13

Prezados Sempre que posso testo os CMRs postados no grupo. Esta é uma forma que encontrei para aprender mais sobre o R. Estou testando um script que requer o ggmap mas ocorre algum erro que não fui capaz de decifrar. O Sistema Operacional que estou usando é o Ubuntu 12.04 LTS 32-bits. A versão do R instalada é a 2.15.2. Tentei instalar os pacotes através do RStudio e também diretamente pelo terminal. Alguém pode me dar uma ajuda para resolver este problema. Agradeço desde já, Atenciosamente Jobenil Júnior > install.packages('ggmap') Installing package(s) into ‘/home/marcia/R/i486-pc-linux-gnu-library/2.15’ (as ‘lib’ is unspecified) --- Please select a CRAN mirror for use in this session --- Carregando a interface Tcl/Tk ... pronto also installing the dependencies ‘RgoogleMaps’, ‘png’ tentando a URL 'http://cran-r.c3sl.ufpr.br/src/contrib/RgoogleMaps_1.2.0.2.tar.gz' Content type 'application/x-gzip' length 471396 bytes (460 Kb) URL aberta ================================================== downloaded 460 Kb tentando a URL 'http://cran-r.c3sl.ufpr.br/src/contrib/png_0.1-4.tar.gz' Content type 'application/x-gzip' length 21774 bytes (21 Kb) URL aberta ================================================== downloaded 21 Kb tentando a URL 'http://cran-r.c3sl.ufpr.br/src/contrib/ggmap_2.2.tar.gz' Content type 'application/x-gzip' length 2920315 bytes (2.8 Mb) URL aberta ================================================== downloaded 2.8 Mb * installing *source* package ‘png’ ... ** package ‘png’ successfully unpacked and MD5 sums checked ** libs gcc -std=gnu99 -I/usr/share/R/include -DNDEBUG `libpng-config --cflags` -fpic -O3 -pipe -g -c read.c -o read.o /bin/bash: libpng-config: comando não encontrado read.c:2:17: erro fatal: png.h: Arquivo ou diretório não encontrado compilação terminada. make: ** [read.o] Erro 1 ERROR: compilation failed for package ‘png’ * removing ‘/home/marcia/R/i486-pc-linux-gnu-library/2.15/png’ ERROR: dependency ‘png’ is not available for package ‘RgoogleMaps’ * removing ‘/home/marcia/R/i486-pc-linux-gnu-library/2.15/RgoogleMaps’ ERROR: dependencies ‘RgoogleMaps’, ‘png’ are not available for package ‘ggmap’ * removing ‘/home/marcia/R/i486-pc-linux-gnu-library/2.15/ggmap’ The downloaded source packages are in ‘/tmp/RtmpDW0XHB/downloaded_packages’ Mensagens de aviso perdidas: 1: In install.packages("ggmap") : installation of package ‘png’ had non-zero exit status 2: In install.packages("ggmap") : installation of package ‘RgoogleMaps’ had non-zero exit status 3: In install.packages("ggmap") : installation of package ‘ggmap’ had non-zero exit status

4 4

Entrada de dados: GNU/Linux
by Vinícius Lionel Mateus 22 Mar '13

22 Mar '13

Prezados useRs, Alguém poderia me ajudar com uma dica ou bibliografia sobre a entrada de dados no R para GNU/Linux? Cheguei a tentar ler a clipboard, mas não fui feliz. Minha distro é Ubuntu 12.04 LTS (64 bits). Vinícius -- -- Atenciosamente, VINÍCIUS LIONEL MATEUS, M.Sc (http://lattes.cnpq.br/6501001637020665) Bacharel em Química - Doutorando em Química Analítica Laboratório de Química Atmosférica - Departamento de Química Pontifícia Universidade Católica - Rio de Janeiro (PUC - Rio) Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - Brasil CEP.: 22451-900 Telefone: (+55) (21) 3527-1327 (+55) (21) 9358-8051 www.puc-rio.br

6 11

GLM ou GAM com dados de proporção
by Pavel Dodonov 22 Mar '13

22 Mar '13

Saudações, Eu estou tentando rodar um modelo aditivo (mas pode ser GLM também) com dados de proporção (variando de 0 a 1). Pensei em usar distribuição binomial, mas parece que ela só funciona com dados binários (ao menos no pacote mgcv, que estou usando). Alguma sugestão? Muito obrigado! - Pavel ------- Pavel Dodonov Biologist, PhD student in Ecology and Natural Resources, São Carlos Federal University (UFSCar), SP, Brazil http://ufscar-br.academia.edu/pdodonov *"The highest function of ecology is understanding consequences." - Frank Herbert, Dune, 1965*

3 2

aRT: atualização do BD
by Gilbert Queiroz 21 Mar '13

21 Mar '13

Olá! Gostaria de saber como eu posso atualizar uma coluna de um BD aberto via aRT ou mesmo criar uma???? Gilbert.

3 2

Busca em texto R
by Pablo . 21 Mar '13

21 Mar '13

Olá, gostaria de saber se alguém sabe como manipular um texto (.txt) comum, para o R, transformá-lo em string para que eu possa buscar palavras nesse texto.

2 1

disposição do dados para usar o INLA
by Gilbert Queiroz 21 Mar '13

21 Mar '13

Qual seria a disposição dos dados para usarmos o INLA???? No caso, estou trabalhando com dados com distribuição espaço-temporais, e ainda algumas covariáveis. Abs. Gilbert.

3 2

(sem assunto)
by Tiago Souza Marçal 21 Mar '13

21 Mar '13

Boa noite pessoal, gostaria de saber como substituir o sinal de vezes por um espaço? CMR plot(0, main=expression(KNO[3]%*%L[-1])) Att. Tiago.

5 7

[1/2 off] Funções multiplicativas em Estatística
by Listeiro 037 20 Mar '13

20 Mar '13

Bom dia. Peço licença aos presentes para esclarecer uma dúvida: aparecem funções com propriedades multiplicativas em Estatística, sejam de variável discreta ou real? Inclusive convolução em alguma teoria, se aparecer. Lembrando que para uma função f(x) ser multiplicativa, deve haver f(x*y) = f(x)*f(y). Grato a todos.

4 8

Teste para decisão do número ótimo de amostra
by ASANTOS 20 Mar '13

20 Mar '13

Boa tarde pessoal, Gostaria de utilizar um teste para decidir o número ótimo de amostral, baseado na abordagem que estou empregando, se faço: pop.total<-rnorm(10000,mean=20,sd=15) ## População total i<-c(5000,50)## Quantidade de amostras n.sim<-999 ## Número de simulações RES=NULL for(Nsim in 1:n.sim){ for(Nam in 1:length(i)){ amostras<-sample(pop.total,i[Nam]) ## Amostragem com reposição m.amostras<-mean(amostras) ## Media encontrada nas amostras err.amost<-100-((mean(amostras)/mean(pop.total))*100) ### Erro amostral RES=rbind(RES,c(Nsim,i[Nam],err.amost)) }} ####### 50 amostras am.50<-RES[RES[,2]==50,] shapiro.test(am.50[,3]-mean(am.50[,3])) ## Distribuição normal,então inferências sobre os parâmetros média e desvio estão corretos mean(am.50[,3])## Média do erro amostral sd(am.50[,3])/sqrt(length(am.50[,3]))## Erro padrão do erro amostral sd(am.50[,3])## Desvio padrão # ####### 5000 amostras am.5000<-RES[RES[,2]==5000,] shapiro.test(am.5000[,3]-mean(am.5000[,3])) ## Distribuição normal,então inferências sobre os parâmetros média e desvio estão corretos mean(am.5000[,3])## Média do erro amostral sd(am.5000[,3])/sqrt(length(am.50[,3]))## Erro padrão do erro amostral sd(am.5000[,3])## Desvio padrão # Percebo que a média do erro amostral e o erro padrão desta média são baixos (próximos de zero) tanto quando minha amostra tem tamanho 50 quanto 5000, porém o desvio padrão, obviamente, aumenta quando reduzo o número de amostras. Gostaria de saber que teste poderia empregar para escolha de um número ótimo de amostras, uma vez que eu não posso usar como critério o erro amostral, quando simulo este cenário 999 vezes e emprego inúmeros tamanhos amostrais? Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br ======================================================================

2 4

Como listar arquivos faltantes usando R
by Thiago V. dos Santos 20 Mar '13

20 Mar '13

Olá pessoal, Baixei centenas de imagens de satélite e agora devo conferir se tenho todas que realmente vou precisar. A nomenclatura das imagens segue o padrão: MCD12Q1.A2001.h10v08.hdf MCD12Q1.A2001.h10v09.hdf MCD12Q1.A2001.h10v10.hdf MCD12Q1.A2001.h10v11.hdf MCD12Q1.A2001.h11v08.hdf MCD12Q1.A2001.h11v09.hdf MCD12Q1.A2001.h11v10.hdf MCD12Q1.A2001.h11v11.hdf MCD12Q1.A2001.h11v12.hdf MCD12Q1.A2001.h11v13.hdf e assim por diante onde 2001 é o ano e vai de 2001 a 2011. h10 é um componente horizontal da imagem e vai de a 10 a 14. v08 é um componente vertical da imagem e vai de a 08 a 13. O que preciso fazer é o seguinte: varrer uma lista de todas as imagens possíveis, que vão de 2001 a 2011, h10 a h14 e v08 a v13, e testar se tem alguma que está faltando. Para isso eu uso a função "file.exists" dentro de três loops: "ano", "h" e "v". No final da mensagem estou colando o código que usei para fazer isso. Ele até funciona, mas não salva a lista inteira, e sim somente o último nome. Estou fazendo alguma coisa errada no loop. Acredito também que fazer três loops não é a maneira mais otimizada de fazer isso, porém foi a primeira que veio na minha cabeça. Alguém pode ajudar, criticando o código? Obrigado desde já, Thiago. ################################################################## # Gather all files dir <- ('/Users/dossa013/data/MODIS/MCD12Q1/original') files <- list.files(dir, pattern='hdf', full.names=T) # Set beginning and ending years byear <- 2001 eyear <- 2011 # Set tiles h1 <- 10 h2 <- 14 v1 <- 08 v2 <- 13 # Loop over year and then over tiles for (x in byear:eyear) { for (y in h1:h2) { for (z in v1:v2) { name <- paste('/Users/dossa013/data/MODIS/MCD12Q1/original/MCD12Q1.A',x,'.h',y,'v',z,'.hdf', sep='') if (file.exists(name)==F) missing <- name } # end of v } # end of h } # end of year > missing [1] "/Users/dossa013/data/MODIS/MCD12Q1/original/MCD12Q1.A2011.h14v9.hdf"

4 7

Nome das linhas em um data frame
by Luciano F. Sgarbi 19 Mar '13

19 Mar '13

Prezados, Estou com um pequeno problema em importar um conjunto de dados para o R. Eu estou criando um conjunto de dados e salvando ele no formato ".csv" Quando eu importo o conjunto de dados para o R, o nome das linhas se transforma em uma nova variável. No entanto quando eu salvo o arquivo em formato ".txt" esse erro não ocorre. Acredito que a resolução deste problema seja muito simples, no entanto não estou conseguindo encontrar essa tal solução. Se puderem me ajudar ficarei muito grato! O CMR segue abaixo: my.data<-data.frame(v1=c(letters[1:5]),v2=1:5,v3=100:104) rownames(my.data)<-LETTERS[1:5] my.data write.csv(my.data,"my_data.csv") dados<-read.csv("my_data.csv",h=T) dados Abraços, Luciano -- *Luciano F. Sgarbi* Mestrando em Ecologia e Evolução Laboratório de Ecologia de Insetos (sl. 222) Departamento de Ecologia Instituto de Ciências Biológicas - ICB 1 Universidade Federal de Goiás, campus II Goiânia-GO 74001-970 Brazil

6 11

Erro ao rodar função
by Ari Clecius 18 Mar '13

18 Mar '13

Boa tarde, rodei o código abaixo, quando tentei gerar o gráfico de cubo pela função cubePlot gerou o erro abaixo. Alguém poderia me ajudar? Att, Ari Clecius Alves de Lima Engenheiro Químico Me. Engenharia Civil (085)87687786 (085)33669042 library(qualityTools) exp1 = facDesign( k = 3 , centerCube = 0, replicates=3) names ( exp1 ) = c ( " Temperatura " , " Concentração" , "Pressão " ) units(exp1)=c("°C", "mg/L", "mmHg") lows ( exp1 ) = c (60, 1,720) highs ( exp1 ) = c (80,2,760) summary( exp1 ) y=c(1.60, 1.85, 1.40, 1.76, 1.70, 1.54, 0.37, 0.40, 0.46, 1.70, 1.80, 1.90, 0.30, 0.35, 1.74, 1.55, 1.80, 1.75, 1.50, 1.55, 0.60, 1.80, 0.55, 1.52) response(exp1)=data.frame(y) summary( exp1 ) Erro em `contrasts<-`(`*tmp*`, value = contr.funs[1 + isOF[nn]]) : contrasts can be applied only to factors with 2 or more levels

2 1

spgam????
by Gilbert Queiroz 18 Mar '13

18 Mar '13

Vi na página do LEG um link sobre esse assunto...Algo mais pode ser dito???? Abs. Gilbert.

2 1

text mining
by Vinicius Brito Rocha 18 Mar '13

18 Mar '13

Senhores, estou começando a estudar text mining. Quais pacotes, os senhores recomendam? -- *Vinicius Brito Rocha.* *Estatístico e Atuário (IM / UFRJ)** Mestre em Pesquisa Operacional (COPPE / UFRJ)* www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

4 3

SoundexBR
by Fátima Lima Paula 17 Mar '13

17 Mar '13

Olá Daniel Marcelino, precisaria de uma ajuda sua no soundex. Verifiquei que erros de digitação não são captados pelo SoundexBR. Exemplo: ?MARIA e MARIA, PE8DRO e PEDRO. Será que seria tranquilo para você implementar algumas coisas? Estou enviando um script que fiz com teste para alguns caracteres. Desde já muito obrigada. Abs Fátima soundexBR<-function(termo){ termo<-toupper(termo) # 1. Retire toda pontuação da palavra; termo<-gsub("[ÃÁÀÂÄ]","A",termo) termo<-gsub("[ÉÈÊ?Ë]","E",termo) termo<-gsub("[ÍÌÏÎI]","I",termo) termo<-gsub("[ÓÕÒÔÖ]","O",termo) termo<-gsub("[ÚÙÛÜU]","U",termo) termo<-gsub("Ç","C",termo) #1.a Substituir todas as letras duplas por uma única letra #http://www.archives.gov/genealogy/census/soundex.html termo<-gsub("([A-Z])\\1", "\\1", termo) # 2. Fique com a primeira letra; N<-nchar(termo) # Modificação PT-BR (Coeli e Camargo) termo<-ifelse(substr(termo,1,2)=="WA",sub("W","V",termo),termo) termo<-ifelse(substr(termo,1,1)=="H",substr(termo,2,N),termo) termo<-ifelse(substr(termo,1,2)=="KA"|substr(termo,1,2)=="KO"|substr(termo,1,2)=="KU",sub("K","C",termo),termo) termo<-ifelse(substr(termo,1,1)=="Y",sub("Y","I",termo),termo) termo<-ifelse(substr(termo,1,2)=="CE"|substr(termo,1,2)=="CI",sub("C","S",termo),termo) termo<-ifelse(substr(termo,1,2)=="GE"|substr(termo,1,2)=="GI",sub("G","J",termo),termo) ############## Segunda Parte ####################### termo.1<-substr(termo,1,1) termo<-substr(termo,2,N) # 3. Mude todas as ocorrências das letras a seguir por '0' (zero): 'A','E','I','O','U','H','W','Y'; termo<-gsub("[A,E,I,O,U,H,W,Y]",0,termo) # 4. Mude as letras restantes de acordo com a tabela abaixo: # Número Letra # 1 'B','F','P','V' termo<-gsub("[B,F,P,V]",1,termo) # 2 'C','G','J','K', 'Q','S','X','Z' termo<-gsub("[C,G,J,K,Q,S,X,Z]",2,termo) # 3 'D','T' termo<-gsub("[D,T]",3,termo) # 4 'L' termo<-gsub("L",4,termo) # 5 'M','N' termo<-gsub("[M,N]",5,termo) # 6 'R' termo<-gsub("R",6,termo) # 5. Remova todos os zeros da string resultante; termo<-gsub(0,"",termo) # Eliminar dois numeros iguais consecutivos por um único numero #http://www.archives.gov/genealogy/census/soundex.html termo<-gsub("([0-9])\\1", "\\1", termo) # Remontar termo<-paste(termo.1,termo,sep="") # 6. Preencha a string resultante com zeros à direita e retorne desta forma: #<letra maiuscula><digito><digito><digito> termo<-paste(termo,"0000",sep="") termo<-substr(termo,1,4) return(termo) } soundexBR("termo") soundexBR("terno") nomea=data.frame(nome=c("GER?D","HARTMUT","RALF","STEFAN","ROBERT")) nomea nomeb=data.frame(nome=c("HARTMUT","GERD","ROBERT","STEFAN","RALF")) nomeb fonoa=soundexBR(nomea$nome) fonob=soundexBR(nomeb$nome) library(RecordLinkage) pairs.mat = outer(fonoa,fonob, "==") pairs.mat #Ponto de interrogação nomea=data.frame(nome=c("GER?D","HARTMUT","?RALF","STEFAN","ROBERT")) nomea nomeb=data.frame(nome=c("HARTMUT?","GERD","ROBERT","STEFAN","RALF")) nomeb fonoa=soundexBR(nomea$nome) fonob=soundexBR(nomeb$nome) pairs.mat = outer(fonoa,fonob, "==") pairs.mat #Não funciona com o ? no início ############################################################################### #Barras inclinadas / nomea=data.frame(nome=c("GER/D","/HARTMUT","RALF/","STEFAN","ROBERT")) nomea nomeb=data.frame(nome=c("HARTMUT","GERD","ROBERT","STEFAN","RALF")) nomeb fonoa=soundexBR(nomea$nome) fonob=soundexBR(nomeb$nome) pairs.mat = outer(fonoa,fonob, "==") pairs.mat #Não funciona com a barra inclinada #Barras inclinadas \ nomea=data.frame(nome=c("GER\D","\HARTMUT","RALF\","STEFAN","ROBERT")) nomea #O R rejeita a \ #Parênteses ( nomea=data.frame(nome=c("GER(D","(HARTMUT","RALF(","STEFA)N","ROBERT)")) nomea nomeb=data.frame(nome=c("HARTMUT","GERD","ROBERT","STEFAN","RALF")) nomeb fonoa=soundexBR(nomea$nome) fonob=soundexBR(nomeb$nome) pairs.mat = outer(fonoa,fonob, "==") pairs.mat #Só funciona com o parênteses depois do nome #Dois pontos e ponto e vírgula nomea=data.frame(nome=c("GER:D",":HARTMUT","RALF:","STEFA;N","ROBERT;")) nomea nomeb=data.frame(nome=c("HARTMUT","GERD","ROBERT","STEFAN","RALF")) fonob=soundexBR(nomeb$nome) pairs.mat = outer(fonoa,fonob, "==") pairs.mat #Só funciona com o : e ; depois do nome #Dólar e ponto nomea=data.frame(nome=c("GER$D","%HARTMUT","RALF%","STEFA.N","ROBERT.")) nomeb=data.frame(nome=c("HARTMUT","GERD","ROBERT","STEFAN","RALF")) nomeb fonoa=soundexBR(nomea$nome) fonob=soundexBR(nomeb$nome) pairs.mat = outer(fonoa,fonob, "==") pairs.mat #Idem #Números nomea=data.frame(nome=c("GER1D","2HARTMUT","RALF3","S4TEFA.N","RO5BERT6")) nomea nomeb=data.frame(nome=c("HARTMUT","GERD","ROBERT","STEFAN","RALF")) nomeb fonoa=soundexBR(nomea$nome) fonob=soundexBR(nomeb$nome) pairs.mat = outer(fonoa,fonob, "==") pairs.mat #Não funciona com números

2 4

Problema para plotar linha
by ASANTOS 17 Mar '13

17 Mar '13

Boa tarde pessoal, Estou tentando plotar uma linha continua ligando 11 pontos, mas quando uso o comando type="l" da função plot() ele faz um zig zag entre os pontos e na realidade eu quero ligá-los, sendo: # x2<-c(69.19,71.13,63.38,69.19,68.60,68.21,63.17,71.13,67.26,63.38,71.90) y0<-c(39.79699,43.82621,29.81363,39.79699,38.64670,37.90465,29.50401,43.82621,36.15617,29.81363,45.53631) plot(x2,y0, type="l",main="") points(x2,y0) # Então eu pergunto, se existe alguma maneira de forçar a ligação entre os pontos, sem considerar a ordem deles? Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br ======================================================================

3 4

Re: [R-br] Digest R-br, volume 25, assunto 23
by Márcio Mendonça 16 Mar '13

16 Mar '13

Prezados, como faço para inserir os eixos(tipo vetores) em um gráfico de ordenação, por exemplo em nmds? existe algum pacote? ATT Márcio ----- Mensagem original ----- De: "r-br-request(a)listas.c3sl.ufpr.br" <r-br-request(a)listas.c3sl.ufpr.br> Para: r-br(a)listas.c3sl.ufpr.br Cc: Enviadas: Sexta-feira, 15 de Março de 2013 12:00 Assunto: Digest R-br, volume 25, assunto 23 Enviar submissões para a lista de discussão R-br para r-br(a)listas.c3sl.ufpr.br Para se cadastrar ou descadastrar via WWW, visite o endereço https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ou, via email, envie uma mensagem com a palavra 'help' no assunto ou corpo da mensagem para r-br-request(a)listas.c3sl.ufpr.br Você poderá entrar em contato com a pessoa que gerencia a lista pelo endereço r-br-owner(a)listas.c3sl.ufpr.br Quando responder, por favor edite sua linha Assunto assim ela será mais específica que "Re: Contents of R-br digest..." Tópicos de Hoje: 1. Re: Eliminar Vetores Repetidos.... (Benilton Carvalho) 2. Re: Eliminar Vetores Repetidos.... (andrebvs(a)bol.com.br) 3. Re: Eliminar Vetores Repetidos.... (Benilton Carvalho) 4. Re: Eliminar Vetores Repetidos.... (andrebvs(a)bol.com.br) 5. Re: Eliminar Vetores Repetidos.... (Benilton Carvalho) 6. Re: Eliminar Vetores Repetidos.... (Benilton Carvalho) 7. Re: Eliminar Vetores Repetidos.... (andrebvs(a)bol.com.br) 8. Re: Eliminar Vetores Repetidos.... (Elias Teixeira Krainski) 9. Re: Nome das linhas em um data frame (Paulo Justiniano) ---------------------------------------------------------------------- Message: 1 Date: Thu, 14 Mar 2013 23:35:35 -0300 From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Eliminar Vetores Repetidos.... Message-ID: <CAO-arWOhcFpUORJqsfobJvOm-W1fw-a7rH-+MpY4pqJSvWfqtQ(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" As linhas 1 e 5 só são iguais depois de ordenar. Correto? Seu objeto final pode ter as linhas ordenadas? On Mar 14, 2013 11:19 PM, <andrebvs(a)bol.com.br> wrote: > Olá colegas! > > Gostaria de saber como posso eliminar de uma matrix vetores repetidos, no > exemplo abaixo, a linha1 e a linha5 são iguais, gostaria que houvesse uma > outra matrix ou dataframe com apenas as linhas não repetidas. > > Seja a matrix 6x4: > > linha1 <- 1,2,4,5 > linha2 <- 7,2,1,5 > linha3 <- 9,3,1,5 > linha4 <- 6,2,1,5 > linha5 <- 5,2,1,4 > linha6 <- 9,2,1,5 > > desde já agradeço! > > *Att.* > *André* > > _______________________________________________ > R-br mailing list > R-br(a)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. > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130314/8c6945e1/atta…> ------------------------------ Message: 2 Date: Thu, 14 Mar 2013 23:36:46 -0300 From: andrebvs(a)bol.com.br To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Eliminar Vetores Repetidos.... Message-ID: <5142893ea22da_53997e1447815f(a)a4-winter14.tmail> Content-Type: text/plain; charset="iso-8859-1" Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130314/bde11ea4/atta…> ------------------------------ Message: 3 Date: Thu, 14 Mar 2013 23:46:45 -0300 From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Eliminar Vetores Repetidos.... Message-ID: <CAO-arWMAE3Sbbys83y-S=8jfQ7s6DQJz9upKJp-OzxfwdSTR+Q(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" unique(t(apply(x, 1, sort))) On Mar 14, 2013 11:36 PM, <andrebvs(a)bol.com.br> wrote: > Olá Benilton! > > No resultado final, tanto faz, ter linha ordenada ou não. > > ------------------------------ > Em 14/03/2013 23:35, *Benilton Carvalho < beniltoncarvalho(a)gmail.com >*escreveu: > > As linhas 1 e 5 só são iguais depois de ordenar. Correto? Seu objeto final > pode ter as linhas ordenadas? > On Mar 14, 2013 11:19 PM, <andrebvs(a)bol.com.br<http://../undefined/compose?to=andrebvs@bol.com.br>> > wrote: > >> Olá colegas! >> >> Gostaria de saber como posso eliminar de uma matrix vetores repetidos, no >> exemplo abaixo, a linha1 e a linha5 são iguais, gostaria que houvesse uma >> outra matrix ou dataframe com apenas as linhas não repetidas. >> >> Seja a matrix 6x4: >> >> linha1 >> linha2 >> linha3linha4 >> linha5linha6 >> >> desde já agradeço! >> >> *Att.* >> *André* >> >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br<http://../undefined/compose?to=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(a)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. > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130314/5236f633/atta…> ------------------------------ Message: 4 Date: Thu, 14 Mar 2013 23:50:45 -0300 From: andrebvs(a)bol.com.br To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Eliminar Vetores Repetidos.... Message-ID: <51428c8569331_5355d95547c1029(a)a4-winter11.tmail> Content-Type: text/plain; charset="iso-8859-1" Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130314/ac92d8b3/atta…> ------------------------------ Message: 5 Date: Fri, 15 Mar 2013 00:00:57 -0300 From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Eliminar Vetores Repetidos.... Message-ID: <CAO-arWMEHpC353V+3+EUG=uzfwi=iYY1opvEiFvwfLnsp0GsCQ(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" duplicated(t(apply(x, 1, sort)), fromLast=FALSE) | duplicated(t(apply(x, 1, sort)), fromLast=TRUE) On Mar 14, 2013 11:50 PM, <andrebvs(a)bol.com.br> wrote: > Obrigado benilton, tinha acabado de utilizar esse simples comando mesmo: > > L1 <- c(1,2,4,5) > L2 <- c(7,2,1,5) > L3 <- c(9,3,1,5) > L4 <- c(6,2,1,5) > L5 <- c(5,2,1,4) > L6 <- c(9,2,1,5) > > M = matrix(c(L1,L2,L3,L4,L5,L6),,4,byrow=T) > > unique(t(apply(M,1,sort))) > > Agora, teria como verificar, por exemplo, quais linhas foram repetidas na > matrix M? > > *Att.* > *André* > > ------------------------------ > Em 14/03/2013 23:46, *Benilton Carvalho < beniltoncarvalho(a)gmail.com >*escreveu: > > unique(t(apply(x, 1, sort))) > On Mar 14, 2013 11:36 PM, <andrebvs(a)bol.com.br<http://../undefined/compose?to=andrebvs@bol.com.br>> > wrote: > >> Olá Benilton! >> >> No resultado final, tanto faz, ter linha ordenada ou não. >> >> ------------------------------ >> Em 14/03/2013 23:35, *Benilton Carvalho < beniltoncarvalho(a)gmail.com<http://../undefined/compose?to=beniltoncarvalho@gmail.com>> >> * escreveu: >> >> As linhas 1 e 5 só são iguais depois de ordenar. Correto? Seu objeto >> final pode ter as linhas ordenadas? >> On Mar 14, 2013 11:19 PM, <andrebvs(a)bol.com.br<http://../undefined/compose?to=andrebvs@bol.com.br>> >> wrote: >> >>> Olá colegas! >>> >>> Gostaria de saber como posso eliminar de uma matrix vetores repetidos, >>> no exemplo abaixo, a linha1 e a linha5 são iguais, gostaria que houvesse >>> uma outra matrix ou dataframe com apenas as linhas não repetidas. >>> >>> Seja a matrix 6x4: >>> >>> linha1 >>> linha2 >>> linha3linha4 >>> linha5linha6 >>> >>> desde já agradeço! >>> >>> *Att.* >>> *André* >>> >>> _______________________________________________ >>> R-br mailing list >>> R-br(a)listas.c3sl.ufpr.br<http://../undefined/compose?to=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(a)listas.c3sl.ufpr.br<http://../undefined/compose?to=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(a)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. > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130315/068a9c48/atta…> ------------------------------ Message: 6 Date: Fri, 15 Mar 2013 00:01:20 -0300 From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Eliminar Vetores Repetidos.... Message-ID: <CAO-arWMuYJhu+kN58urRXLVxY1e1AgDHyGcy-E-i8fUtBT9B2g(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Tudo em uma linha só On Mar 15, 2013 12:00 AM, "Benilton Carvalho" <beniltoncarvalho(a)gmail.com> wrote: > duplicated(t(apply(x, 1, sort)), fromLast=FALSE) | duplicated(t(apply(x, > 1, sort)), fromLast=TRUE) > On Mar 14, 2013 11:50 PM, <andrebvs(a)bol.com.br> wrote: > >> Obrigado benilton, tinha acabado de utilizar esse simples comando mesmo: >> >> L1 <- c(1,2,4,5) >> L2 <- c(7,2,1,5) >> L3 <- c(9,3,1,5) >> L4 <- c(6,2,1,5) >> L5 <- c(5,2,1,4) >> L6 <- c(9,2,1,5) >> >> M = matrix(c(L1,L2,L3,L4,L5,L6),,4,byrow=T) >> >> unique(t(apply(M,1,sort))) >> >> Agora, teria como verificar, por exemplo, quais linhas foram repetidas na >> matrix M? >> >> *Att.* >> *André* >> >> ------------------------------ >> Em 14/03/2013 23:46, *Benilton Carvalho < beniltoncarvalho(a)gmail.com >*escreveu: >> >> unique(t(apply(x, 1, sort))) >> On Mar 14, 2013 11:36 PM, <andrebvs(a)bol.com.br<http://../undefined/compose?to=andrebvs@bol.com.br>> >> wrote: >> >>> Olá Benilton! >>> >>> No resultado final, tanto faz, ter linha ordenada ou não. >>> >>> ------------------------------ >>> Em 14/03/2013 23:35, *Benilton Carvalho < beniltoncarvalho(a)gmail.com<http://../undefined/compose?to=beniltoncarvalho@gmail.com>> >>> * escreveu: >>> >>> As linhas 1 e 5 só são iguais depois de ordenar. Correto? Seu objeto >>> final pode ter as linhas ordenadas? >>> On Mar 14, 2013 11:19 PM, <andrebvs(a)bol.com.br<http://../undefined/compose?to=andrebvs@bol.com.br>> >>> wrote: >>> >>>> Olá colegas! >>>> >>>> Gostaria de saber como posso eliminar de uma matrix vetores repetidos, >>>> no exemplo abaixo, a linha1 e a linha5 são iguais, gostaria que houvesse >>>> uma outra matrix ou dataframe com apenas as linhas não repetidas. >>>> >>>> Seja a matrix 6x4: >>>> >>>> linha1 >>>> linha2 >>>> linha3linha4 >>>> linha5linha6 >>>> >>>> desde já agradeço! >>>> >>>> *Att.* >>>> *André* >>>> >>>> _______________________________________________ >>>> R-br mailing list >>>> R-br(a)listas.c3sl.ufpr.br<http://../undefined/compose?to=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(a)listas.c3sl.ufpr.br<http://../undefined/compose?to=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(a)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. >> > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130315/4352bc55/atta…> ------------------------------ Message: 7 Date: Fri, 15 Mar 2013 00:06:01 -0300 From: andrebvs(a)bol.com.br To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Eliminar Vetores Repetidos.... Message-ID: <5142901979a1c_7a8c3ab6c6c1e1(a)a4-winter18.tmail> Content-Type: text/plain; charset="iso-8859-1" Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130315/8e99627e/atta…> ------------------------------ Message: 8 Date: Fri, 15 Mar 2013 11:21:03 +0100 From: Elias Teixeira Krainski <eliaskrainski(a)yahoo.com.br> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Eliminar Vetores Repetidos.... Message-ID: <1363342863.2119.2.camel@eliasn2> Content-Type: text/plain; charset="UTF-8" Vc não enviou um CMR. Supondo x sua matriz, use unique(x) ------------------------------ Message: 9 Date: Fri, 15 Mar 2013 09:57:40 -0300 (BRT) From: Paulo Justiniano <paulojus(a)leg.ufpr.br> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Nome das linhas em um data frame Message-ID: <alpine.DEB.2.00.1303150955570.5400(a)pataxo.est.ufpr.br> Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" importe com dados<-read.csv("my_data.csv",h=T, row.names=1) uma vez gravado como texto nao há indicacao para o que é me de linha ou variavel se quiser preservar a estrutura do objeto xporte em outro formado em que tal estrutura seja preservada (dput, save, dump ou similar) On Thu, 14 Mar 2013, Luciano F. Sgarbi wrote: > > Prezados, > Estou com um pequeno problema em importar um conjunto de dados para o R. > Eu estou criando um conjunto de dados e salvando ele no formato ".csv" > Quando eu importo o conjunto de dados para o R, o nome das linhas se transforma em uma nova variável. > No entanto quando eu salvo o arquivo em formato ".txt" esse erro > não ocorre. > Acredito que a resolução deste problema seja muito simples, no entanto não estou conseguindo encontrar essa tal solução. > > Se puderem me ajudar ficarei muito grato! > O CMR segue abaixo: > > my.data<-data.frame(v1=c(letters[1:5]),v2=1:5,v3=100:104) > rownames(my.data)<-LETTERS[1:5] > my.data > write.csv(my.data,"my_data.csv") > dados<-read.csv("my_data.csv",h=T) > dados > > > Abraços, > Luciano > > -- > > Luciano F. Sgarbi > > Mestrando em Ecologia e Evolução > Laboratório de Ecologia de Insetos (sl. 222) > Departamento de Ecologia > Instituto de Ciências Biológicas - ICB 1 > Universidade Federal de Goiás, campus II > Goiânia-GO > 74001-970 > Brazil > > > > ------------------------------ _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Fim da Digest R-br, volume 25, assunto 23 *****************************************

2 3

Eliminar Vetores Repetidos....
by andrebvs＠bol.com.br 15 Mar '13

15 Mar '13

3 10