Fevereiro 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

Re: [R-br] linha da RMLE
by Hélio Gallo Rocha 28 Fev '14

28 Fev '14

Valeu Éder Era uma questão muito simples, mas não tava conseguindo resolver, obrigado! Não saiu na lista não. Acho que quando a gente responde, o endereço selecionado é o da pessoa. Com relação a deixar as linhas do contorno dentro dos limites do poligono, resolvi. É que tinha suprimido a linha: gd$borders <- border. Abraço e bom carnaval Em 28 de fevereiro de 2014 08:37, Éder Comunello <comunello.eder(a)gmail.com>escreveu: > Parece que não postou na lista... > > Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> > Dourados, MS - [22 16.5'S, 54 49'W] > > > ---------- Mensagem encaminhada ---------- > De: Éder Comunello <comunello.eder(a)gmail.com> > Data: 27 de fevereiro de 2014 23:34 > Assunto: Re: [R-br] linha da RMLE > Para: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > > > Hélio, boa noite! > > Supondo o ajuste do variograma por geoR::likfit(), onde vA representa o > objeto variograma e LFA.*, os ajustes (objetos likfit): > > ### Não reproduzível! > plot(vA, type='n') ### plota estrutura do gráfico, omitindo pontos > lines(LFA.exp, col=2, lty=1, lwd=2) ### adiciona linha a partir de objeto > likfit > lines(LFA.sph, col=3, lty=2, lwd=2) ### outra linha no mesmo gráfico > > > Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> > Dourados, MS - [22 16.5'S, 54 49'W] > > > -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

1 0

linha da RMLE
by Hélio Gallo Rocha 28 Fev '14

28 Fev '14

Boa tarde a todos, gostaria de uma ajuda dos listeiros. Teria como plotar apenas a linha estimada por RMLE sem os pontos do semivariograma? -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 2

Kappa por categorias
by branco.sk8＠ig.com.br 27 Fev '14

27 Fev '14

Gentileza, Alguém sabe me dizer como se faz para fazermos os kappas por categorias (em função de outras variáveis)? Ex.: fazer o kappa do teste A com o teste B em função da variavel C (sexo, ocupação etc.). obrigado

1 0

Tempo na Sobrevida
by Fátima Lima Paula 27 Fev '14

27 Fev '14

Olá a todos. Estou trabalhando a análise de sobrevivência com um banco. Utilizei o tempo de forma clássica sem problemas. Agora quero utilizar o método de contagem e para isso preciso criar uma coluna com o tempo final de cada um. Só que o tempo final de uns é o tempo de readmissão, de outros é o tempo até o óbito e de outros é o tempo da saída até o final do estudo, pois não readmitiram nem morreram. Bem, o primeiro problema que estou enfrentando é que quando uso ifelse, a data muda de forma. E eu queria que continuasse como data. Segue um exemplo para que vocês possam me ajudar nome=c("Ricardo","Maria","Tereza","Pedro","José", "Germano","Ricardo") sobrenome=c("Cunha","Andrade","Silva","Soares","Silva","Lima","Cunha") dataint=c(20121004,20120907,20120703,20120301,20120805,20121209,20120508) datasaida=c(20130101,20121006,20120912,20120402,20130104,20130102,20120520) sexo=c("M","F","F","M","M","M","M") tentre=c(20,12,67,04,49,27,08) banco1=cbind(nome,sobrenome,sexo,dataint,datasaida,tentre) banco1=data.frame(banco1) banco1 str(banco1) banco1$tentre=as.character(banco1$tentre) banco1$dataint=as.Date(banco1$dataint,"%Y%m%d") banco1$datasaida=as.Date(banco1$datasaida,"%Y%m%d") banco1$datafim=ifelse(banco1$tentre>30,"0",banco1$datasaida) head(banco1) nome sobrenome sexo dataint datasaida tentre datafim 1 Ricardo Cunha M 2012-10-04 2013-01-01 20 15706 2 Maria Andrade F 2012-09-07 2012-10-06 12 15619 3 Tereza Silva F 2012-07-03 2012-09-12 67 0 4 Pedro Soares M 2012-03-01 2012-04-02 4 0 5 José Silva M 2012-08-05 2013-01-04 49 0 6 Germano Lima M 2012-12-09 2013-01-02 27 15707 7 Ricardo Cunha M 2012-05-08 2012-05-20 8 0 Ao invés de continuar as datas do datasaida ele muda a forma. Alguém poderia me ajudar, por favor? Obrigada Fátima -- "Minha felicidade depende da qualidade dos meus pensamentos"

2 3

Foreach + neuralnet
by Vinicius Brito Rocha 26 Fev '14

26 Fev '14

Senhores, ao tentar utilizar o comando foreach para ajustar redes neurais de diferentes fontes recebo a mensagem de que não existe a função neuralnet. Mesmo tendo o pacotes devidamente instalado e carregado. Alguém pode me ajudar?

2 1

Família ply
by Tropidurus Torquatus 26 Fev '14

26 Fev '14

Caros R-users Dada uma matrix qualquer: seq = c("a","a","a","a","a","a","b","b","b","b","b","c","c","c","c","c","c","c","c","c","c","c","c","c”) var1 = c(1,2,3,4,5,6,0,1,2,3,4,0,1,2,3,4,5,6,7,8,9,10,11,12) matrix.x = data.frame(seq,var1) funcao.x = function(x,d=1,T=3) { i = exp ( ( -x ) / ( d * T ) ) net = sum (i) resultado = 100 * ( i / net ) return (resultado) } Preciso aplicar a função dentro de cada subgrupo (seq) da matrix.x e criar ou um novo data.frame ou uma nova linha da matrix.x. Tentei utilizar o comando ddply do pacote, com o comando abaixo, não sei se essa seria a melhor forma, mas mesmo assim não estou conseguindo: library(“plyr”) R2 <- ddply(matrix.x, 'seq', transform, var=funcao.x(matrix.x$var1)) Error in data.frame(list(seq = c(2L, 2L, 2L, 2L, 2L), var1 = c(0, 1, 2, : arguments imply differing number of rows: 5, 24 Na verdade essa formula deveria dar um valor para cada linha… não sei se escolhi a função correta para o que necessito, alguém poderia me mostrar um caminho a seguir? No mais agradeço, Thales Siqueira

3 3