Abril 2012 - R-br - DINF UFPR

Séries Temporais no R
by Sérgio Henrique almeida da silva ju 25 Abr '12

25 Abr '12

Olá Estou com um conjunto de dados de Hanseníase para os anos de 2001 a 2009 e gostaria de fazer algumas análises de séries temporais. Vocês indicariam referências e tutorial sobre o tema no R? Abraços Em 18 de abril de 2012 12:16, Benilton Carvalho <beniltoncarvalho(a)gmail.com>escreveu: > talvez > > f(dados1$linhaa, "A39|G00") > > mas e' apenas achometro. > > boa sorte. > > 2012/4/18 Sérgio Henrique almeida da silva ju <sergio.edfisica(a)gmail.com>: > > Se por exemplo eu busca além do código A39 querer também o G00 > > > > dados1$causaa = f(dados1$linhaa, "A39") > > dados1$causab = f(dados1$linhab, "A39") > > dados1$causac = f(dados1$linhac, "A39") > > dados1$causad = f(dados1$linhad, "A39") > > dados1$causaii = f(dados1$linhaii, "A39") > > dados1$causaa1 = f(dados1$linhaa, "G00") > > dados1$causab1 = f(dados1$linhab, "G00") > > dados1$causac1 = f(dados1$linhac, "G00") > > dados1$causad1 = f(dados1$linhad, "G00") > > dados1$causaii1 = f(dados1$linhaii, "G00") > > E depois agregar essas colunas em uma só! > > > > Tem como eu buscar o código A39 e G00 em um comando só? como por exemplo: > > > > dados1$causaa = f(dados1$linhaa, "A39" | "G00") > > > > Abraços > > > > > > > > Em 18 de abril de 2012 11:56, Benilton Carvalho < > beniltoncarvalho(a)gmail.com> > > escreveu: > > > >> Como sugerido no guia de postagem, sinta-se a vontade para colocar seu > >> arquivo em algum servico como http://www.datafilehost.com/ , assim > >> quem esta' tentando te ajudar tem a chance de entender mais > >> rapidamente qual eh o problema. > >> > >> O que vc quer dizer com "mais de um codigo"? > >> > >> Em teoria, vc precisa apenas de um wrapper para: > >> > >> with(dados1, gsub('NA', '', > >> > >> > paste(dados1$causaa,dados1$causab,dados1$causac,dados1$causad,dados1$causaii, > >> sep=''))) > >> > >> mas sem entender o que exatamente vc quer dizer, minha ajuda acaba por > >> aqui. > >> > >> b > >> > >> 2012/4/18 Sérgio Henrique almeida da silva ju < > sergio.edfisica(a)gmail.com>: > >> > > >> > > >> > ---------- Mensagem encaminhada ---------- > >> > De: Sérgio Henrique almeida da silva ju <sergio.edfisica(a)gmail.com> > >> > Data: 18 de abril de 2012 11:22 > >> > Assunto: Re: [R-br] Duvida String > >> > Para: r-br(a)listas.c3sl.ufpr.br > >> > > >> > > >> > Uma alternativa que encontrei foi > >> > > >> > > >> > f = function(v, dat, v1=sapply(gregexpr(dat,v),"[",i=1)){ > >> > n = ifelse(v1 > 1, nchar(v), NA) > >> > substr(v, v1,n) > >> > } > >> > dados1$causaa = f(dados1$linhaa, "A39") > >> > dados1$causab = f(dados1$linhab, "A39") > >> > dados1$causac = f(dados1$linhac, "A39") > >> > dados1$causad = f(dados1$linhad, "A39") > >> > dados1$causaii = f(dados1$linhaii, "A39") > >> > > >> > dados1$causa<- with(dados1, gsub('NA', '', > >> > > >> > > paste(dados1$causaa,dados1$causab,dados1$causac,dados1$causad,dados1$causaii, > >> > sep=''))) > >> > > >> > > >> > Como faço para buscar mais de um código? Tem que fazer um monte de > >> > comando > >> > como esses? > >> > Existe uma forma mais simples de fazer isso? > >> > > >> > Segue o banco anexo > >> > > >> > > >> > Desculpe enviar o arquivo anexo, mas é uma tentativa para vocês me > >> > ajudarem, > >> > já que o dput não é uma boa saida! > >> > > >> > Abraços > >> > > >> > > >> > > >> > > >> > -- > >> > Sérgio Henrique Almeida da Silva Junior > >> > Doutorando em Epidemiologia em Saúde Pública > >> > Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ > >> > http://lattes.cnpq.br/1611345552843383 > >> > Tel: (21) 94429486/78101651 id: 123*20942 > >> > > >> > > >> > > >> > > >> > > >> > -- > >> > Sérgio Henrique Almeida da Silva Junior > >> > Doutorando em Epidemiologia em Saúde Pública > >> > Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ > >> > http://lattes.cnpq.br/1611345552843383 > >> > Tel: (21) 94429486/78101651 id: 123*20942 > >> > > >> > > >> > > >> > _______________________________________________ > >> > 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. > >> _______________________________________________ > >> 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. > > > > > > > > > > -- > > Sérgio Henrique Almeida da Silva Junior > > Doutorando em Epidemiologia em Saúde Pública > > Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ > > http://lattes.cnpq.br/1611345552843383 > > Tel: (21) 94429486/78101651 id: 123*20942 > > > > > > > > _______________________________________________ > > 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. > _______________________________________________ > 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. > -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

2 1

Divergencia entre o valor do 1 quartil
by Daniel Marcelino 25 Abr '12

25 Abr '12

Ola caros, um colega me perguntou no exemplo abaixo porque o R calcula o primeiro quartil como sendo 1.75 e no boxplot ele apresenta como sendo 1.5. q8<-c(1,1,3,2,3,2,3,2) summary(q8) boxplot(q8) summary(q8) Min. 1st Qu. Median Mean 3rd Qu. Max. 1.000 1.750 2.000 2.125 3.000 3.000

4 3

Re: [R-br] Divergencia entre o valor do 1 quartil
by Ivan Bezerra Allaman 25 Abr '12

25 Abr '12

Complementando... quartil1 <- NULL for(i in 1:9){ quartil1[i] <- quantile(q8,0.25,type=i) } quartil1 quartil3 <- NULL for(i in 1:9){ quartil3[i] <- quantile(q8,0.75,type=i) } quartil3 (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

Informação arquivo xls - Detalhes
by Eder David Borges da Silva 25 Abr '12

25 Abr '12

Boa noite pessoal, Venho mais uma vez recorrer as suas sabias ajudas, em: http://www.leg.ufpr.br/~eder/teste.xls deixei um arquivo xls, ###Caso queiram carregar ele: download.file('http://www.leg.ufpr.br/~eder/teste.xls','teste.xls') file.info('teste.xls') # Minha duvida é a seguinte, se clicarmos com o direito e ir nas abas: # detalhes -> descrição # Veremos que tem a opção de ter um titulo, assunto, marca etc. # minha missão é modificar este titulo (esta como MEU TITULO) para outro nome. # via comando... # Alguma sugestão? # condições de contorno: SO Windows, R 2.15. Obrigado pessoal

2 2

instalação de pacotes
by André Figueiras Rutz 25 Abr '12

25 Abr '12

Caros, Estou tendo dificuldades com o R na instalação de pacotes alguem poderia me ajudar? Parece que o pacote vai para o micro mas de alguma forma na hora de carregar nao funciona sugestões? Grato André Rutz > install.packages("rJava") --- Please select a CRAN mirror for use in this session --- tentando a URL ' http://cran-r.c3sl.ufpr.br/bin/windows/contrib/2.15/rJava_0.9-3.zip' Content type 'application/zip' length 746108 bytes (728 Kb) URL aberta downloaded 728 Kb package ‘rJava’ successfully unpacked and MD5 sums checked The downloaded binary packages are in C:\Users\user\AppData\Local\Temp\RtmpoBqDjS\downloaded_packages > require(rJava) Carregando pacotes exigidos: rJava Error : .onLoad falhou em loadNamespace() para 'rJava', detalhes: chamada: fun(libname, pkgname) erro: JAVA_HOME cannot be determined from the Registry > require("rJava") Carregando pacotes exigidos: rJava Error : .onLoad falhou em loadNamespace() para 'rJava', detalhes: chamada: fun(libname, pkgname) erro: JAVA_HOME cannot be determined from the Registry

3 4

Como saber se é um 'bug' ou uma 'feature', no comportamento de 'plot'?
by Cleber N.Borges 24 Abr '12

24 Abr '12

olá a todos da lista, estava fazendo uns testes e me deparei com um comportamento do qual *acho* que pode ser incoerente. Vamos ao problema: > par(mfrow=c(2,2)) > plot(1:10) # grafico 01 > plot(1:10, cex.axis=0.5, las=1) # grafico 02 > plot(table(1:10)) # grafico 03 > plot(table(1:10), cex.axis=0.5, las=1) # grafico 04 a mudança dos parametros: 'cex.axis' e 'las' funcionam muito bem na mudança do gráfico 1 para 2, mas não há efeito na mudança de 3 para 4. Segunda pergunta: Como saber qual é a verdadeira função 'plot.ALGUMACOISAAQUI' que será chamada para fazer o grafico de um objeto da classe 'table'? dado que não há 'plot.table' como função... Grato pela atenção de todos. Cleber

2 3

Dúvida sobre ANOVA tipo III
by Olympio Neto 24 Abr '12

24 Abr '12

Pessoal, Tenho lido vários livros sobre R, muitos blogs (R-bloggers) e as excelentes páginas do Wiki-Leg da UFPR e as famosas Ridículas do Walmes e em muitas delas vejo que em muitos dos tutoriais, é necessária, para a obtenção da Anova Marginal (Tipo III), setar o contraste como "contr.sum" e outros que nada falam sobre isso, ficando setado o padrão do R, o "contr. treatment". A explicação para a alteração do contraste é que a soma deste deve ser nula, o que não acontece com o "contr.treatment". Caso não seja feita a alteração, os resultados podem estar errados. O artigo que tenta por luz neste assunto é aqui indicado: http://goanna.cs.rmit.edu.au/~fscholer/anova.php Normalmente eu uso o pacote "car" para cálculos que análises de variância do tipo III, mas nos exemplos que vejo, dados que simulo etc, não vejo diferenças. Acredito que nem poderiam haver, mesmo. A ANOVA trivial (tipo I) é invariante à parametrização, mas fico muito em dúvida se a tipo III também seria. O que vocês pensam? Já observaram diferenças importantes em suas análises? Abs e muito obrigado. Olympio

4 10

Paste
by Ricardo Fernandes Campos Junior 24 Abr '12

24 Abr '12

Olá, pessoal. Eu sei que essa pergunta é provavelmente muito simples, mas eu não consigo encontrar como fazê-lo. Eu tenho um objeto de caracteres de tamanho genérico, para fim de exemplo, vou usar 3. Eu quero criar um outro objeto com estes mesmos caracteres, mas concatenados. Porém não consigo. Tentei isso: exemplo = c("A", "B", "C") novoobjeto = paste(exemplo, " ") O que me é retornado é isso: "A" "B" "C" Mas o que eu quero é isso: "A B C" Alguém pode ajudar? -- Abraços, Ricardo Fernandes Campos Junior Departamento de Genética e Biologia Evolutiva Universidade de São Paulo

3 3

Notação
by Edson Lira 24 Abr '12

24 Abr '12

Tenho esse script abaixo: require(PASWR) trat<-rep(LETTERS[1:4],5) med<-c(25,31,22,33,26,25,26,29,20,28,28,31,23,27,25,34,21,24,29,28) exem1_4<-data.frame(med,trat) exem1_4 m1<-lm(exem1_4$med~exem1_4$trat) anova(m1) summary(m1) erro<-resid(m1) erro est<-predict(m1) med-est Quando gero o erro ele me dá essa saída. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 2.00e+00 4.00e+00 -4.00e+00 2.00e+00 3.00e+00 -2.00e+00 -3.33e-16 -2.00e+00 -3.00e+00 1.00e+00 2.00e+00 -1.11e-16 8.88e-16 1.11e-16 -1.00e+00 3.00e+00 -2.00e+00 18 19 20 -3.00e+00 3.00e+00 -3.00e+00 Gostaria de transformar estes números para a notação comum, tipo 2, 3, -4 .... Me lembro que certa vez alguém postou esse problema, mas faz muito tempo. Alguém tem alguma sugestão? Edson Lira Estatístico Manaus-Amazonas

3 2

[Dúvida] Números aleatórios.
by Pedro Rafael 24 Abr '12

24 Abr '12

Formula em LaTeX: F_{T}(t;\alpha,\beta) = \Phi(\frac{1}{\alpha}(\frac{t}{\beta}^(1/2) - \frac{\beta}{t}^(1/2))) Pessoal com gerar números aleatórios com distribuição Birbaum-Saunder em C? Tenho a distribuição acumulada da Birbaum-Saunder que é escrita como acumulada de uma normal e tenho a densidade da Birbaum. A formula da acumulada está escrita acima para t e beta maiores que zero. Pensei primeiramente em pegar números aleatórios normais pela biblioteca GSL e fazer ser igual à t na formula acumulada acima. Daria certo? Se for pela transformada inversa eu teria que gerar números uniformes[0,1] e fazer F(t) = U. Depois eu tinha que ver qual o t em F me daria U, sendo U o número uniforme gerado. Outra complicação é que tenho uma acumulada da normal que é uma integral de -inf a t. Existe outra forma mais fácil? -- Saudações, Pedro Rafael Diniz Marinho.

2 1