- R-br - DINF UFPR

Tabelas em sql para verficar.
by Edson Lira 25 Jul '12

25 Jul '12

Benilton e Leonardo, só agora pude voltar aquele problema do sql. Abaixo o script usado para ler o sql. rm(list = ls()) require(RODBC) options(digits=3) options(max.print=800000) getOption("max.print") channel<-odbcConnect("saddw",uid="saddw",case="tolower") #DOAÇÃO GERAL doacao.dados<- sqlQuery(channel, "select doacao_dw.numerobolsa, doacao_dw.datahoradoacao,doacao_dw.tipodoacao, doacao_dw.tipodoacaodirigida,doacao_dw.localcoleta,doacao_dw.apto, doacao_dw.codigoinaptidao,doador_dw.tipagemsanguinea,doador_dw.rh, doador_dw.idade, doador_dw.sexo, doador_dw.data_primeira_doacao, doador_dw.data_ultima_doacao,bairroresidencia, bairro_nao_tabelado,doador_dw.qtde_doacoes, doador_dw.numerogente from saddw.doador_dw, saddw.doacao_dw where doador_dw.numerogente = doacao_dw.numerogente and to_char(doacao_dw.datahoradoacao,'DD/MM/YYYY') >= to_date('27/08/1992','DD/MM/YYYY') and to_char(doacao_dw.datahoradoacao,'DD/MM/YYYY') <= to_date('28/08/1992','DD/MM/YYYY')" ) No link abaixo o script que gera os dados em sql. http://www.datafilehost.com/download-054b9e31.html Abaixo erro ao rodar no R. + and to_char(doacao_dw.datahoradoacao,'DD/MM/YYYY') >= to_date('28/08/1992','DD/MM/YYYY') + and to_char(doacao_dw.datahoradoacao,'DD/MM/YYYY') <= to_date('28/08/1992','DD/MM/YYYY')" + ) Erro em as.POSIXlt.character(x, tz, ...) : character string is not in a standard unambiguous format [ ]'s. Edson Lira Estatístico Manaus-Amazonas

2 2

Utilização dos pacotes snow e rmpi em cluster
by Junior Beleti 24 Jul '12

24 Jul '12

Boa noite.. Gostaria de saber se alguém possui algum material sobre utilização dos pacotes snow e rmpi em cluster, além é claro dos próprios manuais dos pacotes. Att, Junior.

1 0

Geoprocessamento
by Edmar Caldas 24 Jul '12

24 Jul '12

Olá pessoal, gostaria de saber se alguém tem algum material sobre geo pode ser o básico, grato. estou querendo aprender. Edmar

3 2

Usando o R no Ubuntu
by Vinícius Lionel Mateus 24 Jul '12

24 Jul '12

Prezados, Alguém poderia me indicar algum material para ajudar no uso do R na plataforma Ubuntu 12.04: já instalei o software, mas gostaria de tentar usar o GUI também, assim como no Windows. -- -- Atenciosamente, VINÍCIUS LIONEL MATEUS, B.Sc (http://lattes.cnpq.br/6501001637020665) Bacharel em Química - Doutorando em Química Analítica Laboratório de Química Atmosférica - Departamento de Química Pontifícia Universidade Católica - Rio de Janeiro (PUC - Rio) Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - Brasil CEP.: 22451-900 Telefone: (+55) (21) 3527-1327 (+55) (21) 9358-8051 www.puc-rio.br

7 8

Re: [R-br] Sugestões de interfaces gráficas para o R no UBUNTU
by Gilenio Borges Fernandes 24 Jul '12

24 Jul '12

Eu uso e gosto o Rcomander (Rcmdr) e o RKWARD. A instalação é direta do gerenciador de pacotes. Os icones ficam disponiveis no grupo "Ciencia". Boa sorte Gilenio Fernandes Em 24 de julho de 2012 11:04, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: > 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: Extrair palavras com pontuações! (Ivan Bezerra Allaman) > 2. Variogram model 3D (Tiago Vieira) > 3. comparações multiplas utilizando Kruskall- walis > (Marcelo Claro de Souza) > 4. Re: comparações multiplas utilizando Kruskall- walis > (Walmes Zeviani) > 5. Re: comparações multiplas utilizando Kruskall- walis > (Ivan Bezerra Allaman) > 6. Re: Variogram model 3D (Luis Iván Ortiz Valencia) > 7. Re: Variogram model 3D (Tiago Vieira) > 8. duvida #Include (nani.design(a)bol.com.br) > 9. Re: duvida #Include (Edson Lira) > 10. Re: duvida #Include (nani.design(a)bol.com.br) > 11. Re: duvida #Include (nani.design(a)bol.com.br) > 12. Re: duvida #Include (Daniel C Bezerra) > 13. Usando o R no Ubuntu (Vinícius Lionel Mateus) > 14. Re: Usando o R no Ubuntu (Daniel C Bezerra) > 15. Re: Usando o R no Ubuntu (Pedro Rafael) > 16. Re: Usando o R no Ubuntu (Fabio Mathias Corrêa) > 17. Re: duvida #Include (nani.design(a)bol.com.br) > 18. Re: Usando o R no Ubuntu (Leonard de Assis) > 19. Re: Usando o R no Ubuntu (FHRB Toledo) > 20. Re: Usando o R no Ubuntu (Zhu Sha Zang) > 21. Re: Usando o R no Ubuntu (Vinícius Lionel Mateus) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Mon, 23 Jul 2012 08:03:28 -0700 (PDT) > From: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> > To: Jakson Alves de Aquino <jalvesaq(a)gmail.com> > Cc: R Brasil <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Extrair palavras com pontuações! > Message-ID: > <1343055808.55795.YahooMailNeo(a)web161805.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > > > Valeu!! > > > \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} > > > ________________________________ > De: Jakson Alves de Aquino <jalvesaq(a)gmail.com> > Para: r-br(a)listas.c3sl.ufpr.br; Ivan Bezerra Allaman < > ivanalaman(a)yahoo.com.br> > Enviadas: Segunda-feira, 23 de Julho de 2012 11:50 > Assunto: Re: [R-br] Extrair palavras com pontuações! > > 2012/7/23 Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br>: > > Suponhamos que eu tenha o seguinte vetor: > > > > nfa <- c('A','JA','K','K:P','P:K','K:P:N') > > [...] > > > No entanto, como eu faço para extrair as palavras que tem as duas > > pontuações, ou seja, no caso acima 'K:P:N'? > > nfa[grep(':.:', nfa)] > > ou > > nfa[grep('[[:punct:]].[[:punct:]]',nfa)] > > -- > Jakson > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120723/d8204928/atta… > > > > ------------------------------ > > Message: 2 > Date: Mon, 23 Jul 2012 10:04:14 -0700 (PDT) > From: Tiago Vieira <tiagovieira7(a)yahoo.com.br> > To: R-br <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Variogram model 3D > Message-ID: > <1343063054.7351.YahooMailNeo(a)web122106.mail.ne1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Olá a todos, > > Minha dúvida desta vez é: Como posso alterar o Alcance Vertical na função > vgm do pacote gstat? Pois estou realizando uma krigagem 3D e consigo fazer > isso no gslib, mas ainda não consegui no gstat! > > Grato, > > Tiago Vieira. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120723/5cb66ac1/atta… > > > > ------------------------------ > > Message: 3 > Date: Mon, 23 Jul 2012 10:14:20 -0700 (PDT) > From: Marcelo Claro de Souza <marcelo_claro(a)yahoo.com.br> > To: Rbr <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] comparações multiplas utilizando Kruskall- walis > Message-ID: > <1343063660.70314.YahooMailNeo(a)web160103.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Não consigo instalar o pacote pgirmess para rodar comparações multiplas > pelo Kruskall-walis pois é necessário instalar tb instalar o pacote gdal e > o mesmo não está disponível para a versão 2.14.1 do linux. > Já tentei instalar esse pacote pelo terminal usando comando sudo apt-get > build-dep r-cran-rgdal e apareceu a mesma mensagem de erro. > Existe algum outro pacote que eu possa utilizar no lugar pgirmess que faça > a mesma coisa que a função kruskalmc? > Muito obrigado pela ajuda. > > > Marcelo Claro de Souza > > Biologist, PhD student in Plant Biology > Institute of Bioscience - UNESP, Brazil > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120723/4f19b888/atta… > > > > ------------------------------ > > Message: 4 > Date: Mon, 23 Jul 2012 14:18:22 -0300 > From: Walmes Zeviani <walmeszeviani(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, Marcelo Claro de Souza > <marcelo_claro(a)yahoo.com.br> > Subject: Re: [R-br] comparações multiplas utilizando Kruskall- walis > Message-ID: > <CAFU=EkY=BpY89E-cjK_8tdSHZkPkhdao= > gymJ96a2Fr_C+sO0A(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Porque essas comparações múltiplas tem que ser por esse teste? Seus dados > não atendem suposições de modelos paramétricos? Pensou em bootstrap? > > À 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 > 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/20120723/aa61da6c/atta… > > > > ------------------------------ > > Message: 5 > Date: Mon, 23 Jul 2012 11:25:35 -0700 (PDT) > From: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> > To: R Brasil <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] comparações multiplas utilizando Kruskall- walis > Message-ID: > <1343067935.10095.YahooMailNeo(a)web161801.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Segundo informações no site da cran, o pacote pgirmess depende do pacote > rgdal, e este está disponível no cran. > > > \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/20120723/335722aa/atta… > > > > ------------------------------ > > Message: 6 > Date: Mon, 23 Jul 2012 17:00:38 -0300 > From: Luis Iván Ortiz Valencia <liov2067(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, Tiago Vieira <tiagovieira7(a)yahoo.com.br> > Subject: Re: [R-br] Variogram model 3D > Message-ID: > <CAJ+S-MFkjWWq63joXsZ74_qROeVfwp6cs=UyLuo-DFR= > ZpqhMg(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > take a look at vgm function..from gstat <) > > Em 23 de julho de 2012 14:04, Tiago Vieira <tiagovieira7(a)yahoo.com.br > >escreveu: > > > Olá a todos, > > > > Minha dúvida desta vez é: Como posso alterar o Alcance Vertical na função > > vgm do pacote gstat? Pois estou realizando uma krigagem 3D e consigo > fazer > > isso no gslib, mas ainda não consegui no gstat! > > > > Grato, > > > > Tiago Vieira. > > > > _______________________________________________ > > 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. > > > > > > -- > Luis Iván Ortiz Valencia > Doutorando Saúde Pública - Epidemiologia, IESC, UFRJ > Estatístico Msc. > Spatial Analyst Msc. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120723/053517fa/atta… > > > > ------------------------------ > > Message: 7 > Date: Mon, 23 Jul 2012 13:31:10 -0700 (PDT) > From: Tiago Vieira <tiagovieira7(a)yahoo.com.br> > To: Luis Iván Ortiz Valencia <liov2067(a)gmail.com>, > "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Variogram model 3D > Message-ID: > <1343075470.75580.YahooMailNeo(a)web122101.mail.ne1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Isto foi a primeira coisa que eu fiz, mas a descrição da função não fala > nada a respeito. > > > De: Luis Iván Ortiz Valencia <liov2067(a)gmail.com> > Para: r-br(a)listas.c3sl.ufpr.br; Tiago Vieira <tiagovieira7(a)yahoo.com.br> > Enviadas: Segunda-feira, 23 de Julho de 2012 17:00 > Assunto: Re: [R-br] Variogram model 3D > > > take a look at vgm function..from gstat <) > > > Em 23 de julho de 2012 14:04, Tiago Vieira <tiagovieira7(a)yahoo.com.br> > escreveu: > > Olá a todos, > > > >Minha dúvida desta vez é: Como posso alterar o Alcance Vertical na função > vgm do pacote gstat? Pois estou realizando uma krigagem 3D e consigo fazer > isso no gslib, mas ainda não consegui no gstat! > > > >Grato, > > > >Tiago Vieira. > >_______________________________________________ > >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. > > > > > -- > Luis Iván Ortiz Valencia > Doutorando Saúde Pública - Epidemiologia, IESC, UFRJ > Estatístico Msc. > Spatial Analyst Msc. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120723/90347f89/atta… > > > > ------------------------------ > > Message: 8 > Date: Mon, 23 Jul 2012 22:58:37 -0300 > From: nani.design(a)bol.com.br > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] duvida #Include > Message-ID: <500e014d8a336_12fde2e2c74455(a)a4-winter20.tmail> > Content-Type: text/plain; charset="iso-8859-1" > > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120723/22b2b66c/atta… > > > > ------------------------------ > > Message: 9 > Date: Mon, 23 Jul 2012 22:07:17 -0400 > From: Edson Lira <edinhoestat(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] duvida #Include > Message-ID: <CCD47E10-5663-420F-BAC2-9ABC812188C2(a)yahoo.com.br> > Content-Type: text/plain; charset=utf-8 > > Mande um CMR para que possa ser verificado. > > [. ]'s. > Edson Lira > Estatístico > Ma-Am > > Em 23/07/2012, às 21:58, nani.design(a)bol.com.br escreveu: > > > Boa noite, > > > > Estou com uma pequena dúvida em meu programa. > > Tenho 3 arquivos de funções em R e em um 4º arquivo eu utilizo funções > desses 3 arquivos anteriores. > > Estou fazendo um #Include pra cada arquivo, por exemplo: > > > > #Include(func1) > > #Include(func2) > > #Include(func3) > > > > Porém o problema é que se eu chamo funções do arquivo func1 e func2, > funciona bem. Mas se chamo a função do arquivo func3, dá erro falando que > não encontra o objeto, como se não tivesse inserido. > > > > Alguém tem alguma sugestão do que possa ser? > > Obs.: se rodo o arquivo func3 separadamente, funciona bem. > > > > Obrigada. > > Ariane > > > > _______________________________________________ > > 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. > > > ------------------------------ > > Message: 10 > Date: Mon, 23 Jul 2012 23:15:01 -0300 > From: nani.design(a)bol.com.br > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] duvida #Include > Message-ID: <500e05258ce34_2eba3f5447c1052(a)a4-winter13.tmail> > Content-Type: text/plain; charset="iso-8859-1" > > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120723/1f7ea60e/atta… > > > > ------------------------------ > > Message: 11 > Date: Mon, 23 Jul 2012 23:34:49 -0300 > From: nani.design(a)bol.com.br > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] duvida #Include > Message-ID: <500e09c9b9850_7a80c719c7c24d(a)a4-winter12.tmail> > Content-Type: text/plain; charset="iso-8859-1" > > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120723/8d03de8b/atta… > > > > ------------------------------ > > Message: 12 > Date: Tue, 24 Jul 2012 00:03:36 -0300 > From: Daniel C Bezerra <danielcbezerra(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] duvida #Include > Message-ID: > <CAJGh7uwSSmnfBEPw0dirU+qBB+qr=vJ2Y-avfK= > 1upf8FBBiuA(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > source("nome.do.arquivo") > > > 2012/7/23 <nani.design(a)bol.com.br> > > > Existe alguma outra função para incluir arquivos de funções externos? > > > > Abr, > > > > > > > > ------------------------------ > > Em 23/07/2012 23:07, *Edson Lira < edinhoestat(a)yahoo.com.br >* escreveu: > > Mande um CMR para que possa ser verificado. > > > > [. ]'s. > > Edson Lira > > Estatístico > > Ma-Am > > > > Em 23/07/2012, às 21:58, nani.design(a)bol.com.br escreveu: > > > > > Boa noite, > > > > > > Estou com uma pequena dúvida em meu programa. > > > Tenho 3 arquivos de funções em R e em um 4º arquivo eu utilizo funções > > desses 3 arquivos anteriores. > > > Estou fazendo um #Include pra cada arquivo, por exemplo: > > > > > > #Include(func1) > > > #Include(func2) > > > #Include(func3) > > > > > > Porém o problema é que se eu chamo funções do arquivo func1 e func2, > > funciona bem. Mas se chamo a função do arquivo func3, dá erro falando que > > não encontra o objeto, como se não tivesse inserido. > > > > > > Alguém tem alguma sugestão do que possa ser? > > > Obs.: se rodo o arquivo func3 separadament e, funciona bem. > > > > > > > > Obrigada. > > > Ariane > > > > > > _______________________________________________ > > > 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. > > > > > > _______________________________________________ > > 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/20120724/9d881c81/atta… > > > > ------------------------------ > > Message: 13 > Date: Tue, 24 Jul 2012 00:13:14 -0300 > From: Vinícius Lionel Mateus <vinynegrelli(a)gmail.com> > To: R-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Usando o R no Ubuntu > Message-ID: > <CAMJ4BpD__TXBM7s9ecZ65Jdfg1amL2yeAimi2JTY3= > 2KGTuUaw(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Prezados, > > Alguém poderia me indicar algum material para ajudar no uso do R na > plataforma Ubuntu 12.04: já instalei o software, mas gostaria de tentar > usar o GUI também, assim como no Windows. > > -- > -- > Atenciosamente, > > VINÍCIUS LIONEL MATEUS, B.Sc (http://lattes.cnpq.br/6501001637020665) > Bacharel em Química - Doutorando em Química Analítica > Laboratório de Química Atmosférica - Departamento de Química > Pontifícia Universidade Católica - Rio de Janeiro (PUC - Rio) > Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - Brasil CEP.: > 22451-900 > Telefone: (+55) (21) 3527-1327 > (+55) (21) 9358-8051 > www.puc-rio.br > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120724/8869846c/atta… > > > > ------------------------------ > > Message: 14 > Date: Tue, 24 Jul 2012 00:21:33 -0300 > From: Daniel C Bezerra <danielcbezerra(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Usando o R no Ubuntu > Message-ID: > < > CAJGh7ux-KYr3E06tBdJERMnVqBSzBngfm152VPck_DCwZKN5ew(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > O R no Ubuntu não tem um GUI como no Windows. > > Vc pode executá-lo a partir do terminal (basta escrever 'R' - sem aspas - > no terminal, apertar enter e pronto). > > Algumas alternativas de GUI estão sugeridas aqui: > > http://www.linuxlinks.com/article/20110306113701179/GUIsforR.html > > Uma boa sugestão é o R Studio que vc pode baixar aqui: > > http://rstudio.org/ > > Abs, > > D > > > 2012/7/24 Vinícius Lionel Mateus <vinynegrelli(a)gmail.com> > > > Prezados, > > > > Alguém poderia me indicar algum material para ajudar no uso do R na > > plataforma Ubuntu 12.04: já instalei o software, mas gostaria de tentar > > usar o GUI também, assim como no Windows. > > > > -- > > -- > > Atenciosamente, > > > > VINÍCIUS LIONEL MATEUS, B.Sc (http://lattes.cnpq.br/6501001637020665) > > Bacharel em Química - Doutorando em Química Analítica > > Laboratório de Química Atmosférica - Departamento de Química > > Pontifícia Universidade Católica - Rio de Janeiro (PUC - Rio) > > Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - Brasil > CEP.: > > 22451-900 > > Telefone: (+55) (21) 3527-1327 > > (+55) (21) 9358-8051 > > www.puc-rio.br > > > > _______________________________________________ > > 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/20120724/0bf5e0d3/atta… > > > > ------------------------------ > > Message: 15 > Date: Tue, 24 Jul 2012 01:01:36 -0300 > From: Pedro Rafael <pedro.rafael.marinho(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Usando o R no Ubuntu > Message-ID: > < > CAKwavrkompYOX6oC_rsva_Db0iLCtRXV5dE0+M5YP-6J1jh+5A(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Instale o Rstudio ou utilize algum editor como o Vim ou Emacs para escrever > seus códigos em R. > [ ], > Pedro Rafael Diniz Marinho. > > > Em 24 de julho de 2012 00:22, Daniel C Bezerra [via R-br] < > ml-node+s2285057n4655667h54(a)n4.nabble.com> escreveu: > > > O R no Ubuntu não tem um GUI como no Windows. > > > > Vc pode executá-lo a partir do terminal (basta escrever 'R' - sem aspas - > > no terminal, apertar enter e pronto). > > > > Algumas alternativas de GUI estão sugeridas aqui: > > > > http://www.linuxlinks.com/article/20110306113701179/GUIsforR.html > > > > Uma boa sugestão é o R Studio que vc pode baixar aqui: > > > > http://rstudio.org/ > > > > Abs, > > > > D > > > > > > 2012/7/24 Vinícius Lionel Mateus <[hidden email]< > http://user/SendEmail.jtp?type=node&node=4655667&i=0> > > > > > > >> Prezados, > >> > >> Alguém poderia me indicar algum material para ajudar no uso do R na > >> plataforma Ubuntu 12.04: já instalei o software, mas gostaria de tentar > >> usar o GUI também, assim como no Windows. > >> > >> > >> -- > >> -- > >> Atenciosamente, > >> > >> VINÍCIUS LIONEL MATEUS, B.Sc (http://lattes.cnpq.br/6501001637020665) > >> Bacharel em Química - Doutorando em Química Analítica > >> Laboratório de Química Atmosférica - Departamento de Química > >> Pontifícia Universidade Católica - Rio de Janeiro (PUC - Rio) > >> Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - > >> Brasil CEP.: 22451-900 > >> Telefone: <a href="tel:%28%2B55%29%20%2821%29%203527-1327" value=" > >> +552135271327" target="_blank">(+55) (21) 3527-1327 > >> <a href="tel:%28%2B55%29%20%2821%29%209358-8051" value=" > >> +552193588051" target="_blank">(+55) (21) 9358-8051 > >> www.puc-rio.br > >> > >> _______________________________________________ > >> R-br mailing list > >> [hidden email] <http://user/SendEmail.jtp?type=node&node=4655667&i=1> > >> 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 > > [hidden email] <http://user/SendEmail.jtp?type=node&node=4655667&i=2> > > 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. > > > > ------------------------------ > > If you reply to this email, your message will be added to the discussion > > below: > > > > > http://r-br.2285057.n4.nabble.com/R-br-Usando-o-R-no-Ubuntu-tp4655666p46556… > > To unsubscribe from R-br, click here< > http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib… > > > > . > > NAML< > http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_view… > > > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120724/73b55b72/atta… > > > > ------------------------------ > > Message: 16 > Date: Tue, 24 Jul 2012 04:11:37 -0700 (PDT) > From: Fabio Mathias Corrêa <fabio.ufla(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Usando o R no Ubuntu > Message-ID: > <1343128297.3402.YahooMailNeo(a)web120802.mail.ne1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Olá Vinícius, > > No Ubuntu vc tem diversas opções: Gvim+Vim-r-plugin + Emacs+ESS, > Gedit+R-Plugin, RStudio e outros. > > Eu particularmente uso o Gvim+Vim-r-plugin. > > Dê uma pesquisada e veja qual lhe agrada mais! > > Valeu! > > > 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 > > > ________________________________ > De: Vinícius Lionel Mateus <vinynegrelli(a)gmail.com> > Para: R-br(a)listas.c3sl.ufpr.br > Enviadas: Terça-feira, 24 de Julho de 2012 0:13 > Assunto: [R-br] Usando o R no Ubuntu > > > Prezados, > > Alguém poderia me indicar algum material para ajudar no uso do R na > plataforma Ubuntu 12.04: já instalei o software, mas gostaria de tentar > usar o GUI também, assim como no Windows. > > -- > -- > > Atenciosamente, > VINÍCIUS LIONEL MATEUS, B.Sc (http://lattes.cnpq.br/6501001637020665) > Bacharel em Química - Doutorando em Química Analítica > Laboratório de Química Atmosférica - Departamento de Química > Pontifícia Universidade Católica - Rio de Janeiro (PUC - Rio) > Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - Brasil CEP.: > 22451-900 > Telefone: (+55) (21) 3527-1327 > (+55) (21) 9358-8051 > www.puc-rio.br > > _______________________________________________ > 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/20120724/66f6d331/atta… > > > > ------------------------------ > > Message: 17 > Date: Tue, 24 Jul 2012 09:47:38 -0300 > From: nani.design(a)bol.com.br > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] duvida #Include > Message-ID: <500e996a936b7_21d65d8fc7059c(a)a4-winter16.tmail> > Content-Type: text/plain; charset="iso-8859-1" > > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120724/2c814b37/atta… > > > > ------------------------------ > > Message: 18 > Date: Tue, 24 Jul 2012 10:02:52 -0300 > From: Leonard de Assis <assis.leonard(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Usando o R no Ubuntu > Message-ID: <500E9CFC.4000606(a)gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > to gostando (acho q acostumei com ele) do rstudio > > []s > Leonard de Assis > assis <dot> leonard <at> gmail <dot> com > > Em 24/07/2012 08:11, Fabio Mathias Corrêa escreveu: > > Olá Vinícius, > > > > No Ubuntu vc tem diversas opções: Gvim+Vim-r-plugin + Emacs+ESS, > > Gedit+R-Plugin, RStudio e outros. > > > > Eu particularmente uso o Gvim+Vim-r-plugin. > > > > Dê uma pesquisada e veja qual lhe agrada mais! > > > > Valeu! > > > > > > 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 > > ------------------------------------------------------------------------ > > *De:* Vinícius Lionel Mateus <vinynegrelli(a)gmail.com> > > *Para:* R-br(a)listas.c3sl.ufpr.br > > *Enviadas:* Terça-feira, 24 de Julho de 2012 0:13 > > *Assunto:* [R-br] Usando o R no Ubuntu > > > > Prezados, > > > > Alguém poderia me indicar algum material para ajudar no uso do R na > > plataforma Ubuntu 12.04: já instalei o software, mas gostaria de > > tentar usar o GUI também, assim como no Windows. > > > > -- > > -- > > Atenciosamente, > > > > VINÍCIUS LIONEL MATEUS, B.Sc (http://lattes.cnpq.br/6501001637020665) > > Bacharel em Química - Doutorando em Química Analítica > > Laboratório de Química Atmosférica - Departamento de Química > > Pontifícia Universidade Católica - Rio de Janeiro (PUC - Rio) > > Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - > > Brasil CEP.: 22451-900 > > Telefone: (+55) (21) 3527-1327 > > (+55) (21) 9358-8051 > > www.puc-rio.br <http://www.puc-rio.br/> > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br <mailto:R-br@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. > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120724/7f04590c/atta… > > > > ------------------------------ > > Message: 19 > Date: Tue, 24 Jul 2012 10:05:35 -0300 > From: FHRB Toledo <fernandohtoledo(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, assis.leonard(a)gmail.com > Subject: Re: [R-br] Usando o R no Ubuntu > Message-ID: > < > CAN55XP62bMU4Di4x_7rTt-kPEF5SSofLjUCsC4pjARCQtUsH3Q(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Uma boa sugestão é o Emacs+ESS, que já foi dita... > > Ele entre outras, permite muitas customizações ao seu bel prazer, veja > algumas possibilidades: > http://ridiculas.wordpress.com/2012/05/10/emacs-personalizado/ > > 2012/7/24 Leonard de Assis <assis.leonard(a)gmail.com> > > > to gostando (acho q acostumei com ele) do rstudio > > > > []s > > Leonard de Assis > > assis <dot> leonard <at> gmail <dot> com > > > > Em 24/07/2012 08:11, Fabio Mathias Corrêa escreveu: > > > > Olá Vinícius, > > > > No Ubuntu vc tem diversas opções: Gvim+Vim-r-plugin + Emacs+ESS, > > Gedit+R-Plugin, RStudio e outros. > > > > Eu particularmente uso o Gvim+Vim-r-plugin. > > > > Dê uma pesquisada e veja qual lhe agrada mais! > > > > Valeu! > > > > > > 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 > > ------------------------------ > > *De:* Vinícius Lionel Mateus <vinynegrelli(a)gmail.com>< > vinynegrelli(a)gmail.com> > > *Para:* R-br(a)listas.c3sl.ufpr.br > > *Enviadas:* Terça-feira, 24 de Julho de 2012 0:13 > > *Assunto:* [R-br] Usando o R no Ubuntu > > > > Prezados, > > > > Alguém poderia me indicar algum material para ajudar no uso do R na > > plataforma Ubuntu 12.04: já instalei o software, mas gostaria de tentar > > usar o GUI também, assim como no Windows. > > > > -- > > -- > > Atenciosamente, > > > > VINÍCIUS LIONEL MATEUS, B.Sc (http://lattes.cnpq.br/6501001637020665) > > Bacharel em Química - Doutorando em Química Analítica > > Laboratório de Química Atmosférica - Departamento de Química > > Pontifícia Universidade Católica - Rio de Janeiro (PUC - Rio) > > Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - Brasil > CEP.: > > 22451-900 > > Telefone: (+55) (21) 3527-1327 > > (+55) (21) 9358-8051 > > www.puc-rio.br > > > > _______________________________________________ > > 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. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120724/2a3be7ff/atta… > > > > ------------------------------ > > Message: 20 > Date: Tue, 24 Jul 2012 10:34:14 -0300 > From: Zhu Sha Zang <zhushazang(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Usando o R no Ubuntu > Message-ID: <500EA456.7050109(a)gmail.com> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > Tem esse também: http://rkward.sourceforge.net/ > > Att > > Em 24-07-2012 10:05, FHRB Toledo escreveu: > > Uma boa sugestão é o Emacs+ESS, que já foi dita... > > > > Ele entre outras, permite muitas customizações ao seu bel prazer, veja > > algumas possibilidades: > > http://ridiculas.wordpress.com/2012/05/10/emacs-personalizado/ > > > > 2012/7/24 Leonard de Assis <assis.leonard(a)gmail.com > > <mailto:assis.leonard@gmail.com>> > > > > to gostando (acho q acostumei com ele) do rstudio > > > > []s > > Leonard de Assis > > assis <dot> leonard <at> gmail <dot> com > > > > Em 24/07/2012 08:11, Fabio Mathias Corrêa escreveu: > >> Olá Vinícius, > >> > >> No Ubuntu vc tem diversas opções: Gvim+Vim-r-plugin + Emacs+ESS, > >> Gedit+R-Plugin, RStudio e outros. > >> > >> Eu particularmente uso o Gvim+Vim-r-plugin. > >> > >> Dê uma pesquisada e veja qual lhe agrada mais! > >> > >> Valeu! > >> > >> 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 <tel:73-3680-5076> > >> > ------------------------------------------------------------------------ > >> *De:* Vinícius Lionel Mateus <vinynegrelli(a)gmail.com> > >> <mailto:vinynegrelli@gmail.com> > >> *Para:* R-br(a)listas.c3sl.ufpr.br <mailto:R-br@listas.c3sl.ufpr.br> > >> *Enviadas:* Terça-feira, 24 de Julho de 2012 0:13 > >> *Assunto:* [R-br] Usando o R no Ubuntu > >> > >> Prezados, > >> > >> Alguém poderia me indicar algum material para ajudar no uso do R > >> na plataforma Ubuntu 12.04: já instalei o software, mas gostaria > >> de tentar usar o GUI também, assim como no Windows. > >> > >> -- > >> -- > >> Atenciosamente, > >> > >> VINÍCIUS LIONEL MATEUS, B.Sc ( > http://lattes.cnpq.br/6501001637020665) > >> Bacharel em Química - Doutorando em Química Analítica > >> Laboratório de Química Atmosférica - Departamento de Química > >> Pontifícia Universidade Católica - Rio de Janeiro (PUC - Rio) > >> Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - > >> Brasil CEP.: 22451-900 > >> Telefone: (+55) (21) 3527-1327 > >> <tel:%28%2B55%29%20%2821%29%203527-1327> > >> (+55) (21) 9358-8051 <tel:%28%2B55%29%20%2821%29%209358-8051> > >> www.puc-rio.br <http://www.puc-rio.br/> > >> > >> _______________________________________________ > >> R-br mailing list > >> R-br(a)listas.c3sl.ufpr.br <mailto:R-br@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 <mailto:R-br@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 <mailto:R-br@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. > > > > > -- > > --- > Rodolfo Timóteo da Silva > Linux Counter: 359362 > msn: zhushazang(a)gmail.com > skype: zhushazang > > Ribeirão Preto - SP > > > ------------------------------ > > Message: 21 > Date: Tue, 24 Jul 2012 11:03:59 -0300 > From: Vinícius Lionel Mateus <vinynegrelli(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Usando o R no Ubuntu > Message-ID: > <CAMJ4BpBfXnR-9jv8NQ1J6ASxax5LMWkcCt2Enxp+= > 6oEuab0MQ(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Obrigado a todos. > Eu optei pelo RStudio. > > O JGR não rodou na minha versão, agora preciso desintalar: alguém poderia > me dar um 'help', novamente. > > Em 24 de julho de 2012 10:34, Zhu Sha Zang <zhushazang(a)gmail.com> > escreveu: > > > Tem esse também: http://rkward.sourceforge.net/ > > > > Att > > > > Em 24-07-2012 10:05, FHRB Toledo escreveu: > > > >> Uma boa sugestão é o Emacs+ESS, que já foi dita... > >> > >> Ele entre outras, permite muitas customizações ao seu bel prazer, veja > >> algumas possibilidades: > >> http://ridiculas.wordpress.**com/2012/05/10/emacs-**personalizado/< > http://ridiculas.wordpress.com/2012/05/10/emacs-personalizado/> > >> > >> 2012/7/24 Leonard de Assis <assis.leonard(a)gmail.com > >> <mailto:assis.leonard@gmail.**com <assis.leonard(a)gmail.com>>> > >> > >> > >> to gostando (acho q acostumei com ele) do rstudio > >> > >> []s > >> Leonard de Assis > >> assis <dot> leonard <at> gmail <dot> com > >> > >> Em 24/07/2012 08:11, Fabio Mathias Corrêa escreveu: > >> > >>> Olá Vinícius, > >>> > >>> No Ubuntu vc tem diversas opções: Gvim+Vim-r-plugin + Emacs+ESS, > >>> Gedit+R-Plugin, RStudio e outros. > >>> > >>> Eu particularmente uso o Gvim+Vim-r-plugin. > >>> > >>> Dê uma pesquisada e veja qual lhe agrada mais! > >>> > >>> Valeu! > >>> > >>> 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 <tel:73-3680-5076> > >>> ------------------------------**------------------------------** > >>> ------------ > >>> *De:* Vinícius Lionel Mateus <vinynegrelli(a)gmail.com> > >>> <mailto:vinynegrelli@gmail.com**> > >>> *Para:* R-br(a)listas.c3sl.ufpr.br <mailto:R-br@listas.c3sl.ufpr. > **br<R-br(a)listas.c3sl.ufpr.br> > >>> > > >>> *Enviadas:* Terça-feira, 24 de Julho de 2012 0:13 > >>> *Assunto:* [R-br] Usando o R no Ubuntu > >>> > >>> > >>> Prezados, > >>> > >>> Alguém poderia me indicar algum material para ajudar no uso do R > >>> na plataforma Ubuntu 12.04: já instalei o software, mas gostaria > >>> de tentar usar o GUI também, assim como no Windows. > >>> > >>> -- > >>> -- > >>> Atenciosamente, > >>> > >>> VINÍCIUS LIONEL MATEUS, B.Sc (http://lattes.cnpq.br/** > >>> 6501001637020665 <http://lattes.cnpq.br/6501001637020665>) > >>> Bacharel em Química - Doutorando em Química Analítica > >>> Laboratório de Química Atmosférica - Departamento de Química > >>> Pontifícia Universidade Católica - Rio de Janeiro (PUC - Rio) > >>> Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - > >>> Brasil CEP.: 22451-900 > >>> Telefone: (+55) (21) 3527-1327 > >>> <tel:%28%2B55%29%20%2821%29%**203527-1327> > >>> (+55) (21) 9358-8051 <tel:%28%2B55%29%20%2821%29%**209358-8051> > >>> www.puc-rio.br <http://www.puc-rio.br/> > >>> > >>> ______________________________**_________________ > >>> R-br mailing list > >>> R-br(a)listas.c3sl.ufpr.br <mailto:R-br@listas.c3sl.ufpr.**br< > R-br(a)listas.c3sl.ufpr.br> > >>> > > >>> > >>> https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-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< > 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 <mailto:R-br@listas.c3sl.ufpr.**br< > R-br(a)listas.c3sl.ufpr.br> > >>> > > >>> > >>> https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-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< > 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 <mailto:R-br@listas.c3sl.ufpr.**br< > R-br(a)listas.c3sl.ufpr.br> > >> > > >> > >> https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-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< > 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< > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br> > >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-**guia< > http://www.leg.ufpr.br/r-br-guia>) > >> e forneça código mínimo reproduzível. > >> > >> > > > > -- > > > > --- > > Rodolfo Timóteo da Silva > > Linux Counter: 359362 > > msn: zhushazang(a)gmail.com > > skype: zhushazang > > > > Ribeirão Preto - SP > > > > ______________________________**_________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br > > https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-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< > http://www.leg.ufpr.br/r-br-guia>) > > e forneça código mínimo reproduzível. > > > > > > -- > -- > Atenciosamente, > > VINÍCIUS LIONEL MATEUS, B.Sc (http://lattes.cnpq.br/6501001637020665) > Bacharel em Química - Doutorando em Química Analítica > Laboratório de Química Atmosférica - Departamento de Química > Pontifícia Universidade Católica - Rio de Janeiro (PUC - Rio) > Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - Brasil CEP.: > 22451-900 > Telefone: (+55) (21) 3527-1327 > (+55) (21) 9358-8051 > www.puc-rio.br > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120724/5b7c8009/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 17, assunto 24 > ***************************************** > -- Gilenio Borges Fernandes Universidade Federal da Bahia Instituto de Matemática Departamento de Estatística Av. Adhemar de Barros, s/n – Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6280/6336 Fax: (071)3283-6276 URL: http://lattes.cnpq.br/6764860618464860

1 0

duvida #Include
by nani.design＠bol.com.br 24 Jul '12

24 Jul '12

3 5

Variogram model 3D
by Tiago Vieira 23 Jul '12

23 Jul '12

Olá a todos, Minha dúvida desta vez é: Como posso alterar o Alcance Vertical na função vgm do pacote gstat? Pois estou realizando uma krigagem 3D e consigo fazer isso no gslib, mas ainda não consegui no gstat! Grato, Tiago Vieira.

2 2

Re: [R-br] comparações multiplas utilizando Kruskall- walis
by Ivan Bezerra Allaman 23 Jul '12

23 Jul '12

Segundo informações no site da cran, o pacote pgirmess depende do pacote rgdal, e este está disponível no cran. \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

comparações multiplas utilizando Kruskall- walis
by Marcelo Claro de Souza 23 Jul '12

23 Jul '12

Não consigo instalar o pacote pgirmess para rodar comparações multiplas pelo Kruskall-walis pois é necessário instalar tb instalar o pacote gdal e o mesmo não está disponível para a versão 2.14.1 do linux. Já tentei instalar esse pacote pelo terminal usando comando sudo apt-get build-dep r-cran-rgdal e apareceu a mesma mensagem de erro. Existe algum outro pacote que eu possa utilizar no lugar pgirmess que faça a mesma coisa que a função kruskalmc? Muito obrigado pela ajuda. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

2 1

Extrair palavras com pontuações!
by Ivan Bezerra Allaman 23 Jul '12

23 Jul '12

Bom dia senhores! Suponhamos que eu tenha o seguinte vetor: nfa <- c('A','JA','K','K:P','P:K','K:P:N') Se eu quero extrair as palavras com a pontuação : então: nfa[grep('[[:punct:]]',nfa)] ou nfa[grep(':',nfa)] No entanto, como eu faço para extrair as palavras que tem as duas pontuações, ou seja, no caso acima 'K:P:N'? Desde já grato! (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}

2 2

Reduzindo a fonte em party::plot.BinaryTree
by Pedro Emmanuel Alvarenga Americano do Brasil 23 Jul '12

23 Jul '12

Amigos de R, Utilizando o pacote party para fazer uma ávore de decisão utilizei os seguintes comandos. ct.BF <- ctree(...) plot(ct.BF,drop_terminal=F,terminal_panel=node_barplot(ct.BF)) A minha questão é que essa árvore possui 15 nos e 8 são terminais. Então o boxplot fica com diversos textos automaticos da função sobrepostos com o tamanho padrão da fonte. Fucei e fucei mas não encontrei um jeito de reduzir o tamanho da fonte. Aguem saberia como faze-lo? plot(ct.BF,drop_terminal=F,terminal_panel=node_barplot(ct.BF,beside=T)) No exemplo acima, com o beside=T, o problema melhora um pouco, mas um dos labels das barras não aparece. Uma alternativa seria colocar o texto das barras na vertical. Alguem saberia como colocar o texto na vertical? Obrigado. Dr. Pedro Emmanuel A. A. do Brasil Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365, CEP 21040-360, Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

1 0

Base de Coordenadas Geográficas
by Listeiro 037 20 Jul '12

20 Jul '12

Olá a todos. Há algum lugar que disponibilize gratuitamente base de dados de coordenadas geográficas de cidades do Brasil e do mundo? No caso seria Cidade, Latitude e Longitude. Para cruzar com outros dados. Pode ser em qualquer formato de arquivo que seja exportável para outros. Formato SQL, CVS, TXT, XLS, DBF, MDB, ODF. Menos PDF, que pode ter imagens embutidas ao invés de texto. Neste caso seria péssimo usar um OCR prá tentar capturar o texto. Desde já agradeço qualquer informação.

2 1

Combinações de elementos em grupos distintos
by Leonardo Aguirre 20 Jul '12

20 Jul '12

Pessoal, Preciso criar grupos de indivíduos dessa forma: Tenho 64 indivíduos, para serem divididos em 8 grupos do mesmo tamanho. Portanto, ficam 8 indivíduos por grupo. Preciso replicar a criação desses 8 grupos 8 vezes, de tal forma que o indivíduo i apareça no mesmo grupo do indivíduo j apenas uma vez. Ou seja, se na primeira repetição o João saiu no mesmo grupo que a Maria, eles não podem ficar no mesmo grupo em nenhuma outra das 7 repetições. Sendo assim, suponha que na primeira etapa eu tenha esses grupos definidos: Grupo 1: 1,2,3,4,5,6,7,8 Grupo 2: 9,10,11,12,13,14,15,16 ... Grupo 8: 57,58,59,60,61,62,63,64 Então, nas próximas 7 replicações o indivíduo 1 só poderia sair no mesmo grupo dos indivíduos de 9 a 64. E assim sucessivamente, pra todos os outros indivíduos. Tentei montar esses grupos no R mas não consegui nem sair do zero, por isso não estou enviando um exemplo do programa em R. Alguém teria alguma idéia de como implementar isso no R? Muito obrigado, Leonardo.

7 7

Erro ao ler SQL
by Edson Lira 20 Jul '12

20 Jul '12

Caros amigos, estou lendo um banco de dados de doadores em SQL, uma das variá veis é a datahoradoação, quando leio o mês de agosto, tenho um erro, para identificar o dia fui fazendo a leitura por período, identifiquei que o erro está no dia 08 de agosto, cuja mensagem de erro mostro abaixo: + and to_char(doacao_dw.datahoradoacao,'DD/MM/YYYY') >= to_date('28/08/1992','DD/MM/YYYY') + and to_char(doacao_dw.datahoradoacao,'DD/MM/YYYY') <= to_date('28/08/1992','DD/MM/YYYY')" + ) Erro em as.POSIXlt.character(x, tz, ...) : character string is not in a standard unambiguous format Já exportei os dados para csv e ainda não identifiquei qual o erro, alguém tem alguma sugestão? Abaixo o link com o arquivo em csv. No sql consigo rodar o dia, mas no R me dá a mensagem citada acima. coloquei o dia 27 que não tem problema e o dia 28 que apresenta o problema. http://www.datafilehost.com/download-c0c28ad5.html [ ]'s. Edson Lira Estatístico Manaus-Amazonas

3 10

simulação binomial
by Natalia Martins 19 Jul '12

19 Jul '12

Prezados, Estou simulando dados diários de chuva tentei simular utilizando as funções fitdistr (para obtenção das estimativas dos parâmetros) e rgamma, no entanto tenho bastante zeros (dias sem chuva), com isso não é possível a obtenção dos parâmetros. Erro apresentado no R: Erro em optim(x = c(0, 29.9, 0, 0.7, 0, 0, 8.7, 0, 0, 26.8, 0, 19.4, : valor inicial em vmmin não é finito Posto isso resolvi transformar meus dados 0=não choveu e 1=choveu e simular uma binomial e ver a partir dai se choveu eu simulo o quanto choveu. No entanto para a função fitdistr não é definido a distribuição binomial. Sendo assim alguém saberia como eu posso obter o valor do parâmetro p? Antecipadamente meu obrigada. 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 4

Arquivo metafile à partir do linux
by Marcelo Claro de Souza 19 Jul '12

19 Jul '12

Tem como salvar figuras como metafile em R à partir do Linux? Muito obrigado. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

2 1

Diferença de resultados entre 64 e 32 bits
by Rodrigo Coster 19 Jul '12

19 Jul '12

Caros, programei uma rotina para estimar por máxima verossimilhança os parâmetros de uma cópula e para ver se estava certo comparei os resultados com o comando do pacote copula. Encontrei diferenças apenas na 5a casa decimal em diante, que considerei como sendo por causa do método numérico utilizado. Só que, ao fazer a mesma comparação num computador 64 bits os resultados são bastante divergentes (o meu código muda o valor estimado, enquanto o pacote copula mantem), mudando na 2a casa decimal (no caso o parâmetro é a correlação, então a 2a casa decimal é bem importante). Dai me bateu a seguinte dúvida: em qual confiar? Código: require(copula) require(MASS) set.seed(31415) normCop <- function(param,data) { n <- nrow(data) if (length(param) != n) { param <- rep(param,nrow(data)) } cop <- mapply(normalCopula,param=param,MoreArgs=list(dim=2)) datalist <- apply(data,1,list) for (i in 1:n) { datalist[[i]] <- datalist[[i]][[1]] } out <- -sum(log(mapply(dcopula,copula=cop,u=datalist))) if (out == Inf) { out = exp(100) } return(out) } a <- matrix(0,20,2) n <- 100 Sigma <- matrix(c(10,3,3,2),2,2) for (j in 1:20) { data <- mvrnorm(n=n, rep(0, 2), Sigma) data <- apply(data,2,rank)/(n+1) fitNormCop <- function(data) { optim(cor(data)[2],normCop,data=data, lower = 0, upper = .9999,method="L-BFGS-B") } a[j,1] <- fitNormCop(data)$par # Meu a[j,2] <- fitCopula(normalCopula(.2,2), data, method="ml")@estimate # Pacote copula } E aqui um gráfico de dispersão comparando todos: http://img411.imageshack.us/img411/7527/92588280.png []'s

4 11

Como precificar um treinamento de R básico?
by Paulo Nogueira 18 Jul '12

18 Jul '12

Ola a todos do grupo. Faz um tempo que não recorro a ajuda de vocês, mas dessa vez tenho uma questão um pouco inusitada. Fui questionado sobre a viabilidade de preparar um treinamento básico em R para um grupo de profissionais da área de saúde que, a grosso modo, trabalham com genética. O escopo está razoavelmente bem definido, mas ainda em discussão, e o material em progresso. A grande questão é que não sei precificar um produto desses... Gostaria da ajuda de quem já teve alguma experiência semelhante e possa ajudar em um benchmarking. Seria uma turma de mais ou menos 15 pessoas, com aulas expositivas semanais e material para prática fora de classe. Duração esperada de 3 meses. O cliente é uma empresa da área de saúde. Tenho experiência com outros treinamentos, por exemplo SAS, que é bastante comercial e procurado por grandes instituições. Mas precificar R é duro... ainda mais para quem aprendeu sozinho. Alguém já trabalhou com treinamento R e pode ajudar? Abraços Paulo

6 8

Problemas ao instalar pacote EBImage em windows.
by Augusto Ribas 18 Jul '12

18 Jul '12

Boa tarde a todos. Eu estou com problemas para instalar o pacote EBImage aqui. Eu instalo ele com os comandos: source("http://bioconductor.org/biocLite.R") biocLite("EBImage") Mas qd vou abrir ele eu tenho: > library(EBImage) Carregando pacotes exigidos: abind Error in inDL(x, as.logical(local), as.logical(now), ...) : impossível carregar objeto compartilhado 'C:/Program Files/R/R-2.15.1/library/EBImage/libs/i386/EBImage.dll': LoadLibrary failure: Não foi possível encontrar o módulo especificado. Erro: package/namespace load failed for ‘EBImage’ Eu segui o manual aqui: http://www.bioconductor.org/packages/release/bioc/vignettes/EBImage/inst/do… E instalei o o gtk antes (gtk+-bundle_2.24.10-20120208_win32), e mudei no meu computador->propriedades->Configurações Avançadas do sistema-> Variaveis de ambiente->Variavel do sistema Novo coloquei o caminho pro gtk E também instalei o ImageMagick, marcando pra colocar as opções `Update exe- cutable search path' e `Install developement headers and libraries' . Que foram a segunda e a quinta caixinha em portugues. Eu instalei a versão 6.7.8.... Q16 E por ultimo ainda instalei o Microsoft VC++ 2008 Redistributable Package Ainda sim não vai. Abrindo pelo Tinn-R da ainda a msg: O programa não pode ser iniciado porque está faltando o libgdk-win32-2.0-0.dll no seu computador. Tente reinstalá-lo para resolver o problema. Mas esse arquivo ta la no gtk. Alguem tem alguma ideia? Segue meu sessionInfo se for de alguma ajuda: > sessionInfo() R version 2.15.1 (2012-06-22) Platform: i386-pc-mingw32/i386 (32-bit) locale: [1] LC_COLLATE=Portuguese_Brazil.1252 LC_CTYPE=Portuguese_Brazil.1252 [3] LC_MONETARY=Portuguese_Brazil.1252 LC_NUMERIC=C [5] LC_TIME=Portuguese_Brazil.1252 attached base packages: [1] stats graphics grDevices datasets utils methods base other attached packages: [1] abind_1.4-0 loaded via a namespace (and not attached): [1] tools_2.15.1 -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

3 7

pacote easyanova no CRAN
by Emmanuel Arnhold 18 Jul '12

18 Jul '12

Colegas da lista. A quem interessar esta disponível no cran o pacote easyanova. O pacote possui funções que automatizam anova e cia em vários delineamentos, inclusive em casos de desbalanceamento. http://cran.r-project.org/web/packages/easyanova/index.html intall.packages("easyanova") require(easyanova) ?ea1 ?ea2 ?ea3 ?ea4 ?ea5 ?ea6 ?ea7 ?ea8

2 1

Colocar elipse de confiança num gráfico
by Mauro Sznelwar 17 Jul '12

17 Jul '12

Desculpem minha insistência, até tentei como podem ver, mas não consegui, alguém descobre faz para inserir uma elipse de confiança neste gráfico acima de Análise discriminante: plot(marajoara.norm.ldakm1,dimen=2,ylim=c(-5,5),xlim=c(20,35), cex=1,col=marajoara.norm.km1$cluster) library(ellipse) lines(ellipse(marajoara.norm.km1$cluster, centre = c(27,0)), type = "l") Abaixo tem o script que gerou este gráfico, e estou enviando também o arquivo de dados. http://www.datafilehost.com/download-1eccadf2.html marajoara=read.csv('Marajoara Mario.csv',head=T,dec=',',sep=';') marajoara.norm=sapply(data.frame(marajoara[,2:15]),scale) #marajoara.norm=(data.frame(marajoara[,2:15])) #summary(marajoara.norm) #library(fBasics) #basicStats(marajoara.norm) marajoara.norm.h=hclust(dist(marajoara.norm),method='average') #marajoara.norm.h=hclust(dist(marajoara.norm)^2,method='ward') plot(marajoara.norm.h,main='Marajoara original-Ward Linkage\n Square Euclidian Distance\n Hclust Dendograma') #pdf('figuras/Marajoara-original-dendograma-hclust-average.pdf') library(cluster) marajoara.norm.d=diana(dist(marajoara.norm)) plot(marajoara.norm.d,main='Marajoara original-Divisive Dendograma') #pdf('figuras/Marajoara-original-dendograma-divisive.pdf') #Componentes Principais initial=tapply(marajoara.norm,list(rep(cutree(marajoara.norm.h,3), ncol(marajoara.norm)),col(marajoara.norm)),mean) marajoara.norm.km1=kmeans(marajoara.norm,initial) marajoara.norm.pca=princomp(marajoara.norm) marajoara.norm.px=predict(marajoara.norm.pca) library(MASS) eqscplot(marajoara.norm.px[,1:2],type='n',xlab='component 1',ylab='component 2') text(marajoara.norm.px[,1:2],labels=as.character(marajoara.norm.km1$cluster), col=as.numeric(marajoara.norm.km1$cluster)) title('PCA Marajoara original - clusters kmeans\n centros iniciais corte dendograma hclust 3') #identify(marajoara.norm.px,n=5) #pdf('figuras/PCA-Marajoara-original-centros-kmeans\n centros inicias de corte dendograma hclust 3') initial=tapply(marajoara.norm,list(rep(cutree(marajoara.norm.d,3), ncol(marajoara.norm)),col(marajoara.norm)),mean) marajoara.norm.km2=kmeans(marajoara.norm,initial) marajoara.norm.pca=princomp(marajoara.norm) marajoara.norm.px=predict(marajoara.norm.pca) eqscplot(marajoara.norm.px[,1:2],type='n',xlab='component 1',ylab='component 2') text(marajoara.norm.px[,1:2],labels=as.character(marajoara.norm.km2$cluster), col=as.numeric(marajoara.norm.km2$cluster)) title('PCA Marajoara original - clusters kmeans\n centros iniciais corte dendograma diana 3') #identify(marajoara.norm.px,n=5) #pdf('figuras/PCA-Marajoara-original-centros-kmeans\n centros inicias de corte dendodiana-3.pdf') marajoara.norm.km=kmeans(marajoara.norm,3) marajoara.norm.pca=princomp(marajoara.norm) marajoara.norm.px=predict(marajoara.norm.pca) eqscplot(marajoara.norm.px[,1:2],type='n',xlab='component 1',ylab='component 2') text(marajoara.norm.px[,1:2],labels=as.character(marajoara.norm.km$cluster), col=marajoara.norm.km$cluster) title('PCA Marajoara original - clusters kmeans 3 centros') #identify(marajoara.norm.px,n=5) #pdf('figuras/PCA-Marajoara-original-clusters-kmeans-3-centros.pdf') #Análise Discriminante Linear marajoara.norm.ldakm1=lda(marajoara.norm,as.character(marajoara.norm.km1$cluster)) plot(marajoara.norm.ldakm1,dimen=2,cex=1,col=marajoara.norm.km1$cluster) title('LDA Marajoara-Grupo do K-Means - corte dendograma hclust em 3') #pdf('figuras/LDA-Marajoara-grupos-kmeans-dendohclust-3.pdf') #fazendo um zoom no gráfico corta um ponto plot(marajoara.norm.ldakm1,dimen=2,ylim=c(-5,5),xlim=c(20,35), cex=1,col=marajoara.norm.km1$cluster) title('LDA Marajoara-parte do gráfico para visualização\n Grupos do K-Means-corte dendograma hclust em 3') #identify(marajoara.norm.km1$cluster,n=5) #pdf('figuras/LDA-Marajoara-grupos-kmeans-dendodiana-3.pdf') library(ellipse) lines(ellipse(marajoara.norm.km1$cluster, centre = c(27,0)), type = "l")

1 0

duvida controle de qualidade
by Natalia Martins 17 Jul '12

17 Jul '12

Por favor estou com um dúvida em controle de qualidade e ficaria muito grata se alguem pudesse tirá-la, ou contribuir de alguma forma. Bem é o seguinte, tenho uma máquina que é programada para colocar uma quantidade exata de adubo em um ponto especificado no campo via GPS, ou seja, a maquina vai passando pelo campo e derramando o adubo na quantidade exigida para aquela coordenada (x,y), as quantidades vão sendo mudadas a cada ponto. Eu sei a quantidade que a maquina *deve* derrubar em todo o campo e as quantidades que a maquina *realmente derruba.* Quando a maquina termina o seu serviço um equipamento instalado nela me informa o que ela realmente derrubou de adubo em cada coordenada, assim eu tenho que verificar se esta maquina está sob controle estatístico de qualidade. - tenho o que a maquina deve derramar por coordenada no campo todo - tenho o que a maquina realmente derramou por coordenada no campo todo Minhas duvidas em relação a esse processo são: Para verificar se meu processo esta sob controle eu fiz a diferença entre o que a maquina realmente aplicou e o que deveria ser aplicado, eu posso aplicar a teoria de controle de qualidade a esta diferença? Tenho todos os pontos (2700 pontos) e não uma amostra, a teoria também pode ser aplicada a todos os dados? Quando faço o teste de normalidade de shapiro wilks este mostra me que meus dados não apresentam normalidade, devo transformá-los? Muito obrigada, -- 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

coeficiente de concordância kappa
by Edson Lira 17 Jul '12

17 Jul '12

Caros amigos, em quais pacotes agora posso obter o coeficiente de kappa? Visto que o DiagnosisMed que era o que eu usava não está mais disponível. [ ]'s. Edson Lira Estatístico Manaus-Amazonas

4 3

Combinar linhas entre conjuntos...
by andrebvs 17 Jul '12

17 Jul '12

3 6

Análise Fatorial com dummy
by Felipe Buchbinder 16 Jul '12

16 Jul '12

Professores, Eu tenho uma base de dados com variáveis dummies e gostaria de fazer uma redução de dimensionalidade. Inicialmente pensei em uma análise fatorial com correlação tetracórica. Porém, a matriz de correlação tinha correlações unitárias fora da diagonal principal sem razão aparente (olhando o gráfico, havia pontos nos quatro quadrantes, mas mesmo assim a correlação dava 1!) Como minhas variáveis não são discretizações de variáveis contínuas, não sei se o uso de correlação tetracórica é adequado. Pensei em usar o *phi* de Pearson como medida de correlação, mas ele dá idêntico à correlação de Pearson. Outra idéia foi fazer a extração usando análie por componentes principais para dados categóricos (CATPCA) mas, pelo que andei pesquisando, o resultado não dá muito diferente de uma análise fatorial convencional (com correlação de Pearson e extração por PCA) Eu queria poder contar com uma opinião sobre o que fazer: Usar correlação tetracórica (neste caso, como eu resolvo o fato da matriz de correlação não ser positiva semidefinida?) Usar o phi de Pearson como medida de correlaçao Usar a correlação de Pearson, pura se simplesmente Extrair os fatores por CATPCA ou outra alternativa, que é sempre bem-vinda. Desde já agradeço a atenção e a colaboração, Felipe Buchbinder

2 2

Re: [R-br] Ajuda programação!
by Fernando Mayer 16 Jul '12

16 Jul '12

(encaminhando para a lista para manter arquivado) O tab ali é o objeto criado pelo Ivan na sua mensagem inicial. A única modificação que fiz foi no nome das linhas: k1 <- c(1.8850,2.1575,2.4775,2.7975) k2 <- c(2.9375,2.2475,2.7300,2.9950) tab <- matrix(c(k1,k2),ncol=2) rownames(tab) <- c('p1n1','p1n2','p2n1','p2n2') []s, --- Fernando Mayer Universidade Federal de Santa Catarina - UFSC Departamento de Ecologia e Zoologia - ECZ/CCB URL: http://sites.google.com/site/fernandomayer e-mail: fernandomayer [@] gmail.com 2012/7/16 Mauro Sznelwar <sznelwar(a)uol.com.br>: > Estava olhando o vosso script, e não consegui reproduzir onde tem o tab, > alega que não reconhece este objeto. Precisa de alguma biblioteca pra isto? > > > Também não sei se entendi a pergunta direito, mas se vc quer apenas um > "atalho" para chamar determinada linha usando apenas os valores de "a" > e "b", pode fazer algo como > > f <- function(a, b){ > x <- paste("p", a, "n", b, sep = "") > return(x) > } > >> f(a = 1, b = 2) > [1] "p1n2" >> f(a = 1, b = 1) > [1] "p1n1" >> tab[f(1,1), ] > k1 k2 > 1.8850 2.9375 >> as.vector(tab[f(1,1),]) > [1] 1.8850 2.9375 > > > []s, > > --- > Fernando Mayer > Universidade Federal de Santa Catarina - UFSC > Departamento de Ecologia e Zoologia - ECZ/CCB > URL: http://sites.google.com/site/fernandomayer > e-mail: fernandomayer [@] gmail.com > >

1 0

Duvidas em análise cluster
by edlayne moraes 16 Jul '12

16 Jul '12

Boa tarde, Inseri uma planilha com os dados em .csv, coloquei na linha 1º o nome dos meus dados. Ao importar a tabela eu mandei ler essa linha como rótulo das variáveis. Para fazer o tratamento eu instalei o pacote "vegan" e ao mandar plotar o gráfico o programa informa um erro de "rowsums (x, na.rm=TRUE): x deve ser numérico". Alguém sabe me dizer o que há de erro e se existe outra maneira de fazer o cluster no programa? Desde já agradeço a atenção Att, -- Edlayne Pimentel de Moraes**** Graduanda em Geologia-UFES**** (28) 81117-0585//(28) 3552-0194

4 3

Ajuda programação!
by Ivan Bezerra Allaman 16 Jul '12

16 Jul '12

Bom dia senhores! Suponhamos que tem-se dois níveis de um fator "p" e dois níveis de um fator "n", no qual os níveis do fator "p" irei representar pela letra "a" e os níveis do fator "n" pela letra "b". Então tem-se a seguinte configuração: k1 <- c(1.8850,2.1575,2.4775,2.7975) k2 <- c(2.9375,2.2475,2.7300,2.9950) tab <- matrix(c(k1,k2),ncol=2) rownames(tab) <- c('panb','panb,'panb','panb') colnames(tab) <- c('k1','k2') tab k1 k2 panb 1,8850 2,9375 panb 2,1575 2,2475 panb 2,4775 2,7300 panb 2,7975 2,9950 Vamos supor que eu queira os valores do nível 1 do fator "p" e nível 1 do fator "n". Então, a=1 e b=1; p1n1 1,8850 2,9375 Agora nível 1 do fator "p" e nível 2 do fator "n". Então, a=1, b=2; p1n2 2,1575 2,2475 Tem como automatizar isto em função destes dois parâmetros? algo como: medias = as.vector(tab[ab, ])???? A coluna não interessa, podem ficar tranquilo! Apenas as linhas interessam de acordo com os parâmetros. Se alguém puder me dar uma ajuda ficarei grato! Abraços! 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}

4 3

Colocar elipse de confiança num gráfico
by Mauro Sznelwar 16 Jul '12

16 Jul '12

Como eu faço para inserir uma elipse de confiança neste gráfico acima de Análise discriminante, não sei como usar o comando: plot(marajoara.norm.ldakm1,dimen=2,ylim=c(-5,5),xlim=c(20,35), cex=1,col=marajoara.norm.km1$cluster) library(ellipse) plot(ellipse(S, centre = , t = ), type = "l") Abaixo tem o script que gerou este gráfico, e estou enviando também o arquivo de dados. http://www.datafilehost.com/download-1eccadf2.html marajoara=read.csv('Marajoara Mario.csv',head=T,dec=',',sep=';') marajoara.norm=sapply(data.frame(marajoara[,2:15]),scale) #marajoara.norm=(data.frame(marajoara[,2:15])) #summary(marajoara.norm) #library(fBasics) #basicStats(marajoara.norm) marajoara.norm.h=hclust(dist(marajoara.norm),method='average') #marajoara.norm.h=hclust(dist(marajoara.norm)^2,method='ward') plot(marajoara.norm.h,main='Marajoara original-Ward Linkage\n Square Euclidian Distance\n Hclust Dendograma') #pdf('figuras/Marajoara-original-dendograma-hclust-average.pdf') library(cluster) marajoara.norm.d=diana(dist(marajoara.norm)) plot(marajoara.norm.d,main='Marajoara original-Divisive Dendograma') #pdf('figuras/Marajoara-original-dendograma-divisive.pdf') #Componentes Principais initial=tapply(marajoara.norm,list(rep(cutree(marajoara.norm.h,3), ncol(marajoara.norm)),col(marajoara.norm)),mean) marajoara.norm.km1=kmeans(marajoara.norm,initial) marajoara.norm.pca=princomp(marajoara.norm) marajoara.norm.px=predict(marajoara.norm.pca) library(MASS) eqscplot(marajoara.norm.px[,1:2],type='n',xlab='component 1',ylab='component 2') text(marajoara.norm.px[,1:2],labels=as.character(marajoara.norm.km1$cluster), col=as.numeric(marajoara.norm.km1$cluster)) title('PCA Marajoara original - clusters kmeans\n centros iniciais corte dendograma hclust 3') #identify(marajoara.norm.px,n=5) #pdf('figuras/PCA-Marajoara-original-centros-kmeans\n centros inicias de corte dendograma hclust 3') initial=tapply(marajoara.norm,list(rep(cutree(marajoara.norm.d,3), ncol(marajoara.norm)),col(marajoara.norm)),mean) marajoara.norm.km2=kmeans(marajoara.norm,initial) marajoara.norm.pca=princomp(marajoara.norm) marajoara.norm.px=predict(marajoara.norm.pca) eqscplot(marajoara.norm.px[,1:2],type='n',xlab='component 1',ylab='component 2') text(marajoara.norm.px[,1:2],labels=as.character(marajoara.norm.km2$cluster), col=as.numeric(marajoara.norm.km2$cluster)) title('PCA Marajoara original - clusters kmeans\n centros iniciais corte dendograma diana 3') #identify(marajoara.norm.px,n=5) #pdf('figuras/PCA-Marajoara-original-centros-kmeans\n centros inicias de corte dendodiana-3.pdf') marajoara.norm.km=kmeans(marajoara.norm,3) marajoara.norm.pca=princomp(marajoara.norm) marajoara.norm.px=predict(marajoara.norm.pca) eqscplot(marajoara.norm.px[,1:2],type='n',xlab='component 1',ylab='component 2') text(marajoara.norm.px[,1:2],labels=as.character(marajoara.norm.km$cluster), col=marajoara.norm.km$cluster) title('PCA Marajoara original - clusters kmeans 3 centros') #identify(marajoara.norm.px,n=5) #pdf('figuras/PCA-Marajoara-original-clusters-kmeans-3-centros.pdf') #Análise Discriminante Linear marajoara.norm.ldakm1=lda(marajoara.norm,as.character(marajoara.norm.km1$cluster)) plot(marajoara.norm.ldakm1,dimen=2,cex=1,col=marajoara.norm.km1$cluster) title('LDA Marajoara-Grupo do K-Means - corte dendograma hclust em 3') #pdf('figuras/LDA-Marajoara-grupos-kmeans-dendohclust-3.pdf') #fazendo um zoom no gráfico corta um ponto plot(marajoara.norm.ldakm1,dimen=2,ylim=c(-5,5),xlim=c(20,35), cex=1,col=marajoara.norm.km1$cluster) title('LDA Marajoara-parte do gráfico para visualização\n Grupos do K-Means-corte dendograma hclust em 3') #identify(marajoara.norm.km1$cluster,n=5) #pdf('figuras/LDA-Marajoara-grupos-kmeans-dendodiana-3.pdf') library(ellipse) plot(ellipse(S, centre = , t = ), type = "l")

1 0

Re: [R-br] Digest R-br, volume 17, assunto 14
by Gilenio Borges Fernandes 16 Jul '12

16 Jul '12

Muito obrigado ao Elias T. Krainski. A sugestão foi fundamental para a simulação que pretendo fazer Gilenio Fernandes Em 16 de julho de 2012 12:00, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: > 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. Ajuda programação! (Ivan Bezerra Allaman) > 2. duvida controle de qualidade (Natalia Martins) > 3. Re: duvida controle de qualidade (Eder David Borges da Silva) > 4. Re: Ajuda programação! (Lucas Barbosa) > 5. Re: duvida controle de qualidade (Natalia Martins) > 6. Re: duvida controle de qualidade > (henrique jose de paula alves paula alves) > 7. Re: duvida controle de qualidade (Elton Gean Araújo Araújo) > 8. Re: Ajuda programação! (Walmes Zeviani) > 9. Re: Ajuda para Simulação de variáveis binárias > correlacionadas no R (Elias T. Krainski) > 10. Duvidas em análise cluster (edlayne moraes) > 11. Re: Duvidas em análise cluster (Natalia Martins) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sun, 15 Jul 2012 09:07:18 -0700 (PDT) > From: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> > To: R Brasil <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Ajuda programação! > Message-ID: > <1342368438.19845.YahooMailNeo(a)web161801.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Bom dia senhores! > > Suponhamos que tem-se dois níveis de um fator "p" e dois níveis de um > fator "n", no qual os níveis do fator "p" irei representar pela letra "a" e > os níveis do fator "n" pela letra "b". Então tem-se a seguinte configuração: > > k1 <- c(1.8850,2.1575,2.4775,2.7975) > k2 <- c(2.9375,2.2475,2.7300,2.9950) > tab <- matrix(c(k1,k2),ncol=2) > rownames(tab) <- c('panb','panb,'panb','panb') > colnames(tab) <- c('k1','k2') > tab > k1 k2 > panb 1,8850 2,9375 > panb 2,1575 2,2475 > panb 2,4775 2,7300 > panb 2,7975 2,9950 > > > Vamos supor que eu queira os valores do nível 1 do fator "p" e nível 1 do > fator "n". Então, a=1 e b=1; > > p1n1 1,8850 2,9375 > > Agora nível 1 do fator "p" e nível 2 do fator "n". Então, a=1, b=2; > > p1n2 2,1575 2,2475 > > Tem como automatizar isto em função destes dois parâmetros? > > algo como: > > medias = as.vector(tab[ab, ])???? > > A coluna não interessa, podem ficar tranquilo! Apenas as linhas interessam > de acordo com os parâmetros. > > Se alguém puder me dar uma ajuda ficarei grato! > > Abraços! > > 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/20120715/e65d4cd4/atta… > > > > ------------------------------ > > Message: 2 > Date: Sun, 15 Jul 2012 17:36:31 -0300 > From: Natalia Martins <nsmbarreto(a)gmail.com> > To: Lista R <R-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] duvida controle de qualidade > Message-ID: > < > CAGPGL2FkzWBmCnkL96JZpG9-VQ06EyU-x-3OyanyM34xR+1DBA(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Por favor estou com um dúvida em controle de qualidade e ficaria muito > grata se alguem pudesse tirá-la, ou contribuir de alguma forma. > Bem é o seguinte, tenho uma máquina que é programada para colocar uma > quantidade exata de adubo em um ponto especificado no campo via GPS, ou > seja, a maquina vai passando pelo campo e derramando o adubo na quantidade > exigida para aquela coordenada (x,y), as quantidades vão sendo mudadas a > cada ponto. > Eu sei a quantidade que a maquina *deve* derrubar em todo o campo e as > quantidades que a maquina *realmente derruba.* > Quando a maquina termina o seu serviço um equipamento instalado nela me > informa o que ela realmente derrubou de adubo em cada coordenada, assim eu > tenho que verificar se esta maquina está sob controle estatístico de > qualidade. > - tenho o que a maquina deve derramar por coordenada no campo todo > - tenho o que a maquina realmente derramou por coordenada no campo todo > > Minhas duvidas em relação a esse processo são: > > Para verificar se meu processo esta sob controle eu fiz a diferença entre o > que a maquina realmente aplicou e o que deveria ser aplicado, eu posso > aplicar a teoria de controle de qualidade a esta diferença? > Tenho todos os pontos (2700 pontos) e não uma amostra, a teoria também pode > ser aplicada a todos os dados? > Quando faço o teste de normalidade de shapiro wilks este mostra me que meus > dados não apresentam normalidade, devo transformá-los? > > Muito obrigada, > > -- > > 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 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120715/2a5a2916/atta… > > > > ------------------------------ > > Message: 3 > Date: Sun, 15 Jul 2012 19:50:19 -0300 > From: Eder David Borges da Silva <eder(a)leg.ufpr.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] duvida controle de qualidade > Message-ID: > <CALKkXXqsbB7ACDufTkRL= > mHngB06gtx-196ptJHKHVxpBhsTYw(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Natalia, > Estes dados são legais....eu seguiria uma outro caminho de analise, como a > distribuição é desenvolvida no campo por faixas (A maquina vai passando) e > muito possivelmente em regiões distintas da lavoura por questão de > topografia o erro deve ser diferente, uma analise espacial de seus dados > seria muito legal, pois supostamente o erro não > é distribuído aleatoriamente dentro do talhão, de uma olhada neles por este > ponto de vista podem aparecer informações legais. > Att > > Em 15 de julho de 2012 17:36, Natalia Martins <nsmbarreto(a)gmail.com > >escreveu: > > > Por favor estou com um dúvida em controle de qualidade e ficaria muito > > grata se alguem pudesse tirá-la, ou contribuir de alguma forma. > > Bem é o seguinte, tenho uma máquina que é programada para colocar uma > > quantidade exata de adubo em um ponto especificado no campo via GPS, ou > > seja, a maquina vai passando pelo campo e derramando o adubo na > quantidade > > exigida para aquela coordenada (x,y), as quantidades vão sendo mudadas a > > cada ponto. > > Eu sei a quantidade que a maquina *deve* derrubar em todo o campo e as > > quantidades que a maquina *realmente derruba.* > > Quando a maquina termina o seu serviço um equipamento instalado nela me > > informa o que ela realmente derrubou de adubo em cada coordenada, assim > eu > > tenho que verificar se esta maquina está sob controle estatístico de > > qualidade. > > - tenho o que a maquina deve derramar por coordenada no campo todo > > - tenho o que a maquina realmente derramou por coordenada no campo todo > > > > Minhas duvidas em relação a esse processo são: > > > > Para verificar se meu processo esta sob controle eu fiz a diferença entre > > o que a maquina realmente aplicou e o que deveria ser aplicado, eu posso > > aplicar a teoria de controle de qualidade a esta diferença? > > Tenho todos os pontos (2700 pontos) e não uma amostra, a teoria também > > pode ser aplicada a todos os dados? > > Quando faço o teste de normalidade de shapiro wilks este mostra me que > > meus dados não apresentam normalidade, devo transformá-los? > > > > Muito obrigada, > > > > -- > > > > 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 > > > > > > > > _______________________________________________ > > 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/20120715/7975b619/atta… > > > > ------------------------------ > > Message: 4 > Date: Sun, 15 Jul 2012 21:36:24 -0300 > From: Lucas Barbosa <lucasbcr(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, Ivan Bezerra Allaman > <ivanalaman(a)yahoo.com.br> > Subject: Re: [R-br] Ajuda programação! > Message-ID: > <CAPhSkDFxTn-e8=ZtoKmctr6MqOqYK-ff0cak1= > MVdyENYF9tjg(a)mail.gmail.com> > Content-Type: text/plain; charset=ISO-8859-1 > > Vê aí se é isso mesmo o que vc quer: > > > expand.grid(1:2, 1:2) > Var1 Var2 > 1 1 1 > 2 2 1 > 3 1 2 > 4 2 2 > > apply(expand.grid(1:2, 1:2), 1, function(x){paste('p', x[1], 'n', x[2], > sep='')} ) > [1] "p1n1" "p2n1" "p1n2" "p2n2" > > > Em 15 de julho de 2012 13:07, Ivan Bezerra Allaman > <ivanalaman(a)yahoo.com.br> escreveu: > > Bom dia senhores! > > > > Suponhamos que tem-se dois níveis de um fator "p" e dois níveis de um > fator > > "n", no qual os níveis do fator "p" irei representar pela letra "a" e os > > níveis do fator "n" pela letra "b". Então tem-se a seguinte configuração: > > > > k1 <- c(1.8850,2.1575,2.4775,2.7975) > > k2 <- c(2.9375,2.2475,2.7300,2.9950) > > tab <- matrix(c(k1,k2),ncol=2) > > rownames(tab) <- c('panb','panb,'panb','panb') > > colnames(tab) <- c('k1','k2') > > tab > > k1 k2 > > panb 1,8850 2,9375 > > panb 2,1575 2,2475 > > panb 2,4775 2,7300 > > panb 2,7975 2,9950 > > > > Vamos supor que eu queira os valores do nível 1 do fator "p" e nível 1 do > > fator "n". Então, a=1 e b=1; > > > > p1n1 1,8850 2,9375 > > > > Agora nível 1 do fator "p" e nível 2 do fator "n". Então, a=1, b=2; > > > > p1n2 2,1575 2,2475 > > > > Tem como automatizar isto em função destes dois parâmetros? > > > > algo como: > > > > medias = as.vector(tab[ab, ])???? > > > > A coluna não interessa, podem ficar tranquilo! Apenas as linhas > interessam > > de acordo com os parâmetros. > > > > Se alguém puder me dar uma ajuda ficarei grato! > > > > Abraços! > > > > 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} > > > > _______________________________________________ > > 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. > > > ------------------------------ > > Message: 5 > Date: Mon, 16 Jul 2012 08:33:58 -0300 > From: Natalia Martins <nsmbarreto(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] duvida controle de qualidade > Message-ID: > <CAGPGL2FJvm5SxhU8p4RUJWEAu5v3MphJ-crxoBH8= > y-QXiTFmw(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Oi Eder, bom dia! > muito obrigada pela sugestão, mas uma análise espacial já foi realizada > sim, agora o objetivo era verificar realmente a eficiência da maquina, e > como análises preliminares já fiz um teste de comparação de médias o qual > resultou que a média de aplicação real é diferente da média da aplicação > que deveria ser realizada. > > Att, > > > > > Em 15 de julho de 2012 19:50, Eder David Borges da Silva > <eder(a)leg.ufpr.br>escreveu: > > > Natalia, > > Estes dados são legais....eu seguiria uma outro caminho de analise, como > a > > distribuição é desenvolvida no campo por faixas (A maquina vai passando) > e > > muito possivelmente em regiões distintas da lavoura por questão de > > topografia o erro deve ser diferente, uma analise espacial de seus dados > > seria muito legal, pois supostamente o erro não > > é distribuído aleatoriamente dentro do talhão, de uma olhada neles por > este > > ponto de vista podem aparecer informações legais. > > Att > > > > Em 15 de julho de 2012 17:36, Natalia Martins <nsmbarreto(a)gmail.com > >escreveu: > > > >> Por favor estou com um dúvida em controle de qualidade e ficaria muito > >> grata se alguem pudesse tirá-la, ou contribuir de alguma forma. > >> Bem é o seguinte, tenho uma máquina que é programada para colocar uma > >> quantidade exata de adubo em um ponto especificado no campo via GPS, ou > >> seja, a maquina vai passando pelo campo e derramando o adubo na > quantidade > >> exigida para aquela coordenada (x,y), as quantidades vão sendo mudadas a > >> cada ponto. > >> Eu sei a quantidade que a maquina *deve* derrubar em todo o campo e as > >> quantidades que a maquina *realmente derruba.* > >> Quando a maquina termina o seu serviço um equipamento instalado nela me > >> informa o que ela realmente derrubou de adubo em cada coordenada, assim > eu > >> tenho que verificar se esta maquina está sob controle estatístico de > >> qualidade. > >> - tenho o que a maquina deve derramar por coordenada no campo todo > >> - tenho o que a maquina realmente derramou por coordenada no campo todo > >> > >> Minhas duvidas em relação a esse processo são: > >> > >> Para verificar se meu processo esta sob controle eu fiz a diferença > entre > >> o que a maquina realmente aplicou e o que deveria ser aplicado, eu posso > >> aplicar a teoria de controle de qualidade a esta diferença? > >> Tenho todos os pontos (2700 pontos) e não uma amostra, a teoria também > >> pode ser aplicada a todos os dados? > >> Quando faço o teste de normalidade de shapiro wilks este mostra me que > >> meus dados não apresentam normalidade, devo transformá-los? > >> > >> Muito obrigada, > >> > >> -- > >> > >> 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 > >> > >> > >> > >> _______________________________________________ > >> 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. > > > > > > -- > > 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 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120716/d4006b68/atta… > > > > ------------------------------ > > Message: 6 > Date: Mon, 16 Jul 2012 09:12:52 -0300 > From: henrique jose de paula alves paula alves <jpahenrique(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] duvida controle de qualidade > Message-ID: > < > CAEfjKJTXP8R3BEH_C_WkxuLDieL+qjDZX9_6WxwgBNPTTzuxOA(a)mail.gmail.com> > Content-Type: text/plain; charset=ISO-8859-1 > > escreva pra este professor aqui da ufop, ele pode lhe ajudar melhor. o > cara é fera nesta área. > > > fernandoest(a)ig.com.br > > > Em 16 de julho de 2012 08:33, Natalia Martins <nsmbarreto(a)gmail.com> > escreveu: > > Oi Eder, bom dia! > > muito obrigada pela sugestão, mas uma análise espacial já foi realizada > sim, > > agora o objetivo era verificar realmente a eficiência da maquina, e como > > análises preliminares já fiz um teste de comparação de médias o qual > > resultou que a média de aplicação real é diferente da média da aplicação > que > > deveria ser realizada. > > > > Att, > > > > > > > > > > Em 15 de julho de 2012 19:50, Eder David Borges da Silva < > eder(a)leg.ufpr.br> > > escreveu: > > > >> Natalia, > >> Estes dados são legais....eu seguiria uma outro caminho de analise, > como a > >> distribuição é desenvolvida no campo por faixas (A maquina vai > passando) e > >> muito possivelmente em regiões distintas da lavoura por questão de > >> topografia o erro deve ser diferente, uma analise espacial de seus dados > >> seria muito legal, pois supostamente o erro não é distribuído > aleatoriamente > >> dentro do talhão, de uma olhada neles por este ponto de vista podem > aparecer > >> informações legais. > >> Att > >> > >> Em 15 de julho de 2012 17:36, Natalia Martins <nsmbarreto(a)gmail.com> > >> escreveu: > >>> > >>> Por favor estou com um dúvida em controle de qualidade e ficaria muito > >>> grata se alguem pudesse tirá-la, ou contribuir de alguma forma. > >>> Bem é o seguinte, tenho uma máquina que é programada para colocar uma > >>> quantidade exata de adubo em um ponto especificado no campo via GPS, ou > >>> seja, a maquina vai passando pelo campo e derramando o adubo na > quantidade > >>> exigida para aquela coordenada (x,y), as quantidades vão sendo mudadas > a > >>> cada ponto. > >>> Eu sei a quantidade que a maquina deve derrubar em todo o campo e as > >>> quantidades que a maquina realmente derruba. > >>> Quando a maquina termina o seu serviço um equipamento instalado nela me > >>> informa o que ela realmente derrubou de adubo em cada coordenada, > assim eu > >>> tenho que verificar se esta maquina está sob controle estatístico de > >>> qualidade. > >>> - tenho o que a maquina deve derramar por coordenada no campo todo > >>> - tenho o que a maquina realmente derramou por coordenada no campo todo > >>> > >>> Minhas duvidas em relação a esse processo são: > >>> > >>> Para verificar se meu processo esta sob controle eu fiz a diferença > entre > >>> o que a maquina realmente aplicou e o que deveria ser aplicado, eu > posso > >>> aplicar a teoria de controle de qualidade a esta diferença? > >>> Tenho todos os pontos (2700 pontos) e não uma amostra, a teoria também > >>> pode ser aplicada a todos os dados? > >>> Quando faço o teste de normalidade de shapiro wilks este mostra me que > >>> meus dados não apresentam normalidade, devo transformá-los? > >>> > >>> Muito obrigada, > >>> > >>> -- > >>> > >>> 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 > >>> > >>> > >>> > >>> _______________________________________________ > >>> 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. > > > > > > > > > > -- > > > > 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 > > > > > > > > _______________________________________________ > > 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. > > > > -- > Henrique José de Paula Alves > Graduando em Estatística > UFOP - MG > > > ------------------------------ > > Message: 7 > Date: Mon, 16 Jul 2012 05:20:05 -0700 (PDT) > From: Elton Gean Araújo Araújo <egarauj(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] duvida controle de qualidade > Message-ID: > <1342441205.79483.YahooMailNeo(a)web126004.mail.ne1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Bom dia Natália! > Dá uma olhada no pacote qcc do R, esse pacote realiza analises de controle > de qualidade e pode gerar alguns gráficos para verificar se a máquina > adubadeira está bem regulada ou não. > Elton. > > > ________________________________ > De: henrique jose de paula alves paula alves <jpahenrique(a)gmail.com> > Para: r-br(a)listas.c3sl.ufpr.br > Enviadas: Segunda-feira, 16 de Julho de 2012 9:12 > Assunto: Re: [R-br] duvida controle de qualidade > > escreva pra este professor aqui da ufop, ele pode lhe ajudar melhor. o > cara é fera nesta área. > > > fernandoest(a)ig.com.br > > > Em 16 de julho de 2012 08:33, Natalia Martins <nsmbarreto(a)gmail.com> > escreveu: > > Oi Eder, bom dia! > > muito obrigada pela sugestão, mas uma análise espacial já foi realizada > sim, > > agora o objetivo era verificar realmente a eficiência da maquina, e como > > análises preliminares já fiz um teste de comparação de médias o qual > > resultou que a média de aplicação real é diferente da média da aplicação > que > > deveria ser realizada. > > > > Att, > > > > > > > > > > Em 15 de julho de 2012 19:50, Eder David Borges da Silva < > eder(a)leg.ufpr.br> > > escreveu: > > > >> Natalia, > >> Estes dados são legais....eu seguiria uma outro caminho de analise, > como a > >> distribuição é desenvolvida no campo por faixas (A maquina vai > passando) e > >> muito possivelmente em regiões distintas da lavoura por questão de > >> topografia o erro deve ser diferente, uma analise espacial de seus dados > >> seria muito legal, pois supostamente o erro não é distribuído > aleatoriamente > >> dentro do talhão, de uma olhada neles por este ponto de vista podem > aparecer > >> informações legais. > >> Att > >> > >> Em 15 de julho de 2012 17:36, Natalia Martins <nsmbarreto(a)gmail.com> > >> escreveu: > >>> > >>> Por favor estou com um dúvida em controle de qualidade e ficaria muito > >>> grata se alguem pudesse tirá-la, ou contribuir de alguma forma. > >>> Bem é o seguinte, tenho uma máquina que é programada para colocar uma > >>> quantidade exata de adubo em um ponto especificado no campo via GPS, ou > >>> seja, a maquina vai passando pelo campo e derramando o adubo na > quantidade > >>> exigida para aquela coordenada (x,y), as quantidades vão sendo mudadas > a > >>> cada ponto. > >>> Eu sei a quantidade que a maquina deve derrubar em todo o campo e as > >>> quantidades que a maquina realmente derruba. > >>> Quando a maquina termina o seu serviço um equipamento instalado nela me > >>> informa o que ela realmente derrubou de adubo em cada coordenada, > assim eu > >>> tenho que verificar se esta maquina está sob controle estatístico de > >>> qualidade. > >>> - tenho o que a maquina deve derramar por coordenada no campo todo > >>> - tenho o que a maquina realmente derramou por coordenada no campo todo > >>> > >>> Minhas duvidas em relação a esse processo são: > >>> > >>> Para verificar se meu processo esta sob controle eu fiz a diferença > entre > >>> o que a maquina realmente aplicou e o que deveria ser aplicado, eu > posso > >>> aplicar a teoria de controle de qualidade a esta diferença? > >>> Tenho todos os pontos (2700 pontos) e não uma amostra, a teoria também > >>> pode ser aplicada a todos os dados? > >>> Quando faço o teste de normalidade de shapiro wilks este mostra me que > >>> meus dados não apresentam normalidade, devo transformá-los? > >>> > >>> Muito obrigada, > >>> > >>> -- > >>> > >>> 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 > >>> > >>> > >>> > >>> _______________________________________________ > >>> 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. > > > > > > > > > > -- > > > > 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 > > > > > > > > _______________________________________________ > > 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. > > > > -- > Henrique José de Paula Alves > Graduando em Estatística > UFOP - MG > _______________________________________________ > 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/20120716/2a471d52/atta… > > > > ------------------------------ > > Message: 8 > Date: Mon, 16 Jul 2012 10:48:16 -0300 > From: Walmes Zeviani <walmeszeviani(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Ajuda programação! > Message-ID: > <CAFU=EkY3OHAz8KJCLv3QZ+0ZSRd2= > Ak3KGJNSHKFVDoi7beWHQ(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Coloca as colunas que identificam os fatores junto com as respostas num > data.frame() e usa tapply(), aggregate(), by(), etc para obter resumos por > grupos. Será que entendi a pergunta? > > À 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 > 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/20120716/b10f32cd/atta… > > > > ------------------------------ > > Message: 9 > Date: Mon, 16 Jul 2012 07:26:43 -0700 (PDT) > From: "Elias T. Krainski" <eliaskrainski(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Ajuda para Simulação de variáveis binárias > correlacionadas no R > Message-ID: > <1342448803.39401.YahooMailNeo(a)web163403.mail.gq1.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > Olá Gilenio, seria algo assim: > > ### simula dados > > n <- 1000; beta <- c(-3,1) > y <- integer(n) > y[1] <- rbinom(1,1,.5) > for (i in 2:n) > y[i] <- rbinom(1, 1, 1/(1+exp(-beta[1]-beta[2]*y[i-1]))) > > ### estima modelo > > m <- glm(y ~ y2, family=binomial, data=data.frame(y=y[2:n], y2=y[2:n-1])) > summary(m) > > > Elias T. Krainski > > > >________________________________ > > De: Gilenio Borges Fernandes <gilenio(a)ufba.br> > >Para: r-br(a)listas.c3sl.ufpr.br > >Enviadas: Sábado, 14 de Julho de 2012 20:52 > >Assunto: [R-br] Ajuda para Simulação de variáveis binárias > correlacionadas no R > > > > > >Simulação de variáveis binárias > correlacionadas no R é bem popular, usando os pacotes bindata, bild, > binarySimCLF, > BinNor e binGroup. Mas eu estou precisando simular variáveis binárias > correlacionadas num modelo auto regressivo de primeira ordem de regressão > logística, > a partir dos parâmetros do preditor linear, incluindo o referente à > variável dependente > defasada (Yi-1). Alguém poderia me indicar algum pacote ou função no > R para este fim? > >Desde já, agradeço-- > >Gilenio Borges Fernandes > >Universidade Federal da Bahia > >Instituto de Matemática > >Departamento de Estatística > >Av. Adhemar de Barros, s/n ? Ondina. > >40.170-110 - Salvador - BA, Brasil > >Tel.: (071)3283-6280/6336 Fax: (071)3283-6276 > >URL: http://lattes.cnpq.br/6764860618464860 > > > > > >_______________________________________________ > >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/20120716/0532d17a/atta… > > > > ------------------------------ > > Message: 10 > Date: Mon, 16 Jul 2012 11:45:41 -0300 > From: edlayne moraes <moraes.epgeo(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Duvidas em análise cluster > Message-ID: > <CAPbDHAgFkwcrt= > BzyHsgeQRhwYgBWOaWBG-aUzjGgqYbey1-Tg(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Boa tarde, > Inseri uma planilha com os dados em .csv, coloquei na linha 1º o nome dos > meus dados. Ao importar a tabela eu mandei ler essa linha como rótulo das > variáveis. > Para fazer o tratamento eu instalei o pacote "vegan" > e ao mandar plotar o gráfico o programa informa um erro de "rowsums (x, > na.rm=TRUE): x deve ser numérico". > Alguém sabe me dizer o que há de erro e se existe outra maneira de fazer o > cluster no programa? > Desde já agradeço a atenção > Att, > > -- > > Edlayne Pimentel de Moraes**** > > Graduanda em Geologia-UFES**** > (28) 81117-0585//(28) 3552-0194 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120716/1476e8f7/atta… > > > > ------------------------------ > > Message: 11 > Date: Mon, 16 Jul 2012 11:56:22 -0300 > From: Natalia Martins <nsmbarreto(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Duvidas em análise cluster > Message-ID: > <CAGPGL2EzM27XQ-EpcfKRvT+YzuxX+PD=Ena+Cob= > 8jHSJMJApg(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Oie, tenta na hora de ler os dados colocar o comando head=T, assim ele lerá > a 1 linha como o nome das suas variaveis, ai na 1 coluna só haverá valores > numericos. > > Att, > Em 16 de julho de 2012 11:45, edlayne moraes <moraes.epgeo(a)gmail.com > >escreveu: > > > Boa tarde, > > Inseri uma planilha com os dados em .csv, coloquei na linha 1º o nome dos > > meus dados. Ao importar a tabela eu mandei ler essa linha como rótulo das > > variáveis. > > Para fazer o tratamento eu instalei o pacote "vegan" > > e ao mandar plotar o gráfico o programa informa um erro de "rowsums (x, > > na.rm=TRUE): x deve ser numérico". > > Alguém sabe me dizer o que há de erro e se existe outra maneira de fazer > o > > cluster no programa? > > Desde já agradeço a atenção > > Att, > > > > -- > > > > Edlayne Pimentel de Moraes**** > > > > Graduanda em Geologia-UFES**** > > (28) 81117-0585//(28) 3552-0194 > > > > > > _______________________________________________ > > 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. > > > > > > -- > > 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 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120716/429ebf6b/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 17, assunto 14 > ***************************************** > -- Gilenio Borges Fernandes Universidade Federal da Bahia Instituto de Matemática Departamento de Estatística Av. Adhemar de Barros, s/n – Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6280/6336 Fax: (071)3283-6276 URL: http://lattes.cnpq.br/6764860618464860

1 0

Ajuda para Simulação de variáveis binárias correlacionadas no R
by Gilenio Borges Fernandes 16 Jul '12

16 Jul '12

Simulação de variáveis binárias correlacionadas no R é bem popular, usando os pacotes bindata, bild, binarySimCLF, BinNor e binGroup. Mas eu estou precisando simular variáveis binárias correlacionadas num modelo auto regressivo de primeira ordem de regressão logística, a partir dos parâmetros do preditor linear, incluindo o referente à variável dependente defasada (Yi-1). Alguém poderia me indicar algum pacote ou função no R para este fim? Desde já, agradeço -- Gilenio Borges Fernandes Universidade Federal da Bahia Instituto de Matemática Departamento de Estatística Av. Adhemar de Barros, s/n – Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6280/6336 Fax: (071)3283-6276 URL: http://lattes.cnpq.br/6764860618464860

3 2

Colocar elipse de confiança num gráfico
by Mauro Sznelwar 14 Jul '12

14 Jul '12

Como eu faço para inserir uma elipse de confiança neste gráfico acima de Análise discriminante, não sei como usar o comando: plot(marajoara.norm.ldakm1,dimen=2,ylim=c(-5,5),xlim=c(20,35), cex=1,col=marajoara.norm.km1$cluster) library(ellipse) plot(ellipse(S, centre = , t = ), type = "l") Abaixo tem o script que gerou este gráfico, e estou enviando também o arquivo de dados. http://www.datafilehost.com/download-1eccadf2.html marajoara=read.csv('Marajoara Mario.csv',head=T,dec=',',sep=';') marajoara.norm=sapply(data.frame(marajoara[,2:15]),scale) #marajoara.norm=(data.frame(marajoara[,2:15])) #summary(marajoara.norm) #library(fBasics) #basicStats(marajoara.norm) marajoara.norm.h=hclust(dist(marajoara.norm),method='average') #marajoara.norm.h=hclust(dist(marajoara.norm)^2,method='ward') plot(marajoara.norm.h,main='Marajoara original-Ward Linkage\n Square Euclidian Distance\n Hclust Dendograma') #pdf('figuras/Marajoara-original-dendograma-hclust-average.pdf') library(cluster) marajoara.norm.d=diana(dist(marajoara.norm)) plot(marajoara.norm.d,main='Marajoara original-Divisive Dendograma') #pdf('figuras/Marajoara-original-dendograma-divisive.pdf') #Componentes Principais initial=tapply(marajoara.norm,list(rep(cutree(marajoara.norm.h,3), ncol(marajoara.norm)),col(marajoara.norm)),mean) marajoara.norm.km1=kmeans(marajoara.norm,initial) marajoara.norm.pca=princomp(marajoara.norm) marajoara.norm.px=predict(marajoara.norm.pca) library(MASS) eqscplot(marajoara.norm.px[,1:2],type='n',xlab='component 1',ylab='component 2') text(marajoara.norm.px[,1:2],labels=as.character(marajoara.norm.km1$cluster), col=as.numeric(marajoara.norm.km1$cluster)) title('PCA Marajoara original - clusters kmeans\n centros iniciais corte dendograma hclust 3') #identify(marajoara.norm.px,n=5) #pdf('figuras/PCA-Marajoara-original-centros-kmeans\n centros inicias de corte dendograma hclust 3') initial=tapply(marajoara.norm,list(rep(cutree(marajoara.norm.d,3), ncol(marajoara.norm)),col(marajoara.norm)),mean) marajoara.norm.km2=kmeans(marajoara.norm,initial) marajoara.norm.pca=princomp(marajoara.norm) marajoara.norm.px=predict(marajoara.norm.pca) eqscplot(marajoara.norm.px[,1:2],type='n',xlab='component 1',ylab='component 2') text(marajoara.norm.px[,1:2],labels=as.character(marajoara.norm.km2$cluster), col=as.numeric(marajoara.norm.km2$cluster)) title('PCA Marajoara original - clusters kmeans\n centros iniciais corte dendograma diana 3') #identify(marajoara.norm.px,n=5) #pdf('figuras/PCA-Marajoara-original-centros-kmeans\n centros inicias de corte dendodiana-3.pdf') marajoara.norm.km=kmeans(marajoara.norm,3) marajoara.norm.pca=princomp(marajoara.norm) marajoara.norm.px=predict(marajoara.norm.pca) eqscplot(marajoara.norm.px[,1:2],type='n',xlab='component 1',ylab='component 2') text(marajoara.norm.px[,1:2],labels=as.character(marajoara.norm.km$cluster), col=marajoara.norm.km$cluster) title('PCA Marajoara original - clusters kmeans 3 centros') #identify(marajoara.norm.px,n=5) #pdf('figuras/PCA-Marajoara-original-clusters-kmeans-3-centros.pdf') #Análise Discriminante Linear marajoara.norm.ldakm1=lda(marajoara.norm,as.character(marajoara.norm.km1$cluster)) plot(marajoara.norm.ldakm1,dimen=2,cex=1,col=marajoara.norm.km1$cluster) title('LDA Marajoara-Grupo do K-Means - corte dendograma hclust em 3') #pdf('figuras/LDA-Marajoara-grupos-kmeans-dendohclust-3.pdf') #fazendo um zoom no gráfico corta um ponto plot(marajoara.norm.ldakm1,dimen=2,ylim=c(-5,5),xlim=c(20,35), cex=1,col=marajoara.norm.km1$cluster) title('LDA Marajoara-parte do gráfico para visualização\n Grupos do K-Means-corte dendograma hclust em 3') #identify(marajoara.norm.km1$cluster,n=5) #pdf('figuras/LDA-Marajoara-grupos-kmeans-dendodiana-3.pdf') library(ellipse) plot(ellipse(S, centre = , t = ), type = "l")

1 0

Dúvida sobre o ANOSIM
by Débora Rodrigues de Souza 13 Jul '12

13 Jul '12

Boa tarde a todos, Acabei de encontrar essa página de discussões, e achei muito interessante a ideia, parabéns aos organizadores. Bom, gostaria de auxílio para executar uma análise de similaridade do tipo ANOSIM. Sei que o teste é feito a posteriori o NMDS, ou alguma análise de ordenação, mas tenho dúvida quando ao script do mesmo. Já tentei seguindo os passos indicados no manul do programa R, mas fiquei em dúvida em qual comando devo alterar, para então inserir os meus dados. Estou mexendo nesse programa a pouco tempo, por isso a dificuldade. Ao que me parece, devo entrar com dois tipos de planilhas. Mas mesmo colando meus dados, o resultado que aparece é o da memória do programa. Como alterar isso??? Esse é o script que vejo no programa, data(dune) data(dune.env) dune.dist <- vegdist(dune) attach(dune.env) dune.ano <- anosim(dune.dist, Management) summary(dune.ano) plot(dune.ano) Alguém poderia me ajudar??? Desde já agradeço a atenção, e aguardo resposta. Att, Débora

2 5

Importar tabela com células vazias
by Luciano F. Sgarbi 13 Jul '12

13 Jul '12

Boa tarde, Estou com um problema para importar uma tabela em formato ".txt" com várias células vazias. Para tal, estou usando o comando "read.table("name.file.txt",h=F,fill=T)". O que está acontecendo é que quando eu importo, o arquivo só é importado com as duas primeiras colunas. Alguém poe me ajudar? A tabela está no formato do exemplo abaixo: 1 125 2 151 3 -59,8184 4 -14,009 5 0,00167 6 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 Abraços, Luciano F. Sgarbi -- Luciano F. Sgarbi Mestrando em Ecologia e Evolução - UFG Laboratório de Ecologia de Insetos Cel. (62)8174-2262 Lab. (62)3521-1732

3 2

Acessar banco de dados
by Emerson Cotta Bodevan 12 Jul '12

12 Jul '12

Prezados membros da lista. Tenho um arquivo ".MDB" gerado pelo Epi Info. Estou tentando acessá-lo no R. Com o comando abaixo eu consegui fazer uma conexao com o arquivo. automed<-odbcConnectAccess("Estudo de automedicação em crianças.mdb") O que eu gostaria agoraera poder acessar as variáveis individualmente (como se fosse um data frame), para que eu possa analisar as variáveis separadamente. Alguém sabe como fazê-lo? Estou com o pacote RODBC instalado. Agradeço qualquer ajuda, Emerson. -- "A educação não muda o mundo. A educação muda as pessoas. As pessoas mudam o mundo." * Paulo Freire.*

3 2

Séries temporais [OFF-TOPIC]
by ASANTOS 12 Jul '12

12 Jul '12

Bom dia pessoal, Alguém tem conhecimento de algum livro que tenha em seu conteúdo algo sobre séries temporais inflacionadas com zeros usando o R. Procurei no Times series analysis do Shumway e Stoffer que tenho como referência e não encontrei nada sobre abordagens inflacionada com zeros. Obrigado, -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

2 1

Re: [R-br] Mapa Manaus por zonas e bairros:acentuação e legendas
by Benilton Carvalho 11 Jul '12

11 Jul '12

Copiando a lista novamente.... plot(mao_bz, col = mao_bz$ZONA) legend('bottomright', unique(as.character(mao_bz$ZONA)), fill=unique(mao_bz$ZONA), bty='n') b 2012/7/12 Edson Lira <edinhoestat(a)yahoo.com.br>: > Valeu Benilton!!! > E as legendas? Como faço? São 6 zonas e em torno de 60 bairros. > > Edson Lira > Estatístico > Manaus-Amazonas > ________________________________ > De: Benilton Carvalho <beniltoncarvalho(a)gmail.com> > Para: Edson Lira <edinhoestat(a)yahoo.com.br> > Enviadas: Quarta-feira, 11 de Julho de 2012 20:54 > > Assunto: Re: [R-br] Mapa Manaus por zonas e bairros:acentuação e legendas > > o encoding original dos dados a q vc se refere e' ISO-8859-1. > > use > > toUTF8(mao_bz$NOME, 'ISO-8859-1') > > ou qq um dos compativeis (CP819, CP822... talvez ate' o 1252) > > b > > 2012/7/12 Edson Lira <edinhoestat(a)yahoo.com.br>: >> Link para o shx >> >> http://www.datafilehost.com/download-3ef3a7f9.html >> >> >> sessionInfo() >> R version 2.14.2 (2012-02-29) >> Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit) >> >> locale: >> [1] pt_BR.UTF-8/pt_BR.UTF-8/pt_BR.UTF-8/C/pt_BR.UTF-8/pt_BR.UTF-8 >> >> attached base packages: >> [1] stats graphics grDevices utils datasets methods base >> >> other attached packages: >> [1] classInt_0.1-17 e1071_1.6 class_7.3-3 descr_0.9.7 >> shapefiles_0.6 maptools_0.8-14 lattice_0.20-0 sp_0.9-98 >> foreign_0.8-49 >> >> loaded via a namespace (and not attached): >> [1] grid_2.14.2 tools_2.14.2 xtable_1.7-0 >> >> >> [ ]'s. >> >> Edson Lira >> Estatístico >> Manaus-Amazonas >> ________________________________ >> De: Benilton Carvalho <beniltoncarvalho(a)gmail.com> >> Para: r-br(a)listas.c3sl.ufpr.br; Edson Lira <edinhoestat(a)yahoo.com.br> >> Enviadas: Quarta-feira, 11 de Julho de 2012 20:25 >> Assunto: Re: [R-br] Mapa Manaus por zonas e bairros:acentuação e legendas >> >> e tenho a impressao de q faltou o o link pro shx? >> >> > >

1 0

Mapa Manaus por zonas e bairros:acentuação e legendas
by Edson Lira 11 Jul '12

11 Jul '12

Caros amigos, estou gerando o mapa de Manaus por zona, minha máquina é um Mac i7, 4gb ram. Estou com o problema da acentuação, a linha: mao_bz$NOME <- toUTF8(mao_bz$NOME, "IBM850") me dá a saida abixo: > mao_bz$NOME [1] PONTA NEGRA [2] LĂ˘Â\u0088Â\u009aĂ\u0083ÂŽĂ˘Â\u0088Â\u009aĂ˘Â\u0089Â¤Ă˘Â\u0088Â\u009aĂ\u0083Â˘RIO DO VALE [3] COL TERRA NOVA [4] TARUMĂ˘Â\u0088Â\u009aĂ\u0083ÂŽĂ˘Â\u0088Â\u009aĂ˘Â\u0088Â\u0082Ă\u0082ÂŹĂ\u0082ÂŁ [5] SANTA ETELVINA [6] MONTE DAS OLIVEIRAS [7] FLORES [8] CIDADE NOVA [9] NOVO ISRAEL [10] COL STO ANTĂ˘Â\u0080Â\u009aĂ\u0083ÂŽĂ\u0083ÂşĂ˘Â\u0088Â\u009aĂ˘Â\u0084Â˘NIO [11] REDENĂ˘Â\u0080Â\u009aĂ\u0083ÂŽĂ\u0083ÂşĂ˘Â\u0088Â\u009aĂ\u0082Â˘Ă˘Â\u0088Â\u009aĂ\u0083ÂŽĂ˘Â\u0088Â\u009aĂ˘Â\u0088Â\u0082Ă\u0082ÂŹĂ\u0082ÂŁO . . . O script que estou usando é o abaixo: mao_bz<-readShapePoly("/Users/edsonlira/Edinho/Edinho/Bases_MAO/Setores_MAO/Bairros_e_Zonas.shp") summary(mao_bz) plot(mao_bz) library(descr) mao_bz$NOME <- toUTF8(mao_bz$NOME, "IBM850") mao_bz$NOME # palette(c("#ccddff", "#ffddcc", "#ccffcc", "#ffffcc", "#ccffff")) palette(rainbow(6)) plot(mao_bz, col = mao_bz$ZONA) title("Mapa de Manaus por bairros e zonas geográficas") Os links abaixo tem as base de dados usadas. http://www.datafilehost.com/download-00405c50.html #bairros e zonas.shp http://www.datafilehost.com/download-c36f359a.html #bairros e zonas.dbf http://www.datafilehost.com/delete-a54c1f16-16120426.html #bairros e zonas shx Outra dúvida: Como faço para criar uma legenda com os nomes das zonas e/ou com o nome dos bairros? [ ]'s. Edson Lira Estatístico Manaus-Amazonas

2 3

Como escrever o modelode Michaelis-Menten na linguagem BUGS - R2WinBUGS
by Augusto Ribas 11 Jul '12

11 Jul '12

Opa. Alguem usa o winbugs com o R? Eu estava tentando ajustar um Modelo de michelis-menten e não consegui. Se alguem puder dar uma olhadinha e ver onde estou errado. Segue um exemplo do que estou fazendo: library(R2WinBUGS) ################################################################# #Relação Riqueza~estudos ################################################################# #Modelo do que eu acho que acontece #riqueza~(c*estudos)/(d+estudos) #Fazendo um levantamento do numero de especies de parasitas que parasitam algumas especies de sapos, eu supuz que qd uma especie tinha poucos estudos #tava sempre bem subestimado a riqueza de parasitas associadas aquela especie. #dai como a unica informação que fico é com o número de estudos, pensei que deveria ter uma relação onde qt mais estudos, mais especies eu devo encontrar #até um limite, o limite é o parametro c nesse caso, e d representa como os estudos vão contribuindo pra eu chegar ao valor limite de especies. #existem alguns problemas acredito, como especies podem ser mais faceis ou mais dificies de estudar, e o limite pra cada especie depende de mais coisas e tipos de estudos. #mas ficando so nesse caso, eu fiz esse exemplo set.seed(13) #gerando dados exemplo estudos<- rpois(48,5) #parametros originais #então o limite de especies de parasitos seria 10 c<-10 #e eu preciso de mais ou menos 5 estudos pra ver legal a riqueza de especies, d<-2.5 #resposta riqueza<-round(((c*estudos)/(d+estudos))+rnorm(length(estudos),0,1),0) plot(riqueza~estudos,xlim=c(0,15),ylim=c(0,15),xlab="Riqueza de espécies",ylab="Número de estudos") curve(c*x/(d+x),0,15,add=T,col="black") #ajustando com modelo não linear modelo.estudos<-nls(riqueza~((a*estudos)/(b+estudos)),start=list(a=8,b=1.5)) summary(modelo.estudos) points(predict(modelo.estudos)[order(predict(modelo.estudos))]~estudos[order(estudos)],type="l",col="red",lty=2) #eu li que é possivel lineariza essa função e ajusta com glm da seguinte forma: #ajustando com glm modelo.menten<- glm(riqueza~c(1/estudos),family=gaussian(link = "inverse")) summary(modelo.menten) a1<-1/modelo.menten$coef[1] b1<-a1*modelo.menten$coef[2] a1 b1 curve(a1*x/(b1+x),0,15,add=T,col="blue",lty=2) ### Ajuste com winbugs # Definir modelo sink("michaelis.menten.test.txt") cat(" model { # Priors alpha ~ dnorm(0,0.001) beta ~ dnorm(0,0.001) sigma ~ dunif(0, 100) # Likelihood for (i in 1:n) { riqueza[i] ~ dnorm(mu[i], tau) mu[i] <- alpha + beta * estudos[i] } # Derived quantities tau <- 1/ (sigma * sigma) } ",fill=TRUE) sink() # Juntar os dados - Note que eu não entendi como fazer a função de ligação dentro do modelo winbugs #dai eu registrei os dados como 1/riqueza e 1/estudos aqui, acho que isso que o pepino deve ta nisso win.data <- list(riqueza=c(1/riqueza),estudos=c(1/estudos),n=length(estudos)) # Valores iniciais inits <- function(){ list(alpha=rlnorm(1), beta=rlnorm(1),sigma = rlnorm(1))} # Paramstros a estimar params <- c("alpha","beta","tau") # MCMC settings nc <- 3 ni <- 3000 nb <- 1000 nt <- 2 # Começa o winbugs out <- bugs(data=win.data, inits=inits, parameters.to.save=params, model.file="michaelis.menten.test.txt", n.thin=nt, n.chains=nc, n.burnin=nb, n.iter=ni, debug = TRUE) ###Olhar o modelo print(out, dig = 3) #retornoar os parametros como no glm acima a2<-1/out$summary[1,1] b2<-a*out$summary[2,1] a2 b2 curve(a2*x/(b2+x),0,15,add=T,col="green",lty=2) legend("topright",legend=c("Original","Ajustado com nls","Ajustado com glm","Ajustado com WinBUGS"), lty=c(1,2,2,2),col=c("black","red","blue","green")) #Minhas duvidas são: #01. como montar o modelo na linguagem bugs pra fazer esse tipo de analise? #procurando pela internet eu não achei nenhum exemplo pra esse caso especifico. #02. como eu colocaria dentro do modelo esse função de ligação inversa. #eu tentei usar a função pow(mu[i],-1) mas não rolou, #03. e o que estou fazendo no modelo winBUGS tem muito problemas? entrar com os dados daquele jeito, ao invez de especificar no modelo? #não consigo imaginar a implicação de não definir no modelo a função de ligação. Bom se alguem tiver esperiencia com winbugs e puder dar uma olhadinha no que estou fazendo errado Ou alguma dica se o que estou fazendo tem algum problema critico -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

1 0

Integração de C com R
by Jodavid Ferreira 10 Jul '12

10 Jul '12

Alguém conhece alguma apostila sobre integração de C com R? estou precisando estudar sobre isso, estou procurando nos sites de buscas, mas não estou encontrando resultados muito satisfatórios! -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br 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] Como gerar dados de um experimento
by Walmes Zeviani 10 Jul '12

10 Jul '12

CAR e SAR já não são minha praia mais. Mas você pode usar a função spdep::cell2nb() para obter essa matriz de proximidade. À 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 twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

1 0

help para constrastes
by Clenilso Sehnen Mota 10 Jul '12

10 Jul '12

Prezados (as) Estou analisando os resultado de um trabalho. Minha intensão é verificar a comparação conforme a tabela 2. A tabela 1 mostra os tratamentos e as combinações dos 13 componentes. Já discuti pessoalmente com algumas pessoas, se a forma com estimei esses contrastes esta correta ou não. E não chegamos a um consenso. Gostaria de saber a resposta de você se esta correto ou não. Desde já agradeço pela ajuda. Att Clenilso Tabela 1- Composição dos tratamentos: Trat Combinações 1 A 2 B 3 B+C (3:1) 4 B+C (1:1) 5 B+C (1:3) 6 B+D+E(3:1:6) 7 B+D+C(2:1:2) 8 B+D+C (1:1:1) 9 F 10 G+C(7:3) 11 D+C(1:1) 12 D+H (1:1) 13 I 14 D+J+K(2:2:1) 15 D+K(3:1) 16 L Tabela 2 Composição dos contrastes. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 Contraste 1 (A e L vs Demais) -7.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 -7.0 Contraste 2 (C vs s/C) 0.0 1.0 -1.0 -1.0 -1.0 1.0 -1.0 -1.0 1.0 -1.0 -1.0 1.0 1.0 1.0 1.0 0.0 Contraste 3 (B vs s/B) 0.0 -1.0 -1.0 -1.0 -1.0 -1.0 -1.0 -1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 0.0 Contraste 4 (E, F, H e K vs s/E, F, H e K) 0.0 1.0 1.0 1.0 1.0 -1.8 1.0 1.0 -1.8 1.0 1.0 -1.8 1.0 -1.8 -1.8 0.0 Contraste 5 (D vs s/D) 0.0 1.0 1.0 1.0 1.0 -1.0 -1.0 -1.0 1.0 1.0 -1.0 -1.0 1.0 -1.0 -1.0 0.0 Clenilso Sehnen Mota Eng. Agronômo - Udesc DSc Fisiologia Vegetal - UFV MSc produção Vegetal Udesc/Embrapa E-mail: <mailto:csm.sehnen@gmail.com> csm.sehnen(a)gmail.com Celular: (64) 8150 8808

2 2

[Dúvida] Linguagem R e Julia.
by Pedro Rafael 08 Jul '12

08 Jul '12

Pessoal, está existindo esforços para tornar a linguagem R mais eficiente no que diz respeito eu desempenho computacional de determinadas tarefas, como por exemplo nos laços que em R são mais demorados que algumas outras linguagens interpretadas como Lua, Python, JavaScript, etc? Atualmente vejo a linguagem Julia se promovendo com comparativos de tempo de processamento com outras linguagens como Matlab, Python, Octave e inclusive o R. Acho que daqui a uns 5 ou 10 anos poderemos utilizar essa linguagem como uma boa ferramenta para computação científica geral. Contudo, acredito que ela irá precisar de muitos outros fatores para que seja uma alternativa para estatistica como o R. Posso estar falando cedo demais, mas acredito que Julia não só vislumbra a estatística. Vejo que será uma linguagem alternativa compilada para Scilab, Matlab, Octave e demais outros softwares de computação numérica. Dessa forma, se esse cenário permanecer vejo que o R ainda será o ícone de software livre para estatística. Sem falar na existência de um local (CRAN) que reúne boa parte das contribuições de colaboradores de toda parte do mundo que é de extrema importância para evolução de um software livre. [ ], Pedro Rafael Diniz Marinho.

1 0

Re: [R-br] funcao histogram
by Walmes Zeviani 05 Jul '12

05 Jul '12

Uma sugestão seria my.densityplot <- function(x, cols, ...){ panel.densityplot(x, ...) dx <- density(x) breaks <- sort(c(range(dx$x), quantile(x, c(1,3)/4))) fx <- approxfun(dx$x, dx$y) do.polygon <- function(x){ y <- fx(x) return(list(x=c(min(x), x, max(x)), y=c(0, y, 0))) } seqs <- lapply(1:(length(breaks)-1), function(i){ x <- seq(breaks[i], breaks[i+1], l=30) do.polygon(x) }) for(i in 1:length(seqs)){ seqs[[i]]$col <- cols[i] } lapply(seqs, function(i) do.call(panel.polygon, i)) panel.mathdensity(dmath = dnorm, col = "black", args = list(mean=mean(x),sd=sd(x)), lty=3) } require(RColorBrewer) cols <- brewer.pal(3, "Greens") densityplot(~height|voice.part, data=singer, cols=cols, panel=my.densityplot) À 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 twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== 2012/7/5 Samuel Carvalho <samukajm(a)yahoo.com.br> > Falae Walmes, > Cara é basicamente isso mesmo que queria quando voce gera o ultimo gráfico > com as hachuras nos diferentes quantis. > Mas voce teria alguma dica de como implementar isso dentro da função > xyplot? Talvez a saída seja pelo panel.function() e depois pelo > panel.histogram() > Consegue dar um help? > Valeu > > *====================================* > *Samuel P. C. Carvalho > *Mestre em Ciências Florestais [UFLA] > Doutorando em Recursos Florestais [ESALQ/USP] > ============================================= > >

1 0

Help please
by Danielle Steffen 05 Jul '12

05 Jul '12

Ola, Estou batendo cabeca com o R. Estou com um dos arquivos txt da amostra do censo 2010, e apesar de saber importar com mascara de leitura em SPSS e SAS, estou me sentindo um pre-escolar no R. como faco importacao com mascara de leitura? sabe, aquela coisa de colocar caracter 1 a 4 - var 1 caracter 5 a 7 - var 2 e assim por diante. Bati cabeca com a funcao scan, mas como ela nao me ama e sao quase 3h da manha por aqui, ela nao quis importar, a sintaxe foi ridicula, tentei enfiar a mascara numa lista dentro do what, com os comprimentos de cada campo do arquivo: teste <- scan("C:/Users/Public/novo case/amostra.txt", what=list( length(2),length(5),length(13),length(8),length(3),length(1),length(2),length(3),length(2),length(1),length(2),length(2),length(1),length(3),length(3),length(2),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(4),length(1),length(7),length(7),length(3),length(3),length(1),length(7),length(7),length(7),length(1),length(7),length(7),length(7),length(1),length(1),length(2),length(2),length(1),length(1),length(2),length(1),length(1),length(1),length(3),length(3),length(3),length(1),length(7),length(7),length(7),length(1),length(2),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(4),length(5),length(1),length(1),length(1),length(1),length(6),length(6),length(4),length(1),length(6),length(4),length(7),length(4),length(7),length(4),length(7),length(5),length(6),length(4),length(3),length(1),length(1),length(1),length(1),length(1),length(1),length(6),length(1),length(7),length (7),length(7),length(1),length(1),length(1),length(2),length(2),length(2),length(1),length(2),length(2),length(2),length(1),length(3),length(1),length(1),length(1),length(2),length(4),length(1),length(2),length(2),length(2),length(2),length(1),length(2),length(1),length(1),length(1),length(1),length(1),length(3),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),len gth(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1),length(1) )) Mas esta obviamente errado e claro, nao sei pq. Alguem me ajuda? Nao precisa ser essa funcao, mas tentei achar uma que fosse a mais flexivel possivel, pq as outras que uso (read.table, read.csv, etc) elas supoe que o arquivo tem algum tipo de separador, nem que seja um espacinho qualquer pra quebrar o galho. Agradeco qq ajuda, e nao me zoem, estou morrendo de sono e odio do R. Abs, Danielle

4 3

Quasi AIC
by ASANTOS 05 Jul '12

05 Jul '12

Bom dia pessoal, Estava tentando comparar modelos baseados em distribuição quasi Poisson com Poisson, dentre outros critérios, também realizar a comparação através de Quasi AIC com AIC (Não sei a aceitabilidade dos senhores quanto ao uso do quasi AIC), porém dentre os adeptos, gostaria de saber porque a função não consegue calcular o quasi AIC quando não há sobredispersão, como em: > QAIC(m.2,chat=deviance(m.2)/df.residual(m.2))## Calculo do quasi AIC [1] NA Mensagens de aviso perdidas: In QAIC(m.2, chat = deviance(m.2)/df.residual(m.2)) : 'chat' given is < 1, increased to 1 Segue CRM completo: require(MuMIn) ##Pacote que calcula quasi AIC set.seed(122) y<- c(mapply(rpois, lambda=c(5,20), MoreArgs=list(n=100)))##Criação da variável resposta trat<- gl(2,100) ##Criação dos tratamentos m.1<- glm(y~trat,family="poisson")### Distribuição Poisson summary(m.1) m.2<- glm(y~trat,family="quasipoisson")### Distribuição quasi-Poisson summary(m.2) QAIC(m.2,chat=deviance(m.2)/df.residual(m.2))## Calculo do quasi AIC # Obrigado, -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

3 3

Re: [R-br] funcao histogram
by Walmes Zeviani 04 Jul '12

04 Jul '12

Samuel, tente adaptar isso aqui set.seed(123) x <- rnorm(1000, 5, 2) plot(density(x)) breaks <- sort(c(range(x), quantile(x, c(1,3)/4))) dx <- density(x) fx <- approxfun(dx$x, dx$y) fx(0.5) do.polygon <- function(x){ y <- fx(x) return(list(x=c(min(x), x, max(x)), y=c(0, y, 0))) } seqs <- lapply(1:(length(breaks)-1), function(i){ x <- seq(breaks[i], breaks[i+1], l=30) do.polygon(x) }) str(seqs) plot(density(x)) polygon(seqs[[1]], density=2) seqs[[1]]$density <- 3 seqs[[2]]$density <- 5 seqs[[3]]$density <- 9 plot(density(x)) lapply(seqs, function(i) do.call(polygon, i)) À disposição. Walmes.

1 0

funcao histogram
by Samuel Carvalho 04 Jul '12

04 Jul '12

Bom dia pessoal, Estou precisando de uma ajuda para gerar um histograma de frequencias para dados agrupados e estou reproduzindo o exemplo da função ?histogram A seguir o #RCMR require(lattice) ?histogram histogram( ~ height | voice.part, prob=T, data = singer, xlab = "Height (inches)", type = "density", panel = function(x, ...) { panel.histogram(x, ...) panel.mathdensity(dmath = dnorm, col = "black", args = list(mean=mean(x),sd=sd(x))) panel.abline(v=qnorm(0.25, mean=mean(x), sd=sd(x)), lty=2, col='blue') panel.abline(v=qnorm(0.75, mean=mean(x), sd=sd(x)), lty=2, col='blue') } ) Fiz uma adaptação no exemplo e inseri as duas ultimas linhas para que pudesse ser retornado uma linha para os quantis da normal que me retornam 25% e 75% dos dados porem queria que ao invés de linhas eu possa fazer hachuras divindo a distribuição em tres pontos, sendo que fossem feitas tres hachuras, ou seja, uma primeira para a area abaixo de 0.25, uma segunda entre 0.25 e 0.75 e uma terceira para a area acima de 0.75 Espero que tenha conseguido expressar o que preciso. Obrigado a todos Samuel ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] =============================================

2 1

Ler dados do excell
by Edson Lira 03 Jul '12

03 Jul '12

Boa noite R-Br lista, estou lendo um arquivo em excell com a rotina abaixo: require(xlsx) dat<-read.xlsx("/folder/dbp.xlsx",1,encoding="macroman",row.names=1) names(dat) Mostra a saida abaixo > names(dat) [1] "TRT" "DBP1" "DBP2" "DBP3" "DBP4" "DBP5" "Age" "Sex" Falta a variável da primeira coluna que é "Subject" O arquivo em excell está no link abaixo. http://www.datafilehost.com/download-19ff08a4.html O que está faltando para incluir a primeira coluna? Edson Lira Estatístico Manaus-Amazonas

2 1

Quasi AIC representa a verossimilhança?
by ASANTOS 03 Jul '12

03 Jul '12

Boa noite pessoal, Estou tentando calcular o quasi AIC para distribuições quasi, achei um material do Ben Bolker, que supõem uma alternativa ajustando um modelo para uma família Poisson ou Binomial e outra para a Quasi. Então como solução ele propõem utilizar a log-verossimilhança do primeiro e o parâmetro de dispersão do modelo quasi ajusado. Onde eu pergunto: esses cálculos de Quasi AIC representam a verossimilhança do modelo e servem de comparação direta com o AIC das distribuições Poisson ou Binomial? Obrigado, -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

1 0