Maio 2016 - R-br - DINF UFPR

[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

[OFF-TOPIC] R-consortium
by Felipe 30 Mai '16

30 Mai '16

Prezados colegas, o que acham deste post: https://www.r-consortium.org/news/blogs/2015/08/best-practices-using-r-secu… -- Atenciosamente Felipe E. Barletta Mendes Estatístico(UFPR) - Conre3 9766-A Mestrando em Bioestatística(UEM) +55 (41)-92077191 +55 (41)-33287216

2 1

Output de um loop em uma matriz vazia
by Yury Duarte 30 Mai '16

30 Mai '16

Boa tarde colegas programadores! Estou tendo dificuldades em acertar o formato que devo organizar os outputs de um looping que fiz. Tenho um banco de dados com datas e algumas variáveis climáticas. Dessas datas, preciso simular alguns cálculos para dias específicos. São 36 simulações por ano, durante um período de 34 anos, totalizando 1224 simulações. Para alcançar esse objetivo estou trabalhando com 3 comandos for, um dentro do outro. O primeiro seleciona as datas/posições que devem ser feitos os cálculos. O segundo executa os devidos cálculos para as datas selecionadas e o ultimo integra os valores calculados, gerando a informação de interesse. Acredito que a logica de cálculos dentro de cada for() está ok, mas a saída não está puxando os valores que o looping calcula. Segue o script no corpo do email e, em anexo, um exemplo de um ano do banco de dados que está sendo utilizado. Desde já, agradeço a colaboração dos colegas. Abs #--------------------Remover objetos salvos na memória do R--------------------# rm(list = ls()) #--------------------Seleção de diretorios e arquivos de interesse--------------------# a = read.table("x.txt", sep = "\t", header = T) #--------------------Calculo das Variaves Astronomicas--------------------# head(a) a$DATE = as.Date(levels(a$DATE))[a$DATE] ano = format(a$DATE, trim = T, '%Y'); a$ano=ano mes = format(a$DATE, trim = T, '%m') ; a$mes=mes dia = format(a$DATE, trim = T, '%d') ; a$dia=dia dece = ifelse(dia<11, 1,ifelse(dia<21,2,3)); a$decendio=dece lat = -7.53 corr = pi/180 decl = 23.45*sin(corr*((a$DiaJuliano-80)*360/365)); decl hn = 1/corr*acos(-tan(corr*lat)*tan(corr*decl)) a$N = 2*hn/15 s = (1+0.033*cos(a$DiaJuliano*360/365)) t = corr*hn*sin(corr*lat)*sin(corr*decl) u = cos(corr*lat)*cos(corr*decl)*sin(corr*hn) a$Qo = 37.6*s*(t+u) head(a) x = 0.25 b = 0.50 n.est = (a$SRAD/a$Qo-x)*a$N/b; n.est a$n<- ifelse((n.est)<0,1, n.est); head(a) a$n.N = a$n/a$N; head(a) #--------------------Transformar MJ em caloria--------------------# a$Srad.cal = round((a$SRAD*1000000)/41800,1); head(a) #round é para delimitar o numero de casas apos a virgula a$sem = paste0(dia, '-', mes) a$TMED = (((a$TMAX+a$TMIN))/2) head(a) #-------------------- Determinacao das datas de simulacao--------------------# ndc = 140 datas = 1:366 # referentes aos dia juliano dim(a) anos = table(a$ano); anos potencial = matrix(NA, nrow = 366, ncol = 34 ) rownames(potencial) = 1:366 colnames(potencial) = c(1980:2013) head(potencial); dim(potencial) head(a) for (i in (a$DiaJuliano)){ epocas = a[which( a$sem == '01-01'|a$sem == '11-01' |a$sem =='21-01'| a$sem == '01-02'|a$sem == '11-02' |a$sem =='21-02'| a$sem == '01-03'|a$sem == '11-03' |a$sem =='21-03'| a$sem == '01-04'|a$sem == '11-04' |a$sem =='21-04'| a$sem == '01-05'|a$sem == '11-05' |a$sem =='21-05'| a$sem == '01-06'|a$sem == '11-06' |a$sem =='21-06'| a$sem == '01-07'|a$sem == '11-07' |a$sem =='21-07'| a$sem == '01-08'|a$sem == '11-08' |a$sem =='21-08'| a$sem == '01-09'|a$sem == '11-09' |a$sem =='21-09'| a$sem == '01-10'|a$sem == '11-10' |a$sem =='21-10'| a$sem == '01-11'|a$sem == '11-11' |a$sem =='21-11'| a$sem == '01-12'|a$sem == '11-12' |a$sem =='21-12'& a$ano>=1980),] semear = length(epocas[,1]) for (j in length(semear)){ ini = as.numeric(rownames(epocas[j])) fim = (ini+ndc)-1 pac = a[(ini:fim),] #--------------------Inserindo as Produtividades Parciais--------------------# #----------Inserindo influencia da nebulosidade----------# pac$ctc[j] = -4.16+(0.4325*pac$TMED[j])-(0.00725*((a$TMED[j])^2)) pac$ctn[j] = -1.064+(0.173*a$TMED[j])-(0.0029*((a$TMED[j])^2)) pac$PPBc[j] = (107.2+(0.36*pac$Qo[j]*1000000)/41800)*pac$ctc[j]*pac$n.N[j] pac$PPBn[j] = (31.7+(0.219*pac$Qo[j]*1000000)/41800)*pac$ctn[j]*(1-(pac$n.N[j])) pac$PPB[j] = round((pac$PPBc[j])+(pac$PPBn[j]),1); head(pac) #--------------------Calculo da Produtividade - base seca --------------------# pac$c.resp[j] = ifelse(a$TMED[j] < 20, 0.6, 0.5) pac$c.colheita[j] = 0.5 umidade = 13 c.um = (1/(1-0.01*umidade)) pac$iaf[j] = 5 pac$c.iaf[j] = (0.0093+0.185*pac$iaf[j])-(0.0175*pac$iaf[j]^2) pac$PP.Liquida_seca = NA pac$PP.Liquida_seca[1] = round(pac$PPB[1]*pac$c.resp[1]*pac$c.colheita[1]*pac$c.iaf[1],1); head(pac) # em kg/ha for (w in 2:ndc){ pac$PP.Liquida_seca[w] = round(pac$PPB[w]*pac$c.resp[w]*pac$c.colheita[w]*pac$c.iaf[w]+ pac$PP.Liquida_seca[w-1],1) } #--------------------Soma da PP--------------------# pac$PP_final = round(pac$PP.Liquida_seca[w]*c.um,1) dim(pac) PP = pac[ndc,32] potencial[i,j] = PP } print(i) } View(pac) View(potencial) Yury Duarte Engenheiro Agrônomo - ESALQ/USP

1 0

aceitação
by Edimeire Alexandra Pinto 30 Mai '16

30 Mai '16

Bom dia. Obrigada pela aceitação

1 0

Contar casos em cada variável e em todas as variáveis
by Amiko Bh 28 Mai '16

28 Mai '16

Pessoal, Dado o banco: *ATO_1* *ATO_2* *ATO_3* ART 121 ART 155 ART 180 ART 33 ART 28 ART 10 ART 121 ART 15 ART 155 ART 155 ART 33 ART 121 ART 28 ART 33 ART 33 ART 121 ART 121 ART 121 ART 15 ART 121 ART 180 a) Preciso contar cada elemento em cada variável, para que vislumbre o seguinte resultado: *ATO_1* ART 121 = 3 ART 33 = 1 ART 155= 1 ART 28= 1 ART 15 = 1 e assim com *ATO_2 e **ATO_3.* b) preciso saber o total de cada elemento em todos as colunas, assim: ART 121 = 7 ART 33 = 4 ART 155= 3 ART 28= 2 ART 15 = 2 ART 180=2 ART 10=1 Muito obrigado.

4 5

Colocar legenda sem borda
by Andre Oliveira 27 Mai '16

27 Mai '16

Boa tarde,como posso tirar borda desta legenda? Tentei bty="n" mas não deu certo! B<-matrix(c(578,1289, 439,1061),ncol=2) colnames(B)=c("Feminino","Masculino") rownames(B)=c("População","IBGE") barplot(B, beside=T,legend.text=TRUE,col=c(4,"limegreen"), main="",ylim=c(0,2000)) 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 6

Como colocar barras justapostas com o barplot
by Andre Oliveira 27 Mai '16

27 Mai '16

Pessoal já tentei, mas não consegui! Queria algo do tipo! barplot(dadosXX,beside=T,col = heat.colors(21),space=1) dadosXX= data.frame(Faixa, TFD,TFDMG)dadosXX Faixa TFD TFDMG 1 0-4 anos 3.35 6.52 2 5-9 anos 4.61 7.35 3 10-14 anos 4.09 8.42 4 15-19 anos 5.03 8.37 5 20-24 anos 2.73 7.52 6 25-29 anos 4.19 7.36 7 30-34 anos 6.08 7.53 8 35-39 anos 4.51 7.06 9 40-44 anos 9.22 6.80 10 45-49 anos 7.13 6.72 11 50-54 anos 7.65 5.93 12 55-59 anos 11.11 5.13 13 60-64 anos 7.76 4.21 14 65-69 anos 8.07 3.41 15 70-74 anos 6.60 2.90 16 75-79 anos 3.67 2.22 17 80-84 anos 2.41 1.43 18 85-89 anos 1.15 0.65 19 90-94 anos 0.42 0.31 20 95-99 anos 0.10 0.14 21 100 anos ou mais 0.10 0.03 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

2 2

Moving Average, Modelos mistos
by Fernando Antonio de souza 27 Mai '16

27 Mai '16

Boa tarde, colegas Estou trabalhando com modelos mistos, com dados de medidas repetidas no tempo. Ao definir a estrutura da matriz de covariância, testei a estrutura autoregressiva de ordem 1 AR1(), Moving average ARMA(p=0,q=?), onde testei para ordem de 1 a 6 (acima de 6 o modelo não convergiu). Acontece que com a ordem 6 obtive o melhor modelo (segundo anova()), embora nem todas as correlações foram resolvidas, conforme o gráfico (anexo) Ainda não me sinto confortal em trabalhar com estas estruturas. Faz sentido uma ordem de 6 em ARMA? Vocês teriam algum material que me explique o uso destas estruturas? Atenciosamente ======================================================================= 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/ ========================================================================

3 3

Erro ao fazer download de PDFs: failed, no valid url links detected
by Elias Carvalho 27 Mai '16

27 Mai '16

Bom dia Pessoal Estou baixando tentando baixar alguns pdf de artigos científicos pelo DOI usando o pacote metagear, porém a maioria apresenta a mensagem "failed, no valid url links detected", ou "cannot open: HTTP status was '404 Not Found' PDF download... skipped", mesmo ligando o proxy da universidade. De 1800 artigos que preciso consegui por volta de uns 200 download apenas. Alguém poderia me dar uma ajuda em como resolver este problema ? -- Best regards... 8^) “The mind that is open to new ideas never come back to its original size” *Albert Einstein* -- Obrigado Elias

3 6