- R-br - DINF UFPR

Formatar grafico hexbin
by Bernardo Rangel Tura 28 Mai '12

28 Mai '12

Boa tarde a todos! Estou precisando de uma ajuda na formatação de um gráfico do tipo de hexbin. O problema é o seguinte o eixo Y tem valores muito altos na faixa de 2e8, 4e8, 6e8 etc Eu gostaria de substituir os dizeres para '200 Mi', '400 Mi', '600 Mi' etc para facilitar estou enviando um arquivo 'cmd.csv' que contém apenas 500 observações - o arquivo original tem mais de 100.000 mesmo assim vai ser suficiente O código que estou usando é: rm(list=ls(all=TRUE)) require(hexbin) base<-read.csv('cmd.csv') cores.ATS <- colorRampPalette(c('#2400D8','#181CF7','#2857FF','#3D87FF','#56B0FF','#75D3FF','#99EAFF','#BCF9FF','#EAFFFF','#FFFFEA','#FFF1BC','#FFD699','#FFAC75','#FF7856','#FF3D3D','#F72735','#D8152F','#A50021')) bin<-with(base,hexbin(Angio.sim, delta_modelos, xbins=25)) p <- plot(bin,colramp=cores.ATS,legend=.7,xlab='',ylab='') Alguém pode ajudar ? []s Tura

1 0

Re: [R-br] Digest R-br, volume 15, assunto 34
by Ary Casado 28 Mai '12

28 Mai '12

O meu dá o mesmo resultado. Alguma coisa nebulosa nessa PA. Ary ________________________________ De: "r-br-request(a)listas.c3sl.ufpr.br" <r-br-request(a)listas.c3sl.ufpr.br> Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Domingo, 27 de Maio de 2012 12:00 Assunto: Digest R-br, volume 15, assunto 34 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: problemas com função "seq" (Alexandre Ribeiro Leichsenring) 2. Re: problemas com função "seq" (Fred) 3. Re: problemas com função "seq" (Diego Leal) 4. Re: problemas com função "seq" (Cesar Rabak) 5. Re: problemas com função "seq" (Paulo Justiniano) ---------------------------------------------------------------------- Message: 1 Date: Sat, 26 May 2012 14:18:18 -0300 From: Alexandre Ribeiro Leichsenring <aleleich(a)uol.com.br> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] problemas com função "seq" Message-ID: <CAOZy71uykKy-bTatqiA97eLDWGFjcB_5DejrXObxM9G6GqP1mw(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Ivan, Talvez não tenha ficado claro. O problema não é de formatação: veja que onde deveria dar 0 (zero), a minha sequência está retornando 1.110223e-16! > seq(-.6,.8,by=.2)[1] -6.000000e-01 -4.000000e-01 -2.000000e-01 1.110223e-16 2.000000e-01 [6] 4.000000e-01 6.000000e-01 8.000000e-01 Aliás, formatando a sequência com scientific = FALSE, veja que resultado "curioso": > format(seq(-.6,.8,by=.2),scientific=FALSE)[1] "-0.5999999999999999777955" "-0.3999999999999999666933" [3] "-0.1999999999999999555911" " 0.0000000000000001110223" [5] " 0.2000000000000000666134" " 0.4000000000000000222045" [7] " 0.6000000000000001998401" " 0.8000000000000000444089" Imagino que isso seja um problema de alguma configuração do meu sistema/máquina. Ou alguém tem o mesmo problema? Saudações, Alexandre From: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> To: R Brasil <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] problemas com função "seq" Message-ID: <1337965493.37364.YahooMailNeo(a)web161806.mail.bf1.yahoo.com> Content-Type: text/plain; charset="iso-8859-1" Não há nada de errado! Para retornar exatamente o que você quer existem várias opções. Dentre elas, Não há nada de errado! Para retornar exatamente o que você quer existem várias opções. Dentre elas, a <- seq(-.6,.8,by=.2) round(a,2) format(a,scientific=TRUE) format(a,scientific=FALSE) format(a,digits=3,scientific=FALSE) formatC(a,format='f',digits=2,drop0trailing=TRUE) formatC(a,format='f',decimal.mark=',') as.numeric(formatC(a,format='f',decimal.mark='.')) (S,f,P) Allaman -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120526/5439bb51/atta…> ------------------------------ Message: 2 Date: Sat, 26 May 2012 17:38:25 -0300 From: Fred <eng.fredericoalvarez(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] problemas com função "seq" Message-ID: <4FC13F41.6030008(a)gmail.com> Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" Boa noite. Aparentemente isso está normal, no meu R também acontece isso. Lembre-se que os computadores trabalham em base binária, não decimal. Por esse motivo, é impossível representar números reais com 100% de exatidão. Saudações. Em 26-05-2012 14:18, Alexandre Ribeiro Leichsenring escreveu: > Ivan, > > Talvez não tenha ficado claro. O problema não é de formatação: veja > que onde deveria dar 0 (zero), a minha sequência está retornando > 1.110223e-16! > > > seq(-.6,.8,by=.2) > [1] -6.000000e-01 -4.000000e-01 -2.000000e-01 1.110223e-16 2.000000e-01 > [6] 4.000000e-01 6.000000e-01 8.000000e-01 > > Aliás, formatando a sequência com scientific = FALSE, veja que > resultado "curioso": > > > format(seq(-.6,.8,by=.2),scientific=FALSE) > [1] "-0.5999999999999999777955" "-0.3999999999999999666933" > [3] "-0.1999999999999999555911" " 0.0000000000000001110223" > [5] " 0.2000000000000000666134" " 0.4000000000000000222045" > [7] " 0.6000000000000001998401" " 0.8000000000000000444089" > > Imagino que isso seja um problema de alguma configuração do meu > sistema/máquina. Ou alguém tem o mesmo problema? > > Saudações, > > Alexandre > > > > > > > > From: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br > <mailto:ivanalaman@yahoo.com.br>> > To: R Brasil <r-br(a)listas.c3sl.ufpr.br <mailto:r-br@listas.c3sl.ufpr.br>> > Subject: Re: [R-br] problemas com função "seq" > Message-ID: > <1337965493.37364.YahooMailNeo(a)web161806.mail.bf1.yahoo.com > <mailto:1337965493.37364.YahooMailNeo@web161806.mail.bf1.yahoo.com>> > Content-Type: text/plain; charset="iso-8859-1" > > Não há nada de errado! Para retornar exatamente o que você quer > existem várias opções. Dentre elas, > > Não há nada de errado! Para retornar exatamente o que você quer > existem várias opções. Dentre elas, > > a <- seq(-.6,.8,by=.2) > round(a,2) > format(a,scientific=TRUE) > format(a,scientific=FALSE) > format(a,digits=3,scientific=FALSE) > formatC(a,format='f',digits=2,drop0trailing=TRUE) > formatC(a,format='f',decimal.mark=',') > as.numeric(formatC(a,format='f',decimal.mark='.')) > > (S,f,P) > Allaman > > > _______________________________________________ > 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. -- Frederico Steinmetz Alvarez Engenheiro Eletricista e de Produção Mestrando em Engenharia Eólica Laboratório de Fluidos Grupo de Mecânica dos Fluidos Ambientais da UFPE Home page ( www.ufpe.br/mecfluamb/ ) Universidade Federal de Pernambuco - UFPE Departamento de Engenharia Mecânica Fones: (48) 9143-5141 (VIVO) (81) 8885-1105 (OI) Laboratório Didático Pedagógico - Mecânica dos Flúidos LDP-FLU Universidade Federal de Pernambuco - UFPE -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120526/72c8181c/atta…> ------------------------------ Message: 3 Date: Sat, 26 May 2012 17:39:45 -0300 From: Diego Leal <diego.leight(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] problemas com função "seq" Message-ID: <4FC13F91.6050805(a)gmail.com> Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" Galera, como faço pra parar de receber emails do grupo. Apenas parar de receber emails, quero continuar fazendo parte do grupo... Att, Diego Leal Togni Em 26/05/2012 17:38, Fred escreveu: > Boa noite. > > Aparentemente isso está normal, no meu R também acontece isso. > Lembre-se que os computadores trabalham em base binária, não decimal. > Por esse motivo, é impossível representar números reais com 100% de > exatidão. > > Saudações. > > Em 26-05-2012 14:18, Alexandre Ribeiro Leichsenring escreveu: >> Ivan, >> >> Talvez não tenha ficado claro. O problema não é de formatação: veja >> que onde deveria dar 0 (zero), a minha sequência está retornando >> 1.110223e-16! >> >> > seq(-.6,.8,by=.2) >> [1] -6.000000e-01 -4.000000e-01 -2.000000e-01 1.110223e-16 2.000000e-01 >> [6] 4.000000e-01 6.000000e-01 8.000000e-01 >> >> Aliás, formatando a sequência com scientific = FALSE, veja que >> resultado "curioso": >> >> > format(seq(-.6,.8,by=.2),scientific=FALSE) >> [1] "-0.5999999999999999777955" "-0.3999999999999999666933" >> [3] "-0.1999999999999999555911" " 0.0000000000000001110223" >> [5] " 0.2000000000000000666134" " 0.4000000000000000222045" >> [7] " 0.6000000000000001998401" " 0.8000000000000000444089" >> >> Imagino que isso seja um problema de alguma configuração do meu >> sistema/máquina. Ou alguém tem o mesmo problema? >> >> Saudações, >> >> Alexandre >> >> >> >> >> >> >> >> From: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br >> <mailto:ivanalaman@yahoo.com.br>> >> To: R Brasil <r-br(a)listas.c3sl.ufpr.br <mailto:r-br@listas.c3sl.ufpr.br>> >> Subject: Re: [R-br] problemas com função "seq" >> Message-ID: >> <1337965493.37364.YahooMailNeo(a)web161806.mail.bf1.yahoo.com >> <mailto:1337965493.37364.YahooMailNeo@web161806.mail.bf1.yahoo.com>> >> Content-Type: text/plain; charset="iso-8859-1" >> >> Não há nada de errado! Para retornar exatamente o que você quer >> existem várias opções. Dentre elas, >> >> Não há nada de errado! Para retornar exatamente o que você quer >> existem várias opções. Dentre elas, >> >> a <- seq(-.6,.8,by=.2) >> round(a,2) >> format(a,scientific=TRUE) >> format(a,scientific=FALSE) >> format(a,digits=3,scientific=FALSE) >> formatC(a,format='f',digits=2,drop0trailing=TRUE) >> formatC(a,format='f',decimal.mark=',') >> as.numeric(formatC(a,format='f',decimal.mark='.')) >> >> (S,f,P) >> Allaman >> >> >> _______________________________________________ >> 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. > > > -- > Frederico Steinmetz Alvarez Engenheiro Eletricista e de Produção > Mestrando em Engenharia Eólica Laboratório de Fluidos Grupo de > Mecânica dos Fluidos Ambientais da UFPE Home page ( > www.ufpe.br/mecfluamb/ ) Universidade Federal de Pernambuco - UFPE > Departamento de Engenharia Mecânica Fones: (48) 9143-5141 (VIVO) (81) > 8885-1105 (OI) Laboratório Didático Pedagógico - Mecânica dos Flúidos > LDP-FLU Universidade Federal de Pernambuco - UFPE > > > _______________________________________________ > 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/20120526/c6298e04/atta…> ------------------------------ Message: 4 Date: Sat, 26 May 2012 20:47:02 -0300 From: Cesar Rabak <cesar.rabak(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] problemas com função "seq" Message-ID: <4FC16B76.8010800(a)gmail.com> Content-Type: text/plain; charset=ISO-8859-1; format=flowed Em 26/5/2012 14:18, Alexandre Ribeiro Leichsenring escreveu: > Ivan, > > Talvez não tenha ficado claro. O problema não é de formatação: veja que > onde deveria dar 0 (zero), a minha sequência está retornando 1.110223e-16! > >> seq(-.6,.8,by=.2) > [1] -6.000000e-01 -4.000000e-01 -2.000000e-01 1.110223e-16 2.000000e-01 > > [6] 4.000000e-01 6.000000e-01 8.000000e-01 > Alexandre, Talvez a gente não esteja sendo claro: o seu problema é de [falta] de entendimento mais basal: Por favor, leia o FAQ 7.31 e os documentos referenciados nele. []s -- Cesar Rabak GNU/Linux User 52247. Get counted: http://counter.li.org/ ------------------------------ Message: 5 Date: Sat, 26 May 2012 21:41:22 -0300 (BRT) From: Paulo Justiniano <paulojus(a)leg.ufpr.br> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] problemas com função "seq" Message-ID: <alpine.DEB.2.00.1205262140350.19183(a)pataxo.est.ufpr.br> Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" Vá na pagina da lista com sue usuário cadansgtrado onde voce pode definir as confiugurações de recebimento (e interrupcao se desejado) On Sat, 26 May 2012, Diego Leal wrote: > Galera, > como faço pra parar de receber emails do grupo. Apenas parar de receber emails, quero continuar fazendo parte do > grupo... > > Att, > Diego Leal Togni > > Em 26/05/2012 17:38, Fred escreveu: > Boa noite. > > Aparentemente isso está normal, no meu R também acontece isso. Lembre-se que os computadores trabalham > em base binária, não decimal. Por esse motivo, é impossível representar números reais com 100% de > exatidão. > > Saudações. > > Em 26-05-2012 14:18, Alexandre Ribeiro Leichsenring escreveu: > Ivan, > > Talvez não tenha ficado claro. O problema não é de formatação: veja que onde deveria dar 0 (zero), a > minha sequência está retornando 1.110223e-16! > > > seq(-.6,.8,by=.2) > [1] -6.000000e-01 -4.000000e-01 -2.000000e-01 1.110223e-16 2.000000e-01 > [6] 4.000000e-01 6.000000e-01 8.000000e-01 > > Aliás, formatando a sequência com scientific = FALSE, veja que resultado "curioso": > > > format(seq(-.6,.8,by=.2),scientific=FALSE) > [1] "-0.5999999999999999777955" "-0.3999999999999999666933" > [3] "-0.1999999999999999555911" " 0.0000000000000001110223" > [5] " 0.2000000000000000666134" " 0.4000000000000000222045" > [7] " 0.6000000000000001998401" " 0.8000000000000000444089" > > Imagino que isso seja um problema de alguma configuração do meu sistema/máquina. Ou alguém tem o mesmo > problema? > > Saudações, > > Alexandre > > > > > > > > From: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> > To: R Brasil <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] problemas com função "seq" > Message-ID: > <1337965493.37364.YahooMailNeo(a)web161806.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Não há nada de errado! Para retornar exatamente o que você quer existem várias opções. Dentre elas, > > Não há nada de errado! Para retornar exatamente o que você quer existem várias opções. Dentre elas, > > a <- seq(-.6,.8,by=.2) > round(a,2) > format(a,scientific=TRUE) > format(a,scientific=FALSE) > format(a,digits=3,scientific=FALSE) > formatC(a,format='f',digits=2,drop0trailing=TRUE) > formatC(a,format='f',decimal.mark=',') > as.numeric(formatC(a,format='f',decimal.mark='.')) > > (S,f,P) > Allaman > > > _______________________________________________ > 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. > > > > -- > Frederico Steinmetz Alvarez Engenheiro Eletricista e de Produção Mestrando em Engenharia Eólica Laboratório de > Fluidos Grupo de Mecânica dos Fluidos Ambientais da UFPE Home page ( www.ufpe.br/mecfluamb/ ) Universidade > Federal de Pernambuco - UFPE Departamento de Engenharia Mecânica Fones: (48) 9143-5141 (VIVO) (81) 8885-1105 > (OI) Laboratório Didático Pedagógico - Mecânica dos Flúidos LDP-FLU Universidade Federal de Pernambuco - UFPE > > > _______________________________________________ > 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 Fim da Digest R-br, volume 15, assunto 34 *****************************************

1 0

problemas com função "seq"
by Alexandre Ribeiro Leichsenring 28 Mai '12

28 Mai '12

Caros, Volta e meia me deparo com um comportamento estranho da função "seq", compartilho agora um exemplo: > seq(-.6,.8,by=.2) ao invés de resultar nisso: -.6 -.4 -.2 0 .2 .4 .6 .8 resulta nisso: -6.000000e-01 -4.000000e-01 -2.000000e-01 1.110223e-16 2.000000e-01 4.000000e-01 6.000000e-01 8.000000e-01 ! Outras sequências ao redor de zero resultam no mesmo "erro" (pelo menos na minha máquina)... que tipo de problema será esse? Obrigado, Alexandre

4 3

Re: [R-br] problemas com função "seq"
by Alexandre Ribeiro Leichsenring 26 Mai '12

26 Mai '12

Ivan, Talvez não tenha ficado claro. O problema não é de formatação: veja que onde deveria dar 0 (zero), a minha sequência está retornando 1.110223e-16! > seq(-.6,.8,by=.2)[1] -6.000000e-01 -4.000000e-01 -2.000000e-01 1.110223e-16 2.000000e-01 [6] 4.000000e-01 6.000000e-01 8.000000e-01 Aliás, formatando a sequência com scientific = FALSE, veja que resultado "curioso": > format(seq(-.6,.8,by=.2),scientific=FALSE)[1] "-0.5999999999999999777955" "-0.3999999999999999666933" [3] "-0.1999999999999999555911" " 0.0000000000000001110223" [5] " 0.2000000000000000666134" " 0.4000000000000000222045" [7] " 0.6000000000000001998401" " 0.8000000000000000444089" Imagino que isso seja um problema de alguma configuração do meu sistema/máquina. Ou alguém tem o mesmo problema? Saudações, Alexandre From: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> To: R Brasil <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] problemas com função "seq" Message-ID: <1337965493.37364.YahooMailNeo(a)web161806.mail.bf1.yahoo.com> Content-Type: text/plain; charset="iso-8859-1" Não há nada de errado! Para retornar exatamente o que você quer existem várias opções. Dentre elas, Não há nada de errado! Para retornar exatamente o que você quer existem várias opções. Dentre elas, a <- seq(-.6,.8,by=.2) round(a,2) format(a,scientific=TRUE) format(a,scientific=FALSE) format(a,digits=3,scientific=FALSE) formatC(a,format='f',digits=2,drop0trailing=TRUE) formatC(a,format='f',decimal.mark=',') as.numeric(formatC(a,format='f',decimal.mark='.')) (S,f,P) Allaman

5 4

Re: [R-br] problemas com função "seq"
by Ivan Bezerra Allaman 25 Mai '12

25 Mai '12

Não há nada de errado! Para retornar exatamente o que você quer existem várias opções. Dentre elas, a <- seq(-.6,.8,by=.2) round(a,2) format(a,scientific=TRUE) format(a,scientific=FALSE) format(a,digits=3,scientific=FALSE) formatC(a,format='f',digits=2,drop0trailing=TRUE) formatC(a,format='f',decimal.mark=',') as.numeric(formatC(a,format='f',decimal.mark='.')) (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

Dificuldade em maximização/minimização de funções! (Repassando help Adailton Alcantara)
by Ivan Bezerra Allaman 25 Mai '12

25 Mai '12

Adailton, Procure seguir o guia de postagem e evite anexar arquivos! Repassando então a dúvido de Adailton Alcantara para os nobres amigos do R. OBS: Segue o link no final da postagem com as funções que o Adailton criou! Prezados, Nome: Adeilton Pedro de Alcântara Venho pedir ajuda. Definição do problema: Eu tenho uma superfície (mistura de cinco distribui��es normal bivariada) lambda_23(x,y) e uma reta lambda_1(t). Considerando a reta e a superf�cie como conhcidas ( distribui��es verdadeiras) o meu objetivo � recuperar essas distribui��es maximuzando a equa��o de log verossimilhan�a ( no programa que enviei denominada matricaalg1). Ou seja o meu objetivo � recuperar a estrutura do polin�mio que gerou um determinado conjunto de dados utilizando uma equa��o de log-verossimilhan�a sem a necessidade de escrever um modelo param�trico. A minha dificuldade est� no comando nlminb() do R que encontra-se na fun��o metricaalg1. Eu n�o consigo lidar com estes comandos de maximiza��o/minimiza��o de fun��es. Quando altero o chute inicial, embora tenha converg�ncia, os valores encontrados s�o completamente diferentes do valor que eu quero encontrar. Ou seja dependendo do chute inicial a fun��o me retorna coeficientes completamente diferentes. Eu n�o sei se o problema est� na op��o control.list() ou no chute inicial. Provavelmente n�o. O programa que gera o PPNH( Processo de Poisson n�o homog�neo) est� funcionando. A dificuldade est� no uso do comando nlminb(). J� utilizei nlm(), optim(), e DEoptim() para maximizar a equa��o de log-verossimilhan�a por�m n�o consigo enxergar o que est� errado . O problema deve ser de m�ltiplos m�ximos na verossimilhan�a, pois tenho muitos coeficientes. Já estou a quatro anos e não consigo sair disso. Atenciosamente, Adeilton Pedro de Alcântara. Estou enviando um c�digo na linguagem R. O objetivO � recuperar uma reta (lambda_1(t) )e uma mistura de 5 distribui��es normal bivariada (lambda_23(x,y) ). A minha dificuldade est� no chute inicial, ou seja os resultados dependem muito do chute inicial. Al�m disso, cada vez que eu modifici o chute inicial tenho um resultado completamente diferente. O senhor verifica a fun��o onde eu gero o chute inicial (CHUTEINICIALalg1) e a fun��o metricaalg1 (fun��o onde tenho a verossimilhan�a). Nesta fun��o estou fazendo a maximiza��o num�rica da equa��o de log-verossimilhan�a utilizaddo o comando "nlminb() do R. Onde tem CC<-(as.single(CC)>0) # Coloquei positivo para n�o ficar muito tempo rodando. O senhor pode deixar CC<-(as.single(CC)). Para rodar a fun��o basta seguir os comandos a partir de ################# #Rodando a fun��o# ################# O arquivo com as rotinas do R est� no anexo deste e-mail Me ajude, por favor. Atenciosamente, Adeilton Pedro de Alcântara http://dl.dropbox.com/u/33619290/codigoAdailton.txt \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

odfWeave problemas com acentos
by Walmes Zeviani 24 Mai '12

24 Mai '12

Saudações, Estou montando um micro tutorial de odfWeave e observei problemas na geração de documentos. Eu tenho um arquivo odf, gerado pelo OpenOffice Writer com o seguinte conteúdo #---- conteúdo do arquivo odf, walmes_teste1.odf ---- Isso é um teste <<teste, echo=TRUE>> x <- rnorm(100,50,3) summary(x) @ #---- fim do arquivo walmes_teste1.odf --- Abro uma sessão R para compilar por meio dos comandos que apresentam a mensagem de erro abaixo. *> odfWeave("walmes_teste1.odt", "walmes_saida1.odt")* Copying walmes_teste1.odt Setting wd to /tmp/RtmpCE5KJb/odfWeave24115817941 Unzipping ODF file using unzip -o walmes_teste1.odt Archive: walmes_teste1.odt extracting: mimetype creating: Configurations2/statusbar/ inflating: Configurations2/accelerator/current.xml creating: Configurations2/floater/ creating: Configurations2/popupmenu/ creating: Configurations2/progressbar/ creating: Configurations2/menubar/ creating: Configurations2/toolbar/ creating: Configurations2/images/Bitmaps/ inflating: content.xml inflating: manifest.rdf inflating: styles.xml extracting: meta.xml inflating: Thumbnails/thumbnail.png inflating: settings.xml inflating: META-INF/manifest.xml Removing walmes_teste1.odt Creating a Pictures directory Pre-processing the contents Sweaving content.Rnw *Erro: ‘content.Rnw’ is not ASCII and does not declare an encoding* Se eu removo o e acentuado da frame "isso é um teste" para "isso e um teste", a compilação tem sucesso. Googlei e não encontrei alternativas para superar o problema. Minha sessão R é *> sessionInfo()* R version 2.15.0 (2012-03-30) Platform: i486-pc-linux-gnu (32-bit) locale: [1] LC_CTYPE=pt_BR.utf8 LC_NUMERIC=C [3] LC_TIME=pt_BR.utf8 LC_COLLATE=pt_BR.utf8 [5] LC_MONETARY=pt_BR.utf8 LC_MESSAGES=pt_BR.utf8 [7] LC_PAPER=C LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C [11] LC_MEASUREMENT=pt_BR.utf8 LC_IDENTIFICATION=C attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] odfWeave_0.7.17 XML_3.6-1 lattice_0.20-6 loaded via a namespace (and not attached): [1] grid_2.15.0 tools_2.15.0 Estou no Ubuntu 10.04 com o seguinte resultado de locale walmes@galibi:~$ locale LANG=pt_BR.utf8 LANGUAGE=pt_BR:pt:en LC_CTYPE="pt_BR.utf8" LC_NUMERIC="pt_BR.utf8" LC_TIME="pt_BR.utf8" LC_COLLATE="pt_BR.utf8" LC_MONETARY="pt_BR.utf8" LC_MESSAGES="pt_BR.utf8" LC_PAPER="pt_BR.utf8" LC_NAME="pt_BR.utf8" LC_ADDRESS="pt_BR.utf8" LC_TELEPHONE="pt_BR.utf8" LC_MEASUREMENT="pt_BR.utf8" LC_IDENTIFICATION="pt_BR.utf8" LC_ALL= walmes@galibi:~$ Alguém tem alguma sugestão? Grato. 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 ==========================================================================

3 4

Help
by Adeilton Alcantara 24 Mai '12

24 Mai '12

Nome: Adeilton Pedro de Alcântara

1 0

pergunta simples - regular expressions
by Vitor Aguiar 24 Mai '12

24 Mai '12

Caros, faz um tempo que o Benílton me passou um código muito útil, mas como não entendo de regular expressions não sei como exatamente ele funciona. Será que alguém poderia me ajuda a entender? É algo muito simples. Eu tenho um dataframe com várias colunas com nomes como: "PentaE.1" "PentaE.2" "D18S51.1" "D18S51.2" "D21S11.1" "D21S11.2"... Então Benilton me passou o seguinte código: Names = unique(gsub("(.*)\\.\\d{1}$", "\\1", names(x))) que me gera o seguinte resultado: "PentaE" "D18S51" "D21S11" Eu sei que gsub está fazendo um matching do primeiro argumento no meu daframe e substituindo pelo segundo argumento, retirando tudo que está após o "." (ponto) e criando apenas 1 nome onde havia 2. Mas não entendo tudo, como o código produz esse resultado? Quais as funções dessas \\ e d{}... Obrigado, Vitor

1 1

eliminar valores de data.frame
by Vitor Aguiar 24 Mai '12

24 Mai '12

Caros, tenho um data frame com 30 colunas e 100 mil linhas, composto de valores numéricos de 1 a 65. Gostaria de eliminar todos os valores que aparecem menos de 5 vezes, em cada 2 colunas consecutivas separadamente. Por exemplo, se o valor "18" aparece 4 vezes entre os 200.000 valores das minhas 2 primeiras colunas, eu gostaria de eliminar todos os 4. Eu sei que no final eu teria colunas de lengths diferentes no meu data frame, mas eu não preciso recriar o data frame, eu só preciso usar os pares de colunas separadamente. Como fazer isso? Desde já, muito obrigado aos que puderem ajudar. Vitor Aguiar

2 6