Março 2015 - R-br

[OFF -Tópic] Cotas em concursos para negros e índios!
by Andre Oliveira 12 Jun '49

12 Jun '49

Pessoal, tem um colegar meu que vou fazer um concurso ele é descendente de negro ele é pardo. Ele veio com a pergunta como provar que sou negro, pois, tem vagas reservada para negros e índios? Eu sinceramente não sei e não tenho ideia pois nunca tive esta curiosidade. O que de fato comprova sua raça? Algum documento? O que são usados em concursos para provar? obrigado André Oliveira Souza

5 6

Legenda eixo X
by Luciano F. Sgarbi 10 Jun '49

10 Jun '49

Prezados, Estou tentando centralizar a legenda do eixo X em uma figura com vários gráficos. Estou conseguindo centralizar a legenda da linha inferior como no exemplo, mas as linhas do meio não. par(mfrow=c(2,2)) plot(rnorm(100),xlab="",ylab="y1") plot(rnorm(100),xlab="",ylab="y2") mtext(side=1,cex=2,line=-1.8,text="Variável Explanatória 1",outer=TRUE) ### o problema está aqui, pois a legenda desce para a parte inferior da figura. plot(rnorm(100),xlab="",ylab="y1") plot(rnorm(100),xlab="",ylab="y2") mtext(side=1,cex=2,line=-1.8,text="Variável Explanatória 2",outer=TRUE)## Esse dá certo Alguém podia dar uma luz de como eu posso fazer isso? Abraços -- *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

2 2

Estatísticas do R
by Andre Oliveira 26 Abr '16

26 Abr '16

Pessoal,preciso de estatísticas do R tais como nº de bibliotecas oficiais, pacotes que o R já vem na instalação. André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

8 8

Re: [R-br] Render {rmarkdown} não funciona em emacs
by Rafael Tieppo 13 Fev '16

13 Fev '16

Fernando, por acaso testei ontem o "render" com 'pandoc'. Aqui funcionou 100% para gerar o pdf. Só que rodei o comando 'render(xxx.Rmd , pdf)' de uma janela aberta no modo terminal fora do emacs. não testei de outro modo. Rafael Tieppo State University of Mato Grosso - Department of Agricultural Engineering site: http://www2.unemat.br/rafaeltieppo blog: https://sistemasagricolas.wordpress.com "Evite o desperdício: antes de imprimir pense na sua responsabilidade com o ambiente". On Monday, March 16, 2015 11:00 AM, "r-br-request(a)listas.c3sl.ufpr.br" <r-br-request(a)listas.c3sl.ufpr.br> wrote: 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. Render {rmarkdown} não funciona em emacs (Fernando Antonio de souza) 2. loop lento (Roney Fraga Souza) 3. RES: [R-brUso do laço for para impessão de mais de um arquivo gerado. (Mauro Sznelwar) 4. Re: RES: [R-brUso do laço for para impessão de mais de um arquivo gerado. (izi) 5. Re: loop lento (Daniel Batista Lemes) 6. Re: loop lento (Roney Fraga Souza) 7. Re: loop lento (Daniel Batista Lemes) 8. Chamada de Trabalhos - VIII Reunião da ABAVE (Leandro Marino) 9. RES: RES: [R-brUso do laço for para impessão de mais de um arquivo gerado. (Mauro Sznelwar) ---------------------------------------------------------------------- Message: 1 Date: Sun, 15 Mar 2015 12:57:45 -0300 From: Fernando Antonio de souza <nandodesouza(a)gmail.com> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Render {rmarkdown} não funciona em emacs Message-ID: <CAFrWFskiv2-XGy0jo=4FdyqzRsh3xqFzeaQ-y1Gj1m8HPjgjEg(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Caros amigos Preciso rodar um arquivo *.rmd e gostaria de utilizar a função Render do pacote rmarkdown. Acontece que ao utilizá-la ela retorna a seguinte mensagem "Error: pandoc version 1.12.3 or higher is required and was not found." Eu já atualizee a versão do pandoc mas a mensagem continau. Fiz uma busca na internte e encontrei esta postagem do Walmes que teve o mesmo problema http://r.789695.n4.nabble.com/Apply-rmarkdown-render-outside-the-RStudio-do… Segundo o Walmes ele solucionou criando um soft-link. Sou usuário novo em linux. Como posso fazer isso? abraços -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150315/59c3c8e2/atta…> ------------------------------ Message: 2 Date: Sun, 15 Mar 2015 19:25:03 +0000 From: Roney Fraga Souza <roneyfraga(a)gmail.com> To: R-br <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] loop lento Message-ID: <10C667AD-0A7A-425E-9D05-2CA56ACFE7D1(a)gmail.com> Content-Type: text/plain; charset="utf-8" Olá pessoal, Eu tenho uma lista (ver2) composta por com 14 elementos, cada elemento é um data.frame() contendo as mesmas variáveis, o que varia é apenas a quantidade de linhas. Cada elemento da lista é um ano específico do data.frame(), exemplo, ver2[[1]] é o ano de 1977, ver2[[2]] é 1978, ver2[[3]] é 1979 até o último ano ver2[[14]] que é 1990. Nos data.frame() tenho informações sobre artigos científicos: nome (código que não se repete), grupo (extraído via processo de clusterização), PY (ano que o artigo foi publicado) e one (variável apenas com valor 1 que uso em outros procedimentos). Os data.frame() são cumulativos, no ano de 1978 estão todos os artigos publicados até 1978, em 1979 estão todos os artigos publicados até 1978 mais os artigos publicados em 1979, no em de 1990 estão todos os artigos publicados até 1990. O procedimento de clusterização é realizado de modo independente para cada ano, logo o artigo 'Baker, 1976, V4, P5' pode estar em grupos diferentes para cada ano. Quero saber a origem dos artigos que estão em um determinado grupo. Exemplo, grupo 1 de 1990, é composto por 100 artigos, desses 25 faziam parte do grupo 1 em 1989, 30 do grupo 2 em 1989, 35 do grupo 3 em 1989, e os 10 restantes eu não tenho interesse, pois, foram publicados no ano de 1990 logo não existim nos anos anteriores. Segue exemplo do resultado que eu obtenho o código atual: label qtde tm1.grupo tm1.ano t.grupo t.ano 1-2 6 1 1989 2 1990 2-1 5 2 1989 1 1990 3-1 3 3 1989 1 1990 4-5 3 4 1989 5 1990 5-4 2 5 1989 4 1990 6-9 2 6 1989 9 1990 7-6 2 7 1989 6 1990 8-8 2 8 1989 8 1990 9-10 2 9 1989 10 1990 qtde = quantidade de artigos em cada grupo tm1.grupo = ano t menos 1 determinado grupo tm1.ano = ano t menos 1 t.grupo = grupo t.ano = ano Para chegar a tal resultado eu utilizo o seguinte código: # ------------------------------ # início do código # ------------------------------ # baixar o arquivo 'ver2.rds' nesse link: # https://db.tt/13dT9XQI <https://db.tt/13dT9XQI> ver2 <- readRDS("~/Downloads/ver2.rds") grupo <- list() for(k in 2:length(ver2)){ grp <- sort(unique((ver2[[k-1]]$grupo))) ### qtde de grupos no ano t-1 grp2 <- sort(unique((ver2[[k]]$grupo))) ### qtde de grupos no ano t RES <- LAB <- list() for (i in grp) { for (j in grp2) { RES <- append(RES,list(length(intersect(ver2[[k-1]]$name[ver2[[k-1]]$grupo==i], ver2[[k]]$name[ver2[[k]]$grupo==j])))) LAB <- append(LAB, paste(i,j, sep='-')) } } grupo[[k]] <- data.frame(label=sapply(LAB, "["), qtde=sapply(RES, sum)) grupo[[k]] <- subset(grupo[[k]],qtde>0) grupo[[k]]$tm1.grupo <- as.numeric(gsub('-.*','',grupo[[k]]$label)) grupo[[k]]$tm1.ano <- as.numeric(max(ver2[[k]]$PY)-1) grupo[[k]]$t.grupo <- as.numeric(gsub('^.-','',grupo[[k]]$label)) grupo[[k]]$t.ano <- as.numeric(max(ver2[[k]]$PY)) } # grupo # ------------------------------ # fim do código # ------------------------------ Meu problema é: esse código funciona mas é muito lento, muito mesmo, inviabilizando sua aplicação para o volume de dados que trabalho no dia a dia. Alguém conhece uma forma mais rápida de fazer isso. Obs.: parte desse código foi contribuição de Eder Comunello através dessa lista, mas naquela altura o problema estava estruturado de outra forma. Abraço Roney -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150315/a40bddcd/atta…> ------------------------------ Message: 3 Date: Sun, 15 Mar 2015 20:43:09 -0300 From: "Mauro Sznelwar" <sznelwar(a)uol.com.br> To: <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] RES: [R-brUso do laço for para impessão de mais de um arquivo gerado. Message-ID: <00e001d05f79$cb0b2090$612161b0$(a)com.br> Content-Type: text/plain; charset="utf-8" Por que não estou conseguindo rodar? symbol<-c("XOM","AAPL", "DIS") getCSV <- function(symbol){ + URL <- paste0("http://chartapi.finance.yahoo.com/instrument/1.0/", symbol, "/chartdata;type=quote;range=1d/csv") + tab <- read.table(URL, sep = ",", dec = ".", skip = 17) + colnames(tab) <- paste(c("timestamp","close","high","low","open","volume"), 1:6, sep = "-") + rownames(tab) <- paste(symbol, 1:nrow(tab)) + SAIDA <- paste0("Dados",symbol,"yahoo20150312.txt") + write.table(tab, file = SAIDA, sep = ";", dec = ".", header = T) + } sapply(symbol, getCSV) Error in write.table(tab, file = SAIDA, sep = ";", dec = ".", header = T) : unused argument (header = T) Oi Luis, Assim fica bem mais elegante =) Muito obrigada pela dica! Queria aproveitar para perguntar se você saberia um jeito de obter do yahoo usando o R todos os tickers das companhias pertencentes ao NIKKEI, porque ao invés de escrever os 225 tickers no symbol, eu puxaria direto. Muito obrigada mais uma vez! Em 13 de março de 2015 09:50, Luis G. S. e Silva <lgsilvaesilva(a)gmail.com> escreveu: Michelle, Essa tarefa pode ser realizada eliminado os loops (for) e algumas linha também. Dá uma olhadinha no script abaixo. Abraço symbol<-c("XOM","AAPL", "DIS") getCSV <- function(symbol){ URL <- paste0("http://chartapi.finance.yahoo.com/instrument/1.0/", symbol, "/chartdata;type=quote;range=1d/csv") tab <- read.table(URL, sep = ",", dec = ".", skip = 17) colnames(tab) <- paste(c("timestamp","close","high","low","open","volume"), 1:6, sep = "-") rownames(tab) <- paste(symbol, 1:nrow(tab)) SAIDA <- paste0("Dados",symbol,"yahoo20150312.txt") write.table(tab, file = SAIDA, sep = ";", dec = ".", header = T) } sapply(symbol, getCSV) --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150315/88a0159d/atta…> ------------------------------ Message: 4 Date: Sun, 15 Mar 2015 21:27:02 -0300 From: izi <salah3.1416(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] RES: [R-brUso do laço for para impessão de mais de um arquivo gerado. Message-ID: <1426465622.4713.2.camel(a)gmail.com> Content-Type: text/plain; charset="UTF-8" Retire o argumento <header = T> da função write.table ficando assim: write.table(tab, file = SAIDA, sep = ";", dec = ".") -- Salah Saudações! ------------------------------ Message: 5 Date: Sun, 15 Mar 2015 21:28:50 -0300 From: Daniel Batista Lemes <dlemes(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] loop lento Message-ID: <CAPa99wrNeQ8JZ_=d9Rn0B7GkauicANPgxKjehDmU3=7n_fL2Dw(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Roney, Não sei teu conhecimento de programação(não uso R, mas sou programador) o que tu está fazendo é 3 loop, onde, me parece, que tu percorre toda a matriz todas as vezes, uma forma de otimizar o que tu está fazendo é não fazer aquele sort ali, tu pode organizar os dados forra do loop, porque aquilo diminui muito a performance. Não consegui ver o código, mas quem sabe tu posta ele aqui http://pastebin.com/ assim como código fica melhor de analisar. Me parece que tu esta usando esse código apenas para organizar os dados, se for isso, existe linguagens para eficientes para fazer isso(python é um exemplo) e ai tu usa o R só para fazer as análises estatísticas. Caso tu não conheça python e queira manter o R existem implementações com melhor performance como o pqR(http://www.pqr-project.org/) ou o fastr ( https://github.com/allr/fastr) Sei que não te ajudei muito, mas espero ter te dado um norte :D Daniel Lemes Em 15 de março de 2015 16:25, Roney Fraga Souza <roneyfraga(a)gmail.com> escreveu: > Olá pessoal, > > Eu tenho uma lista (ver2) composta por com 14 elementos, cada elemento é > um data.frame() contendo as mesmas variáveis, o que varia é apenas a > quantidade de linhas. Cada elemento da lista é um ano específico do > data.frame(), exemplo, ver2[[1]] é o ano de 1977, ver2[[2]] é 1978, > ver2[[3]] é 1979 até o último ano ver2[[14]] que é 1990. Nos data.frame() > tenho informações sobre artigos científicos: nome (código que não se > repete), grupo (extraído via processo de clusterização), PY (ano que o > artigo foi publicado) e one (variável apenas com valor 1 que uso em outros > procedimentos). Os data.frame() são cumulativos, no ano de 1978 estão todos > os artigos publicados até 1978, em 1979 estão todos os artigos publicados > até 1978 mais os artigos publicados em 1979, no em de 1990 estão todos os > artigos publicados até 1990. O procedimento de clusterização é realizado de > modo independente para cada ano, logo o artigo 'Baker, 1976, V4, P5' pode > estar em grupos diferentes para cada ano. > > Quero saber a origem dos artigos que estão em um determinado grupo. > Exemplo, grupo 1 de 1990, é composto por 100 artigos, desses 25 faziam > parte do grupo 1 em 1989, 30 do grupo 2 em 1989, 35 do grupo 3 em 1989, e > os 10 restantes eu não tenho interesse, pois, foram publicados no ano de > 1990 logo não existim nos anos anteriores. Segue exemplo do resultado que > eu obtenho o código atual: > > label qtde tm1.grupo tm1.ano t.grupo t.ano > 1-2 6 1 1989 2 1990 > 2-1 5 2 1989 1 1990 > 3-1 3 3 1989 1 1990 > 4-5 3 4 1989 5 1990 > 5-4 2 5 1989 4 1990 > 6-9 2 6 1989 9 1990 > 7-6 2 7 1989 6 1990 > 8-8 2 8 1989 8 1990 > 9-10 2 9 1989 10 1990 > > qtde = quantidade de artigos em cada grupo > tm1.grupo = ano t menos 1 determinado grupo > tm1.ano = ano t menos 1 > t.grupo = grupo > t.ano = ano > > Para chegar a tal resultado eu utilizo o seguinte código: > # ------------------------------ > # início do código > # ------------------------------ > > # baixar o arquivo 'ver2.rds' nesse link: > # https://db.tt/13dT9XQI > > ver2 <- readRDS("~/Downloads/ver2.rds") > > grupo <- list() > > for(k in 2:length(ver2)){ > grp <- sort(unique((ver2[[k-1]]$grupo))) ### qtde de grupos no ano t-1 > grp2 <- sort(unique((ver2[[k]]$grupo))) ### qtde de grupos no ano t > RES <- LAB <- list() > for (i in grp) { > for (j in grp2) { > RES <- > append(RES,list(length(intersect(ver2[[k-1]]$name[ver2[[k-1]]$grupo==i], > ver2[[k]]$name[ver2[[k]]$grupo==j])))) > LAB <- append(LAB, paste(i,j, sep='-')) > } > } > grupo[[k]] <- data.frame(label=sapply(LAB, "["), qtde=sapply(RES, sum)) > grupo[[k]] <- subset(grupo[[k]],qtde>0) > grupo[[k]]$tm1.grupo <- as.numeric(gsub('-.*','',grupo[[k]]$label)) > grupo[[k]]$tm1.ano <- as.numeric(max(ver2[[k]]$PY)-1) > grupo[[k]]$t.grupo <- as.numeric(gsub('^.-','',grupo[[k]]$label)) > grupo[[k]]$t.ano <- as.numeric(max(ver2[[k]]$PY)) > } > # grupo > # ------------------------------ > # fim do código > # ------------------------------ > > Meu problema é: esse código funciona mas é muito lento, muito mesmo, > inviabilizando sua aplicação para o volume de dados que trabalho no dia a > dia. Alguém conhece uma forma mais rápida de fazer isso. > > Obs.: parte desse código foi contribuição de Eder Comunello através dessa > lista, mas naquela altura o problema estava estruturado de outra forma. > > Abraço > Roney > > _______________________________________________ > 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. > -- @lemes_daniel -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150315/301f5b32/atta…> ------------------------------ Message: 6 Date: Mon, 16 Mar 2015 01:09:48 +0000 From: Roney Fraga Souza <roneyfraga(a)gmail.com> To: R-br <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] loop lento Message-ID: <3DA46750-883C-49E3-8F43-3D3D79625304(a)gmail.com> Content-Type: text/plain; charset="windows-1252" Olá Daniel, obrigado pela ajuda! Sou economista sem estudo formal de algoritmo, daí a chance de fazer cagada é grande. Segue link do código: https://gist.github.com/roneyfraga/debb242d919d9fdb3412 <https://gist.github.com/roneyfraga/debb242d919d9fdb3412> Quanto ao sort(), eu preciso fazer ele dentro do primeiro loop, pois e feito um para cada ano e/ou data.frame(). Conversando no IRC freenode #R recomendaram eu substituir a ?append? pelo ?outer?, sendo a segunda uma função mais rápida. Vou tentar ir nesse caminho no primeiro momento, caso não tenho sucesso estou pensando em me aventurar na mundo de Julia. Abraço Roney -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150316/d6e9d207/atta…> ------------------------------ Message: 7 Date: Sun, 15 Mar 2015 23:54:10 -0300 From: Daniel Batista Lemes <dlemes(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] loop lento Message-ID: <CAPa99wqAvLd9dLr7GsxQRT+AXRiq7fPYFSt5-piG5bxMNDYfuA(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Bom a primeira coisa que eu faria era gerar um data.frame com os dados organizados, ai tu teria um loop separado para organizar os dados e depois poderia usar os outros, já te daria uma certa performance Em 15 de março de 2015 22:09, Roney Fraga Souza <roneyfraga(a)gmail.com> escreveu: > Olá Daniel, obrigado pela ajuda! > > Sou economista sem estudo formal de algoritmo, daí a chance de fazer > cagada é grande. Segue link do código: > https://gist.github.com/roneyfraga/debb242d919d9fdb3412 > > Quanto ao sort(), eu preciso fazer ele dentro do primeiro loop, pois e > feito um para cada ano e/ou data.frame(). > > Conversando no IRC freenode #R recomendaram eu substituir a ?append? pelo > ?outer?, sendo a segunda uma função mais rápida. Vou tentar ir nesse > caminho no primeiro momento, caso não tenho sucesso estou pensando em me > aventurar na mundo de Julia. > > Abraço > Roney > > > _______________________________________________ > 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. > -- @lemes_daniel -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150315/f354ba5e/atta…> ------------------------------ Message: 8 Date: Wed, 11 Mar 2015 16:11:40 -0300 From: Leandro Marino <leandromarino(a)leandromarino.com.br> To: undisclosed-recipients:; Subject: [R-br] Chamada de Trabalhos - VIII Reunião da ABAVE Message-ID: <CAKSaaFnKY5K+u0mp8eqDUD3qZWWK3uF3Dpdr-_5OQVP8Ww3WJw(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" *VIII Reunião da ABAVE * *Avaliação de Larga Escala no Brasil: * *ensinamentos, aprendizagens e tendências* Florianópolis/SC ? 19 a 21 de agosto de 2015 A Diretoria e os membros da Diretoria Científica da ABAVE convocam os interessados a enviar os trabalhos originais, contendo resultados de suas atividades de pesquisa científica ou de execução de projetos de avaliação educacional tanto da educação básica, quanto da educação superior, para apresentação durante a reunião. O formato da reunião foi pensado para favorecer a discussão e troca de experiências. Os autores deverão submeter seus artigos em duas categorias: - *Apresentação oral:* Serão aceitos textos que apresentem pesquisas com conclusão parciais ou finais, abordando temáticas novas ou já estabelecidas na área de avaliação educacional, preferencialmente utilizando os dados disponíveis das avaliações educacionais em larga escala, que evidenciem elaboração teórica e rigor conceitual e metodológico na análise. - *Apresentação pôster: *Serão aceitos textos que apresentem projetos de pesquisa em andamento ou relatos de experiências na área de avaliação educacional. Os trabalhos deverão ser submetidos eletronicamente até 30 de março de 2015. Em breve serão divulgados os procedimentos para a submissão eletrônica dos trabalhos. No ato da submissão, os próprios autores deverão classificar o trabalho em uma das categorias: apresentação oral ou pôster. *Regras para Submissão de Trabalhos:* Os artigos que serão submetidos à análise devem conter as seguintes especificações: - Deve estar no formato PDF, na fonte Arial, corpo 12. - Os artigos podem ser escritos em português, espanhol ou inglês. - A página deverá ser configurada em papel A4, utilizando os seguintes parâmetros: margem superior 2,5 cm; inferior 2,5 cm; lateral esquerda 3,5 cm; lateral direita 2,5 cm. - O trabalho deve conter apenas o título no corpo do texto. NENHUMA informação de autoria deve ser inserida no texto do artigo ou do seu resumo. - A primeira página do artigo deverá conter APENAS as seguintes informações de identificação: 1. Título do texto 2. Nome do autor ? sigla da Instituição 3. Nome(s) do co/autor(es) ? sigla da(s) instituição(es) (quando houver) 4. Agência ou Instituição Financiadora (quando houver) *Para Apresentação Ora**l* ?O texto deve conter no máximo 20 páginas, incluindo resumo, figuras, diagramas,bibliografias e anexos. Colocar resumo em português com, no máximo, 200 palavras, terminado pelas palavras-chave, em negrito. ?As notas de rodapé, em caso de serem indispensáveis, devem ser colocadas no final dotrabalho, antes das referências bibliográficas. ?Os títulos das seções internas devem estar em negrito e posicionado à margemesquerda. ?As figuras e tabelas devem ter legendas e estas, como as referências, devem seguir opadrão ABNT. *Para Apresentação de Pôster* ? O texto deve conter no máximo 10 páginas, incluindo resumo, figuras, diagramas, bibliografias e anexos. ? Colocar resumo em português com, no máximo, 200 palavras, terminado pelas palavras-chave, em negrito. ? As notas de rodapé, em caso de serem indispensáveis, devem ser colocadas no final do trabalho, antes das referências bibliográficas. ? Os títulos das seções internas devem estar em negrito e posicionado à margem esquerda. ? As figuras e tabelas devem ter legendas e estas, como as referências, devem seguir o padrão ABNT. *Avaliação dos Trabalhos* ? Os artigos submetidos serão avaliados pelos membros da Diretoria Científica e por pareceristas ad hoc indicados pela Diretoria. ? Os trabalhos serão avaliados pelo método do double blind review, que possibilita a análise inominada dos artigos, garantindo a imparcialidade da avaliação. ? Serão selecionados, para a apresentação na VIII Reunião Anual da ABAVE, 50 trabalhos na categoria Apresentação Oral e 70 na categoria Pôster. ? Cabe à Comissão Científica a possibilidade de reclassificar o trabalho, a partir de critérios da natureza da pesquisa desenvolvida, como também por disponibilidade de espaço para alocação das sessões apresentação oral. *Datas Relevantes* ? 30 de março: data final para submissão de trabalhos completos, nas categorias apresentação oral e apresentação de pôster; ? 01 de abril a 15 de maio: julgamento dos trabalhos pela Comissão Científica; ? 20 de maio: Divulgação do resultado final sobre os trabalhos aprovados. *Outras Informações* ? Os artigos aprovados e apresentados serão publicados eletronicamente nos anais do evento no sítio da ABAVE. ? O autor que apresentará o artigo submetido deverá estar inscrito no evento para que o artigo seja considerado para apresentação. Outras Informações devem ser solicitadas à: Secretaria da ABAVE (21) 2103-9635 abave(a)abave.org.br -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150311/7d692d4c/atta…> ------------------------------ Message: 9 Date: Mon, 16 Mar 2015 11:13:46 -0300 From: "Mauro Sznelwar" <sznelwar(a)uol.com.br> To: <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] RES: RES: [R-brUso do laço for para impessão de mais de um arquivo gerado. Message-ID: <001a01d05ff3$6c6abad0$45403070$(a)com.br> Content-Type: text/plain; charset="iso-8859-1" Muito obrigado pelo retorno,mas poderia me enviar o código fonte otimizado de novo, eu sem querer acabei acabei apagando. Retire o argumento <header = T> da função write.table ficando assim: write.table(tab, file = SAIDA, sep = ";", dec = ".") -- Salah Saudações! --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com ------------------------------ Subject: Legenda do Digest _______________________________________________ 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 51, assunto 17 *****************************************

3 3

Fwd: Enc: model.matrix()
by Ze Henrique 23 Abr '15

23 Abr '15

> Assunto: model.matrix() > > Bom tarde pessoal, > > gostaria de saber se é possível transformar um nível de um fator em numérico no próprio vetor, de forma que ele não seja codificado na função model.matrix(). > > Segue abaixo um CMR: > > k = rep(0:3,each=3) > k = factor(k) > k = model.matrix(~-1+k) > > # Essa é a matriz gerada. > > k0 k1 k2 k3 > 1 1 0 0 0 > 2 1 0 0 0 > 3 1 0 0 0 > 4 0 1 0 0 > 5 0 1 0 0 > 6 0 1 0 0 > 7 0 0 1 0 > 8 0 0 1 0 > 9 0 0 1 0 > 10 0 0 0 1 > 11 0 0 0 1 > 12 0 0 0 1 > > # Entretanto gostaria de obter a seguinte matriz sem a necessidade de eliminar as colunas posteriormente. > k0 k3 > 1 1 0 > 2 1 0 > 3 1 0 > 4 0 0 > 5 0 0 > 6 0 0 > 7 0 0 > 8 0 0 > 9 0 0 > 10 0 1 > 11 0 1 > 12 0 1 > > Att. > > Tiago. > > ################################################################# > > Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES > > Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas > > #################################################################

7 22

mudar rank de um vetor
by Odirley Campos 18 Abr '15

18 Abr '15

Caros listeiros, bom dia!!Recorro a vcs para resolver um problema talves simples mas que vem me incomodando a alguns meses.Tenho o seguinte vetor retirado de outro data frame: trt<-c(0,30,60,120,180,0,30,60,180,120,0,30,60,120,180) Mas cujo rank é: rank(trt) [1] 2 11 14 5 8 2 11 14 8 5 2 11 14 5 8 Como eu poderia modificar o rank do vetor? Talves colocar o rank com o valor do proprio elemento do vetor. Antecipadamente agradeço. Odirley R. Campos Engenheiro Agrônomo UFV/MG

3 6

Probit/CL50
by Flavia Batista Gomes - Embrapa Amazonia Ocidental - CPAA 10 Abr '15

10 Abr '15

Olá pessoal, Preciso de orientação para fazer a análise de Probit, estimar a CL50 e fazer o gráfico. Meu experimento avaliou o efeito inseticida de doses de uma substância. Se puderem...será muito válido! Flávia Batista ____________________________________________________________________________ Aviso de confidencialidade Esta mensagem da Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), empresa publica federal regida pelo disposto na Lei Federal no. 5.851, de 7 de dezembro de 1972, e enviada exclusivamente a seu destinatario e pode conter informacoes confidenciais, protegidas por sigilo profissional. Sua utilizacao desautorizada e ilegal e sujeita o infrator as penas da lei. Se voce a recebeu indevidamente, queira, por gentileza, reenvia-la ao emitente, esclarecendo o equivoco. Confidentiality note This message from Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), a government company established under Brazilian law (5.851/72), is directed exclusively to its addressee and may contain confidential data, protected under professional secrecy rules. Its unauthorized use is illegal and may subject the transgressor to the law's penalties. If you are not the addressee, please send it back, elucidating the failure.

5 11

Como fazer esse grafico no ggplot2
by Pedro Emmanuel Alvarenga Americano do Brasil 06 Abr '15

06 Abr '15

Amigos de R, Essa semana eu decidi que iria aprender a usar o ggplot2. Logo na primeira tarefa empaquei uns 40 min. Como fazer esse grafico que faço no graphics com o ggplot2? > rec.cum b2 r3 r24 r25 1 1 2014-04-09 Sim 2014-04-09 10 2 2014-04-15 Não 2014-04-15 12 3 2014-04-16 Sim 2014-04-16 21 4 2014-05-16 Sim 2014-05-16 28 5 2014-05-28 Sim 2014-05-28 35 6 2014-07-29 Sim 2014-07-29 42 7 2014-08-27 Sim 2014-08-27 49 8 2014-10-07 Sim 2014-10-08 55 9 2014-10-22 Sim 2014-10-22 61 10 2014-10-27 Não 2014-10-27 78 11 2014-10-29 Sim 2014-10-29 84 12 2014-11-10 Sim 2014-11-10 89 13 2014-11-10 Sim 2014-11-10 93 14 2014-11-19 Não 2014-11-19 94 15 2014-12-02 Sim 2014-12-02 95 16 2014-12-08 Sim 2014-12-08 99 17 2015-02-25 Sim 2015-02-25 100 18 2015-03-03 Sim 2015-03-03 101 19 2015-03-10 Sim 2015-03-10 102 20 2015-03-17 Sim 2015-03-17 plot.recrut <- function(){ plot(rec.cum$r3,1:nrow(rec.cum),type='b',xlab='',ylab='Número de pacientes rastreados/recrutados',lty=2,col=2,pch=6,cex=1,lwd=3,xlim=c(min(rec.cum$r3,na.rm=T),max(rec.cum$r25,na.rm=T)),axes=F);grid();box();axis(2) tmp <- rec.cum[order(rec.cum$r25),] tmp <- tmp[which(tmp$r24=='Sim'),] lines(tmp$r25,1:nrow(tmp),type='b',lty=3,col=4,pch=2,cex=1,lwd=3) rm(tmp) axis.Date(1,rec.cum$r3,format='%d') # axis.Date(1,rec.cum$r3,format='%b',line=1.5,tick=F) axis.Date(1,rec.cum$r3,format='%Y',line=2.5,tick=F) legend('bottomrigh',c(paste('Primeiro rastreado:',min(rec.cum$r3,na.rm=T)),paste('Último rastreado:',max(rec.cum$r3,na.rm=T))),bty='n') legend('topleft',c(paste('Primeiro recrutado:',min(rec.cum$r25,na.rm=T)),paste('Último recrutado:',max(rec.cum$r25,na.rm=T))),bty='n') legend('top',c('Rastreamento','Recrutamento'),lty=c(2,3),col=c(2,4),pch=c(6,2),pt.cex=1.5,lwd=3,bty='n',xpd=T,inset=-.15,horiz=T) } plot.recrut() Abraço forte, Pedro Brasil

3 5

Mudar formatação de uma timestamp separando em colunas.
by Michelle Bau Graczyk 31 Mar '15

31 Mar '15

Boa tarde, Eu tenho um arquivo nesse formato e queria que a coluna de timestamp fosse separada em uma de data e outra de hora pois preciso calcular a media do volume por horário. Para isso eu tentei usar o programa abaixo mas a saida foi NA: > INTC<-read.table("/Users/bau/Documents/Mestrado/Dadosbaixadosdothebonnotgang/INTC_1m.txt", header=TRUE, sep=";", dec=",") > INTC[1:25,] symbol nr timestamp open high low close volume 1 INTC 1 2011-09-09 09:30:59 19.7500 19.86 19.75 19.8400 615100 2 INTC 2 2011-09-09 09:31:59 19.8350 19.86 19.79 19.8600 242900 3 INTC 3 2011-09-09 09:32:58 19.8501 19.86 19.80 19.8001 314400 4 INTC 4 2011-09-09 09:33:56 19.8068 19.81 19.77 19.7800 183500 5 INTC 5 2011-09-09 09:34:59 19.7900 19.81 19.78 19.7901 150900 6 INTC 6 2011-09-09 09:35:59 19.8000 19.81 19.74 19.7700 264600 7 INTC 7 2011-09-09 09:36:57 19.7700 19.81 19.77 19.8001 174300 8 INTC 8 2011-09-09 09:37:57 19.8100 19.81 19.76 19.7700 146200 9 INTC 9 2011-09-09 09:38:59 19.7600 19.83 19.76 19.8300 192100 10 INTC 10 2011-09-09 09:39:59 19.8400 19.87 19.83 19.8600 172500 11 INTC 11 2011-09-09 09:40:59 19.8700 19.87 19.76 19.7700 275000 12 INTC 12 2011-09-09 09:41:58 19.7700 19.77 19.72 19.7500 269200 13 INTC 13 2011-09-09 09:42:59 19.7500 19.82 19.75 19.8100 209600 14 INTC 14 2011-09-09 09:43:57 19.8197 19.82 19.79 19.7900 303200 15 INTC 15 2011-09-09 09:44:57 19.7900 19.82 19.78 19.7900 148000 16 INTC 16 2011-09-09 09:45:58 19.7900 19.79 19.75 19.7550 168700 17 INTC 17 2011-09-09 09:46:59 19.7600 19.85 19.75 19.8500 224200 18 INTC 18 2011-09-09 09:47:57 19.8594 19.88 19.85 19.8700 158600 19 INTC 19 2011-09-09 09:48:01 19.8800 19.91 19.87 19.9000 219500 20 INTC 20 2011-09-09 09:49:59 19.8900 19.92 19.89 19.9200 121000 21 INTC 21 2011-09-09 09:50:59 19.9200 19.93 19.90 19.9200 241400 22 INTC 22 2011-09-09 09:51:59 19.9100 19.95 19.90 19.9300 214000 23 INTC 23 2011-09-09 09:52:59 19.9300 19.94 19.92 19.9300 74200 24 INTC 24 2011-09-09 09:53:59 19.9292 19.93 19.88 19.9000 274400 25 INTC 25 2011-09-09 09:54:59 19.9000 19.92 19.87 19.9200 75800 > strptime(INTC$timestamp," %H:%M:%S") [1] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [31] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [61] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [91] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [121] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [151] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [181] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [211] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [241] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [271] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [ reached getOption("max.print") -- omitted 335637 entries ] Muito obrigada!!!

2 3

Manipulação matriz com laço
by Michelle Bau Graczyk 31 Mar '15

31 Mar '15

Boa noite, Eu estou normalizando autovetores de uma matriz e colocando eles em outra. Para isso eu montei o programa a seguir. A primeira forma, ele só calcula o ultimo elemento da coluna 4 da matriz, e da segunda forma os valores do resultado estão errados. Primeira forma: > a<-read.table("/Users/bau/Documents/Mestrado/teste/eigen1.txt", header=TRUE) > b<-as.matrix(a) > matriznorm<-matrix(0,30,4) > matriznorm[,1]<-b[,1]# só vou usar este valores do arquivo acima > show(matriznorm[,1]) [1] 10.373173554 6.125265914 3.148949951 2.722401273 2.045041946 1.498078581 1.350122780 [8] 1.335457110 1.190051106 0.999610112 0.850212256 0.716720583 0.600704889 0.474433269 [15] 0.431826856 0.388576986 0.266709936 0.176960029 0.045935755 0.020126910 -0.003897135 [22] -0.030136936 -0.074815367 -0.145410630 -0.294553778 -0.453450244 -0.551343912 -0.710794973 [29] -1.222537293 -1.273419528 > for( i in 1:30){ + matriznorm[i,2]<-matriznorm[i,1]*matriznorm[i,1] + } > > somaQuadrados<-sum(matriznorm[,2]) > raiz<-sqrt(somaQuadrados) > matriznorm[,3]<-raiz > > for(j in i:30){ + matriznorm[j,4]<-matriznorm[j,1]/raiz + } > > write.table(matriznorm,"/Users/bau/Documents/Mestrado/teste/eigenValuesNorm1.txt") > show(matriznorm) [,1] [,2] [,3] [,4] [1,] 10.373173554 1.076027e+02 13.46876 0.00000000 [2,] 6.125265914 3.751888e+01 13.46876 0.00000000 [3,] 3.148949951 9.915886e+00 13.46876 0.00000000 [4,] 2.722401273 7.411469e+00 13.46876 0.00000000 [5,] 2.045041946 4.182197e+00 13.46876 0.00000000 [6,] 1.498078581 2.244239e+00 13.46876 0.00000000 [7,] 1.350122780 1.822832e+00 13.46876 0.00000000 [8,] 1.335457110 1.783446e+00 13.46876 0.00000000 [9,] 1.190051106 1.416222e+00 13.46876 0.00000000 [10,] 0.999610112 9.992204e-01 13.46876 0.00000000 [11,] 0.850212256 7.228609e-01 13.46876 0.00000000 [12,] 0.716720583 5.136884e-01 13.46876 0.00000000 [13,] 0.600704889 3.608464e-01 13.46876 0.00000000 [14,] 0.474433269 2.250869e-01 13.46876 0.00000000 [15,] 0.431826856 1.864744e-01 13.46876 0.00000000 [16,] 0.388576986 1.509921e-01 13.46876 0.00000000 [17,] 0.266709936 7.113419e-02 13.46876 0.00000000 [18,] 0.176960029 3.131485e-02 13.46876 0.00000000 [19,] 0.045935755 2.110094e-03 13.46876 0.00000000 [20,] 0.020126910 4.050925e-04 13.46876 0.00000000 [21,] -0.003897135 1.518766e-05 13.46876 0.00000000 [22,] -0.030136936 9.082349e-04 13.46876 0.00000000 [23,] -0.074815367 5.597339e-03 13.46876 0.00000000 [24,] -0.145410630 2.114425e-02 13.46876 0.00000000 [25,] -0.294553778 8.676193e-02 13.46876 0.00000000 [26,] -0.453450244 2.056171e-01 13.46876 0.00000000 [27,] -0.551343912 3.039801e-01 13.46876 0.00000000 [28,] -0.710794973 5.052295e-01 13.46876 0.00000000 [29,] -1.222537293 1.494597e+00 13.46876 0.00000000 [30,] -1.273419528 1.621597e+00 13.46876 -0.09454616 > Da segunda forma: > a<-read.table("/Users/bau/Documents/Mestrado/teste/eigen1.txt", header=TRUE) > b<-as.matrix(a) > matriznorm<-matrix(0,30,4) > matriznorm[,1]<-b[,1] > show(matriznorm[,1]) [1] 10.373173554 6.125265914 3.148949951 2.722401273 2.045041946 1.498078581 1.350122780 [8] 1.335457110 1.190051106 0.999610112 0.850212256 0.716720583 0.600704889 0.474433269 [15] 0.431826856 0.388576986 0.266709936 0.176960029 0.045935755 0.020126910 -0.003897135 [22] -0.030136936 -0.074815367 -0.145410630 -0.294553778 -0.453450244 -0.551343912 -0.710794973 [29] -1.222537293 -1.273419528 > for( i in 1:30){ + matriznorm[i,2]<-matriznorm[i,1]*matriznorm[i,1] + + somaQuadrados<-sum(matriznorm[,2]) + raiz<-sqrt(somaQuadrados) + matriznorm[,3]<-raiz + + for(j in i:30){ + matriznorm[j,4]<-matriznorm[j,1]/raiz + } + } > write.table(matriznorm,"/Users/bau/Documents/Mestrado/teste/eigenValuesNorm1.txt") > show(matriznorm) [,1] [,2] [,3] [,4] [1,] 10.373173554 1.076027e+02 13.46876 1.0000000000 [2,] 6.125265914 3.751888e+01 13.46876 0.5084624670 [3,] 3.148949951 9.915886e+00 13.46876 0.2528991576 [4,] 2.722401273 7.411469e+00 13.46876 0.2135962692 [5,] 2.045041946 4.182197e+00 13.46876 0.1584251428 [6,] 1.498078581 2.244239e+00 13.46876 0.1152793173 [7,] 1.350122780 1.822832e+00 13.46876 0.1033376904 [8,] 1.335457110 1.783446e+00 13.46876 0.1016853664 [9,] 1.190051106 1.416222e+00 13.46876 0.0902440218 [10,] 0.999610112 9.992204e-01 13.46876 0.0755856426 [11,] 0.850212256 7.228609e-01 13.46876 0.0641564605 [12,] 0.716720583 5.136884e-01 13.46876 0.0540043425 [13,] 0.600704889 3.608464e-01 13.46876 0.0452163558 [14,] 0.474433269 2.250869e-01 13.46876 0.0356888678 [15,] 0.431826856 1.864744e-01 13.46876 0.0324667092 [16,] 0.388576986 1.509921e-01 13.46876 0.0292025272 [17,] 0.266709936 7.113419e-02 13.46876 0.0200398900 [18,] 0.176960029 3.131485e-02 13.46876 0.0132951406 [19,] 0.045935755 2.110094e-03 13.46876 0.0034511674 [20,] 0.020126910 4.050925e-04 13.46876 0.0015121392 [21,] -0.003897135 1.518766e-05 13.46876 -0.0002927926 [22,] -0.030136936 9.082349e-04 13.46876 -0.0022641887 [23,] -0.074815367 5.597339e-03 13.46876 -0.0056207917 [24,] -0.145410630 2.114425e-02 13.46876 -0.0109238800 [25,] -0.294553778 8.676193e-02 13.46876 -0.0221227473 [26,] -0.453450244 2.056171e-01 13.46876 -0.0340370871 [27,] -0.551343912 3.039801e-01 13.46876 -0.0413498303 [28,] -0.710794973 5.052295e-01 13.46876 -0.0532327969 [29,] -1.222537293 1.494597e+00 13.46876 -0.0911767942 [30,] -1.273419528 1.621597e+00 13.46876 -0.0945461625 Muito obrigada!

2 4