- R-br - DINF UFPR

tidyverse e vroom: Leitura de informações alvo em um arquivo *txt e também computo de informação ausente
by ASANTOS 29 Set '20

29 Set '20

Prezados membros, Surgiu um novo problema para a extração de informações em um arquivo *txt quando eu tenho que computar também quando eu não tenho determinada informação (lca) e preciso considerar isso no meu data.frame final. Em meu exemplo: #Pacotes library(vroom) library(tidyverse) #Download dos dados myfile <- vroom("https://raw.githubusercontent.com/Leprechault/yolov3ants/master/result_imgv…") As informações alvo são: 1) No caso quando ocorreu o evento lca o padrão das linhas é: imgval/sample_21915.jpg: Predicted in 7337.147000 milli-seconds. lca: 65.648048% left_x: 0 top_y: 18 width: 13 height: 27 lca: 45.191948% left_x: 71 top_y: 31 width: 7 height: 12 End Então meu data.frame precisa receber a informação: image lca sample_21915.jpg 1 sample_21915.jpg 1 2) Quando não ocorrer o evento lca, o padrão das linhas será diferente: Start processing imgval/sample_3908.jpg imgval/sample_3908.jpg: Predicted in 7282.433000 milli-seconds. End E nesse caso o meu data.frame precisa receber a seguinte informação: image lca sample_3908.jpg 0 Gostaria de saber ideias de como filtrar essa informação e especialmente do que fazer para gerar os zeros, se o meu data.frame final que preciso, após a análise da informação de 100 imagens deve ser: image lca sample_14715 1 sample_3101 0 sample_5769 0 sample_21884 1 sample_4163 0 sample_20368 0 sample_20371 0 sample_3370 0 sample_4434 0 sample_20888 1 sample_19295 1 sample_22196 0 ... sample_20314.jpg 0 Antecipadamente agradeço, Alexandre -- Alexandre dos Santos Geotechnologies and Spatial Statistics applied to Forest Entomology Instituto Federal de Mato Grosso (IFMT) - Campus Caceres Caixa Postal 244 (PO Box) Avenida dos Ramires, s/n - Vila Real Caceres - MT - CEP 78201-380 (ZIP code) Phone: (+55) 65 99686-6970 / (+55) 65 3221-2674 Lattes CV: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 ResearchGate: www.researchgate.net/profile/Alexandre_Santos10 Publons: https://publons.com/researcher/3085587/alexandre-dos-santos/ --

2 4

Dúvida para conversão de um output em *txt para data.frame usando tidyverse [NÃO FORMATADO]
by ASANTOS 29 Set '20

29 Set '20

Estou com um problema para criar um dataframe a partir de um conjunto de duas linhas de interesse, sendo: library(tidyverse) myfile<-read_lines("https://raw.githubusercontent.com/Leprechault/yolov3ants/master/yolov3-lca_…") head(myfile) Em algumas parte do meu arquivo eu tenho trechos que quero extrair que são: #Last accuracy mAP(a)0.5 = 73.10 % #1585: 0.299288, 0.423193 avg loss, 0.001000 rate, 3.927486 seconds, 101440 images Nesses pontos eu quero remover as informações: mAP(a)0.5 = 73.10 % (na primeira linha), 1585: e 0.423193 avg loss (da linha seguinte) e criar um dataframe final que seria: # iteration avg_loss mAP # 1585 0.423193 73.10 mas quando tento remover a informação de interesse e por ser beginner no tidyverse, não tenho sucesso, sendo a minha tentativa no sentido de: names_col2 <- c("iteration", "avg_loss", "mAP") mydf <- myfile %>% str_subset("Last accuracy mAP(a)0.5 =") %>% enframe(name = NULL) %>% separate(col = value, into = names_col2, sep = "[=,]") %>% mutate_all(parse_number) head(as.data.frame(mydf)) Alguém poderia me dar uma luz? Obrigado e abraços, Alexandre -- Alexandre dos Santos Geotechnologies and Spatial Statistics applied to Forest Entomology Instituto Federal de Mato Grosso (IFMT) - Campus Caceres Caixa Postal 244 (PO Box) Avenida dos Ramires, s/n - Vila Real Caceres - MT - CEP 78201-380 (ZIP code) Phone: (+55) 65 99686-6970 / (+55) 65 3221-2674 Lattes CV: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 ResearchGate: www.researchgate.net/profile/Alexandre_Santos10 Publons: https://publons.com/researcher/3085587/alexandre-dos-santos/ --

2 2

Dúvida para conversão de um output em *txt para data.frame usando tidyverse
by ASANTOS 29 Set '20

29 Set '20

Estoucomumproblemaparacriarumdataframeapartirdeumconjuntodeduaslinhasdeinteresse, sendo: library(tidyverse) myfile<-read_lines("https://raw.githubusercontent.com/Leprechault/yolov3ants/master/yolov3-lca_…") head(myfile) Emalgumaspartedomeuarquivoeutenhotrechosquequeroextrairquesão: #Last accuracy mAP(a)0.5 = 73.10 % #1585: 0.299288, 0.423193 avg loss, 0.001000 rate, 3.927486 seconds, 101440 images Nessespontoseuqueroremoverasinformações:mAP@0.5 = 73.10 % (naprimeiralinha), 1585:e0.423193avgloss (dalinhaseguinte) ecriarumdata.framefinalqueseria: # iteration avg_loss mAP # 1585 0.423193 73.10 # ... masquandotentoremoverainformaçãodeinteresseeporserbeginnernotidyverse, nãotenhosucesso, sendoaminhatentativana seguinte direção: names_col2 <- c("iteration", "avg_loss", "mAP") mydf <- myfile %>% str_subset("Last accuracy mAP(a)0.5 =") %>% enframe(name = NULL) %>% separate(col = value, into = names_col2, sep = "[=,]") %>% mutate_all(parse_number) head(as.data.frame(mydf)) Alguémpoderiamedarumaluz? Obrigado, Alexandre -- Alexandre dos Santos Geotechnologies and Spatial Statistics applied to Forest Entomology Instituto Federal de Mato Grosso (IFMT) - Campus Caceres Caixa Postal 244 (PO Box) Avenida dos Ramires, s/n - Vila Real Caceres - MT - CEP 78201-380 (ZIP code) Phone: (+55) 65 99686-6970 / (+55) 65 3221-2674 Lattes CV: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 ResearchGate: www.researchgate.net/profile/Alexandre_Santos10 Publons: https://publons.com/researcher/3085587/alexandre-dos-santos/ --

2 1

Assimetria
by ana paula coelho madeira 24 Set '20

24 Set '20

Prezados, Bom dia! Gostaria de sugestões de como plotar uma densidade com assimetria negativa. Desde já agradeço! At.te Ana Paula

2 1

Diferença entre o teste de normalidade normalidade sobre a variável resposta e sobre a aov() dos dados em função da variável explicativa
by ASANTOS 22 Ago '20

22 Ago '20

Prezados Membro do R-br, Gostaria de saber a diferença entre realizar: 1) Análise de normalidade sobre a variável resposta x <- rep(seq(from=10, to=50, by=0.5),6) y1 <- rnorm(length(x), mean=10, sd=1.5) y1.res<-y1-mean(y1) shapiro.test(y1.res) # Shapiro-Wilk normality test # #data: y1.res #W = 0.99688, p-value = 0.4787 2) Sobre a análise de variância do modelo em função da variável explicativa y1.av<-aov(y1 ~ x) shapiro.test(y1.av$res) # Shapiro-Wilk normality test # #data: y1.av$res #W = 0.99736, p-value = 0.6338 Obrigado, Alexandre -- Alexandre dos Santos Geotechnologies and Spatial Statistics applied to Forest Entomology Instituto Federal de Mato Grosso (IFMT) - Campus Caceres Caixa Postal 244 (PO Box) Avenida dos Ramires, s/n - Vila Real Caceres - MT - CEP 78201-380 (ZIP code) Phone: (+55) 65 99686-6970 / (+55) 65 3221-2674 Lattes CV: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 ResearchGate: www.researchgate.net/profile/Alexandre_Santos10 Publons: https://publons.com/researcher/3085587/alexandre-dos-santos/ --

4 3

Mudar eixo numérico para labels
by Diogo Jerônimo 10 Ago '20

10 Ago '20

Gente, mais uma dúvida, para isso o exemplo abaixo (em vermelho): Item=c("Lote1","Lote2","Lote3","Lote4","Lote5") R1=c(8.30,5.69,3.55,1.29,0.01)R2=c(7.70,5.62,3.36,1.28,0.01)Título="Análise de Linearidade"CR=seq(1:mean(length(R1),length(R2))) dados=data.frame(Item,R1,R2,CR)i=1 for (i in (1:nrow(dados))){dados$MRES[i]=mean(c(dados$R1[i],dados$R2[i]))i=i+1}dados$CR2=(dados$CR)^2 dados$CR3=(dados$CR)^3 dados$reglin=as.numeric(matrix(predict(lm(MRES~CR,data=dados)))) valmax=trunc(max(max(dados$MRES),max(dados$reglin)),2)+1 minesc=0maxesc=valmaxplot(c(min(dados$CR),max(dados$CR)),c(minesc,maxesc),type="n",main=Título,xlab="",ylab="") points(dados$CR,dados$MRES,col="blue",pch=16,cex=1.8)lines(dados$CR,dados$MRES,col="blue",lwd=2)points(dados$CR,dados$reglin,col="green",pch=16,cex=1.8)lines(dados$CR,dados$reglin,col="green",lwd=2) No exemplo, com esses dados, fiz um gráfico de pontos e linhas. sendo que como era de esperar, ele vai apresentar valores numéricos nas duas colunas. O que queria saber é se é possível substituir no eixo "x" os números 1,...5, pelos labels de texto que aparecem na variável "Item" ("Lote 1",..., Lote5). Tentei mudar transformando o "class" de CR de "numeric" para "factor", mas não deu certo... Obrigado novamente pela atenção!!! Diogo JerônimoBacharel em Ciências Estatísticas - ENCE/IBGEMestre em Metrologia - PUC-Rio/PósMQICONRE: 8514 - SÉRIE Ahttp://lattes.cnpq.br/8996149312896520

2 2

Artigo sobre p.value
by Luiz Roberto Martins Pinto 10 Ago '20

10 Ago '20

Bom dia. Há algum tempo alguém encaminhou para a lista um artigo sobre p.value. O autor se referiu ao artigo como "o artigo que ele sempre quiz escrever", ou algo parecido. Eu tive problema no computador e perdi. Peço ajuda para acessar este artigo. Outros artigos temáticos sobre p.value também serão bem-vindos. Grato. Luiz Roberto. Luiz Roberto Martins Pinto Prof. Pleno/DCET/UESC Laboratório de Estatística Computacional Universidade Estadual de Santa Cruz Ilhéus-Bahia-Brasil luizroberto.uesc(a)gmail.com skype: lrmpinto http://lattes.cnpq.br/2732314327604831 "*The s**cience exists because there are patterns. * * The patterns exist because God created them*. * The statistic exists to research the patterns that God created.*"

4 10

Diferença de todos os valores
by Diogo Jerônimo 07 Ago '20

07 Ago '20

Boa tarde a todos. Seja por exemplo o seguinte vetor: x=(19,23,27,29) É possível se obter no R um vetor com TODAS as diferenças absolutas entre TODOS os pares de valores combinados de um vetor? Ou seja: vetor=(|19-23|,|19-27|,|19-29|,|23-27|,|23-29|,|27-29|) Obrigado por qualquer ajuda!!! Diogo JerônimoBacharel em Ciências Estatísticas - ENCE/IBGEMestre em Metrologia - PUC-Rio/PósMQICONRE: 8514 - SÉRIE Ahttp://lattes.cnpq.br/8996149312896520

4 8

Erro ao rodar análise bifatorial - ExpDes.pt
by Tiago Camponogara Tomazetti 07 Ago '20

07 Ago '20

Prezados, gostaria de saber se alguém pode me tirar uma dúvida quanto a um "possível" bug que tenho tido quando uso a função fat2.dic() do pacote ExpDes.pt. Sempre que vou rodar a análise bifatorial, onde há interação, a função da o seguinte erro: ######################################################################## ... Interacao significativa: desdobrando a interacao ------------------------------------------------------------------------ Desdobrando F1 dentro de cada nivel de F2 ------------------------------------------------------------------------ ------------------------------------------------------------------------ Quadro da analise de variancia ------------------------------------------------------------------------ GL SQ QM Fc Pr.Fc F2 4 194102.2 48525.55 4.768 0.0016 F1:F2 5 NA NA <NA> <NA> F1:F2 D0 5 NA NA <NA> <NA> F1:F2 D3 5 NA NA <NA> <NA> F1:F2 D5 5 NA NA <NA> <NA> F1:F2 D7 5 NA NA <NA> <NA> Residuo 87 885431.5 10177.37 Total 107 3341005.8 31224.35 ------------------------------------------------------------------------ Error in if (des1.tab[(i + 2), 5] <= sigF) { : valor ausente onde TRUE/FALSE necessário ######################################################################## Aparentemente é um erro dentro do pacote, contudo, se alguém conseguiu contornar este bug de alguma forma, ou conhecer outro modo de rodar uma análise bifatorial no R, seria grato por receber suas opiniões. A título de informação, estou chamando a função deste modo: "fat2.dic(data$TratA, data$TratB, data$Resp, quali=c(TRUE, TRUE), mcomp='tukey')" Desde já, agradeço pela atenção, Tiago Dr. Tiago Camponogara Tomazetti Plant Breeding professor - Federal University of Santa Catarina +55 (48) 9-9681-3848 Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4408476D6 ResearchGate https://www.researchgate.net/profile/Tiago_Tomazetti

2 2

Ponto de corte Cook's Distance
by Fernando Souza 07 Ago '20

07 Ago '20

Prezados colegas, Gostaria de saber como vocês definem o ponto de corte para a distância de Cook? Na literatura e listas de discussão tenho vistos várias sugestões, como valores de cook > 1 , valores maiores que 4/n , valor de cook maior que 3x a média, entre outros. Como eu devo proceder para definir o melhor ponto de corte. Eu estou analisando uma variável onde do total de 35 observações , 6 observações etão acima de 0.10 e dessas duas, as mais distantes estão entre 0.25 e 0.30. Eu deveria considerar as 6 observações como pontos de influencia? Se eu utilizar a regra de 4/35 = 0.11 então às 6 observações seriam removidas; Se eu considerar D>1, então nenhuma observação será removida Gostaria de saber sobre a experiencia de vocês e como procedem ATt

2 1