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

RCurl: passando parâmetros e executando função javascript
by Éder Comunello 04 Dez '14

04 Dez '14

Colegas, boa tarde! Estou precisando de ajuda com relação ao pacote *RCurl*. Tomando por exemplo o endereço < http://hidroweb.ana.gov.br/Estacao.asp?Codigo=2254001>, ao acessá-lo via *browser *pode-se notar dois botões que executam trechos de código em *javascript *e geram links para o *download *de arquivos de dados: * Botão [Arquivo Access] remete ao comando "javascript:criarArq(2254001,2)" * Botão [Arquivo Texto] remete ao comando "javascript:criarArq(2254001,1)" Estou tentando usar o *RCurl *para passar diretamente os parâmetros para criação do *link *para o "Arquivo Texto". Isto feito, posteriormente pretendo usar o *link *gerado para *download* dos dados. Contudo, não estou tendo sucesso e ficarei muito grato por qualquer ajuda. Abaixo segue o código que tentei aplicar: ### <BEGIN> ### library(RCurl) myURL <- 'http://hidroweb.ana.gov.br/Estacao.asp?Codigo=2254001' ### URL resultante da operação via browser: ### myURL2 <- ' http://hidroweb.ana.gov.br/Estacao.asp?Codigo=2254001&CriaArq=true&TipoArq=1 ' myCurl <- getCurlHandle(cookiefile = "myCookies.txt", useragent="Mozilla/5.0", verbose = T, autoreferer = T, followlocation = T) myParams <- list(Codigo=2254001, CriaArq='true', TipoArq=1) ### Codigo: ID da estação ### CriaArq: true (executa função?), ### TipoArq: 1 (tipo texto) e 2 (Access) postForm(myURL, .params=myParams, curl=myCurl) getURL(myURL, curl=myCurl) ### Operação falha, retornando página original! ### Em caso de sucesso embutiria no código um link no formato: ### <a href="ARQ/A20131008-164720-14/CHUVAS.ZIP"> ### <END> ### -- Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]

4 3

Como faço para descobrir o poder de um teste de hipótese
by Jodavid Ferreira 15 Abr '14

15 Abr '14

Olá, Estou gerando duas amostras aleatórias de poisson: amostra1<-rpois(500,lambda=2.5) amostra2<-rpois(500,lambda=3.5) e aplicando o teste de wilcoxon e o teste t para dados pareados. wilcox.test(amostra1,amostra2, paired=T) t.test(amostra1,amostra2, paired=T) e gostaria de obter o poder de cada um dos testes acima. -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://rccpbjovem.com Membro da equipe Estadual do Ministério de Comunicação da RCCPB - http://rccpb.com.br Responsável pelo Ministério de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com Twitter - @JodavidFerreira <http://twitter.com/JodavidFerreira>

3 2

Criar função utilizando "for"
by Rodrigo Muniz 25 Mar '14

25 Mar '14

Caros membros da lista, Me desculpem a ignorância, porém sou novo como usuário e estou tentando criar uma função para facilitar uma operação no meu banco de dados. considerem o exemplo: dados = as.data.frame(matrix(rnorm(200), ncol=10, nrow=50)) cor.multipla <- function(x, y){ n <- ncol(y) cor.mult <- numeric(0) + for(i in 1:n){cor.mult<- cor.test(x, y[,i]) + return(cor.mult$p.value)} } cor.multipla(dados$V1, dados[2:10]) Nesta função, eu gostaria fazer a correlação da coluna 1 com as colunas de 2 a 10, e extrair um vetor de p.value para cada correlação, no entanto, a função está me retornando apenas um valor de p.value. Alguém poderia me dizer onde eu estou errando ou mesmo se estou indo pelo caminho errado? Desde já agradeço... -- *Rodrigo A. Muniz* Eng. Agrônomo. Ms Produção Vegetal (UENF) Doutorando em Engenharia de Sistemas Agrícolas (ESALQ/USP) Cel (19) 8300-4333 (Pessoal) Tel (19) 3375-1196 (Residêncial) Tel (19) 3447-8537 (Trabalho)

4 5

rodando código R na web
by Paulo Nogueira Starzynski 25 Mar '14

25 Mar '14

Caros colegas de lista, algum tempo atrás deparei-me com uma ferramenta online onde era possível executar códigos em R. Na época não guardei nos favoritos e agora não consigo voltar a encontrar. Lembra um pouco essa outra, para Python: http://www.codeskulptor.org/ Alguém sabe qual é e pode me passar a localização? Abraços, Paulo Nogueira Starzynski

4 4

Grárico de controle
by Edson Lira 25 Mar '14

25 Mar '14

Caros amigos da lista. Tenho 3507 medidas de tempo de doação distribuídas por dias conforme a tabela abaixo. 01 02 03 04 05 07 08 09 10 11 14 15 16 17 18 19 21 22 23 25 26 28 29 30 31 92 165 154 113 237 146 139 156 92 149 147 142 120 123 142 125 111 94 84 85 73 91 140 130 105 Gostaria de plotar os gráficos xbarra e o da amplitude para esses tempos. Usei a rotina abaixo como o Eder me indicou: amostra = sample(rep(1:631,5),3155) #rever esses parâmetros, pois só postei uma amostra require(qcc) tem <- qcc.groups(tempo,amostra) Plotei o gráfico e funcionou, me mostrando os tempos e os grupos de amostra. Funcioniou, só que me deparei com outro problema, essa rotina me dá o tempo e o número da amostra, como faço para ver o dia em que o tempo ou amostra de tamanho 5 foi registrada. Ou seja, necessito mostrar o gráfico com o tempo e os dias do registro do tempo em amostras de tamanho 5. No link abaixo tem um link para o arquivo que estou usando com 300 observações. http://www.datafilehost.com/d/e38f36e0 [ ]'s. Edson Lira Estatístico Manaus-Amazonas

2 12

Tirar maiuscula apos o ponto
by Bernardo Rangel Tura 25 Mar '14

25 Mar '14

Boa noite a todos Primeiro de tudo quero pedir desculpas por uma duvida tao basica considerem o seguinte exemplo: texto <- 'Deus abencoou Noe e seus filhos, e lhes disse: "Sede fecundos, multiplicai-vos, enchei a terra. Sede o medo e o pavor de todos os animais da terra e de todas as aves do ceu, como de tudo o que se move na terra e de todos os peixes do mar: eles sao entregues nas vossas maos. Tudo o que se move e possui a vida vos servira de alimento, tudo isso eu vos dou, como vos dei a verdura das plantas. Mas nao comereis a carne com sua alma, isto e, o sangue. Pedirei contas porem, do sangue de cada um de vos. Pedirei contas a todos os animais e ao homem, aos homens entre si, eu pedirei contas da alma do homem. Quem derrama o sangue do homem pelo homem terra seu sangue derramado. Pois a imagem de Deus o homem foi feito. Quanto a vos, sede fecundos, multiplicai-vos, povoai a terra e dominai-a." Agora desejo transformar as maiusculas apos o ponto em minusuclas, por exemplo: ". Sede o medo" em ". sede o medo" sei que posso fazer com gsub e pensei que gsub('. [A-Z]','. [a-z]',texto) seria suficiente, mas nao foi. Onde estou errando? Antecipadamente agradeco aos R-mestres []s Tura

6 13

Re: [R-br] Criar função utilizando "for"
by Robert Iquiapaza 25 Mar '14

25 Mar '14

Caso se sinta mais confortável com código explícito: dados = data.frame(matrix(rnorm(500), ncol=10, nrow=50)) cor.multipla <- function(x, y){ x=as.numeric(x) n <- ncol(y) y=as.matrix(y,ncol=n) cor.mult <- vector() for(i in 1:n){cor.mult[i]<- cor.test(x, y[,i])$p.value} return(cor.mult) } dim(dados) cor.multipla(dados$X1, dados[2:10]) Sds Robert -----Mensagem Original----- From: Paulo Justiniano Sent: Monday, March 24, 2014 10:08 PM To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Criar função utilizando "for" O seu loop for tá sobre escrevendo resultados em cor.mult<- cor.test(x, y[,i]) e portanto somente o ultimo está sendo guardado (note tb que sua matrz de dados esá reciclando valores) Uma sugestão: x <- 1:10 y <- matrix(rnorm(200), nrow=10) apply(y, 2, function(y, x=x) cor.test(x,y)$p.value On Mon, 24 Mar 2014, Rodrigo Muniz wrote: > Caros membros da lista, Me desculpem a ignorância, porém sou novo como > usuário e estou tentando criar uma função para facilitar uma operação no > meu banco de dados. > considerem o exemplo: > > dados = as.data.frame(matrix(rnorm(200), ncol=10, nrow=50)) > cor.multipla <- function(x, y){ > n <- ncol(y) > cor.mult <- numeric(0) > + for(i in 1:n){cor.mult<- cor.test(x, y[,i]) > + return(cor.mult$p.value)} > } > > cor.multipla(dados$V1, dados[2:10]) > > Nesta função, eu gostaria fazer a correlação da coluna 1 com as colunas de > 2 a 10, e extrair um vetor de p.value para cada correlação, no entanto, > a função está me retornando apenas um valor de p.value. > Alguém poderia me dizer onde eu estou errando ou mesmo se estou indo pelo > caminho errado? > Desde já agradeço... > -- > Rodrigo A. MunizEng. Agrônomo. Ms Produção Vegetal (UENF) > Doutorando em Engenharia de Sistemas Agrícolas (ESALQ/USP) > Cel (19) 8300-4333 (Pessoal) > Tel (19) 3375-1196 (Residêncial) > Tel (19) 3447-8537 (Trabalho) > > _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

1 0

problemas com comparações múltiplas
by Maurício Lordêlo 25 Mar '14

25 Mar '14

Saudações a todos! Estou realizando inferência simultânea para um conjunto de dados correlacionados onde desejo comparar o efeito de quatro tratamentos. Para isto estou utilizando a função "simplesimint" do pacote "BSagri". Consigo fazer utilizando Dunnet porém quando mudo para Tukey ocorre um problema. Abaixo tenho um CMR onde tenho o vetor de coeficientes estimados (coefs), a matriz de variância e covariâncias (VarCov), o número de obervações em cada grupo (n) e a matriz de contrastes (cmat). O problema principal surge quando faço utilizando todas as linhas (são seis) da matriz de contrastes (cmat2). Já para duas, três e até quatro linhas da matriz de contrastes não há problemas. Grato caso alguém possa ajudar. Abraço library(BSagri) coefs=c("TratamentoA"=-1.0618,"TratamentoB"=1.1470,"TratamentoC"=0.3068,"TratamentoC"=0.2501) VarCov=matrix(data = NA, nrow = 4, ncol = 4, byrow = FALSE,dimnames=NULL) VarCov[1,]=c(0.3133311,-0.1618691,-0.1736806,-0.1712345) VarCov[2,]=c(-0.1618697,0.4510964,0.1595845,0.1622442) VarCov[3,]=c(-0.1736809,0.1595843,0.2729564,0.1617069) VarCov[4,]=c(-0.1712371,0.1622446,0.1617080,0.2757233) VarCov n=c("TratamentoA"=40,"TratamentoB"=20,"TratamentoC"=60,"TratamentoD"=60) cmat1<-contrMat(n=n, type="Dunnett") cmat1 test1<-simplesimint(coef=coefs, vcov=VarCov, df=NULL, cmat=cmat1, alternative="two.sided") plotCI(test1) cmat2<-contrMat(n=n, type="Tukey") cmat2 test2<-simplesimint(coef=coefs, vcov=VarCov, df=NULL, cmat=cmat2, alternative="two.sided") ##ERRO!! #Erro em uniroot(pfct, interval = interval) : # f() values at end points not of opposite sign test3<-simplesimint(coef=coefs, vcov=VarCov, df=114, cmat=cmat2[1:2,1:4], alternative="two.sided") summary(test3) plotCI(test3) test4<-simplesimint(coef=coefs, vcov=VarCov, df=114, cmat=cmat2[1:3,1:4], alternative="two.sided") summary(test4) plotCI(test4) test5<-simplesimint(coef=coefs, vcov=VarCov, df=114, cmat=cmat2[1:4,1:4], alternative="two.sided") summary(test5) plotCI(test5) test6<-simplesimint(coef=coefs, vcov=VarCov, df=114, cmat=cmat2[1:5,1:4], alternative="two.sided") ##ERRO!! #Erro em uniroot(pfct, interval = interval) : # f() values at end points not of opposite sign

1 0