- R-br - DINF UFPR

Teste
by Rodrigo Coster 04 Jun '16

04 Jun '16

Alguém recebeu?

4 3

[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

ajuda com implementação de Programação Linear com o R - constrOptim
by Cleber N.Borges 26 Mai '16

26 Mai '16

Prezados, Bom dia! Como na lista tem pessoas com conhecimento e experiência, peço algumas orientações... :-) Imagino que o método de "Programação Linear (LP)" deve servir para resolver meu problema de achar os parâmetros de uma dada função objetivo. Sei que no R tem a opção nativa "constrOptim" (poderia ser outra, posso instalar pacotes) Entretanto eu não sei como elaborar minhas condições de contorno em formato matricial que é o modo como o comando constrOptim trabalha (até onde consegui entender) Agradeceria muito se alguém puder me ajudar, respondendo se eu realmente estou certo em pensar que a LP é a ferramenta certa para meu caso e /(ou que) me indique um material básico (introdutório). Abaixo segue um exemplo simplificado do tipos de função que pretendo obter parâmetros. Agradeço antecipadamente pela ajuda e por qualquer dica Cleber Borges ##################### # fobjetivo <- function( x = 0.1 ){ k1 <- 0.007585776 k2 <- 6.16595e-08 k3 <- 2.137962e-13 k4 <- 1e-14 k <- k1 * k2 * k3 * k4 # Determinar as 6 variáveis: a, b, c, d, e, f ################### # 6 Condições de contorno: ( Igualdades ) ################### k1 == a * d / c # 1 k2 == a * e / d # 2 k3 == a * f / e # 3 k4 == a * b # 4 4 * a == b + d + e + f # 5 x == c + d + e + f # 6 return( c( a, b, c, d, e, f ) ) } --- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus

3 4