- R-br - DINF UFPR

Bando de dado .txt
by Andre Oliveira 02 Out '16

02 Out '16

Boa noite,não estou conseguindo identificar o motivo de não abrir este arquivo .txt corretamente, ele tem em torno de 6.3 GBAlguém teria uma dica para compartilhar? Acredito que não estou trabalhando corretamente! library(data.table) readLines("DADOS_ENEM_2011.TXT", n = 1) [1] "3000000000012011 2803548708SAO BERNARDO DO CAMPO SP110000000000000000000011 . . ...3548708SAO BERNARDO DO CAMPO SP1111 577.10 576.00 575.50 737.10CEBDCACBEAABBAABAEBDCABDCDEAAADCBBEBADDBEBECDDEECABABEBDDEBAEADCACABCBBECADDAADBADDAABABCAEEAEAEDDDAAEEBDECADCEABDDBAECAAABDAEBCCACBAECBAEBECCAAEBABECECCDDAACDADADBACCBAADCAECCCEDD1241201251290CECDDACBCAEBEBBACEEDCAEBDAEDADBEBBDDADCBCBACDDEECEBABAADDEBAECDCACEDCBBEBCDCAADBADDDBAEACAEEDBDBDDCAEDDDACDEBBECAACDABADBAECCEADDCEBCBACBAEEBAEEECCEEEBCBECEBADCBDBEADCDBDCCCBADCAECCCCDDP 140.00 140.00 140.00 140.00 60.00 620.000 . . . require(bit64) enem2011<- fread("DADOS_ENEM_2011.TXT", nrow = 1, header = TRUE, sep=";") head(enem2011) 3000000000012011 2803548708SAO BERNARDO DO CAMPO SP110000000000000000000011 . . ...3548708SAO BERNARDO DO CAMPO SP1111 577.10 576.00 575.50 737.10CEBDCACBEAABBAABAEBDCABDCDEAAADCBBEBADDBEBECDDEECABABEBDDEBAEADCACABCBBECADDAADBADDAABABCAEEAEAEDDDAAEEBDECADCEABDDBAECAAABDAEBCCACBAECBAEBECCAAEBABECECCDDAACDADADBACCBAADCAECCCEDD1241201251290CECDDACBCAEBEBBACEEDCAEBDAEDADBEBBDDADCBCBACDDEECEBABAADDEBAECDCACEDCBBEBCDCAADBADDDBAEACAEEDBDBDDCAEDDDACDEBBECAACDABADBAECCEADDCEBCBACBAEEBAEEECCEEEBCBECEBADCBDBEADCDBDCCCBADCAECCCCDDP 140.00 140.00 140.00 140.00 60.00 620.000 . . . ... 1: 3000000000022011 2214125506SAO JOSE DOS PINHAIS PR110000000000000000000011 . . ...4125506SAO JOSE DOS PINHAIS PR0000 . . . . 1211171251290DCEACDBECBCAEBACBEAECEDEDAABDBEDBBDDDCABCBDCAEDEBCEAEDBDAACACECDBBCACEDCDCBAEADADBBDDEAABAEEDBDBDDCAEDDDACDEBBECAACDABADBAECCEADDCEBCBACBAEEBAEEECCEEEBCBECEBADCBDBEADCDBDCCCBADCAECCCCDDF 0.00 0.00 0.00 0.00 0.00 0.000 . . . obrigado. 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

3 8

Re: [R-br] Digest R-br, volume 70, assunto 3rrrregghtgrgfffdf
by Josmar Mazucheli 02 Out '16

02 Out '16

8ttgdftfese Em 02/10/2016 12:00, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: > 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: Bando de dado .txt (Andre Oliveira) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sun, 2 Oct 2016 12:03:27 +0000 (UTC) > From: Andre Oliveira <andreolsouza(a)yahoo.com.br> > To: Marcos Silva <marcosfs2006(a)gmail.com> > Cc: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Bando de dado .txt > Message-ID: <1862678262.297135.1475409807537(a)mail.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > Oi Marcos,obrigado pelas sugestões. Eu estava indo muito bem com a > data.table, tentei library(ff) e library(ffbase), mas sem sucesso. > > 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 > > Em Domingo, 2 de Outubro de 2016 9:00, Marcos Silva < > marcosfs2006(a)gmail.com> escreveu: > > > Dá uma olhada nesse post: > http://rcoster.blogspot.com.br/2014/02/lendo-grandes-bancos-de-dados.html > > Em 2 de outubro de 2016 08:36, Marcos Silva <marcosfs2006(a)gmail.com> > escreveu: > > André, > Este tipo de arquivo não tem um caractere delimitador, tal como virgula, > ponto e virgula, pipe ("|"), Você não vai conseguir importar com > read.table() e derivações (read.csv(), read.csv2()) ou mesmo com fread(). > Mas para vc "ler" esse arquivo, vc precisa saber qual o comprimento de > cada variável (quantidade de caracteres que ocupa no arquivo), informação > essa que deveria estar no dicionário de variáveis. > > Este tipo de arquivo vc importa com read.fwf() ou read_fwf() do pacote > readr. Mas como essa base é muito grande, recomendo utilizar algum pacote, > a exemplo do pacote LaF. Existem outras possibilidades talvez melhores. > > > > > > Em 2 de outubro de 2016 03:13, Andre Oliveira <andreolsouza(a)yahoo.com.br> > escreveu: > > Oi Marcos,neste dicionário fala que é formato em formato ASCII. Não tem > mais instruções pelo que percebi. Descreve os nomes da varáveis também! > 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 > > Em Domingo, 2 de Outubro de 2016 0:13, Marcos Silva < > marcosfs2006(a)gmail.com> escreveu: > > > Acho que este e um arquivo de formato fixo. Vc vai precisar do dicionario > de dados para importar. Normalmente este dicionario e disponibilizado com > os dados pelo inep. > Em 01/10/2016 22:31, "Andre Oliveira via R-br" <r-br(a)listas.c3sl.ufpr.br> > escreveu: > > Boa noite,não estou conseguindo identificar o motivo de não abrir este > arquivo .txt corretamente, ele tem em torno de 6.3 GBAlguém teria uma dica > para compartilhar? Acredito que não estou trabalhando corretamente! > > library(data.table) > readLines("DADOS_ENEM_2011. TXT", n = 1) > > [1] "3000000000012011 2803548708SAO BERNARDO DO > CAMPO > > SP110000000000000000000011 . . > > > ...3548708SAO BERNARDO DO CAMPO > > SP1111 577.10 576.00 > 575.50 737. 10CEBDCACBEAABBAABAEBDCABDCDEA AADCBBEBADDBEBECDDEECABABEBDDE > BAEADCACABCBBECADDAADBADDAABAB CAEEAEAEDDDAAEEBDECADCEABDDBAE > CAAABDAEBCCACBAECBAEBECCAAEBAB ECECCDDAACDADADBACCBAADCAECCCE > DD1241201251290CECDDACBCAEBEBB ACEEDCAEBDAEDADBEBBDDADCBCBACD > DEECEBABAADDEBAECDCACEDCBBEBCD CAADBADDDBAEACAEEDBDBDDCAEDDDA > CDEBBECAACDABADBAECCEADDCEBCBA CBAEEBAEEECCEEEBCBECEBADCBDBEA > DCDBDCCCBADCAECCCCDDP 140.00 140.00 140.00 140.00 60.00 > 620.000 . . . > > > require(bit64) > enem2011<- fread("DADOS_ENEM_2011.TXT", nrow = 1, header = TRUE, sep=";") > head(enem2011) > > 3000000000012011 2803548708SAO BERNARDO DO CAMPO > > SP110000000000000000000011 > . . > > ...3548708SAO BERNARDO DO > CAMPO > > SP1111 577.10 576.00 575.50 737. 10CEBDCACBEAABBAABAEBDCABDCDEA > AADCBBEBADDBEBECDDEECABABEBDDE BAEADCACABCBBECADDAADBADDAABAB > CAEEAEAEDDDAAEEBDECADCEABDDBAE CAAABDAEBCCACBAECBAEBECCAAEBAB > ECECCDDAACDADADBACCBAADCAECCCE DD1241201251290CECDDACBCAEBEBB > ACEEDCAEBDAEDADBEBBDDADCBCBACD DEECEBABAADDEBAECDCACEDCBBEBCD > CAADBADDDBAEACAEEDBDBDDCAEDDDA CDEBBECAACDABADBAECCEADDCEBCBA > CBAEEBAEEECCEEEBCBECEBADCBDBEA DCDBDCCCBADCAECCCCDDP 140.00 140.00 > 140.00 140.00 60.00 620.000 . . > . > > ... > 1: 3000000000022011 2214125506SAO JOSE DOS PINHAIS > > > PR110000000000000000000011 . . > > > ...4125506SAO JOSE DOS PINHAIS > > PR0000 . > . . . > > > 1211171251290DCEACDBECBCAEBACB > EAECEDEDAABDBEDBBDDDCABCBDCAED EBCEAEDBDAACACECDBBCACEDCDCBAE > ADADBBDDEAABAEEDBDBDDCAEDDDACD EBBECAACDABADBAECCEADDCEBCBACB > AEEBAEEECCEEEBCBECEBADCBDBEADC DBDCCCBADCAECCCCDDF 0.00 0.00 > 0.00 0.00 0.00 0.000 . . > . obrigado. > > > 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 > ______________________________ _________________ > 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. > > > > > > > > -- > Marcos F. Silva > http://sites.google.com/site/ marcosfs2006 > > > > -- > Marcos F. Silva > http://sites.google.com/site/marcosfs2006 > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20161002/24cbb345/attachment-0001.html> > > ------------------------------ > > 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 70, assunto 3 > **************************************** >

1 0

Título de gráficos em mais de uma linha
by Guilherme Biz 30 Set '16

30 Set '16

Boa tarde. Estou tentando fazer um gráfico com título extenso e gostaria que fosse dividido em duas ou mais linhas de forma automática, sem ser pelo comando \n. Ex: a=rnorm(100) b=rnorm(100) plot(a,b,main="Este título é muito grande para ficar só em uma linha e gostaria de dividi-lo em duas ou mais linhas de forma automática") Sei que usando "\n" eu passo para a linha de baixo e no Rstudio se eu escrever o título em duas linhas também ocorre isso, mas preciso de forma automática, pois são muitos gráficos e estou fazendo de forma iterativa. Alguém tem uma solução? Agradeço desde já. Att, Guilherme Biz

3 3

Re: [R-br] Parcela Subdividida
by Angélica Ricarte 27 Set '16

27 Set '16

Cesar Rabak e demais colegas, boa tarde. Agradeço o retorno e peço desculpas pela demora em responder esta postagem. Respondendo os questionamentos: - Por que então que você constrói (especifica) o seu termo Error como PARC:REP? R: o modelo da anova que eu utilizei [mod<-aov(var~PARC*SubP+Error(PARC:REP))] retirei de outras postagens feitas no fórum e também o vi em outras consultas, portanto, não sei o porquê do "Error como PARC:REP" (e gostaria de entender), mas entendo que seria o erro referente às parcelas ('Erro a' do resultado da anova). Em todas as minhas consultas o modelo para parcela subdividida é sempre este, mas comigo não funcionou (ou melhor, funcionou, mas apareceu uma mensagem de erro, como mostrado anteriormente). Entendo que se eu não inserir o termo do 'Error', a análise será gerada considerando um esquema fatorial. - Ou fazendo uma pergunta mais próxima da teoria da análise que você quer fazer (uma ANOVA com medidas repetidas), qual é o erro "dentro" que você quer levar em conta na sua ANOVA? R: como dito anteriormente, o erro referente à Parcela (Erro a), visto que no resultado da anova devo ter dois Erros (a e b). - Outra perguntinha: você entende que ao especificar o termo erro você reduziu o número de pontos para ajuste para dez e esse é o número de fatores PARC vezes o número de parcelas SubP? R: não consigo compreender isto. Se você, ou alguém aqui do fórum, puder me ajudar, agradeço muito, estou ansiosa para resolver este problema. Se eu estiver errada em algum ponto, por favor, corrijam-me. Se o problema está no modelo, qual seria o que eu poderia utilizar para parcela subdividida (PS) (tenho um experimento em DIC, com parcelas subdividas (2 x 5) e 4 repetições)? Em suma, meu objetivo é gerar o resíduo (que não sei se utilizei o comando correto, visto que é uma PS) para poder realizar o teste de Shapiro-Wilk (SW), pois o quadro da anova que utilizo é aquele do pacote 'ExpDes.pt' (e, infelizmente, para o caso de PS este pacote não apresenta o resultado do teste de SW). Assim, se houver uma solução ainda mais simples para que possa atingir este objetivo, peço, por gentileza, que me apresentem. Muito obrigada. Aguardo retorno. Att., Angélica. ________________________________ De: Cesar Rabak <cesar.rabak(a)gmail.com> Enviado: sexta-feira, 9 de setembro de 2016 16:14 Para: Angélica Ricarte Assunto: Re: [R-br] Parcela Subdividida Olá Angélica, CMR (que está descrito no FAQ deste grupo) é "código mínimo reproduzível". A ideia é permitir que a gente refaça seus passos e acelere a ajuda. Seu código acima não é exatamente reproduzível porque esta linha: > PS<-read.xlsx(file.choose(),sheetName="ParcSub_DIC",h=T) diz respeito a um arquivo que está na sua máquina. Contudo, como você imprimiu o total do dataframe importado co o comando seguinte, está quase! OK, então você tem uma resposta "var" para o fator PARC e aninhado em PARC você tem quatro subparcelas, com quatro repetições por medida. Por que então que você constrói (especifica) o seu termo Error como PARC:REP? Ou fazendo uma pergunta mais próxima da teoria da análise que você quer fazer (uma ANOVA com medidas repetidas), qual é o erro "dentro" que você quer levar em conta na sua ANOVA? Outra perguntinha: você entende que ao especificar o termo erro você reduziu o número de pontos para ajuste para dez e esse é o número de fatores PARC vezes o número de parcelas SubP? HTH -- Cesar Rabak 2016-09-08 14:27 GMT-03:00 Angélica Ricarte <angelica.ricarte(a)hotmail.com<mailto:angelica.ricarte@hotmail.com>>: Boa tarde, Cezar Rabak. Entendi melhor, porém ainda não por completo. Mas, o meu problema agora é que não consigo mais obter o mesmo resultado. Estou importando exatamente os mesmos dados e utilizando exatamente os mesmos comandos, porém não consigo mais o mesmo resultado que enviei no fórum dia 04/09/16. Quando realizo a ANOVA o resultado gerado está sendo diferente (e está errado) daquele que apresentei anteriormente (este está correto). Gostaria de saber como posso obter o "resíduo" para inserir no comando para o teste de Shapiro Willk: "shapiro.test(resíduo)". Seguem abaixo os comandos utilizados (acho que é o 'CMR' a que você se referiu) e o dados, alguém poderia me ajudar a encontrar o meu erro e mostrar uma solução? Já procurei muito, mas não encontro. _____________________________________________________________ > PS<-read.xlsx(file.choose(),sheetName="ParcSub_DIC",h=T);PS PARC SubP REP var 1 C 30 1 12 2 C 30 2 13 3 C 30 3 14 4 C 30 4 12 5 C 45 1 18 6 C 45 2 17 7 C 45 3 16 8 C 45 4 17 9 C 60 1 22 10 C 60 2 21 11 C 60 3 24 12 C 60 4 23 13 C 75 1 24 14 C 75 2 23 15 C 75 3 22 16 C 75 4 24 17 C 90 1 21 18 C 90 2 20 19 C 90 3 21 20 C 90 4 19 21 S 30 1 22 22 S 30 2 21 23 S 30 3 24 24 S 30 4 23 25 S 45 1 21 26 S 45 2 22 27 S 45 3 23 28 S 45 4 21 29 S 60 1 20 30 S 60 2 18 31 S 60 3 19 32 S 60 4 20 33 S 75 1 17 34 S 75 2 15 35 S 75 3 16 36 S 75 4 15 37 S 90 1 12 38 S 90 2 13 39 S 90 3 12 40 S 90 4 11 > attach(PS) > SubP<-as.factor(SubP);is.factor(SubP) [1] TRUE > mod<-aov(var~PARC*SubP+Error(PARC:REP));summary(mod) Error: PARC:REP Df Sum Sq Mean Sq F value Pr(>F) PARC 1 6.453 6.453 80.67 0.0706 . Residuals 1 0.080 0.080 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Error: Within Df Sum Sq Mean Sq F value Pr(>F) PARC 1 1.7 1.67 1.535 0.226 SubP 4 108.9 27.23 25.076 6.64e-09 *** PARC:SubP 4 496.9 124.22 114.418 < 2e-16 *** Residuals 28 30.4 1.09 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 > resíduo<-residuals(mod);resíduo NULL ___________________________________________________________________ Aguardo retorno. Obrigada. Att., Angélica. ________________________________ From: cesar.rabak(a)gmail.com<mailto:cesar.rabak@gmail.com> Date: Tue, 6 Sep 2016 21:49:25 -0300 Subject: Re: [R-br] Parcela Subdividida To: angelica.ricarte(a)hotmail.com<mailto:angelica.ricarte@hotmail.com>; r-br(a)listas.c3sl.ufpr.br<mailto:r-br@listas.c3sl.ufpr.br> Angélica, Como você deu nomes "crípticos" para seu dataframe, e não postou um CMR, dá só para fazer uma 'interpretação' das mensagens. Aparentemente você quer a resposta PS[,4] para a interação entre as variáveis (grupos[?]) PS[,1] e PS[,2] e considera a combinação dos fatores PS[,1] e PS[,3] como medidas repetidas e tem 28 casos ao todo. Além do mais você considera o fator PS[,1] como resposta e Error() ou seja parte da especificação das medidas repetidas. A mensagem sobre o modelo Error() ser singular ocorre quando você tenta ajustar mais parâmetros que seus dados permitem... Se você continuar com dificuldades, você precisa explicar melhor o que deseja fazer ou consultar uma obra sobre ANOVA com medidas repetidas. HTH -- Cesar Rabak 2016-09-05 21:05 GMT-03:00 Angélica Ricarte via R-br <r-br(a)listas.c3sl.ufpr.br<mailto:r-br@listas.c3sl.ufpr.br>>: Prezado Cesar, Não consegui compreender a mensagem de erro e nem o porquê que minha matriz é nula. Você pode ajudar? Agradeço e aguardo retorno. Obrigada. Att., Angélica ________________________________ Date: Mon, 5 Sep 2016 12:52:41 -0700 From: ml-node+s2285057n4666617h13(a)n4.nabble.com<mailto:ml-node%2Bs2285057n4666617h13@n4.nabble.com> To: angelica.ricarte(a)hotmail.com<mailto:angelica.ricarte@hotmail.com> Subject: Re: [R-br] Parcela Subdividida Angélica, Você entendeu a mensagem de erro? Não é possível analisar uma matriz nula.... 2016-09-04 1:22 GMT-03:00 Angélica Ricarte via R-br <[hidden email]>: Prezados, Estou tentando realizar uma análise de variância para um experimentos em parcelas subdivididas em DIC, mas, apesar de ser gerado um resultado, aparece uma mensagem de erro. Quando tento realizar o cálculo dos resíduos também dá erro e, consequentemente, não consigo fazer o teste de normalidade. Alguém pode me ajudar? Vejam como realizei as análises e os seus resultados: #ANOVA: > mod<-aov(PS[,4]~PS[,1]*PS[,2]+Error(PS[,1]:PS[,3]));summary(mod) Warning message: In aov(PS[, 4] ~ PS[, 1] * PS[, 2] + Error(PS[, 1]:PS[, 3])) : modelo Error() é singular Error: PS[, 1]:PS[, 3] Df Sum Sq Mean Sq F value Pr(>F) PS[, 1] 1 8.1 8.100 11.3 0.0152 * Residuals 6 4.3 0.717 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Error: Within Df Sum Sq Mean Sq F value Pr(>F) PS[, 2] 4 108.9 27.23 24.94 3.02e-08 *** PS[, 1]:PS[, 2] 4 496.9 124.23 113.79 3.09e-15 *** Residuals 24 26.2 1.09 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 #Resíduos: > residuals(mod) NULL #Teste de Shapiro-Wilk: > res<-residuals(mod) > shapiro.test(res) Error: is.numeric(x) is not TRUE Obrigada. Att., Angélica. _______________________________________________ R-br mailing list [hidden email] 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 [hidden email] 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-Parcela-Subdividida-tp4666613p466661… 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…> _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br<mailto: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.

1 0

Salvar gráfico no Linux
by Alexandre Loures 26 Set '16

26 Set '16

Boa tarde! Para conseguir visualizar os gráficos preciso maximizar a janela. Porém, não estou sabendo salvar a janela maximizada no Linux. Segue pdf (salvo no windows) em anexo para exemplificar. Alguém poderia me ajudar? 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 Site: www.ccsa.ufpb.br/ppge www.orcid.org/0000-0002-1288-0135

3 2

Precipitação diária baseada em uma rede neural
by Alisson Lucrécio 25 Set '16

25 Set '16

Caro Colegas, Boa tarde. Estou tentando calcular a ocorrência de dias secos e chuvosos e a quantidade de precipitação diária através de uma rede neural. Algum colega teria alguns exemplo que poderia me ajudar? Ou algum livro que poderia recorrer? Obrigado. -- Alisson Lucrecio da Costa

2 1

Extração de dados climáticos
by ana paula coelho madeira 23 Set '16

23 Set '16

________________________________ Prezados colegas, bom dia! Então, é o seguinte, no site http://berkeleyearth.org/ eu tenho dados de temperatura média de toda a terra. Eu preciso extrair de acordo com os meus pontos geográficos, a temperatura média anual para cada ponto. Só que eu tenho os específicos que eu preciso. Comecei o script de forma que eu pudesse extrair para um ponto específico e para todos os anos, para depois arrumar o script para fazer para todos os meus pontos e anos de interesse. Segue o script utilizado! library(raster) b <- brick("D:/Documentos/Desktop/analises_cap2/dados_climaticos/precipt.nc", varname <- "precip") print(b) b idx <- getZ(b) coords <- matrix(c(-55, -10), ncol = 2) vals <- extract(b, coords, df = T) df <- data.frame(idx, t(vals)[-1,]) rownames(df) <- NULL names(df) <- c('date','value') head(df) Seguindo o script no local onde deveria apresentar as variáveis extraídas , aparece NA. Alguém já fez este tipo de extração? Peço ajuda. Desde já, muito obrigada. Ana Paula

2 1

"Append" arquivos csv de diversos subdiretórios
by Alexandre Loures 22 Set '16

22 Set '16

Bom dia! Tenho um diretório com 50 subdiretórios sendo que cada um possui 11 arquivos csv. O que preciso é "empilhar" esses arquivos (em um único arquivo) e para isso estou fazendo o seguinte: filenames <- list.files () mydata <- do.call ("rbind", lapply (filenames, read.csv, header = TRUE)) Entretanto, gostaria de "automatizar" esse script, ou seja, o próprio comando fosse em cada subdiretório e fizesse o "append" (dos 550 arquivos csv) dando o resultado final. Alguém poderia me ajudar? Desde já muito obrigado!

4 4

Salvar uma tabela/matriz como imagem
by Elias Carvalho 20 Set '16

20 Set '16

Boa tarde a todos. Eu tenho uma matriz, chamada "matrix.results", como resultado de algumas análises: Consigo imprimir uma tabela no console com o comando: kable(head(matrix.results ), format = "rst", align = "r") Agora eu gostaria de salvar uma imagem desta tabela para inserir em um documento. No entanto eu gostaria de salvar como imagem uma tabela mais bem formatada, como esta: Alguém poderia ajudar ? -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho*

3 2

Re: [R-br] Distribuição exponencial mista
by Alisson Lucrécio 20 Set '16

20 Set '16

Jaildo Rocha, Muito obrigado por sua ajuda. Alisson 2016-09-20 20:53 GMT-03:00 Alisson Lucrécio < alisson.lucrecio(a)ifgoiano.edu.br>: > Jaildo Rocha, > > Muito obrigado por sua ajuda. > > Alisson > > 2016-09-20 19:27 GMT-03:00 Jaildo Rocha <j.vieirarocha(a)gmail.com>: > >> Prezado Alisson, >> >> Para realizar as estimativas dos parâmetros via o package “Renext” >> utilize a função “EM.mixexp”. Apesar disso, é importante notar dois >> detalhes desta distribuição; >> >> >> - O primeiro é o caso do coeficiente de variação maior ou igual a um >> que tende a gerar “erros” (convergir para uma exponencial simples) quando a >> amostra apresenta tal coeficiente menor que um. >> >> >> >> - Um segundo detalhe é que as estimativas dos parâmetros,no que diz >> respeito a convergência em relação a amostra, é bem mais “lenta” com o >> crescimento da amostra do que distribuições usuals como a gama, weibull, >> gumbel e outras (o segundo detalhe é apenas um feeling particular). >> >> >> Caso não queira usar o “Renext”, outra possibilidade é utilizar o pacote >> “fitdistplus” para a realização das estimativas dos parâmetros, o único >> inconveniente é a necessidade de inserir um valor inicial.Por fim, segue os >> exemplos de estimativas utilizando os dois pacotes diferentes. >> >> >> ##### >> >> # Via Renext >> >> library("Renext") >> mix <- rmixexp2(n = 1000, prob1 = 0.5, rate1 = 1.0, rate2 = 3.0) >> fit_mixexp <- EM.mixexp(x = mix) >> >> # Parâmetros >> fit_mixexp[[1]] >> >> >> # Via fitdistrplus >> >> library("Renext") >> library("fitdistrplus") >> >> fit_mixexp2 <- fitdist(mix, "mixexp2", start = list(prob1 = 0.8,rate1 = >> 1, rate2 = 2)) >> fit_mixexp2[[1]] >> >> >> #### >> >> Obs: No segundo exemplo, apesar de implícito, utilizei as funções >> pmixexp2, dmixexp2, qmixexp2 do pacote “Renext” mas caso seja de sua >> preferência é possível você mesmo defini-las. >> >> Att, >> ------- >> Jaildo Rocha >> MSc. pelo Programa Pós Graduação em Saneamento Meio Ambiente e Recursos >> Hídricos – UFMG >> Área de concentração: Recursos Hídricos >> >> >> Em 20 de set de 2016, à(s) 00:21, Alisson Lucrécio via R-br < >> r-br(a)listas.c3sl.ufpr.br> escreveu: >> >> Cesar, >> >> Preciso solucionar essa equação f(x) = (alpha/beta1*exp(-x/beta1)) + >> ((1-alpha)/beta2*exp(-x/beta2)) para obter os paramentos da >> distribuição. Assim como solucionado pelo fitdistr do pacote MASS para >> as ditribuções gamma, webull e exponencial. Observe que fiz uma tentativa >> com o pacote Renext, mas sem sucesso. >> >> Muito obrigado. >> >> 2016-09-19 16:12 GMT-03:00 Alisson Lucrécio <alisson.lucrecio@ifg >> oiano.edu.br>: >> >>> Como solucionar essa equação f(x) = (alpha/beta1*exp(-x/beta1)) + >>> ((1-alpha)/beta2*exp(-x/beta2)) para obter os paramentos da >>> distribuição. Assim como é feito pelo fitdistr do pacote MASS para as >>> ditribuções gamma, webull e exponencial. Observe que fiz uma tentativa com >>> o pacote Renext, mas sem sucesso. >>> >>> Muito obrigado. >>> >>> >>> On 19 September 2016 at 16:05, Cesar Rabak <cesar.rabak(a)gmail.com> >>> wrote: >>> >>>> . . . e sua dúvida seria? >>>> >>>> On Mon, Sep 19, 2016 at 10:54 AM, Alisson Lucrécio via R-br < >>>> r-br(a)listas.c3sl.ufpr.br> wrote: >>>> >>>>> Caro Colegas, >>>>> >>>>> Bom dia. >>>>> >>>>> Algum colega da poderia me ajudar na solução de uma distribução >>>>> exponencial mista. >>>>> >>>>> # Script >>>>> >>>>> library(MASS) >>>>> >>>>> # Distribuição Gamma >>>>> >>>>> rg <- rgamma(100, shape = 1.5, rate = 1) >>>>> >>>>> fit_gamma <- fitdistr(rg,"gamma") >>>>> >>>>> par(bg="white", las=1, cex=1.1) >>>>> plot(density(rg, bw=0.5, cut=0), las=1, lwd=2, >>>>> xlim=c(0,5),col="steelblue") >>>>> >>>>> # Distribuição Weibull >>>>> >>>>> rw <- rweibull(100, scale=1, shape=1.5) >>>>> >>>>> par(bg="white", las=1, cex=1.1) >>>>> plot(density(rw, bw=0.5, cut=0), las=1, lwd=2, >>>>> xlim=c(0,5),col="steelblue") >>>>> >>>>> fit_weibull <- fitdistr(rw,"weibull") >>>>> >>>>> # Distribuição Exponential >>>>> >>>>> ex <- rexp(100, rate = 1) >>>>> >>>>> par(bg="white", las=1, cex=1.1) >>>>> plot(density(ex, bw=0.5, cut=0), las=1, lwd=2, >>>>> xlim=c(0,5),col="steelblue") >>>>> >>>>> fit_exp <- fitdistr(ex, "exponential") >>>>> >>>>> plot(density(rw, bw=0.5, cut=0), las=1, lwd=2, >>>>> xlim=c(0,5),col="steelblue") >>>>> lines(density(rg, bw=0.5, cut=0),col="green") >>>>> lines(density(ex, bw=0.5, cut=0),col="red") >>>>> >>>>> # Distribuição Exponencial mista >>>>> >>>>> library(Renext) >>>>> >>>>> # f(x) = (alpha/beta1*exp(-x/beta1)) + ((1-alpha)/beta2*exp(-x/beta2)) >>>>> >>>>> mixexp <- rmixexp2(n = 100, prob1 = 0.5, rate1 = 1.0, rate2 = 3.0) >>>>> >>>>> fit_mixexp <- mom.mixexp2(mixexp) >>>>> >>>>> plot(density(mixexp, bw=0.5, cut=1), lwd=2, >>>>> xlim=c(0, 5),col="steelblue") >>>>> lines(density(ex, bw=0.5, cut=0),col="red") >>>>> >>>>> mod <- nls(mixexp ~ prob1 * dexp(x, rate1) + (1 - prob1) * dexp(x, >>>>> rate2), >>>>> start = list(prob1 = 1, rate1 = 1, rate2 = 3.0), control = >>>>> nls.control(maxiter = 100)) >>>>> >>>>> Obrigado. >>>>> -- >>>>> Alisson Lucrecio da Costa >>>>> >>>>> _______________________________________________ >>>>> 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. >>>>> >>>> >>>> >>> >>> >>> -- >>> Alisson Lucrecio da Costa >>> >> >> >> >> -- >> Alisson Lucrecio da Costa >> _______________________________________________ >> 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. >> >> >> > > > -- > Alisson Lucrecio da Costa > -- Alisson Lucrecio da Costa

1 0