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

Estatísticas do R
by Andre Oliveira 26 Abr '16

26 Abr '16

Pessoal,preciso de estatísticas do R tais como nº de bibliotecas oficiais, pacotes que o R já vem na instalação. 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

8 8

Microdados IBGE
by Thiago Silva Pereira 04 Abr '16

04 Abr '16

#Packages > library("dicionariosIBGE", lib.loc="~/R/win-library/3.2") #Carregar dados armazenados no pacote > data(dicPNAD2005) #ver dados > View(dicpes2005) #localizar micro dados > caminho.microdados <- file.path("C:","DADOS","PES2005.TXT") > DADOS <- le.pesquisa(dicionario=dicpes2005, pathname.in =caminho.microdados, + pathname.out="dados.Rdata", + codigos=c("V4703","V4704","V4705","V4707","V4729"), + nomes=c("ANOSEST","CONDATV","CONDOCP","HTRAB","PESO"), + quant=408148) Error in le.pesquisa(dicionario = dicpes2005, pathname.in = caminho.microdados, : unused arguments (pathname.out = "dados.Rdata", nomes = c("ANOSEST", "CONDATV", "CONDOCP", "HTRAB", "PESO"), quant = 408148) Vou mandar o PDF que estou estudando, Obrigado deis de já Obs: Os codigos e explicações estão lá.

7 12

"Unmelt" data frame
by Jefferson Ferreira-Ferreira 31 Mar '16

31 Mar '16

Pessoal, estou com uma questão aberta no stackoverflow. Se alguém puder ajudar, eu ficaria muito grato! http://stackoverflow.com/questions/36271393/r-unmelt-a-melted-dataframe-bac… -- *Jefferson Ferreira-Ferreira* Geógrafo – GEOPROCESSAMENTO IDSM | Coordenadoria de TI Jefferson.ferreira(a)mamiraua.org.br *Instituto de Desenvolvimento Sustentável Mamirauá* Ministério da Ciência, Tecnologia e Inovação Telefone: +55 97 3343-9710 *Google Maps* - Mapas deste e-mail: Exibir mapa ampliado <https://maps.google.com.br/maps?q=-3.355557,-64.731151&ll=-3.355471,-64.731…> *Contatos particulares:* *(55) 9615-0100*

3 5

Nova library ssmn
by Clécio da Silva Ferreira 30 Mar '16

30 Mar '16

Prezados, gostaria de divulgar a library confeccionada pelo Luis Benites Sanchez e por mim, baseada nas distribuições SSMN. ssmn: Skew Scale Mixtures of Normal Distributions Performs the EM algorithm for regression models using Skew Scale Mixtures of Normal Distributions. Version: 1.0 Author: Luis Benites Sanchez and Clecio da Silva Ferreira Atenciosamente, Prof. Clécio da Silva Ferreira Depto de Estatística - UFJF Tel: (32) 2102-3306 Página: www.ufjf.br/clecio_ferreira

2 1

Teste de diferça entre proporções via loop
by Wagner Tassinari 29 Mar '16

29 Mar '16

Olá pessoal, tudo bem ? Estou querendo rodar um teste de diferença entre proporções em várias linhas do meu banco de dados, daí gostaria de saber como poderia rodar de uma vez só, através de um loop, esse comando e gerar uma coluna só com os p-valores. Veja como estõa tentando fazer abaixo: O comando seria esse: prop.test(x = c(x1, x2), n = c(n1, n2), correct = FALSE) O banco seria dessa forma: AREA X1 X2 N p_valor 1 1982 9416 11398 ? 2 6546 4084 10630 ? 3 3418 667 4085 ? 4 1503 4408 5911 ? 5 966 339 305 ? AREA1 <- prop.test(x = c(1982, 9416), n = c(11398, 11398), correct = FALSE) AREA1$p.value AREA2 <- prop.test(x = c(6546, 4084), n = c(10630, 10630), correct = FALSE) AREA2$p.value AREA3 <- prop.test(x = c(3418, 667), n = c(4085, 4085), correct = FALSE) AREA3$p.value ... Obrigado -- Wagner S. Tassinari Departamento de Matemática Universidade Federal Rural do Rio de Janeiro. BR-465, Km 7 - Seropedica, RJ - Brasil CEP: 23890-000 Skype: wagner.tassinari wtassinari(a)gmail.com tassinari(a)ufrrj.br ------------------------------------------------------- "Statistical thinking will one day be as necessary for efficient citizenship as the ability to read and write." (H.G.Wellis)

5 8

Criando um fatorial completo com dados desbalanceados
by ASANTOS 29 Mar '16

29 Mar '16

Prezados Listeiros, Gostaria de criar um fatorial completo entre 4 tratamentos com dados desbalanceados usando a função gen.factorial() com o pacote AlgDesign. Alguém poderia me ajudar, segue CRM abaixo: #Dados artificiais tratamentos_b <- gl(4, 15, labels = paste("t", 1:3, sep=""))## Dados balanceados tratamentos_nb <- gl(1, 30, labels = paste("t", 4, sep="")) ### Dados desbalanceados tratamentos<-as.character(tratamentos_b,tratamentos_nb) sp <- cbind(c(rnorm(10, 5, 0.25), rnorm(50, 2.5, 0.25)), rnorm(60, 2.5, 0.25), c(rnorm(10, 12, 0.25), rnorm(50, 2.5, 0.25)), rnorm(60, 2.5, 0.25)) colnames(sp) <- c("sp1", "sp2", "sp3", "sp4") dados<-data.frame(sp,tratamentos) ## Criando o fatorial completo require(AlgDesign) contr2df2 <-AlgDesign::gen.factorial(dados$tratamentos, TRUE, varNames=c("T1_2", "T1_3", "T1_4", "T2_3", "T2_4", "T3_4")) contr2df2 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 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) 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 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================

2 2

Função para criar agrupamento de variáveis em função do valor de p
by ASANTOS 29 Mar '16

29 Mar '16

Prezados Listeiros, Gostaria de uma ajuda do grupo para criação de uma função para criação de uma variável dados$groups, em função da significância de res$aov.tab$Pr=> 0.05, de maneira que níveis de fatores iguais serão fundidos, como no exemplo abaixo onde trat1 e trat2 são iguais e vão ficar trat1_trat2, e os não significativos ficam do mesmo jeito, sem alteração. O output desejado esta no final do meu CRM, não estou conseguindo chegar a resultados satisfatórios com a função ifelse (res$aov.tab$Pr=> 0.05) { ... Alguém teria alguma sugestão? ### <code r> require(vegan) #Fator trat <- gl(3, 15, labels = paste("t", 1:3, sep="")) # Variáveis resposta set.seed(124) sp <- cbind(c(rnorm(10, 5, 0.25), rnorm(35, 5, 0.25)), rnorm(45, 2.5, 0.25), c(rnorm(10, 4.5, 0.25), rnorm(35, 5, 0.25)), rnorm(45, 5, 0.25)) colnames(sp) <- c("sp1", "sp2", "sp3", "sp4") #Comparações contrasts(trat) <- cbind(c(0,1,0),c(0,0,1)) trat_c <- model.matrix(~ trat)[, -1] trat1_trat2 <- trat_c[, 1] trat1_trat3 <- trat_c[, 2] res <- adonis(sp ~ trat1_trat2 + trat1_trat3, method = "euclidean", perm = 999) res # Tentativa de função ifelse (res$aov.tab$Pr=> 0.05) { dados$groups = ## dificultade aqui para agrupar os fatores não significativos , dados$groups = dados$trat } # Output desejado trat sp1 sp2 sp3 sp4 groups [1,] 1 4.653732 2.490977 4.236323 5.382113 trat1_trat2 [2,] 1 5.009581 2.254713 4.604529 4.842553 trat1_trat2 [3,] 1 4.809242 2.391435 4.675318 4.732977 trat1_trat2 [4,] 1 5.053077 2.483129 4.561690 5.311215 trat1_trat2 [5,] 1 5.356384 2.745474 4.616074 5.114969 trat1_trat2 ... # 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 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) 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 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================

1 1

XLConnect no R-studio
by Luciana Barbosa 28 Mar '16

28 Mar '16

Mudei de computador, e instalei Rx64 3.24 e Ri386 3.2.4, porém quando fui usar a biblioteca XLConnect no Ri386 3.2.4 funciona normalmente, porém para o Rx64 3.24 e o Rstudio aparece sempre o mesmo erro, mesmo fazendo atualizações ou seguindo o que pede para ser feito, parece que é incompatibilidade no JAVA, gostaria de voltar a usar o Rstudio. Por favor alguém pode ajudar-me? Essa é a mensagem que aparece: > library(XLConnect) Carregando pacotes exigidos: XLConnectJars Error : .onLoad falhou em loadNamespace() para 'rJava', detalhes: chamada: fun(libname, pkgname) erro: No CurrentVersion entry in Software/JavaSoft registry! Try re-installing Java and make sure R and Java have matching architectures. Error: package ‘XLConnectJars’ could not be loaded

2 2