Março 2015 - R-br

Michelle Baú Gaczyk compartilhou "AIGparaCalculodaMedia.txt" com você.
by Michelle Baú Gaczyk via Dropbox 18 Mar '15

18 Mar '15

Michelle usou o Dropbox para compartilhar um arquivo com você! Clique aqui para visualizar.: https://www.dropbox.com/l/oqyecuqKWf5ydyT1wHijAg?text=1

1 0

Michelle Baú Gaczyk compartilhou "AAparaCalculodaMedia.txt" com você.
by Michelle Baú Gaczyk via Dropbox 18 Mar '15

18 Mar '15

Michelle usou o Dropbox para compartilhar um arquivo com você! Clique aqui para visualizar.: https://www.dropbox.com/l/c2Ba0YfWF90baTnAarcYzg?text=1

1 0

Calculo de correlcao entre coluas que apresentam NAS de varias matrizes
by Michelle Bau Graczyk 18 Mar '15

18 Mar '15

Boa tarde, Eu tenho tres matrizes de tres empresas diferentes que estao dispostas da seguinte maneira: Empresa AA data 9:31 data 09:31 09:32 09:33 09:34 09:35 09:36 09:37 09:38 09:39 09:40 09:41 09:42 09:43 09:44 09:45 09:46 09:47 2004-07-01 700 200 7000 9400 1300 1500 6300 1800 4300 1300 4800 3900 2500 4300 4300 2004-07-02 NA NA 68900 2800 4900 1200 12300 NA NA 6500 NA 2900 8900 11600 26900 2004-07-06 800 4400 13700 10700 4700 4000 18700 9000 2100 2300 700 8900 10500 3900 8400 2004-07-07 NA 5600 1000 7300 5600 11200 3900 5800 5700 4500 34400 4100 500 15300 500 2004-07-08 NA NA NA NA NA NA NA NA NA NA 154100 42400 16700 4700 5900 2004-07-09 NA 96900 1000 20400 18300 18000 22500 4100 6000 1000 1900 3700 600 2900 14600 2004-07-12 1400 100 4500 3000 1300 4200 6500 4500 1600 1600 5600 7600 1900 4600 8300 2004-07-13 NA 5600 3400 16900 9000 5700 16500 1800 21000 6300 5900 12500 51000 7500 13100 2004-07-14 NA NA NA NA NA NA 164400 35000 41600 56600 27800 36500 55300 29900 13100 2004-07-15 NA 8100 2000 6100 4900 4500 2500 1400 5200 4100 5400 4000 13400 3200 11300 2004-07-16 282300 18100 30400 NA 15900 40800 9000 5300 5200 27600 7600 12100 8400 1000 8400 2004-07-19 28700 1700 2900 9500 400 2800 1100 6200 11500 17900 11100 16000 6800 5200 7100 Onde as datas vao de 2004-07-01 ate 2004-12-31(ao todo 128) e os horarios vao de 09:31 ate 16:04( ao todo 397) Eu quero calcular a correlacao entre essas tres empresas para cada hora fixa,ou seja, para a empresa AA eu calculo com ela mesma e com a AIG e a AXP para cada hora fixa, e por isso montei este programa: A<-read.table(file="/home/usuario/Documentos/AAparaCalculodaMedia.txt", header=TRUE) AA.<-as.matrix(AA) ##write.table(AA.,"AA..txt") AXP<-read.table("/home/usuario/Documentos/AXPparaCalculodaMedia.txt",header=TRUE) AXP.<-as.matrix(AXP) ##write.table(AXP.,"AXP..txt") AIG<-read.table("/home/usuario/Documentos/AIGparaCalculodaMedia.txt",header=TRUE) AIG.<-as.matrix(AIG) Symbols<-c("AA.","AIG.","AXP.") matrizAA<-matrix() j<-1 for (j in 1:length(Symbols)){ empresa<-paste0(Symbols)#,"[",j,"]") for (i in 2:397){ empresa1<-as.numeric(paste0(empresa,"[,",i,"]")) head(empresa1) matrizAA[,i-1]<-cor(as.numeric(AA.[,i]), as.numeric(empresa1), use="pairwise.complete.obs") } } Mas da o erro: Erro em cor(as.numeric(AA.[, i]), as.numeric(empresa1), use = "pairwise.complete.obs") : incompatible dimensions Além disso: Mensagens de aviso perdidas: NAs introduzidos por coerção Eu achava que usando o use = "pairwise.complete.obs" as correlacoes seriam calculadas apenas entre pares de numeros, mas acho que isso nao esta acontecendo. Os tres arquivos vou deixar em dropbox. Muito obrigada, Michelle

1 0

Merge com mais de 2 bancos
by Vanessa Felix 18 Mar '15

18 Mar '15

Pessoal, eu consigo juntar 3 ou 4 bancos de uma vez só? Não queria ter que fazer vários merges... Abç, -- ------------------------------------------------------------------------------------------ Vanessa Felix do Nascimento Sergio Graduada em Estatística - ICE/UFJF Mestranda em Saúde Coletiva - IMS/UERJ Assistente de pesquisa I - Instituto de Pesquisa Econômica Aplicada – IPEA (32) 9142-8428 (21) 97557-3956 -------------------------------------------------------------------------------------------

5 6

Extrair linhas específicas de um data.frame
by Alexandre Loures 17 Mar '15

17 Mar '15

Boa tarde! Tenho um arquivo em que uma das variáveis tem o seguinte formato: data 1958Q3 2001Q1 1989M7 1948M12 1949 2011 . . . E assim por diante. Ou seja, tenho informações quadrimestrais, mensais e anuais. Todavia, gostaria de obter um novo data.frame apenas com os dados anuais. Alguém saberia como fazer isso? Desde já muito obrigado! -- *Alexandre Rodrigues Loures* Doutorando em Economia Aplicada Universidade Federal da Paraíba - UFPB Centro de Ciências Sociais Aplicadas - CCSA Programa de Pós-Graduação em Economia - PPGE e-Mail: alexandre.loures(a)ymail.com Site: www.ccsa.ufpb.br/ppge

6 7

loop lento
by Roney Fraga Souza 17 Mar '15

17 Mar '15

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

4 7

Uso do laço for para impessão de mais de um arquivo gerado.
by Michelle Bau Graczyk 17 Mar '15

17 Mar '15

Boa noite, Com ajuda de colegas do grupo montei um programa para pegar dados do yahoo e organiza-los. Porém, eu preciso pegar para vårias empresas e por isso fiz um laço com o uso da função for para imprimir a saida. Porem, eu não estou conseguindo fazer com que ele imprima todos e sim, somente o ultimo elemento. No caso abaixo, ele só imprime e salva o arquivo referente a empresa DIS enquanto eu também quero das outras. Já tentei colocar o write.table no final da última chaves mas também não imprimiu de todos os arquivos. Muito obrigada! Michelle symbol<-c("XOM","AAPL", "DIS") for(i in seq_along(symbols)){ symbol <- symbols[i] URL <- paste0("http://chartapi.finance.yahoo.com/instrument/1.0/", symbol, "/chartdata;type=quote;range=1d/csv") tmp <- readLines(URL) tabela<-tmp[18:length(tmp)] tabela_split<-strsplit(tabela,",") # matriz para receber os dados dados<-matrix(NA,ncol=length(tabela_split[[1]]),nrow=length(tabela_split)) for(j in 1:length(tabela_split)) { dados[j,]<-as.numeric(tabela_split[[j]]) } # organizar os nomes de colunas linhas dimnames(dados)<-list(paste(symbol,1:length(tabela_split)),paste(c("timestamp","close","high","low","open","volume"),1:length(tabela_split[[1]]))) SAIDA<-paste0("Dados",symbol,"yahoo20150312.txt") write.table(dados,file= SAIDA) }

6 11

Render {rmarkdown} não funciona em emacs
by Fernando Antonio de souza 15 Mar '15

15 Mar '15

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

1 0

Knitr-emacs dúvida
by Fernando Antonio de souza 14 Mar '15

14 Mar '15

Caros amigos Estou criando um relatório utilizando as funcionalidades do knitr+ html+markdown. Eu desejo após gerar o arquivo *.rmd convertê-lo em documento word/odt utilizando a função pandoc{knitr}. O meu problema é o seguinte. ao fazer a conversão para docx/odt o código chunk é impresso de forma literal no documento. No chunk abaixo mesmo a especificação "echo=FALSE" não foi atendida. O meu objetivo é que no documento doc/odt apareça apenas a tabela gerada pela função kable{knitr} e não o código utilizado para gerá-lo Utilizo o editor emacs+ess. Rodei os comandos no RStudio e nele o comando funcionou perfeitamente. Por isso acredito que seja algo do emacs que desconheço. Para trabalhar com knitr instalei o pacote polymode. ```{r,echo=FALSE} opts_chunk$set(echo=FALSE,message=FALSE,results = "asis") ``` ```{r} #install.packages("repmis") #install.packages("knitr") library(pander) library(repmis) #funçao para importar dados do dropbox library(knitr) dados <- source_DropboxData(file = "DQL.csv",key = "xtp04eaa1316kt2", sep = ";", header = TRUE) dados <- transform(dados,VACA=factor(VACA),PERIODO=factor(PERIODO),TRAT=factor(TRAT)) modelo <- aov(LEITE~VACA+PERIODO+TRAT,data=dados) #install.packages("multcomp") #library(multcomp) comp <- glht(modelo,linfct=mcp(TRAT="Tukey")) nometrat <- colnames(summary(comp)$test$pvalues) medias<- aggregate(LEITE~TRAT,data=dados,FUN=mean) comparacao <-unname( cld(summary(comp))$mcletters$Letters) tabela <- cbind(medias,comparacao) colnames(tabela) <- c("TRATAMENTO","LEITE(Kg)","Tukey (<0,05)") ``` ```{r echo=FALSE,results="asis"} kable(tabela,format="markdown",digits=3,align="c") ```

1 0

Chamada de Trabalhos - VIII Reunião da ABAVE
by Leandro Marino 11 Mar '15

11 Mar '15

*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

1 0