Janeiro 2017 - 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

Diferença resultados agricolae::HSD.test e lsmeas
by Fernando Antonio de souza 05 Nov '18

05 Nov '18

Pessoal, poque o resultado dos testes HSD.test e lsmeans diferem neste exemplo? O resultado fornecido pelo lsmeas está em acordo com o material de onde tirei esse exemplo e que havia sido rodado no SAS. Já o HSD.test difere tanto no valor das médias estimadas como na significancia. install.packages(c("agricolae",lsmeans")) library(agricolae) library(lsmeans) COVARIANCIA<-read.csv(" https://www.dropbox.com/s/0uvwqnnvsdsvphc/Covariancia.csv?raw=1",head=TRUE) modeloCOVAR<-lm(GANHO~PESOINIC+DIETA,data=COVARIANCIA) anova(modeloCOVAR) HSD.test(modeloCOVAR,"DIETA",console=TRUE) lsmeans(modeloCOVAR,pairwise~DIETA,adjust="tukey") -- ========================================= Fernando Souza Zootecnista, DSc. Produção e Alimentação Animal Celular: (31)99796-8781 (Vivo) E-mail:nandodesouza@gmail.com <e-mail%3Anandodesouza(a)gmail.com> Lattes: http://lattes.cnpq.br/6519538815038307 Blog: https://producaoanimalcomr.wordpress.com/ ==========================================

3 2

Selecionar a extração de valores dos pixels em linhas a cada número de pixels arbitrário
by ASANTOS 03 Fev '17

03 Fev '17

Prezados Membros, Tenho um raster hipotético r: require(raster); require(sp) ## Criando um raster r <- raster(nc=10, nr=10) r <- setValues(r, round(runif(ncell(r))* 255)) Consigo extrair os valores dos pixels de interesse dadas as coordenadas centrais dos mesmos: x <- c(54,18,54,54) y <- c(27,81,9,63) pontos <- SpatialPoints(cbind(x,y)) cells <- cellFromXY(r, pontos); cells p.cen <- xyFromCell(r, cells); p.cen values <- extract(r,p.cen) dados<-cbind(p.cen,cells,values) No entanto, eu gostaria de encontrar um modo para extrair o valor dos pixels, no qual, eu selecionaria uma coordenada e a função realizaria a extração dos pixels de linhas verticais inteiras do raster a cada dois pixels à partir da coordenada dada, por exemplo? Alguém poderia me dar uma ideia de como chegar a esse resultado? Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 - ResearcherID: A-5790-2016 Researchgate: www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635 Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/ ======================================================================

3 5

Finalizando a lista de blogs brasileiro de Rq
by Marcos Vital 31 Jan '17

31 Jan '17

Oi, pessoal! Finalizei a primeira versão da lista de blogs brasileiros de R, e lancei como um site usando o GitHub, aqui: https://marcosvital.github.io/blogs-de-R-no-Brasil/ Agradeço aos que contribuíram com a lista, que continua aberta para adições de novos blogs (basta me escrever ou contribuir pelo github). Abraços! Marcos -- Marcos Vinícius Carneiro Vital Universidade Federal de Alagoas Instituto de Ciências Biológicas e da Saúde Setor de Biodiversidade e Ecologia

1 0

função glht{multcomp} comparaçãos função lme
by Fernando Souza 31 Jan '17

31 Jan '17

Caros colegas, Estou realizando um exercício de parcela subdividida utilizando a função lme{nlme}. Acontece que quando rodo a função summary.glht as comparações pairwise são feitas e eu recebo as seguintes mensagens warnings: Warning messages: 1: In RET$pfunction("adjusted", ...) : Completion with error > abseps 2: In RET$pfunction("adjusted", ...) : Completion with error > abseps 3: In RET$pfunction("adjusted", ...) : Completion with error > abseps 4: In RET$pfunction("adjusted", ...) : Completion with error > abseps Essas mensagens no entanto não aparecem quando defino outro teste na função summary.glht summary.glht(compara, test=adjusted("fdr")) Porque isso ocorre? Segue o CMR. att #install.packages(c("nlme","multcomp")) library(nlme) library(multcomp) PARCSUB<-read.csv("https://www.dropbox.com/s/y8qmhwau0i8euva/ExemploSub.csv?raw=1",head=TRUE) PARCSUB<-transform(PARCSUB,PARCELA=factor(PARCELA),BLOCO=factor(BLOCO),PASTO=factor(PASTO),MINERAL=factor(MINERAL)) modeloPARSUBlme<-lme(PRODLEIT~BLOCO+PASTO*MINERAL,random=~1|PARCELA,data=PARCSUB) compara<- glht(modeloPARSUBlme, linfct=mcp(PASTO="Tukey")) summary.glht(compara) ================================== Fernando Souza Zootecnista, DSc. Produção e Alimentação Animal celular: (31)99796-8781 (Vivo) / (31)97358-4685 (Tim) e-mail:nandodesouza@gmail.com Lattes: http://lattes.cnpq.br/6519538815038307 blog: https://producaoanimalcomr.wordpress.com/ =================================

2 1

barplot agrupado e empilhado
by Christina Grupioni 31 Jan '17

31 Jan '17

Olá pessoal, Estou precisando fazer um grafico de barras mixto, que mistura agrupado e empilhado. É possível? Qual pacote faz isso? De qualquer forma, sou usuária nova e não consegui fazer nem ao menos o agrupado sem o empilhado. A necessidade é a seguinte: a primeira barra será empilhada (variaveis: colhidos e não colhidos). As outras três barras agrupadas devem ser a proporção de semdanos, dleves e dgraves. Na tabela, o X.semdanos representa a porcentagem de semdanos. e o X.leves e o X.graves representam a porcentagem de danos também. Preciso fazer isso para os 8 tratamentos. Preferencialmente no mesmo plot, para efeito de comparação. Meu data frame é o seguinte: trat colhidos dleves dgraves naocolhidos carga X._leves X._graves cap ef semdanos X.semdanos 1 810 85 81 36 846 10,49 10 332,95 95,74 644 79,51 2 531 65 43 27 558 12,24 8,1 548,63 95,16 423 79,66 3 805 111 34 44 849 13,79 4,22 559,86 94,82 660 81,99 4 1248 192 56 15 1263 15,38 4,49 849,94 98,81 1000 80,13 5 557 49 27 19 576 8,8 4,85 436,82 96,7 481 86,36 6 485 111 85 112 597 22,89 17,53 270,65 81,24 289 59,59 7 530 45 34 25 555 8,49 6,42 577,23 95,5 451 85,09 8 613 194 132 48 661 31,65 21,53 246,54 92,74 287 46,82 Na escuta, -- CHRISTINA GRUPIONI Engenheira Agrícola e Ambiental- CREA BA: 84142 *OCA (Organização Cooperativa de Agroecologia**)* *Tecnologias Ecológicas e Sociais* (31) 8863-0005, (31) 3892-2236 *skype:chrisgrupioni@gmail.com <chrisgrupioni(a)gmail.com>*

2 1

Extrair Período de datas
by Amiko Bh 30 Jan '17

30 Jan '17

Prezados, Estou postando, por não ter conseguido encontrar a solução aqui na lista, a despeito de procurar... Tenho um banco de dados, alimentado diariamente, sendo uma das variáveis o dia de entrada (dd/mm/aaaa). Preciso extrair deste banco todas as entradas que se deram no período do carnaval... Como faço? Uso subset? Atenciosamente,

5 7

problema em rJava
by Christina Grupioni 30 Jan '17

30 Jan '17

Olá pessoal, Estou com problema no rJava. Está instalado, mas sempre que preciso de usar, ele dá a noticia do erro, da forma abaixo. Alguém sabe como resolver isso? Na escuta, chris -- Carregando pacotes exigidos: xlsx Carregando pacotes exigidos: rJava Error : .onLoad falhou em loadNamespace() para 'rJava', detalhes: chamada: fun(libname, pkgname) erro: No CurrentVersion entry in Software/JavaSoft registry! Try re-installing Java and make sure R and Java have matching architectures. In addition: Warning messages: 1: package ‘xlsx’ was built under R version 3.2.5 2: package ‘rJava’ was built under R version 3.2.5 Failed with error: ‘package ‘rJava’ could not be loaded’ CHRISTINA GRUPIONI Engenheira Agrícola e Ambiental- CREA BA: 84142 *OCA (Organização Cooperativa de Agroecologia**)* *Tecnologias Ecológicas e Sociais* (31) 8863-0005, (31) 3892-2236 *skype:chrisgrupioni@gmail.com <chrisgrupioni(a)gmail.com>*

3 4

Extraindo valores de interesse em um dataframe
by Yury Duarte 30 Jan '17

30 Jan '17

Bom dia colegas programadores! Há algum tempo venho trabalhando com o modelo de simulação DSSAT. Agora estou tentando organizar suas saídas de uma forma que facilite as analises. Para que possam entender melhor o formato das saídas, anexei dois links para download do output do modelo. Toda simulação é antecedida de um pequeno cabeçalho e cada uma delas representa uma data de plantio (primeira linha da coluna DOY) em um determinado ano (coluna YEAR). Dessa forma, dentro de cada arquivo eu tenho 36 simulações de plantio por ano e 34 anos de histórico. A informação que eu gostaria de extrair é exatamente a última linha de cada simulação, referente a décima coluna (GWAD). Como cada simulação tem um numero diferente de linhas (dias do ciclo da cultura em questão), fiquei sem ter como resgatar os valores da coluna de interesse em intervalos fixos. Infelizmente não tenho nenhum script decente em andamento que possa ser compartilhado para elucidar melhor minhas intenções, mas espero ter conseguido explicar minimamente a minha dúvida. Desde já, agradeço a todos! link para a saída parcial: https://drive.google.com/file/d/0B6tnaf2fmoCcbUtPQ0NfUmpxdXM/ view?usp=sharing link para a saída completa: https://drive.google.com/file/d/0B6tnaf2fmoCcWjlFclc2R2NkVEU/ view?usp=sharing Yury Duarte Engenheiro Agrônomo - ESALQ/USP

3 6