Março 2019 - 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

Redistribuição das contagens de uma estrutura de dados
by Paulo Eduardo de Mesquita 02 Abr '19

02 Abr '19

Prezado(a)s Colegas, Peço a vossa ajuda para encontrar solução para o seguinte problema: Eu tenho uma estrutura de dados do tipo tabela dinâmica cujos nomes de colunas podem ser "simples" ou "combinações separadas por vírgulas" (ver código anexo). As colunas com nomes combinados são construídas através da "concatenação" dos nomes de colunas simples. As colunas com nomes combinados assumem formas imprevisíveis e ocorrem aleatoriamente conforme a base de dados que lhe dá origem vai crescendo. Eu preciso redistribuir as contagens contidas nas colunas com nomes combinados para colunas de nomes simples, de modo que se, em uma linha qualquer da minha tabela de origem, houver uma coluna AB com valor 3, na minha tabela ajustada, o valor 3 será somado uma vez ao que já havia na mesma linha correspondente das colunas A e B. O exemplo que segue anexo, imagino eu, facilita a compreensão do problema. Muito obrigado pela atenção, -- Paulo Eduardo de Mesquita Disciplina de Infectologia - Faculdade de Medicina Universidade do Oeste Paulista Presidente Prudente - São Paulo - Brasil telefone: 5518 97718261

4 3

two-way anova - violação das pressuposições
by Marcelo Laia 28 Mar '19

28 Mar '19

Colegas, Estamos com um conjunto de dados que viola pressuposições da ANOVA, principalmente a normalidade. Eu já li, por diversas vezes, que essa violação pode não ser tão danosa assim, devido a robustez do teste. Por favor, tire um tempo para ver os resultados e depois peço uma ajuda. Genotipo e Isolado são fatores Area está em centímetros quadrados > leveneTest(Area ~ Genotipo*Isolado, data = cerato.desc) Levene's Test for Homogeneity of Variance (center = median) Df F value Pr(>F) group 41 3.4755 4.718e-08 *** 126 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 > cat("Normality p-values by Factor Genotipo: ") for (i in unique(factor(cerato.desc$Genotipo))){ cat(shapiro.test(cerato.desc[cerato.desc$Genotipo==i, ]$Area)$p.value," ") } 7.074459e-10 3.200422e-06 #Shapiro-Wilk normality tests by Isolado for (i in unique(factor(cerato.desc$Isolado))){ cat(shapiro.test(cerato.desc[cerato.desc$Isolado==i, ]$Area)$p.value," ") } 0.09534117 0.4006495 0.6065291 0.2093362 0.6138097 0.5604402 0.1302976 0.3135567 0.905537 0.7294285 0.0966383 0.1512716 0.8469947 0.1226855 0.2713435 0.9695489 0.2747097 0.5476302 0.0008750702 0.03693436 0.4197769 > bartlett.test(Area~Genotipo,data = cerato.desc ) Bartlett test of homogeneity of variances data: Area by Genotipo Bartlett's K-squared = 19.769, df = 1, p-value = 8.738e-06 > bartlett.test(Area~Isolado,data = cerato.desc ) Bartlett test of homogeneity of variances data: Area by Isolado Bartlett's K-squared = 171.26, df = 20, p-value < 2.2e-16 A pergunta é: mesmo com esses resultados eu poderia afirmar que o teste, neste caso, será robusto o suficiente para essas violações? Obrigado! -- Marcelo

6 14

Análise circular - fenologia
by lucas milanesi 14 Mar '19

14 Mar '19

Necessito encontrar o ângulo médios em graus de uma análise circular, onde os meses são representados em graus e as frequências da fenofase o outro vetor.... Alguém poderia me auxiliar numa função para executar isso? Att.Lucas

2 1

Adicionar barras em um mapa (ggplot2)
by Carlos Augusto Jardim Chiarelli 13 Mar '19

13 Mar '19

Como adicionar barras (bars) em um mapa ggplot2 (geom_polygon) no R 3.5?

2 1

Simular números
by Angélica Ricarte 12 Mar '19

12 Mar '19

Prezados, bom dia. Por gentileza, peço ajuda para simular 4 números de forma que a soma dê 100, como na tabela abaixo, por exemplo. Qual o comando e/ou pacote posso utilizar? Comun. Execul. Planej. Analis. TOTAL A1 21,26 20,45 28,74 29,55 100 A2 A3 A4 A5 A6 A7 A8 A9 A10 Grata, aguardo retorno. Att. Angélica Ricarte da Silva Batista Autarquia Educacional do Belo Jardim/Faculdade do Belo Jardim - AEB/FBJ Rodovia PE-166, Km 05 - BR, Belo Jardim - PE Telefone/Phone: +55 (81) 3726-1800 CV<http://lattes.cnpq.br/0601801978475226> <http://lattes.cnpq.br/0601801978475226>

3 3

Conflito de funções
by Elias Carvalho 10 Mar '19

10 Mar '19

Esto enfrentando o seguinte problema: Estou fazendo uma análise onde uso o pacote bnlearn e depois o pacote lavaan. Porém ambos tem a função "print.bn" e no bnlearn está função está embutida em outras funções. Em certos momentos eu crio uma rede Bayesiana e em seguida um modelo SEM com lavaan e depois volto a criar aa rede Bayesiana em em seguida o modelo SEM novamente, e assim sucessivamente, dai o erro ocorre porque o bnlearn acaba usando a função do lavaan. [image: Screenshot from 2017-10-05 10-57-11.png] Alguem poderia me ajudar a solucionar esse problema? -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho* *"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he who has been able to understand the cause of things"*

2 1

Teste de Qui-quadrado aplicado para a distribuição binomial
by ASANTOS 07 Mar '19

07 Mar '19

Prezados Membros, Estou tentando realizar um teste de aderência (Qui-quadrado) aplicado a dados com distribuição binomial negativa sem sucesso. Em meu CRM: #Pacote MASS library(MASS) #Crio um banco de dados simulado Features<-sort(rep(1:3,100)) Resp1<-rnbinom(200,0.75,0.10) ##60 dados com distribuição binomial negativa Resp2<-rpois(100,10) ##30 dados com distribuição de Poisson Resp<-c(Resp1,Resp2) d<-as.data.frame(cbind(Features,Resp)) #Criando um loop para realizar o teste de Chi quadrado para os três conjuntos de dados - Feature 1, 2 e 3 uniq <- unique(unlist(d$Features)) res<-NULL for (i in 1:length(uniq)){ data_1 <- subset(d, Features == uniq[i]) k <- fitdistr(data_1$Resp,"negative binomial") #Extração dos parâmetros média e mu da binomial negativa par <- k$estimate size <- par[1]#Parâmetro k mu <- par[2]#Média N <- length(data_1$Resp) est <-N*dnbinom(data_1$Resp,size=size,mu=mu) ## Criando valores estimados fecdf <- ecdf(data_1$Resp) ###ecdf- Cria a distribuição cumulativa empírica knotsX <- knots(fecdf) emp <- fecdf(c(knotsX,Inf)) # Distribuição empírica testChi<-chisq.test(table(emp),table(est),correct=TRUE) #Teste de Qui-quadrado sobre a contagem dos dados empíricos versus estimados result<-ifelse(testChi$p.value>0.05,"Binomial Negative","Other distribution") ## Classifica os testes em binomial negativa ou outra distribuição res<-rbind(res,c(uniq[i],result)) colnames(res)<-c("Features","Distribution") } res Meu output é: Features Distribution [1,] "1" "Other distribution" [2,] "2" "Other distribution" [3,] "3" "Binomial Negative" Ou seja, dá o resultado errado, na verdade as Features 1 e 2 é que são binomial negativa e a Feature 3 que são dados com distribuição de Poisson que tem que dar "Other distribution". Alguém teria alguma sugestão para dar? 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 4

Re: [R-br] Modelos Mistos - Duvida sobre transformação/distribuição logistica
by sznelwar＠uol.com.br 05 Mar '19

05 Mar '19

2 1