- R-br - DINF UFPR

Dúvida_experimentação_estimação
by andrebvs＠bol.com.br 02 Jan '13

02 Jan '13

3 2

Re: [R-br] R nas Nuvens
by Tiago Vieira 31 Dez '12

31 Dez '12

Ok Crysttian, Vou tentar realizar esta tarefa usando sua dica. Obrigado! Att, Tiago Vieira. ________________________________ De: Crysttian Arantes Paixão <crysttian(a)gmail.com> Para: tiagovieira7(a)yahoo.com.br Enviadas: Segunda-feira, 31 de Dezembro de 2012 8:33 Assunto: R nas Nuvens Olá Tiago, Do que você precisa, o CENAPAD que eu falei te fornece de graça. É um cluster na Unicamp, que possui o R instalado. Para poder usar, basta preencher um formulário e mandar para os administradores. Vale a pena dar uma olhada no site que te passei. Qualquer dúvida é só falar. Até Crysttian Arantes Paixão Escola de Matemática Aplicada - Fundação Getulio Vargas - FGV/RJ Cel: 021 8038-3085 (Tim) Linux User: 456734 - Debian “Se fiz descobertas valiosas, foi mais por ter paciência do que qualquer outro talento.” [Isaac Newton] "Otimismo é esperar pelo melhor. Confiança é saber lidar com o pior." [Roberto Simonsen] "Mesmo as noites totalmente sem estrelas podem anunciar a aurora de uma grande realização." [Martin Luther King] OUVIR ESTRELAS "Ora (direis) ouvir estrelas! Certo Perdeste o senso!" E eu vos direi, no entanto, Que, para ouvi-las, muita vez desperto E abro as janelas, pálido de espanto... E conversamos toda a noite, enquanto A via-láctea, como um pálio aberto, Cintila. E, ao vir do sol, saudoso e em pranto, Inda as procuro pelo céu deserto. Direis agora: "Tresloucado amigo! Que conversas com elas? Que sentido Tem o que dizem, quando estão contigo?" E eu vos direi: "Amai para entendê-las! Pois só quem ama pode ter ouvido Capaz de ouvir e de entender estrelas." (Olavo Bilac, Poesias, Via-Láctea, 1888.)

1 0

Ajuda com gráfico de pontos
by Julia Maria 29 Dez '12

29 Dez '12

Caros colegas ... Quando vou rodar este script abaixo eu não consigo fazer com que as observações saiam em cada mês ... Esta saaindo apenas um pontinho preto. Será que algum colega poderia me auxiliar neste script? Os dados estão em anexo. Grata dados <- read.table("D:\\Dados.txt",h=F) dadost <- as.matrix(dados[,3]) x <- c(710000,712000,712000,710000) y <- c(7700000,7700000,7900000, 7900000) poly <- cbind(x,y) plot(poly) n <- max(dadost) ############################## ############################## ############################ GRAFICO DE TODOS OS TEMPOS ############################## ############################## ############################## dts <- c("Nov2001","Dez2001"," Jan2002","Fev2002","Mar2002", "Abr2002","Maio2002","Jun2002" ,"Jul2002","Ago2002","Set2002" ) par(mfrow=c(7,3),mar=rep(2,4)) for(i in 1:n){ plot(poly,type="n",xlab="x", ylab="y",main=dts[i]) points(dados[,1][dadost==i], dados[,2][dadost==i],xlab="x", ylab="y",main="")}

2 1

Estrutura interna de arquivos no R
by Daniel Marcelino 27 Dez '12

27 Dez '12

Caros, estou com o seguinte problema. Quero fazer o download de um shape e extraí-lo, mas parece que eu não estou achando o caminho exato dos arquivos. Alguém tem alguma cola pra isso? fl <- "http://metodologiapolitica.com/download/Brasil.zip" download.file(fl, "Brasil.zip") br <- unzip("Brasil.zip") > br [1] "./Brasil.dbf" "./__MACOSX/._Brasil.dbf" "./Brasil.shp" "./__MACOSX/._Brasil.shp" [5] "./Brasil.shx" "./__MACOSX/._Brasil.shx" shp <- readShapePoly(system.file("./Brasil.shp") ) Error in getinfo.shape(filen) : Error opening SHP file > shp <- readShapePoly(system.file("./__MACOSX/._Brasil.shp") ) Error in getinfo.shape(filen) : Error opening SHP file > shp <- readShapePoly("./__MACOSX/._Brasil.shp") Error in getinfo.shape(filen) : Error opening SHP file -- "Small steps toward a much better world" \begin{signature} Daniel Marcelino ☁ dm.silva(a)umontreal.ca ☎ (514) 343 6111 #3799 ✎ 3200 Jean Brillant, Office C5071 Montreal, QC; H3T 1N8 Canada \end{signature}

2 2

Script S100 revisão 1
by Hélio Gallo Rocha 27 Dez '12

27 Dez '12

Caros, Com havia implementado os gráficos(envelopamento, seria isso?) de detecção de dependência espacial, vi que tinha um errinho, ele gerava sempre sobre a mesma distância, ajustei... Também mudei a sequencia do script, coloquei quatro gráficos na mesma tela, com a seguinte ideia: 1. visualizar os semivariogramas e verificar, ao mesmo tempo: 1.1 a existência de dep. espacial 1.2 auxilia a escolha da dist. máxima O ultimo gráfico é com a distância escolhida, dai pra frente todo o script vai rodar com ela... espero ter sido claro na exposição das ideias, estou a disposição. estou afastado para mestrado, assim tenho meu tempo pra isso. gostaria de um retorno de quem já testou o script, ok? Outra coisa, fiz o ajuste das imagens de krigagem, mas espero que alguém tenha uma ideia mais fácil -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

4 18

Ajuda ggmap
by Sérgio Henrique almeida da silva ju 27 Dez '12

27 Dez '12

Olá Estou trabalhando com o ggmap e gostaria de plotar a taxa de uma doença usando o cols[findInterval...]. Isso é possível ou só posso usar o geom_point? Meu script é: stores <- data.frame(name=c("Jaboatão dos Guararapes", "Olinda","Paulista", "Recife", "Maceió"), longitude=c(-35.00329,-34.86707,-34.88437,-34.93294, -35.71009), latitude=c(-8.152749,-7.993280,-7.927770,-8.039684,-9.523036), disease=c(143,226,24,142,14)) location = c(-45.6,-10) nordeste = get_map(location = location, source = "osm", zoom=5) nordesteMap = ggmap(nordeste) Obrigado -- 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) 68463637 / 94429486

1 0

(sem assunto)
by Marcelo Claro de Souza 27 Dez '12

27 Dez '12

http://webdesignersmn.com/wp-content/themes/twentyeleven/images/perfectway.…

1 0

CONCURSO_ESTATISTICA
by NERIO CARDOSO 27 Dez '12

27 Dez '12

Bom dia, Não sei se este é o melhor caminho, mas é um meio. Peço ajuda para divulgar o edital de concurso para Universidade federal de Rondônia, teremos 6 vagas para graduandos em Estatísticas, entre outras vagas. http://www.unir.br/noticias_arquivos/7351_edital_006_final_professor_2012.p… Coordialmente Estatístico Nerio Aparecido Cardoso Professor do Departamento de Matemática e Estatística Fundação Universidade Federal de Rondônia - UNIR Campus de Ji-Paraná-RO Rua Rio Amazonas, Nº351 - Jardins dos Migrantes Cep: 76900-726 - Ji-Paraná-RO, Brasil Fone: +55 69 34213595 +55 69 81106457 +55 69 84546930 Aos que seguem rigorosamente protocolos: "Nos experimentos científicos o que vemos não é a natureza em si, mas a natureza submetida ao nosso modo peculiar de interrogá-la. (Heisenber)" > From: r-br-request(a)listas.c3sl.ufpr.br > Subject: Digest R-br, volume 22, assunto 36 > To: r-br(a)listas.c3sl.ufpr.br > Date: Thu, 27 Dec 2012 12:00:01 -0200 > > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > > Tópicos de Hoje: > > 1. Re: Script S100 revisão 1 (Mauro Sznelwar) > 2. Intervalos de credibilidade e HPD (Edmundo Caetano) > 3. Apostila de modelos mistos (Alisson Lucrecio) > 4. Re: Apostila de modelos mistos (Walmes Zeviani) > 5. Re: Apostila de modelos mistos (Gilbert Queiroz) > 6. Re: Intervalos de credibilidade e HPD (Ivan Bezerra Allaman) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Thu, 27 Dec 2012 03:34:43 -0200 > From: "Mauro Sznelwar" <sznelwar(a)uol.com.br> > To: <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Script S100 revisão 1 > Message-ID: <9902AE9B77264B8BB4BCB829535CE066@ACER4520> > Content-Type: text/plain; charset="iso-8859-1" > > Já procurei e não encontrei. Não o tem para enviar? > > Prezado Mauro, > > > se fizeres uma pesquisa nos últimos posts da lista, verás que o arquivo "recorte.txt" está anexado > > a um dos posts anteriores do Sr. Hélio. Ok? > > > Feliz Natal a você e a todos da lista, > > > Boas festas. > > > cordialmente, > > > > _______________ > Rodrigo Sant'Ana > Oceanógrafo > > > . > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121227/bd8e99cf/atta…> > > ------------------------------ > > Message: 2 > Date: Thu, 27 Dec 2012 04:39:12 -0800 (PST) > From: Edmundo Caetano <eddycaetano(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Intervalos de credibilidade e HPD > Message-ID: > <1356611952.56394.YahooMailNeo(a)web142503.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Prezados, > Estou com uma dúvida pretendo estabelecer intervalos > > de credibilidade ou HPD para alguns parâmetros. Alguém pode > me ajudar. > > Abrços > > Edmundo Caetano > "Ninguém é tão sábio que nada tenha para aprender, nem tão tolo que nada tenha para ensinar." > Blaise Pascal > > > > ________________________________ > De: "r-br-request(a)listas.c3sl.ufpr.br" <r-br-request(a)listas.c3sl.ufpr.br> > Para: r-br(a)listas.c3sl.ufpr.br > Enviadas: Quarta-feira, 26 de Dezembro de 2012 12:00 > Assunto: Digest R-br, volume 22, assunto 35 > > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > > Tópicos de Hoje: > > 1. Re: Script S100 revisão 1 (Rodrigo Sant'Ana) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 26 Dec 2012 00:15:09 -0200 > From: "Rodrigo Sant'Ana" <rodrigo.gringo(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Script S100 revisão 1 > Message-ID: > <CACpsV-wF1SV6EKoBb4A6RgAeTG8MQG3XOBP8LQyA_Vb-4+Pr+w(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Prezado Mauro, > > se fizeres uma pesquisa nos últimos posts da lista, verás que o arquivo > "recorte.txt" está anexado > a um dos posts anteriores do Sr. Hélio. Ok? > > Feliz Natal a você e a todos da lista, > > Boas festas. > > cordialmente, > > _______________ > Rodrigo Sant'Ana > Oceanógrafo > > > Em 24 de dezembro de 2012 04:59, Mauro Sznelwar <sznelwar(a)uol.com.br>escreveu: > > > ** > > > > Alguém teria o arquivo recorte.txt que faz parte do Script? > > > > Pessoal, > > > > seguindo a proposta do Sr. Benilton, criei um gist publico para > > otimização desta discussão. > > Assim podemos maximizar a troca de informações e experiências sobre estas > > análises, e quaisquer > > outras análises vindouras, sem criar uma rotina de emails cheios de anexos > > para a lista da R-BR. > > > > Acho que a ideia do grupo de estudos super válida, e ela se tornará ainda > > melhor se utilizarmos algo > > como o gist para otimizar nosso trabalho. Fica a dica, e o link para > > darmos prossegmento nesta > > discussão. > > > > link: https://gist.github.com/4262669 > > > > Como primeiro passo, a titulo de exemplo, fiz uma pequena proposta de > > alteração no script original > > proposto pelo Sr. Hélio, simplesmente alterei a forma de carregar a base > > de dados s100 prevista no > > pacote geoR. > > > > Quem nunca utilizou o gist, é simples, para inserir novas alterações no > > script, sejam novas análises ou mesmo > > otimizações no código, basta clicar no botão "Revisions", após clicar em > > "edit", realizar as inserções, > > alterações, sugestões, etc, e clicar em "Update Gist". > > Todas as revisões ficaram salvas e para visualiza-las bastará clicar no > > botão "Revisions" novamente. > > > > Peço desculpas caso tenha sido muito invasivo, mas acredito que este seja > > um caminho interessante para > > adotarmos, assim como inúmeros outros exemplos sugeridos anteriormente. > > > > Boa noite a todos, > > > > _______________ > > Rodrigo Sant'Ana > > Oceanógrafo > > > > > > > > . > > > > > > _______________________________________________ > > 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. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121226/a3909095/atta…> > > ------------------------------ > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > Fim da Digest R-br, volume 22, assunto 35 > ***************************************** > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121227/ad94b6c7/atta…> > > ------------------------------ > > Message: 3 > Date: Thu, 27 Dec 2012 05:05:59 -0800 (PST) > From: Alisson Lucrecio <alissonluc(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Apostila de modelos mistos > Message-ID: > <1356613559.28071.YahooMailNeo(a)web160706.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Caros colegas, > > Boa tarde. > > Alguém poderia me sugerir alguma apostila ou livro para modelos mistos? > > Obrigado. > > Alisson Lucrécio da Costa > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121227/751d1693/atta…> > > ------------------------------ > > Message: 4 > Date: Thu, 27 Dec 2012 11:13:09 -0200 > From: Walmes Zeviani <walmeszeviani(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, Alisson Lucrecio > <alissonluc(a)yahoo.com.br> > Subject: Re: [R-br] Apostila de modelos mistos > Message-ID: > <CAFU=EkYsNJ+h3a46dip-JP8fCNHVMoBPWFsO3ghgN8Jt6Z4B=g(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Livro do Pinheiro e Bates, principalmente para quem usa R, é indicado. Uma > boa revisão de literatura você pode fazer no google com termos > relacionados/derivados dos termos "mixed models", "lectures/class notes", > "overview", "introduction". Em geral os resultados são textos curtos e > consultar fontes de diversos autores é enriquecedor. > > À disposição. > Walmes. > > ========================================================================== > Walmes Marques Zeviani > LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) > Departamento de Estatística - Universidade Federal do Paraná > fone: (+55) 41 3361 3573 > VoIP: (3361 3600) 1053 1173 > e-mail: walmes(a)ufpr.br > skype: walmeszeviani > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121227/5087493f/atta…> > > ------------------------------ > > Message: 5 > Date: Thu, 27 Dec 2012 05:38:17 -0800 (PST) > From: Gilbert Queiroz <gilbert_queiroz(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Apostila de modelos mistos > Message-ID: > <1356615497.97469.YahooMailNeo(a)web162704.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Walmes, > > É possível aplicar os modelos mistos em dados espaço-temporais???? > > Abs. > > Gilbert. > > > > > ________________________________ > De: Walmes Zeviani <walmeszeviani(a)gmail.com> > Para: r-br(a)listas.c3sl.ufpr.br; Alisson Lucrecio <alissonluc(a)yahoo.com.br> > Enviadas: Quinta-feira, 27 de Dezembro de 2012 11:13 > Assunto: Re: [R-br] Apostila de modelos mistos > > > Livro do Pinheiro e Bates, principalmente para quem usa R, é indicado. Uma boa revisão de literatura você pode fazer no google com termos relacionados/derivados dos termos "mixed models", "lectures/class notes", "overview", "introduction". Em geral os resultados são textos curtos e consultar fontes de diversos autores é enriquecedor. > > À disposição. > Walmes. > > > ========================================================================== > Walmes Marques Zeviani > LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) > Departamento de Estatística - Universidade Federal do Paraná > fone: (+55) 41 3361 3573 > VoIP: (3361 3600) 1053 1173 > e-mail: walmes(a)ufpr.br > skype: walmeszeviani > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > > _______________________________________________ > 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. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121227/fc5c95b8/atta…> > > ------------------------------ > > Message: 6 > Date: Thu, 27 Dec 2012 05:55:53 -0800 (PST) > From: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> > To: R Brasil <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Intervalos de credibilidade e HPD > Message-ID: > <1356616553.76275.YahooMailNeo(a)web162303.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Embora sua dúvida seja muito ampla, o pacote "languageR" faz intervalos HPD. Sobre o uso da técnica, procure um estatístico local que entenda de inferência bayesiana. > > (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} > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121227/7b98943a/atta…> > > ------------------------------ > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > Fim da Digest R-br, volume 22, assunto 36 > *****************************************

1 0

Apostila de modelos mistos
by Alisson Lucrecio 27 Dez '12

27 Dez '12

Caros colegas, Boa tarde. Alguém poderia me sugerir alguma apostila ou livro para modelos mistos? Obrigado. Alisson Lucrécio da Costa

3 3

Re: [R-br] Intervalos de credibilidade e HPD
by Ivan Bezerra Allaman 27 Dez '12

27 Dez '12

Embora sua dúvida seja muito ampla, o pacote "languageR" faz intervalos HPD. Sobre o uso da técnica, procure um estatístico local que entenda de inferência bayesiana. (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

Intervalos de credibilidade e HPD
by Edmundo Caetano 27 Dez '12

27 Dez '12

Prezados, Estou com uma dúvida pretendo estabelecer intervalos de credibilidade ou HPD para alguns parâmetros. Alguém pode me ajudar. Abrços Edmundo Caetano "Ninguém é tão sábio que nada tenha para aprender, nem tão tolo que nada tenha para ensinar." Blaise Pascal ________________________________ De: "r-br-request(a)listas.c3sl.ufpr.br" <r-br-request(a)listas.c3sl.ufpr.br> Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Quarta-feira, 26 de Dezembro de 2012 12:00 Assunto: Digest R-br, volume 22, assunto 35 Enviar submissões para a lista de discussão R-br para r-br(a)listas.c3sl.ufpr.br Para se cadastrar ou descadastrar via WWW, visite o endereço https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ou, via email, envie uma mensagem com a palavra 'help' no assunto ou corpo da mensagem para r-br-request(a)listas.c3sl.ufpr.br Você poderá entrar em contato com a pessoa que gerencia a lista pelo endereço r-br-owner(a)listas.c3sl.ufpr.br Quando responder, por favor edite sua linha Assunto assim ela será mais específica que "Re: Contents of R-br digest..." Tópicos de Hoje: 1. Re: Script S100 revisão 1 (Rodrigo Sant'Ana) ---------------------------------------------------------------------- Message: 1 Date: Wed, 26 Dec 2012 00:15:09 -0200 From: "Rodrigo Sant'Ana" <rodrigo.gringo(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Script S100 revisão 1 Message-ID: <CACpsV-wF1SV6EKoBb4A6RgAeTG8MQG3XOBP8LQyA_Vb-4+Pr+w(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Prezado Mauro, se fizeres uma pesquisa nos últimos posts da lista, verás que o arquivo "recorte.txt" está anexado a um dos posts anteriores do Sr. Hélio. Ok? Feliz Natal a você e a todos da lista, Boas festas. cordialmente, _______________ Rodrigo Sant'Ana Oceanógrafo Em 24 de dezembro de 2012 04:59, Mauro Sznelwar <sznelwar(a)uol.com.br>escreveu: > ** > > Alguém teria o arquivo recorte.txt que faz parte do Script? > > Pessoal, > > seguindo a proposta do Sr. Benilton, criei um gist publico para > otimização desta discussão. > Assim podemos maximizar a troca de informações e experiências sobre estas > análises, e quaisquer > outras análises vindouras, sem criar uma rotina de emails cheios de anexos > para a lista da R-BR. > > Acho que a ideia do grupo de estudos super válida, e ela se tornará ainda > melhor se utilizarmos algo > como o gist para otimizar nosso trabalho. Fica a dica, e o link para > darmos prossegmento nesta > discussão. > > link: https://gist.github.com/4262669 > > Como primeiro passo, a titulo de exemplo, fiz uma pequena proposta de > alteração no script original > proposto pelo Sr. Hélio, simplesmente alterei a forma de carregar a base > de dados s100 prevista no > pacote geoR. > > Quem nunca utilizou o gist, é simples, para inserir novas alterações no > script, sejam novas análises ou mesmo > otimizações no código, basta clicar no botão "Revisions", após clicar em > "edit", realizar as inserções, > alterações, sugestões, etc, e clicar em "Update Gist". > Todas as revisões ficaram salvas e para visualiza-las bastará clicar no > botão "Revisions" novamente. > > Peço desculpas caso tenha sido muito invasivo, mas acredito que este seja > um caminho interessante para > adotarmos, assim como inúmeros outros exemplos sugeridos anteriormente. > > Boa noite a todos, > > _______________ > Rodrigo Sant'Ana > Oceanógrafo > > > > . > > > _______________________________________________ > 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. > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121226/a3909095/atta…> ------------------------------ _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Fim da Digest R-br, volume 22, assunto 35 *****************************************

1 0

Experimento de mistura
by Ari Clecius 24 Dez '12

24 Dez '12

Boa tarde, estou tentando realizar um planejamento de mistura de três componentes onde um dos componentes não pode passar de 0.05 ou 5%, estou usando o pacote qualityTools, algoém poderia me dizer como configurar os limites? Att, Ari Clecius Alves de Lima Engenheiro Químico Me. Engenharia Civil (085)87687786 (085)33669042

1 1

R para Linux
by Gilbert Queiroz 23 Dez '12

23 Dez '12

Pessoal, R para Linux, roda legal??? Abs.

4 3

Re: [R-br] R para Linux
by Sérgio Henrique 23 Dez '12

23 Dez '12

Eu gosto bastante do emacs! Enviado via iPad Em Dec 23, 2012, às 8:29 AM, Rodolfo Timoteo da Silva <zhushazang(a)gmail.com> escreveu: Além do RStudio (rstudio.org) tem o rkward (rkward.sourceforge.net) que podem te auxiliar no uso. Att Em 22-12-2012 20:14, FHRB Toledo escreveu: Roda! 2012/12/22 Gilbert Queiroz <gilbert_queiroz(a)yahoo.com.br> > Pessoal, > R para Linux, roda legal??? > Abs. > > _______________________________________________ > 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 listR-br@listas.c3sl.ufpr.brhttps://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.

1 0

(sem assunto)
by Raul Carlos de Mello 22 Dez '12

22 Dez '12

Boa tarde a todos a lista. Estou fazendo leitura de dados do site da bovespa e melhores e maiores da exame e tenho o seguinte problema que gostaria de ajuda: - Faço a teste da existência do site com o comando; - Ok, tudo funciona sem problema - Após a verificação existência da página busco pelos dados contidos nas tabelas. Entretanto nem todas a tabelas existem dados. Como posso fazer, ao utilizar "readHTMLTable", para evitar que o programa seja interrompido em decorrência de um erro pela tabela não existir. -- Att, Prof. Adm. Raul Carlos de Mello http://lattes.cnpq.br/3784085036882680

1 0

Re: [R-br] Cálculo do percentage prediction error em ZINB
by Humberto Hazin 21 Dez '12

21 Dez '12

Opa Leonard, Muito obrigado vou dar uma olhada! Sim estou usando o pacote pscl para rodar os modelos ZINB! Se eu conseguir resolver o problema coloco na lista. Cara se não me engano tive esse problema também a um tempo atrás, veja se código abaixo para ver se resolve o teu problema. Só uma observação, com esse código não precisas rodar o pacote pscl. https://www.dropbox.com/s/8l480on83dvm3m7/zeroinflCorreto.R Um abraço, Humberto Enviado do Email do Windows De: Leonard de Assis Enviado: ‎19‎ de ‎dezembro‎ de ‎2012 ‎22‎:‎04 Para: r-br(a)listas.c3sl.ufpr.br Assunto: Re: [R-br] Cálculo do percentage prediction error em ZINB Humberto vc está usando o pscl? Eu estou no momento com um problema envolvendo ZINB, só que o meu é um pouco mais chato, pois meu modelo é multivariado (E pra complicar, está dando que uma das dimensões é singular, hehehe) Sugiro ler os tutoriais, vc acha um monte via rseek. Esse aqui é um deles (De antemão, não tem exatamente o que vc precisa, mas tem muita coisa) http://www.ats.ucla.edu/stat/r/dae/zinbreg.htm []s Leonard de Assis http://about.me/ldeassis Em 19/12/2012 17:44, Humberto Hazin escreveu: Prezados, Estou com uma dúvida que pode parecer besteira, porém gostaria da ajuda de vocês! A dúvida é a seguinte como calcular o percentual de "Bias" (percentage prediction error) por linha de um modelo ZINB (package- pscl)? O meu objetivo é fazer um boxplot com esse percentual por ano Eu iniciei da seguinte forma: dois bancos de dados um de try.csv (10% retirados do banco original) e outro original.csv try<-read.csv(”try.csv”, header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE) https://www.dropbox.com/s/v7vruztus11nnbi/try.csv original<-read.csv(”original.csv”, header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE) https://www.dropbox.com/s/7htqc2bqi45rkku/original.csv # Criar Fatores para o banco original original$area<-rep(1,nrow(original)) original$area[original$lat<=-15&original$lat>=-45]<-2 original$A=as.factor(original$area) summary(original) original$S<-factor(original$id.group.fleet) original$Q<-factor(original$Quarter) original$Y<-factor(original$Year) original$logEf<-log(original$effort1) summary(original) # Criar Fatores para o banco try try$area<-rep(1,nrow(try)) try$area[try$lat<=-15&try$lat>=-45]<-2 try$A=as.factor(try$area) summary(try) try$S<-factor(try$id.group.fleet) try$Q<-factor(try$Quarter) try$Y<-factor(try$Year) try$logEf<-try$effort1*0 summary(try) zeroinfl.fit_original<-zeroinfl(BSH~Y+S+Q+A+offset(logEf)|Y+S+Q+A,dist="negbin",data=original) pred<-predict(zeroinfl.fit_original,newdata=try,type="response") Daqui em diante não tenho mais ideia de como gerar esse percentage prediction error Alguem poderia me ajudar! Desde já Agradeço Humberto Enviado do Email do Windows _______________________________________________ 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.

4 4

Junção de dados no mclapply
by Junior Beleti 21 Dez '12

21 Dez '12

Bom dia a todos. Estou realizando experimentos para execução paralela. Para isso criei uma função que realiza a multiplicação de matrizes na sua forma digamos "primitiva", ou seja, como na matemática mesmo, realizando a multiplicação linha por coluna. Sei que existe o operador %*% que já realiza a multiplicação matricial, porém, como já disse antes, estou apenas realizando experimentos. Segue código: m1 <- matrix(0,3,3) #criando matriz m1[1,1]<-1 m1[1,2]<-2 m1[1,3]<-3 m1[2,1]<-4 m1[2,2]<-1 m1[2,3]<-4 m1[3,1]<-2 m1[3,2]<-3 m1[3,3]<-2 m2 <- matrix(0,3,3) #criando matriz m2[1,1]<-2 m2[1,2]<-1 m2[1,3]<-2 m2[2,1]<-3 m2[2,2]<-4 m2[2,3]<-1 m2[3,1]<-2 m2[3,2]<-3 m2[3,3]<-3 m3 <- matrix(0,3,3) #criando matriz numlinha<-3 numcoluna<-3 N<-3 #funcao que realiza a multiplicação myf = function(N,mat1,mat2,num,x,y){ num + (mat1[x,N] * mat2[N,y]) } #forma sequencial (já testada e que está correta) aux <-0 for(i in 1:numlinha){ for(j in 1:numcoluna){ for(k in 1:N){ aux <-myf(k,m1,m2,aux,i,j) } m3[i,j] <- aux aux <-0 } } ## resultado da execucao: ## matriz m1 multiplicada pela matriz m2 que gera a matriz m3 > m1 [,1] [,2] [,3] [1,] 1 2 3 [2,] 4 1 4 [3,] 2 3 2 > m2 [,1] [,2] [,3] [1,] 2 1 2 [2,] 3 4 1 [3,] 2 3 3 > m3 [,1] [,2] [,3] [1,] 14 18 13 [2,] 19 20 21 [3,] 17 20 13 Se eu fosse utilizar a função mclapply, ao invés das tres repetições (for), bastaria colocar: require(parallel) aux <- 0 for(i in 1:numlinha){ mclapply(1:numcoluna,function(h){ for(k in 1:N){ aux <-myf(k,m1,m2,aux,i,h) } m3[i,h] <- aux aux <-0 },mc.cores=2) } O problema é que não consigo armazenar os valores em cada posição da matriz (m3[i,h] <- aux), ou seja, ele não salva na matriz de destino m3, visto que isso ocorre dentro da chamado ao mclapply. Tentei utilizar o "unlist" mas não obtive sucesso. Alguem pode me auxiliar com isso? Att, Junior Beleti.

2 1

ler arquivo excel 2007 com o R (64bit)
by Alexandro (Yahoo) 20 Dez '12

20 Dez '12

#Saudações a todos da lista R-Br, utilizei o comando abaixo no R (64bit) para abrir um arquivo excel 2007 require(RODBC);cnct = odbcConnectExcel2007(file.choose()) #escolhi o arquivo de interesse e apareceu a mensagem: #Carregando pacotes exigidos: RODBC #Mensagens de aviso perdidas: #1: In odbcDriverConnect(con, tabQuote = c("[", "]"), ...) : # [RODBC] ERROR: state IM002, code 0, message [Microsoft][ODBC Driver Manager] Nome da fonte de dados não encontrado e nenhum driver padrão especificado #2: In odbcDriverConnect(con, tabQuote = c("[", "]"), ...) : # ODBC connection failed #Uso processador i5 4GB de memória e windows 7 (64bit). #Quando executei o mesmo comando no R (32bit), funcionou perfeitamente, isto é o que me deixa intrigado! #o que pode ter ocorrido e o que fazer para ter este problema solucionado? #será que este comando no R (64bit) pode não estar implementado? #Muito Obrigado ... #Alexandro

9 10

Cálculo do percentage prediction error em ZINB
by Humberto Hazin 19 Dez '12

19 Dez '12

Prezados, Estou com uma dúvida que pode parecer besteira, porém gostaria da ajuda de vocês! A dúvida é a seguinte como calcular o percentual de "Bias" (percentage prediction error) por linha de um modelo ZINB (package- pscl)? O meu objetivo é fazer um boxplot com esse percentual por ano Eu iniciei da seguinte forma: dois bancos de dados um de try.csv (10% retirados do banco original) e outro original.csv try<-read.csv(”try.csv”, header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE) https://www.dropbox.com/s/v7vruztus11nnbi/try.csv original<-read.csv(”original.csv”, header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE) https://www.dropbox.com/s/7htqc2bqi45rkku/original.csv # Criar Fatores para o banco original original$area<-rep(1,nrow(original)) original$area[original$lat<=-15&original$lat>=-45]<-2 original$A=as.factor(original$area) summary(original) original$S<-factor(original$id.group.fleet) original$Q<-factor(original$Quarter) original$Y<-factor(original$Year) original$logEf<-log(original$effort1) summary(original) # Criar Fatores para o banco try try$area<-rep(1,nrow(try)) try$area[try$lat<=-15&try$lat>=-45]<-2 try$A=as.factor(try$area) summary(try) try$S<-factor(try$id.group.fleet) try$Q<-factor(try$Quarter) try$Y<-factor(try$Year) try$logEf<-try$effort1*0 summary(try) zeroinfl.fit_original<-zeroinfl(BSH~Y+S+Q+A+offset(logEf)|Y+S+Q+A,dist="negbin",data=original) pred<-predict(zeroinfl.fit_original,newdata=try,type="response") Daqui em diante não tenho mais ideia de como gerar esse percentage prediction error Alguem poderia me ajudar! Desde já Agradeço Humberto Enviado do Email do Windows

2 1

problemas com R-br listas
by Fernando A. Souza 19 Dez '12

19 Dez '12

Caros colegas estou com um problema Gostaria de mudar a forma como eu recebo as postagem na lista. Geralmente recebo vários e-mails de todos os que postam. Recentemente um amigo fez o cadastro na lista e ele geralmente recebe um único e-mail com a listas das mensagens de um dado período. Gostaria de receber desta forma tambem. Pois evitaria a sobrecarca de minha caixa de e-mail. Algúem sabe como se faz isto? att -- Fernando Antônio de Souza Zootecnista, Msc. Nutricão Animal, UFMG lattes: http://lattes.cnpq.br/6519538815038307

3 3

Problema com strsplit()
by Augusto Ribas 19 Dez '12

19 Dez '12

Ola. Eu estou usando o pacote taxize, e estava tentando usar a função iucn_summary() deste pacote, mas ela funciona perfeitamente com os exemplos, mas quando eu vou usar um nome que digito de uma lista minha eu recebo a seguinte mensagem: ##### > iucn_summary("Molossus currentium") Erro em strsplit(pop, split = ":") : argumento modo não caractere ###### Eu olhei no site que a função busca e a espécie esta la. Não entendo o que acontece que ela não funcionou no meu caso, mas funcionou nos exemplos. Os exemplos que digo são esses: # library(taxize) iucn_summary("Lynx lynx") iucn_summary("aaaa") iucn_summary("Molossus currentium") # Eu não tenho ideia do que estou fazendo errado, mesmo quando a espécie não é encontrada no site, deveria ser retornado algo, como no exemplo do "aaaa". Eu estou usando Ubunto aqui e usando o R 2.15.2 do Emacs. -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.wordpress.com/ Lattes: http://lattes.cnpq.br/7355685961127056

2 1

Expression
by Fabio Mathias Corrêa 19 Dez '12

19 Dez '12

Pessoal, Ao utilizar o expression em um plot, o resultado não sai como o esperado. O mesmo argumento é usado no xlab e ylab, mas só funciona corretamente no xlab. CMR plot(1:10, xlab=expression(gamma[2]^"**"), ylab=expression(gamma[2]^"**")) Muito obrigado! Fábio Mathias Corrêa Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e da Terra - DCET Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna CEP 45662-900. Ilhéus-Bahia Tel.: 73-3680-5076

2 1

Chamada de Artigo Revista Mediações
by Ronaldo Baltar 18 Dez '12

18 Dez '12

Gostaríamos de convidar os interessados e também pedir a divulgação para Chamada de artigos – Revista Mediações, com o próximo dossiê sobre "Análise quantitativa e indicadores sociais". A Comissão Editorial de Mediações, Revista de Ciências Sociais, do Departamento e Programa de Pós-Graduação em Ciências Sociais, da Universidade Estadual de Londrina, receberá artigos para o volume 18, número 1, correspondente ao primeiro semestre de 2013. Esse Dossiê se propõe a incentivar a discussão sobre: aplicação de modelos quantitativos e indicadores sociais para análise das diferenciações sociais; procedimentos estatísticos para análise multivariada; agentes e simulação social; análise de redes sociais; entre outros. Serão especialmente bem-vindos artigos que usem R como código fonte para replicação científica. Um dos objetivos deste número é ajudar a disseminação do R como ferramenta para a pesquisa em ciências sociais. Organizadores: Prof. Ronaldo Baltar e Profª Claudia Siqueira Baltar Infosoc - Laboratório de Sociologia Computacional Departamento de Ciências Sociais Universidade Estadual de Londrina Recebimento de contribuições até 28 de fevereiro de 2013. Além do dossiê, Mediações publica trabalhos sobre temas diversos das Ciências Sociais na Seção Artigos, como também Resenhas de livros relevantes em Ciências Sociais. As submissões em português podem ser enviadas através do portal da Revista http://www.uel.br/revistas/uel/index.php/mediacoes ou pelo e-mail mediacoes(a)uel.br MEDIAÇÕES – REVISTA DE CIÊNCIAS SOCIAIS Departamento de Ciências Sociais/Centro de Letras e Ciências Humanas Universidade Estadual de Londrina Campus Universitário - Caixa Postal 6001 – Londrina/PR – 86055-900 Brasil

1 0

glm - variável factor
by Vinicius Brito Rocha 18 Dez '12

18 Dez '12

Senhores, tenho um modelo glm, logistico, onde um dos coef é uma variável factor com as Unidades Federativas. modelo=glm(y ~ UF+VA1+VA2-1) summary(modelo) O problema é que ao olhar o summary dos resultados do ajustes algumas das UF's não tem significância no teste para o parâmetro. OU seja tenho cenário com RJ significativo e SP não. Como faço para excluir a classe SP do meu novo modelo? Abs Vinicius -- *Vinicius Brito Rocha.* *Estatístico e Atuário (IM / UFRJ)** Mestre em Pesquisa Operacional (COPPE / UFRJ)* www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

3 4

Dúvida_experimentação_estimação
by andrebvs＠bol.com.br 18 Dez '12

18 Dez '12

1 0

Trabalhar com pares do RecordLinkage
by Fátima Lima Paula 17 Dez '12

17 Dez '12

Prezados, preciso trabalhar com os pares resultantes do RecordLinkage e não sei como fazer. Enviei email para o autor do pacote, mas ele disse que terei que programar, pois não existe uma maneira automática de fazer isso. Bem, aí vai, meu desafio: Tenho dois bancos ( um retirado do outro) e vou fazer o linkage para obter as readmissões. No resultado sempre haverá um par de cada um, com mesma data de internação, já que o banco menor é retirado do maior (o menor é de internados por fratura de fêmur e o maior é de internados por todas as causas). Aí vai um exemplo criado para você verem como é a saída dos pares: nome=c("Ricardo","Maria","Tereza","Pedro","José", "Germano") sobrenome=c("Cunha","Andrade","Silva","Soares","Silva","Lima") idade=c(67,89,78,65,68,67) ano=c(1945,1923,1934,1947,1944,1945) dataint=c(20120907,20120703,20120301,20120805,20121004,20121209) banco1=cbind(nome,sobrenome,idade,ano,dataint) banco1=data.frame(banco1) banco1 nome=c("Maria","Lúcia","Paulo","Marcos", "Ricardo", "Germanio") sobrenome=c("Andrade","Silva","Soares","Pereira","Cunha","Lima") idade=c(67,88,78,60,68,80) ano=c(1945,1924,1934,1952,1944,1932) dataint=c(20121208,20121103,20120302,20120105,20121004,20121209) banco2=cbind(nome,sobrenome,idade,ano,dataint) banco2=data.frame(banco2) banco2 str(banco1) banco1$nome=as.character(banco1$nome) banco1$sobrenome=as.character(banco1$sobrenome) banco1$idade=as.character(banco1$idade) banco1$ano=as.character(banco1$ano) banco1$dataint=as.character(banco1$dataint) banco2$nome=as.character(banco2$nome) banco2$sobrenome=as.character(banco2$sobrenome) banco2$idade=as.character(banco2$idade) banco2$ano=as.character(banco2$ano) banco2$dataint=as.character(banco2$dataint) library(RecordLinkage) soundexBR<-function(termo){ termo<-toupper(termo) # 1. Retire toda pontuação da palavra; termo<-gsub("[ÃÁÀÂÄ]","A",termo) termo<-gsub("[ÉÈÊ?Ë]","E",termo) termo<-gsub("[ÍÌÏÎI]","I",termo) termo<-gsub("[ÓÕÒÔÖ]","O",termo) termo<-gsub("[ÚÙÛÜU]","U",termo) termo<-gsub("Ç","C",termo) #1.a Substituir todas as letras duplas por uma única letra #http://www.archives.gov/genealogy/census/soundex.html termo<-gsub("([A-Z])\\1", "\\1", termo) # 2. Fique com a primeira letra; N<-nchar(termo) # Modificação PT-BR (Coeli e Camargo) termo<-ifelse(substr(termo,1,2)=="WA",sub("W","V",termo),termo) termo<-ifelse(substr(termo,1,1)=="H",substr(termo,2,N),termo) termo<-ifelse(substr(termo,1,2)=="KA"|substr(termo,1,2)=="KO"|substr(termo,1,2)=="KU",sub("K","C",termo),termo) termo<-ifelse(substr(termo,1,1)=="Y",sub("Y","I",termo),termo) termo<-ifelse(substr(termo,1,2)=="CE"|substr(termo,1,2)=="CI",sub("C","S",termo),termo) termo<-ifelse(substr(termo,1,2)=="GE"|substr(termo,1,2)=="GI",sub("G","J",termo),termo) ############## Segunda Parte ####################### termo.1<-substr(termo,1,1) termo<-substr(termo,2,N) # 3. Mude todas as ocorrências das letras a seguir por '0' (zero): 'A','E','I','O','U','H','W','Y'; termo<-gsub("[A,E,I,O,U,H,W,Y]",0,termo) # 4. Mude as letras restantes de acordo com a tabela abaixo: # Número Letra # 1 'B','F','P','V' termo<-gsub("[B,F,P,V]",1,termo) # 2 'C','G','J','K', 'Q','S','X','Z' termo<-gsub("[C,G,J,K,Q,S,X,Z]",2,termo) # 3 'D','T' termo<-gsub("[D,T]",3,termo) # 4 'L' termo<-gsub("L",4,termo) # 5 'M','N' termo<-gsub("[M,N]",5,termo) # 6 'R' termo<-gsub("R",6,termo) # 5. Remova todos os zeros da string resultante; termo<-gsub(0,"",termo) # Eliminar dois numeros iguais consecutivos por um único numero #http://www.archives.gov/genealogy/census/soundex.html termo<-gsub("([0-9])\\1", "\\1", termo) # Remontar termo<-paste(termo.1,termo,sep="") # 6. Preencha a string resultante com zeros à direita e retorne desta forma: #<letra maiuscula><digito><digito><digito> termo<-paste(termo,"0000",sep="") termo<-substr(termo,1,4) return(termo) } soundexBR("termo") soundexBR("terno") pares=compare.linkage(banco1,banco2,blockfld=list(c(1,2,4),c(1,2)), phonetic=c(1,2),phonfun=soundexBR,strcmp=FALSE, strcmpfun=jarowinkler, exclude=FALSE,identity1=NA,identity2=NA, n_match=NA, n_non_match=NA) editMatch(pares) pares pesos=epiWeights(pares, e=0.01,f=pares$frequencies) hist(pesos$Wdata,plot=F) getPairs(pesos,max.weight=Inf, min.weight=-Inf) Agora, como tenho as datas da internação de todos, preciso verificar, par a par, quem tem data de internação diferente. Estes serão as readmissões. Existe uma forma de eu criar uma coluna para identificar os pares? E depois, testar se as datas de internação são diferentes? Obrigada Fátima

2 1

ifelse composto
by Edson Lira 17 Dez '12

17 Dez '12

Bom dia caros r-listeiros, estou com um probleminha. Tenho um banco com a variável Naturalidade que tem as seguintes respostas: MANAUS, BAHIA, RIO DE JANEIRO, ANAMÃ, BORBA,SALVADOR,CAREIRO... Quero criar uma variável chamda de prodedência na qual terei os seguintes níveis: MANAUS=CAPITAL ANAMÃ,BORBA,CAREIRO=INTERIOR BAHIA,RIO DE JANEIRO,SALVADOR=OUTRO ESTADO. COMO FAZER ISSO COM O ifelse? [ ]'s. Edson Lira Estatístico Manaus-Amazonas

4 5

média acumulada
by Leonard de Assis 17 Dez '12

17 Dez '12

Existe pronto um comando estilo cumsum que retorne ao invés de soma, a média? exemplo: set.seed(1) op <- par(mfrow = c(2, 2)) x <-rnorm(500) y <- cumsum(x) ts.plot(x) ts.plot(y) # para eu fazer o que quero, tenho que fazer: y.cmean <-cumsum(x)/(1:length(x)) ts.plot(y.cmean) par(op) existe algum jeito mais fácil? -- []s Leonard de Assis http://about.me/ldeassis

2 1

Dúvida_experimental_estimação
by andrebvs 16 Dez '12

16 Dez '12

1 0

R-br Dúvida - Windows 7 não consigo ler banco em excel
by Marcela Arruda 15 Dez '12

15 Dez '12

Boa tarde, estou com dificuldade de ler a base em excel no windows 7. A versão do meu R 2.15.2 e o Tinn-R 2.3.7.3. Ele funcionava muito bem antes de mudar a versão. Isso já ocorreu com vocês? Obrigada Marcela SVS/MS <http://about.me/ldeassis>

6 5

SoundexBR
by Fátima Lima Paula 15 Dez '12

15 Dez '12

Prezados, recebi do Daniel Marcelino o soundex para ser utilizado no RecordLinkage. Gostaria de saber se esse soundex é validado e pode ser utilizado sem problemas. Obrigada Fátima

3 4

Mallows - estimador robusto
by Natália Bordin Barbieri 15 Dez '12

15 Dez '12

Olá, preciso ajustar o seguinte modelo de estimação robusta, através do método de "Mallows' para o modelo de regressão logística que segue: require(robust) DM.glmRob3<-glmRob(DM ~ INFLAMACAO + IDADE + SEXOM + RACACOR + HIPERT + IMC + RCQ + TRIGT, family=binomial, data=DM, weights=NULL, method ="mallows", model = TRUE, control = glmRob.control) Porém, ao rodar, o que retorna é o seguinte erro: 'The covariance matrix has become singular during the iterations of the MCD algorithm. There are 284 observations (in the entire dataset of 400 obs.) lying on the hyperplane with equation a_1*(x_i1 - m_1) + ... + a_p*(x_ip - m_p) = 0 with (m_1,...,m_p) the mean of these observations and coefficients a_i from the vector a <- c(1, 0, 0, 0, 0, 0, 0, 0) Erro em solve.default(V) rotina Lapack dgesv: sistema é exatamente singular' Alguém já trabalhou com isso?Sabe me ajudar? Obrigada, Natália

2 7

Analise dados proporção e contagem
by José Francisco Arruda e Silva 14 Dez '12

14 Dez '12

Ola pessoal.Alguem poderia me dar uma dicas de como posso analisar dados de proporcao e contagem, uma vez que necessito saber se existe diferenca entre os tratamentos (dietas) no que se refere ao numero de individuos que nasceram e da razao sexual. segue abaixo uma pequena amostra. dados <- expand.grid(rept=1:15, trat=LETTERS[1:4]) dados$razao_sexual <- c(1,0.5,0.8,1,0.5,0.8,1,1,0.8,0.6,0.8,1,0.5,1,0.8,1,0.4,1,1,0.4,0.8,0.8,1,1,0.6,0.6,0.8,0.8,0.8,1,0.8,1,1,0.6,1,1,0.8,1,0.8,1,0.8,1,0.8, 1,0.8,1,0.6,0.7,0.7,1,1,1,0.8,1,1,1,1,1,0.8,1) dados$nascidos<-c(9,10,9,11,13,8,10,11,10,9,10,10,9,10,9,10,14,10,9,10,10,7,10,9,10,8,10,10,12,10,14,10,9,10,10,11,10,10,9,10,10,12,11,10,15,10,10,10,11, 9,10,10,9,10,9,9,10,10,10,9) Obrigado. Chicao.

3 4

AJUSTE DOS BETAS DE REGRESSÃO
by tiago souza marçal 14 Dez '12

14 Dez '12

Caros colegas, não estou conseguindo estimar os betas dos modelos quadrático e cúbico da regressão através da função reg.poly() do pacote ExpDes. y<-c(13.337250,24.96667,22.34600,24.31983) x<-c(1000000, 10000000, 50000000, 100000000) Glresíduo<- 72 Sqresíduo<- 1159.9 Gltratamento<- 3 Sqtratamento<- 167.6 reg.poly(y,x,72,1159.9,3,167.6) Agradeço pela atenção Att. Tiago.

2 1

Pacote RNetLogo
by Felipe Buchbinder 14 Dez '12

14 Dez '12

Olá, Estou tentando utilizar o pacote RNetLogo. Eu consigo abrir o NetLogo a partir do R mas, quando tento qualquer outro comando, recebo a seguinte mensagem: Error in .jcall(nl.obj, "V", "loadModel", .jnew("java/lang/String", model.path)) : RcallMethod: invalid object parameter Alguém sabe o que está acontecendo? Obrigado, Felipe

1 0

Seleção de modelo misto
by Diego Pujoni 14 Dez '12

14 Dez '12

Olá pessoal, tenho um experimento de medidas repetidas que estou analisando através de um modelo misto. Nove unidades amostrais independentes (frascos com meio de cultura com algas) foram alocadas aleatoriamente em 3 grupos ("c","t1","t2"). Não há a necessidade de inclusão de efeito aleatório de intercepto, pois as nove unidades amostrais são homogêneas entre si (retiradas da mesma amostra mãe). A concentração de algas foi mensurada a cada dois dias durante 10 dias. O objetivo é testar diferenças entre os grupos (i.e. tratamentos). Estimei um modelo somente com o intercepto e com a interação Group*Day, para testar qual é o melhor library(nlme) library(lme4) Day = rep(c(0,2,4,6,8,10),each=9) Group = as.factor(rep(c("c","c","c","t1","t1","t1","t2","t2","t2"),6)) Individual = rep(1:9,6) X = c(0.71,0.72,0.71,0.72,0.72,0.72,0.70,0.69,0.70,0.72,0.72, 0.71,0.72,0.72,0.71,0.71,0.70,0.71,0.73,0.73,0.69,0.74, 0.69,0.73,0.67,0.71,0.69,0.71,0.71,0.72,0.70,0.71,0.70, 0.52,0.64,0.60,0.70,0.73,0.73,0.67,0.66,0.71,0.47,0.56, 0.54,0.65,0.73,0.73,0.67,0.71,0.58,0.44,0.52,0.58) xyplot(X~Day, groups=Group) LME = lme(X ~ 1, random = ~Day|Individual) Erro em lme.formula(X ~ 1, random = ~Day | Individual) : nlminb problem, convergence error code = 1 message = iteration limit reached without convergence (10) LME1 = lme(X ~ Group*Day, random = ~Day|Individual) Erro em lme.formula(X ~ Group * Day, random = ~Day | Individual) : nlminb problem, convergence error code = 1 message = iteration limit reached without convergence (10) LMER = lmer(X ~ 1 + (Day|Individual)) LMER1 = lmer(X ~ Group*Day + (Day|Individual)) AIC(LMER) [1] -179.0302 AIC(LMER1) [1] -151.1938 anova(LMER,LMER1) Data: Models: LMER: X ~ 1 + (Day | Individual) LMER1: X ~ Group * Day + (Day | Individual) Df AIC BIC logLik Chisq Chi Df Pr(>Chisq) LMER 5 -187.2 -177.26 98.602 LMER1 10 -203.2 -183.31 111.600 25.996 5 8.939e-05 *** xyplot(fitted(LMER)~Day, groups=Group) xyplot(fitted(LMER1)~Day, groups=Group) 1a pergunta: Por que a função nlme:lme não converge, mas a lme4:lmer converge? A segunda é melhor que a primeira? 2a pergunta: Por que a função "anova" dá valores distintos de AIC para os dois modelos. Se olharmos o valor do AIC de cada um, o modelo LMER é melhor que o LMER1, mas a "anova" diz o contrário. 3a pergunta: Por que os valores ajustados do modelo somente com o intercepto (LMER) variam ao longo do tempo? Muito obrigado pela ajuda Diego PJ

1 0

mínimos quadrados bivariado
by Luiz Henrique da Conceicao Leal 14 Dez '12

14 Dez '12

Prezados Existe algum pacote no R que realize análise de regressão através dos mínimos quadrados bivariado - Bivariate Least Squares (BLS)?

3 3

[Dúvida] Tratamento de erro.
by Pedro Rafael 13 Dez '12

13 Dez '12

Pessoal como realizar tratamento de erros em R? Estou com uma função em que é executado um loop. Neste laço é executado o comando arima que em determinada interação poderá não convergir dando assim uma mensagem de erro que ocasionará o termino do laço. Gostaria de ignorar esses erros. O código segue: nrep = 50 matriz_result = matrix(,nrow=nrep,ncol=7) for(i in 1:nrep){ par = trunc(runif(5,0,4)) resultado = arima(retirar(serie,3)$Serie,order=c(par[1],2,par[2]), seasonal = list(order = c(par[3], par[4], par[5]))) aic = resultado$aic matriz_result[i,] = c(par[1],1,par[2:5],aic) } Gostaria de que no caso de erro a minha do objeto matriz_result seja completada por NA. Grato [ ], Pedro Rafael Diniz Marinho.

2 2

Qual pacote que tem o comando scripchart
by jairan dionizio 13 Dez '12

13 Dez '12

Bom dia todos e todas! gostaria de uma ajudinha. Qual pacote do contém comando scripchart? Jairan Dionizio Skype:jairan.dionizio Facebook: Jairan Dionizio Twitter: @jairandionizio Tel.:(84) 9423-5286 "A força não provém da capacidade física e sim de uma vontade indomável." (Mahatma Gandhi) Cuidados ao enviar e-mail Ao reenviar apague o meu e-mail e o meu nome, apague também os endereços dos amigos antes de reenviar qualquer e-mail, encaminhe como cópia oculta (Cco ou Bcc) aos SEUS destinatários.Agindo sempre assim dificultaremos a disseminação de vírus, spams, banners e outros.

3 2

Problemas Git-Hub
by Leandro Marino 13 Dez '12

13 Dez '12

Caros, alguém saberia me dizer o motivo do código abaixo não estar funcionando?! classcol = c(rep('character',5),rep('numeric',2),rep('character',2)) url <- 'https://raw.github.com/leandromarino/TBAmoMunic/master/Munic.dat' munic <- read.table(url,colClasses=classcol,sep='\t',header=T) Obrigado. Leandro

5 11

Manipulação de String
by Eder David Borges da Silva 13 Dez '12

13 Dez '12

Pessoal, peço ajuda aos mestres em manipulação de String... eu ainda não me acertei com isso. A duvida é simples: x <- c("Eder David Borges da Silva","João Silva","Pedro Augusto Paulo") gostaria de transformar estes nomes em formato "Americano" ficando: Silva, Eder David Borges da Silva, João Paulo, Pedro Augusto. Em síntese pegar oque esta depois do ultimo espaço, e colar com tudo que esta antes dele, separando isso por uma virgula. Já tentei algumas coisa, mas sempre me confundo nas expressões regulares. Alguma dica? Obrigado

6 7

Re: [R-br] Qual pacote que tem o comando scripchart
by Walmes Zeviani 13 Dez '12

13 Dez '12

Infelizmente não sou capaz de te ajudar porque não entendi o seu problema. Forneça código reproduzível que sem dúvida é esclarecedor. Envie suas mensagens para lista de discussão pois mais pessoas podem colaborar, obter soluções diferentes, ter resposta em menos tempo, tópico ser útil para outras pessoas, ter as mensagens arquivadas para referência futura. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

1 0

Scatter plot de matriz
by Luciano F. Sgarbi 13 Dez '12

13 Dez '12

Prezados, Estou tentando inserir no "scatter plot de matriz" (veja o conjunto dados de exemplo, abaixo) as linhas de correlação entre as variáveis, entretanto não estou conseguindo de um modo mais direto. Seria para como usar para cada um dos gráficos o seguinte exemplo: #------- z <- lm(dist ~ speed, data = cars) plot(cars) abline(z) # equivalent to abline(reg = z) or abline(coef = coef(z)) #--------- Entretanto eu preciso colocar uma reta para cada plot abaixo #DADOS DE EXEMPLO var_1<-c(0.492077114759013 , 0.836016362765804 , 0.481298430124298 , 0.386469675228 , 0.777098514838144 , 0.617958343820646 , 0.219213913194835 , 0.678557483945042 , 0.423462828854099 , 0.829378235153854 , 0.806068887002766 , 0.337134317029268 , 0.661742311669514 , 0.619118098868057 , 0.788589890813455 , 0.279048688709736 , 0.820936637930572 , 0.247459845617414 , 0.752754796063527 , 0.934190914966166) var_2<-c(0.658414233683847 , 0.605210445823305 , 0.801341124716519 , 1.011018222403 , 0.572025150910094 , 0.836177895660953 , 1.03193629108894 , 0.656809315729027 , 0.949461899587905 , 0.553121786467357 , 0.817865236941146 , 0.925669452041572 , 0.578475899228501 , 0.476615221718288 , 0.571342170684858 , 0.867686931852481 , 0.615533868096807 , 0.926396173911956 , 0.667834629397337 , 0.65446214987638) var_3<-c(0.805190744576976 , 0.00357919884845614 , 0.296429808018729 , 0.147386643337086 , 0.224931408418342 , 0.967352513456717 , 0.88843358354643 , 0.991240894189104 , 0.906759183388203 , 0.635862796101719 , 0.90428285067901 , 0.755293083842844 , 0.484490071190521 , 0.0830997943412513 , 0.66727612237446 , 0.314308791654184 , 0.478725346503779 , 0.885625332361087 , 0.567443900275975 , 0.175529750995338) var_4<-c(1.51875793231285 , 2.08332774568708 , 1.85743044603302 , 1.94252530517918 , 1.70102115019736 , 1.4099218798547 , 1.37410255920602 , 1.63332942855045 , 1.59918963127413 , 1.54217901497363 , 1.39090041665641 , 1.56455567202422 , 1.70518441073333 , 1.96772731408659 , 1.76799247833419 , 1.85562083102216 , 1.81231983166867 , 1.58934006350458 , 1.77421167457209 , 2.01928129901513) data<-data.frame(var_1,var_2,var_3,var_4) plot(data,pch=16) Certo de vossa atenção, Atenciosamente, Luciano -- Luciano F. Sgarbi Mestrando em Ecologia e Evolução - UFG Laboratório de Ecologia de Insetos Cel. (62)8174-2262 Lab. (62)3521-1732

2 2

Guardar objetos dentro de uma função
by Leandro Marino 13 Dez '12

13 Dez '12

Caros, boa tarde. Gostaria de tirar um dúvida. Estou preparando uma função básica e ela tem 3 objetos distintos de saída sendo um deles o resumo para o usuário e outros a informação linha por linha do arquivo de dados de entrada. O que gostaria de fazer é colocar estes dois últimos objetos escondidos na função de tal forma que quando o usuário rodasse só seria apresentado o resumo. Seria algo equivalente ao que o lm() faz. Se rodo o lm() só aparece o resumo, mas posso extrair os resíduos com o comando lm()$residuals . y <- rnorm(100) x <- rbinom(100,10,prob=.2) obj <- lm(y~x) names(obj) []s, Leandro

4 3

respostas
by Gilbert Queiroz 13 Dez '12

13 Dez '12

Só queria saber pq minhas respostas nunca vão p/ lista...

3 2

Dúvida strsplit
by Leandro Marino 13 Dez '12

13 Dez '12

Caros colegas, estou com a seguinte situação: vec <- c("001011001110010110000011010100", "101111111011100111110101111101", "101111111111111110111111111101", "111001100010111000100000101010", "100011110011000000100000000000", "110111111010110001111111011111", "001000010100010010001010001111", "000011000000000011001001110101", "011111110101111101111111111111", "110010111111101000010001010000") temp <- strsplit(vec,NULL) temp <- do.call(rbind,temp) no entanto, gostaria de transformar as linhas e colunas do temp para integer. Como fazer isto de forma eficiente? poderia transformar no momento do strsplit?! Abraços,

3 5

Lógica Fuzzy
by marcelo claro de souza 13 Dez '12

13 Dez '12

Olá pessoal, Vocês poderiam me sugerir algum material básico sobre construção de modelos utilizando lógica Fuzzy no R ou um exemplo comentado passo a passo. Muito obrigado. Marcelo

4 4

Derivar função não linear
by Fernando Antonio de souza 12 Dez '12

12 Dez '12

Caros amigos Estou precisando de uma ajuda que demanda conhecimentos matemáticos maiores que possuo. Procurei alguns professores aqui no departamento mas ninguem sabe como fazer, por isto estou novamente recorrendo a comunidade na esperança de encontrar uma solução. Estou necessitando fazer derivadas de algumas equações, conforme uma metodologia que estou seguindo. Mas não lembro mais como se faz isto. Da última vez o Tiago Marçal me ajudou, mas acontece que surgiu outra. possuo o seguinte modelo: log10(Y) = a-b*exp(-c*t) preciso derivá-la em função do tempo(t). onde Y é o conteúdo de energia em dado tempo (t) Na metodologia que estou seguindo após fazer esta derivada chegou-se ao seguinte modelo d(E)/d(t) = E(t)*0,0737*exp(-0,099643*t) onde E(t) é o conteúdo de energia no tempo(t) Necessito saber como ele fez para chegar nisto?E se possível como se chega! Desde já agradeçoa a atenção

2 10

instalação R studio
by Natalia Martins 12 Dez '12

12 Dez '12

Prezados bom dia, estou tentando instalar o R studio mas no momento da instalaçao ha o pedido de uma senha, algume sabe o porque isso ocorre? Att -- Natália da Silva Martins Bacharel em Estatística - Universidade Estadual de Maringá/ UEM Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 8306-4743

5 5

ajuda ifelse
by Sérgio Henrique almeida da silva ju 12 Dez '12

12 Dez '12

Amigos estou trabalhando com uma lista de doenças por municípios e vou plotar as mesmas em um mapa. Tenho por exemplo Doenças: a,b,c,d,e,d,f, categorizadas com 1.5 (o município possui a doença) e 0.5 (não possui), gostaria de categorizar os municípios, segundo a quantidade de doenças e para isso fiz: Municípios com todas as doenças: banco$todas <- ifelse(banco$a=="1.5" & banco$b=="1.5" & banco$c==1.5 & banco$d=="1.5" & banco$e=="1.5" & banco$f=="1.5"),1.5, 0.5) Municípios com pelo menos uma doença: banco$pelomenosuma <- ifelse(banco$a=="1.5" | banco$b=="1.5" | banco$c==1.5 | banco$d=="1.5" | banco$e=="1.5" | banco$f=="1.5"),1.5, 0.5) a partir dai tenho dificuldade de categorizar os municípios com por exemplo: 2 ou mais doenças, três ou mais doenças e assim por diante, eu tentei algumas coisas, porém não consegui. Como posso faze isso? 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) 68463637 / 94429486

5 8