Setembro 2016 - 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

servidor cloud para rodar o R + shiny
by Marcio B 20 Out '16

20 Out '16

Olá pessoal boa tarde. Estou pesquisando um serviço de cloud para utilizar o R + Shiny. Minha duvida esta no custo e por isso gostaria de saber se alguém pode me ajudar para entender melhor: Quem é a empresa (cloud) que esta usando? AWS, Azure, Google ou IBM Qual o tipo de serviço para configurar o R + Shiny? Nestes servidores oferecem vários tipos, o que mais encontro na web é o EC2 - AWS, como vários tutoriais de instalação e configuração Se puder nos informar o uso e o custo deste serviço? Principalmente no uso de Memoria e Storage de dados, algum parâmetro para temos uma noção. Somos uma startup e queremos oferecer nos serviços utilizando o R + Shiny. Porem gostaríamos de um parâmetro para tomar a decisão de configurar este serviço e começar a utilizar. Espero que fui claro e toda a ajuda será bem vinda. Obrigado desde já. Marcio. -- Marcio B. Cel.: 19 99648.3949 tim 99798.0104 vivo skype: maborbaa http://www.linkedin.com/in/marciobar https://www.facebook.com/marciobar Empresa http://www.techagr.com/ ==== PAZ E BEM ====

6 8

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 26 Set '16

26 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 23 Set '16

23 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 21 Set '16

21 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