Maio 2020 - 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

Duvidas Análise de Componentes Principais
by Fernando Souza 09 Dez '20

09 Dez '20

Olá pessoal tudo bem Estou realizando uma análise de componentes principais, com o objetivo de selecionar as melhores variáveis para uma análise multivariada. Gostaria de uma orientação sobre a interpretação. Eu interpretei de forma correta? Variáveis que possuem alta correlação com um CP também estão correlacionadas entre si. Isso significa que posso utilizar somente uma das principais variáveis do CP1 (por exemplo) para incluir no modelo multivariado? No caso a variável que apresentar significado para o trabalho. Os círculos de correlação, screem plot e corrplot geradas pelo Pca estão disponíveis no link https://drive.google.com/open?id=1jS5tCX-2-qHo2Gbej00DkFQBpLztwXks A interpretação que fiz foi a seguinte: Considerei as variáveis PC,IC e Prod leite (Azul no gráfico) como variável suplementar quantitativa. Baseado no ScreenPlot selecionei somente 5 CP para avaliação que explicaram 67,80% da variancia total. Pela minha interpretação cheguei as seguintes conclusões Principais variáveis para cada componente: Consumos (CP1): CMS,CMO,CPB,CEE,CCNF,CFDNcp,CFDA,CNDT,CEB,CED,CEM,CEL Energia (CP2) : EB,ED,EM,EL Peso (CP3): PV,PCJ,PCVZ,PCVZ075 Composição do Leite(CP4) : Solido e DCNF Variaveis sanguíneas (CP5): PT e ALBUMINA Segue uma descrição do trabalho para compreensão. Aceito sugestões sobre a melhor forma trabalhar esses dados O banco de dados é oriundo dois estudos realizados pelo mesmo grupo de pesquisa nas mesmas condições experimentais (animais,instalação,dieta,manejo) e avaliaram 48 variáveis. Os fatores experimentais foram: *Raça*: Estudo 1 avaliou 12 animais (6 raça1 e 6 da raça 2). O estudo 2 avaliou 13 animais (6 raça 1 e 7 raça 2) *Plano nutricional:* O estudo 1 avaliou os níveis (0,15,30) o estudo 2 avaliou os níveis (0,5,10,20) *Dias em lactação-DEL* (medidas repetidas no tempo): Estudo 1, dias dispersos entre os períodos de 30 a 120 dias. No estudo 2 os DEL foram avaliados em dias entre 150 a 250 dias. Eu realizei uma análise gráfica e não há indicações de haver diferenças entre estudos, o que era de certo modo esperado. Raça e DEL parecem ser os fatores que afetam as variáveis respostas. -- ========================================= 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/ ==========================================

6 14

Comparar 2 rasterbricks usando pacote TSdist
by Jackson Rodrigues 28 Mai '20

28 Mai '20

Pessoal Bom Dia! Estou com uma demanda aqui e não estou consigo enxergar onde estou errando. Preciso comparar 2 rasterbricks usando a função de dissimilaridade CorDistance do pacote TSdist. No entanto ela não funciona com NA's. Eu tentei remover ou ignorar os NA's sem sucesso. Alguém poderia me dar uma luz sobre como proceder? Abaixo está uma cópia do código que estou trabalhando. Saudações, Jackson ################### library(raster) library(TSdist) set.seed(12) r <- raster(nrow=10, ncol=10) s <- lapply(1:200, function(i) setValues(r, rnorm(ncell(r), sample.int(5,1), 0.5))) s <- stack(s) s[s < 0] <- NA s2<-s^2 # Visualizando que alguns pixels são NA hasna <- stackApply(s, indices = 1, fun = function(x, na.rm){anyNA(x)}) hasna2 <- stackApply(s2, indices = 1, fun = function(x, na.rm){anyNA(x)}) par(mfrow=c(1,2)) plot(hasna);plot(hasna2) # CorDistancenão funciona com NA CorDistance(as.vector(s[1]),as.vector(s2[1])) [1] 0.2192382 CorDistance(as.vector(s[2]),as.vector(s2[2])) <simpleError in .common.ts.sanity.check(x): NA in the series> [1] NA fun <- function(x,y){ out <- ifelse(anyNA(x) | anyNA(y), #Preciso resolver esta questão yes = NA, no = unname(CorDistance(x,y))) return(out) } # Aqui funciona fun(as.vector(s[1]),as.vector(s2[1])) [1] 0.2192382 # Aqui não funciona fun(as.vector(s[2]),as.vector(s2[2])) [1] NA new_s1 <- fun(s,s2) <simpleError in if (any(is.na(x))) { stop("NA in the series")}: argument is not interpretable as logical> # Veja que não funciona com calc new_s2 <- calc(s,s2, fun = fun) Error in .calcTest(x[1:5], fun, na.rm, forcefun, forceapply) : cannot use this function plot(s_out) ############# -- Jackson M. Rodrigues

1 0

Duvida teste de correlação
by Rodrigo Zucaratto 27 Mai '20

27 Mai '20

Prezados, Estou com uma duvida em relação ao teste de correlação. Quando faço uma correlação, não ha uma relação de causa e efeito (até aqui entendo). No entanto, estou querendo avaliar a inluência de variáveis ambientais (temperatura e pluviosidade) na fenologia de uma palmeira. To fazendo a correlação usando a função "cor.test". Um dos argumentos na função é method=pearson (normalidade) ou method=sperman (dados sem normalidade). Gostaria de saber como testo a normalidade em uma correlação, alguem pode ajudar? Abraços *Rodrigo Zucaratto* Bacharel em Ciências Biológicas Mestre em Ciências - Conservação da Natureza Doutor em Ecologia e Evolução Laboratório de Ecologia de Mamíferos - Universidade do Estado do Rio de Janeiro - UERJ Laboratório de Ecologia e Conservação de Florestas - LECF - UFRuralRJ Currículo: http://lattes.cnpq.br/2138780069263641 BSc - Biology MSc - Nature Conservation PhD - Ecology and Evolution Mammalian Ecology Lab - Rio de Janeiro State University - UERJ Ecology and Conservation Forests Lab - LECF - Federal Rural University of Rio da Janeiro Curriculum: http://lattes.cnpq.br/2138780069263641 https://scholar.google.com/citations?view_op=list_works&hl=en&user=on_Y4WYA…

3 4

Expoente de Lyapunov: horizonte de predição
by Nilson Guiotoku 27 Mai '20

27 Mai '20

Pessoal, quero encontrar o horizonte de previsões, usando o expoente de Lyapunov. O exemplo proposto no R, a serie de lorenz.ts funciona. Porém, quando troco pela minha variável que estou predizendo, da erro. Gostaria de entender qual parâmetro estou colocando errado. Segue o código mínimo reproduzível. #===================================== set.seed(22) y_obs <- 1:100 + rnorm(100,mean=0, sd=5) # supondo essa ser a variavel observada (100 observações temporais) y_estimado <- 1:80 # supondo essa ser a variavel estimada, até o tempo 80 y_pred <- 81:100 + rnorm(20,mean=0, sd=3) # suponha que usando uma rede recorrente chegue nesse resultado, como codigo minimo reproduzido, #simplifiquemos com uma rnorm com tendencia linear, para predição usando 80 tempos de predição plot(y_obs) points(1:80,1:80, type = "l") points(81:100,y_pred, type = "l", col="blue") legend("topleft", c("Observado","Estimado", "Predito"), col=c("black","black", "blue"), pch = c(1,NA, NA), lty = c(NA,1, 1)) # Horizonte de predicao usando o expoente de Lyapunov require("tseriesChaos") y_pred = as.ts(y_pred) output <-lyap_k(y_pred, m=3, d=2, s=200, t=2, ref=2, k=2, eps=4) par = lyap(output, 0.73, 2.47) lambda = par[2] plot(output) #=================================== O erro é esse: > output <-lyap_k(y_pred, m=3, d=2, s=200, t=2, ref=2, k=2, eps=4) Finding nearests Error in find_knearests(series, m = m, d = d, t = t, ref = length(ref), : cannot allocate memory block of size 134217728 Tb OBS: No meu caso, já havia testado e a dimensão m e d, são exatamente 3 e 2. Porém, não sei como controlar os demais parâmetros. Fico no aguardo, se alguém conseguir me explicar o problema. Agradeço e obrigado.

1 0

Re: [R-br] Selecionar bancos específicos de uma pasta
by sznelwar＠uol.com.br 22 Mai '20

22 Mai '20

1 0

Ajuda na função matplot
by Gilenio Borges Fernandes 21 Mai '20

21 Mai '20

Prezados(as) No CMR seguinte, esperava 255 linhas no matplot. Está plotando apenas uma linha e a seguinte mensagem: Error in plot.xy(xy.coords(x, y), type = type, ...) : tipo de gráfico inválido ',' Alguma ideia/sugestão? set.seed(1000) u = 0.25 dp = 0.15 periodo = 1:256 n_sim = 50 total_sim = matrix(0, nrow = 256, ncol = n_sim) p_inicial = 10 for(j in 1:n_sim) { aux = 2 p = p_inicial preco = c(p) dist = rnorm(255, 0, 1) for (i in dist) { P = p+p*(u/255 + dp/sqrt(255)*i) preco[aux] = P p = P aux = aux +1 } total_sim[,j] = preco } matplot(periodo,total_sim, main = "Simulacao Monte Carlo", xlab= "Tempo", ylab = "Valor", type = "l,", col = 1:n_sim) Gilênio Borges Fernandes Professor substituto da UFBA (03/2019 a ...) Professor substituto da UEFS/DEXA (01/2018 a 03/2019) Professor substituto Adjunto A (UFBA - 03/2015 a 03/2017) Professor Associado IV (Aposentado) Universidade Federal da Bahia Instituto de Matemática e Estatística Departamento de Estatística Av. Adhemar de Barros, s/n – Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6340/6341/6337 Fax: (071)3283-6336 ORCID: http://orcid.org/0000-0002-2364-0954 Lattes: http://lattes.cnpq.br/6764860618464860 https://gilenio.github.io/ Skype: gilenio.fernandes

2 1

Re: [R-br] Selecionar bancos específicos de uma pasta
by sznelwar＠uol.com.br 20 Mai '20

20 Mai '20

5 12

Selecionar bancos específicos de uma pasta
by Diogo Jerônimo 18 Mai '20

18 Mai '20

Lista, bom dia e boa semana!!! Criei o "reproduzível" para explicar a dúvida. Nesse exemplo, de uma pasta, identifiquei arquivos "*.csv", e depois apliquei o lapply e o do.call para importar e juntar esses bancos: dir.create('c:/users/diogo/desktop/bancosex/')banco1=data.frame(cbind(c("Diogo","Patrícia"),c(42,40),c(5000,7000))) write.csv2(banco1,file='c:/users/diogo/desktop/bancosex/banco1.csv',row.names=FALSE)banco2=data.frame(cbind(c("João","Alexandre"),c(41,40),c(8000,9000))) write.csv2(banco2,file='c:/users/diogo/desktop/bancosex/banco2.csv',row.names=FALSE)banco3=data.frame(cbind(c("Angélica","Nádia"),c(40,38),c(9500,7500))) write.csv2(banco3,file='c:/users/diogo/desktop/bancosex/banco3.csv',row.names=FALSE) setwd('c:/users/diogo/desktop/bancosex/') lista<-list.files() arquivos<-lapply(lista, function(x) read.csv2(x, header=TRUE, sep=";")) dados<-do.call("rbind", arquivos) Meu problema: esse código vai "pegar" TODOS os bancos da pasta antes de importar e juntar. Eu gostaria de pegar somente PARTE desses bancos (ex: apenas banco1 e banco2, apenas banco2 e banco3...). Alguém saberia como fazer isso? Obrigado!!! Diogo JerônimoBacharel em Ciências Estatísticas - ENCE/IBGEMestre em Metrologia - PUC-Rio/PósMQICONRE: 8514 - SÉRIE A

1 0