- R-br - DINF UFPR

Identação do Código R
by Pedro Rafael 25 Ago '11

25 Ago '11

Pessoal, uma dúvida besta. Existe um padrão para identação de um código R? Várias linguagens usam um padrão específico como Python por exemplo. -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

6 8

Plotar dados em Mapa
by Leandro Marino 25 Ago '11

25 Ago '11

Caros, Mais uma vez venho pedir um caminho... Tenho um indicador por município que gostaria que fosse “plotado” em um mapa do Brasil aonde tivesse uma escala de cores para cada um dos níveis deste indicador. Sei que existe alguma forma de fazer, no entanto, nunca trabalhei com mapas no R e não sei nem por onde começar (Pacote, Função e principalmente MAPA). Agradeço qualquer ajuda!!! Um abraço Leandro Atenciosamente, Leandro Marino http://www.leandromarino.com.br (Fotógrafo) http://est.leandromarino.com.br/Blog (Estatístico) Cel.: + 55 21 9845-7707 Cel.: + 55 21 8777-7907

5 10

Re: [R-br] Digest R-br, volume 6, assunto 41
by Julio Bueno 25 Ago '11

25 Ago '11

No caso específico do exemplo dado, consegui alocar esta matriz da seguinte forma: a <- expand.grid(1:5,1:5,1:5,1:5,1:5,1:5,1:5,1:5,1:5,1:5,1:5) o problema é que ela ocupa 7.1 Gb! Em qual sistema (hardw. e s.o.) você está rodando o R? []s Júlio 2011/8/25 <r-br-request(a)listas.c3sl.ufpr.br> > 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: Problema de memória: MATRIZES GIGANTES (Alexandre Castagna) > 2. Função eigen (Junior Beleti) > 3. Re: Otimizar tempo de processamento (Tiago de Oliveira Vieira) > 4. Re: Otimizar tempo de processamento (Benilton Carvalho) > 5. Representação usando expression() (Alexandre Santos) > 6. Re: Representação usando expression() (Benilton Carvalho) > 7. Re: Representação usando expression() (Alexandre Santos) > 8. Re: Objetos da Classe SpatialPolygonsDataFrame > (Leonel Lucena Rodrigues) > 9. [OFF- Dúvida] Linguagem de programação para trabalhar com > estatística. (Pedro Rafael) > 10. Pedro Rafael convidou você para abrir uma conta de e-mail do > Google (Pedro Rafael) > 11. [OFF TOPIC] - Learn R Webinar (Gilson Sanchez) > 12. Plotar dados em Mapa (Leandro Marino) > 13. [Dúvida] Invocar funções C++ no R (Pedro Rafael) > 14. Re: [Dúvida] Invocar funções C++ no R (Benilton Carvalho) > 15. Re: [Dúvida] Invocar funções C++ no R (Pedro Rafael) > 16. Re: [Dúvida] Invocar funções C++ no R (Benilton Carvalho) > 17. Re: [Dúvida] Invocar funções C++ no R (Pedro Rafael) > 18. Re: Plotar dados em Mapa (Jakson Alves de Aquino) > 19. Re: [Dúvida] Invocar funções C++ no R (Benilton Carvalho) > 20. Agrupamento de Tocher no R (Emmanuel Arnhold) > 21. Re: Plotar dados em Mapa (Leandro Marino) > 22. Re: Plotar dados em Mapa (Rodrigo Plei) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 24 Aug 2011 23:06:56 -0300 > From: Alexandre Castagna <alexandre.castagna(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Problema de memória: MATRIZES GIGANTES > Message-ID: > <CAFiuohLh13nV29aDu=KuPvhqaaUwAtdg_gvk9tsMsoBp7ARQcg(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Se não me engano, há um tempo atrás correu na lista o mesmo problema e me > lembro de terem citado uma solução que fazia o R trabalhar com dados > reescrevendo e lendo do disco, permitindo a análise de bancos maiores, mas > com custo de tempo de processamento. > > Em 24 de agosto de 2011 19:56, Pedro Rafael > <pedro.rafael.marinho(a)gmail.com>escreveu: > > > Sim, o Revolution R é bom mas novas mudanças demoram mais para chegar > nele, > > como novos pacotes, correções, etc. > > > > Pedro > > > > Em 24 de agosto de 2011 18:06, Alexandre Castagna [via R-br] < > > ml-node+3766656-34166862-223914(a)n4.nabble.com> escreveu: > > > >> Caros, > >> > >> Também estou enfrentando um problema e memória e encontrei um software > >> chamado Revolution R Enterprise (http://www.revolutionanalytics.com/) > que > >> é compatível com o R, mas permite análise de grandes bancos de dados; > ele é > >> pago, mas tem licença gratuita para acadêmicos. Ainda não testei (acabei > de > >> instalar), mas segue a dica se alguém se interessar. > >> > >> Abraços, > >> > >> Em 23 de agosto de 2011 20:03, Márcio Luís Moreira de Souza <[hidden > >> email] <http://user/SendEmail.jtp?type=node&node=3766656&i=0>> > escreveu: > >> > >> Olá, Henrique! > >>> > >>> Muito obrigado pela atenção. > >>> Darei uma olhada. > >>> > >>> Abração!! > >>> > >>> > >>> Em 23 de agosto de 2011 18:05, Henrique Leme Felizatti <[hidden email]< > http://user/SendEmail.jtp?type=node&node=3766656&i=1> > >>> > escreveu: > >>> > >>> Existe um projeto em R, com respectivo pacote, denominado BigMemory ( > >>>> www.*bigmemory*.org/) que serve exatamente para isso. Nem todas as > >>>> funções do R estão disponíveis nesse projeto no entanto as principais > >>>> funções de manipulação de dados está implementada. > >>>> > >>>> Atenciosamente > >>>> > >>>> Henrique > >>>> > >>>> > >>>> > >>>> Em 23 de agosto de 2011 17:01, Márcio Luís Moreira de Souza <[hidden > >>>> email] <http://user/SendEmail.jtp?type=node&node=3766656&i=2 > >>escreveu: > >>>> > >>>> Tentei rodar: > >>>>> > >>>>> expand.grid(1:5,1:5,1:5,1:5,1:5,1:5,1:5,1:5,1:5,1:5,1:5)E deu o > resultado > >>>>> Error: cannot allocate vector of size 186.3 Mb > >>>>> > >>>>> O que posso fazer?? Alguma sugestão? > >>>>> > >>>>> > >>>>> p.s.: tal comando geraria um data.frame com números inteiros. > >>>>> Meu PC é HP pavilon (processador I5), 4 Gb de memória... > >>>>> > >>>>> > >>>>> > >>>>> > >>>>> Em 23 de agosto de 2011 16:30, Benilton Carvalho <[hidden email]< > http://user/SendEmail.jtp?type=node&node=3766656&i=3> > >>>>> > escreveu: > >>>>> > >>>>>> defina gigantes. > >>>>>> > >>>>>> eu trabalho com 7M+ linhas e centenas de colunas diariamente e sem > >>>>>> problemas. > >>>>>> > >>>>>> vc precisa ser especifico com o seu problema. matrizes de reais? > >>>>>> inteiros? logicos? > >>>>>> > >>>>>> que OS? qto de RAM? que operacoes precisa na matriz? > >>>>>> > >>>>>> b > >>>>>> > >>>>>> > >>>>>> _______________________________________________ > >>>>>> R-br mailing list > >>>>>> [hidden email] < > http://user/SendEmail.jtp?type=node&node=3766656&i=4> > >>>>>> > >>>>>> 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. > >>>>>> > >>>>> > >>>>> > >>>>> > >>>>> -- > >>>>> > >>>>> Prof. Márcio Souza, M.Sc. (DE/UFJF) > >>>>> www.profmsouza.blogspot.com > >>>>> http://lattes.cnpq.br/4578008002785684 > >>>>> > >>>>> > >>>>> _______________________________________________ > >>>>> R-br mailing list > >>>>> [hidden email] <http://user/SendEmail.jtp?type=node&node=3766656&i=5 > > > >>>>> > >>>>> 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=3766656&i=6> > >>>> > >>>> 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. > >>>> > >>> > >>> > >>> > >>> -- > >>> > >>> Prof. Márcio Souza, M.Sc. (DE/UFJF) > >>> www.profmsouza.blogspot.com > >>> http://lattes.cnpq.br/4578008002785684 > >>> > >>> > >>> _______________________________________________ > >>> R-br mailing list > >>> [hidden email] <http://user/SendEmail.jtp?type=node&node=3766656&i=7> > >>> > >>> 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. > >>> > >> > >> > >> > >> -- > >> Alexandre Castagna Mourão e Lima > >> Laboratório de Radioecologia e Mudanças Globais - UERJ > >> Rua São Francisco Xavier, 524, Maracanã > >> Pav. Haroldo Lisboa da Cunha, Sub-solo > >> > >> _______________________________________________ > >> R-br mailing list > >> [hidden email] <http://user/SendEmail.jtp?type=node&node=3766656&i=8> > >> 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-Problema-de-memoria-MATRIZES-GIGANTE… > >> To unsubscribe from R-br, click here< > http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib… > >. > >> > >> > > > > > > > > -- > > Saudações, > > Pedro Rafael Diniz Marinho - Estatístico SES-PB > > Currículo Lattes: > > http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 > > > > > > _______________________________________________ > > 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. > > > > > > -- > Alexandre Castagna Mourão e Lima > Laboratório de Radioecologia e Mudanças Globais - UERJ > Rua São Francisco Xavier, 524, Maracanã > Pav. Haroldo Lisboa da Cunha, Sub-solo > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110824/17055c17/atta… > > > > ------------------------------ > > Message: 2 > Date: Thu, 25 Aug 2011 00:48:44 -0300 > From: Junior Beleti <beleti.junior(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Função eigen > Message-ID: > <CA+1Tbo9c-c=1sKdj17JireTVMi3ccnABvuQWdKZApqicJ6nwGA(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Boa Noite.. > > Estou analisando um script que realiza um loop e chama a função eigen. > Com o pacote debug constatei que a execução desse script gasta muito tempo > nessa função. > Vi também que ela faz chamada para código fortran. > > Gostaria de saber se existe outra função que realiza o mesmo que a eigen. > > Att, > > -- > Carlos Roberto Beleti Junior > Mestrado em Ciência da Computação > Departamento de Informática > Universidade Estadual de Maringá > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110825/ce329385/atta… > > > > ------------------------------ > > Message: 3 > Date: Thu, 25 Aug 2011 01:10:37 -0300 > From: Tiago de Oliveira Vieira <tiagovieira7(a)yahoo.com.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Otimizar tempo de processamento > Message-ID: <1314245437.2540.5.camel@tiago-Positivo-Mobile> > Content-Type: text/plain; charset="UTF-8" > > Benilton, > > Segui suas dicas e o programa rodou bem, mas quando aumento o número > máximo de filhos "f" ou o número de gerações "n" e programa não roda, > veja: > > > library(doMC) > > registerDoMC() > > programa2 <- function(f,n,m) { > + ## f: número máximo de filhos > + ## n: número de gerações > + ## m: número de iterações > + S <- foreach(p=seq(0,1,by=0.01), .combine=rbind) %dopar% { > + pp <- 0 > + for (g in 1:m) { > + y <- 1 > + for (k in 1:n) y <- rbinom(1, y*f, p) > + if (y > 0) pp <- pp + 1 > + } > + c(p, pp/m) > + } > + return(S) > + } > > system.time(prog2 <- programa2(3,100,1000)) > Erro em { : task 48 failed - "valor ausente onde TRUE/FALSE necessário" > Além disso: Houve 50 ou mais avisos (use warnings() para ver os > primeiros 50) > Timing stopped at: 61.27 0.03 61.58 > > warnings() > Mensagens de aviso: > 1: In rbinom(1, y * f, p) : NAs produzidos > 2: In rbinom(1, y * f, p) : NAs produzidos > 3: In rbinom(1, y * f, p) : NAs produzidos > 4: In rbinom(1, y * f, p) : NAs produzidos > > Em Qui, 2011-08-18 às 17:45 -0300, Leonard Assis escreveu: > > Entendi sua explanação. heehhe > > > > olhei o código bem por alto depois do almoço num tempó livre, única > > coisa "estranha" que notei foi esse loop que pra mim estava sem muito > > sentido. > > > > Mas muito dessa minha opinião vai do fato de não conhecer ao fundo o > > problema, pois para mim, ele estaria interessado em cada uma das > > iterações (eu, quando faço simulações, procuro guardar cada uma delas, > > nem que seja pra brincar de bootstrap depois, hehehe) > > > > Quanto à copa américa, estava meio na cara que daria Uruguai, é o time > > com mais conjunto, Argentina e Brasil estavam em início de trabalho. > > No betclick eu apostei no uruguai e salvei algum, ehhehe > > > > lmassis <at> yahoo <dot> com <dot> br > > assis.leonard <at> gmail <dot> com > > > > > > > > 2011/8/18 Benilton Carvalho <beniltoncarvalho(a)gmail.com> > > Leonard, > > > > para cada 'g', ele determina uma realizacao (aleatoria) de Y_n > > condicionado em Y_{n-1}. Se essa realizacao por nao-nula, ele > > incrementa 'pp'. Como a chance de Y_m ser zero e' bem pequena, > > ele > > repete o laco 'm' vezes, contando quantas vezes (em 'm') ele > > observa > > Y_n > 0... > > > > Estou bem certo de que o codigo esta' correto... mas eu tambem > > estava > > certo que o Brasil ganharia a Copa America... > > > > b > > > > 2011/8/18 Leonard Assis <assis.leonard(a)gmail.com>: > > > > > Benilton, eu nem cheguei a rodar, analisei o que ele faz > > mesmo. > > > > > > Tem que resolver esse problema daquele loop que só faz > > perder tempo senão o > > > próprio estudo de simulação dele está inválido. > > > lmassis <at> yahoo <dot> com <dot> br > > > assis.leonard <at> gmail <dot> com > > > > > > > > > 2011/8/18 Benilton Carvalho <beniltoncarvalho(a)gmail.com> > > >> > > >> Tiago, > > >> > > >> esse e' um problema legal... se vc derivar as distribuicoes > > a mao, > > >> certamente vc consegue ganhos bem maiores... > > >> > > >> FYI: usando propriedades de estatistica, vc chega na > > mudanca que eu > > >> fiz na geracao das suas variaveis aleatorias. > > >> > > >> Adicionalmente, eu nao tive paciencia de esperar o seu > > codigo original > > >> acabar de rodar... eu o cancelei depois de um pouco mais de > > 1h10m. > > >> > > >> Abaixo uma sugestao para voce: > > >> > > >> > > >> library(foreach) > > >> library(doMC) > > >> registerDoMC() > > >> programa2 <- function(f,n,m) { > > >> ## f: número máximo de filhos > > >> ## n: número de gerações > > >> ## m: número de iterações > > >> S <- foreach(p=seq(0,1,by=0.01), .combine=rbind) %dopar% > > { > > >> pp <- 0 > > >> for (g in 1:m) { > > >> y <- 1 > > >> for (k in 1:n) > > >> y <- rbinom(1, y*f, p) > > >> if (y > 0) pp <- pp + 1 > > >> } > > >> c(p, pp/m) > > >> } > > >> return(S) > > >> } > > >> system.time(prog2 <- programa2(2,200,1000)) > > >> plot(prog2, xlab = "p", ylab = "Probabilidade de > > Percolação") > > >> > > >> para matar sua curiosidade, no meu computador, eu obtive: > > >> > > >> > system.time(prog2 <- programa2(2,200,1000)) > > >> user system elapsed > > >> 257.072 2.562 20.835 > > >> _______________________________________________ > > >> 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. > > > > > > > > > > > > > -- > > Successful people ask better questions, and as a result, they > > get > > better answers. (Tony Robbins) > > > > _______________________________________________ > > 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. > > > > > ------------------------------ > > Message: 4 > Date: Thu, 25 Aug 2011 10:08:45 +0100 > From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Otimizar tempo de processamento > Message-ID: > <CAO-arWPqAgQdoTvREqL0YPxpEDycvEPm398=PLdYr=wmsMj33A(a)mail.gmail.com > > > Content-Type: text/plain; charset=ISO-8859-1 > > Tiago, > > isso e' resultado de y sendo igual a NA ou integer(0). > > Vc vai precisar trocar o %dopar% por %do% para debugar. Uma sugestao > e' comecar por adiciionar umas condicoes para teste dentro do for() > > for (k in 1:n) y <- rbinom(1, y*f, p) > > usar algo como > > for (k in 1:n){ > y <- rbinom(1, y*f, p) > if (length(y) == 0){ > browser() > }else{ > if (!is.finite(y)) > browser() > } > } > > b > > > ------------------------------ > > Message: 5 > Date: Thu, 25 Aug 2011 04:02:04 -0700 (PDT) > From: Alexandre Santos <alexandresantosbr(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Representação usando expression() > Message-ID: > <1314270124.69929.YahooMailNeo(a)web38404.mail.mud.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Bom > > Bom dia pessoal, > Gostaria de saber s é possível fazer com que em: > > > text(33,0.9,expression(italic("y")==italic("exp")(-"0.9588"+"0.0841"*italic("*x")-"0.0010"*italic("*x"))/(1+italic("exp")(-"0.9588"+"0.0841"*italic("*x")-"0.0010"*italic("*x")))),cex=0.8) > > A equação aparecesse no gráfico como: > ("exp")(-"0.9588"+"0.0841"*italic("*x")-"0.0010"*italic("*x")) sobre > (1+italic("exp")(-"0.9588"+"0.0841"*italic("*x")-"0.0010"*italic("*x") e não > separado por / , isso é possível? > Obrigado, > > > Alexandre dos Santos > Engenheiro Florestal, MSc. > Universidade Federal de Lavras > Departamento de Entomologia > Laboratório de Entomologia Florestal > Caixa Postal 3037 > 37200-000 - Lavras/MG > Fone: +55 (35) 3829-5122 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110825/eb737094/atta… > > > > ------------------------------ > > Message: 6 > Date: Thu, 25 Aug 2011 12:17:06 +0100 > From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, Alexandre Santos > <alexandresantosbr(a)yahoo.com.br> > Subject: Re: [R-br] Representação usando expression() > Message-ID: > <CAO-arWNKiAoPemWRGgichpb+JkNU2FVQg_gX96RmvpVtDPXJ6w(a)mail.gmail.com > > > Content-Type: text/plain; charset=ISO-8859-1 > > Tomei a liberdade de simplificar um pouco sua expressao (e' > equivalente ao modelo que vc descreve): > > txt = expression(frac(1, 1+e^{0.9588-0.0841*x+0.0010*x^2})) > plot(1, type='n') > legend(1, 1, txt, bty='n') > > b > > > ------------------------------ > > Message: 7 > Date: Thu, 25 Aug 2011 04:29:53 -0700 (PDT) > From: Alexandre Santos <alexandresantosbr(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Representação usando expression() > Message-ID: <1314271793.2934.YahooMailNeo(a)web38406.mail.mud.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Obrigado Benilton, > > Resolveu o meu problema e ficaram muito mais elegante os argumentos > também, > Abraço, > > > Alexandre dos Santos > Engenheiro Florestal, MSc. > Universidade Federal de Lavras > Departamento de Entomologia > Laboratório de Entomologia Florestal > Caixa Postal 3037 > 37200-000 - Lavras/MG > Fone: +55 (35) 3829-5122 > > > > Tomei a liberdade de simplificar um pouco sua expressao (e' > equivalente ao modelo que vc descreve): > > txt = expression(frac(1, 1+e^{0.9588-0.0841*x+0.0010*x^2})) > plot(1, type='n') > legend(1, 1, txt, bty='n') > > b > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110825/5385b15f/atta… > > > > ------------------------------ > > Message: 8 > Date: Thu, 25 Aug 2011 08:29:54 -0300 > From: Leonel Lucena Rodrigues <leonellucena(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Objetos da Classe SpatialPolygonsDataFrame > Message-ID: > <CADovVG3hzHenVknnN2HVuBgEUTS-pwzhsjfQU-4G_6S0rsxZxw(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Caro Prof. Paulo > > Eu esta me referinto a matrix de coordenadas que define cada poligono, > pois > o motivo da obtenção dessas coordenadas é para usar > na função kernel2d que, junto com a função image e um raio especificado, > estimará a desnidade kernel. > Se o senhor poder me ajudar ficaria grato. > > Att, Leonel Lucena > > > > Em 24 de agosto de 2011 16:43, Paulo Justiniano <paulojus(a)leg.ufpr.br > >escreveu: > > > Leonel > > > > certamente mas precisa ficar claro o que voce quer dizzer por > coordenadas. > > Voce deseja: > > 1) a matrix de coordenadas que define cada poligono > > 2) as coordenadas de um ponto que represent o polygono (pr exemplo o > > centroide dado por labpt > > > > se for 2) fica ainda mais direto bastando > > > > coordinates(dados) > > > > onde dados é um SpatialPolygonDataFrame > > (ou SpatialPolygons) > > > > > > tudo isto está nas funcionalidades do pacote sp > > > > > > > > > > On Tue, 23 Aug 2011, Leonel Lucena Rodrigues wrote: > > > > > >> Muito obrigado Prof. pela sua atenção, mas eu não fui muito explicativo > >> ao me referir às coordenadas. > >> > >> Se o Sr. me permitir enviar o shp de Fortaleza e os comandos para que eu > >> possa me explicar melhor ficaria muito grato. > >> > >> Abraço, Leonel. > >> > >> > >> > >> > >> Em 23 de agosto de 2011 15:38, Paulo Justiniano <paulojus(a)leg.ufpr.br> > >> escreveu: > >> 1o as coords estarao disponíveis para cada poligono individualmente > >> > >> Por exemplo, suponha "dados" em que cada polygono seja definido por > >> um unico conjun to de coordenadas > >> > >> dados@polygons[[1]]@Polygons[[**1]]@coords > >> > >> Com isto se quiser extrair dos varios polygonos voce pode montar > ujm > >> lapply() > >> > >> > >> 2o possiveis vieses podem ser dados por: forma de considerar (ou > nao > >> considerar...) bordas e sistemas de coordenadas (por > >> ex usar lat/long que é inadequada paa calculos baseados em > distancias > >> euclidianas) > >> > >> > >> > >> > >> On Tue, 23 Aug 2011, Leonel Lucena Rodrigues wrote: > >> > >> Boa Tarde, > >> > >> Estou com duas dúvidas sobre objetos da classe > >> SpatialPolygonsDataFrame. > >> Tenho arquivos em formato shapefile, que é considerado pelo R como > um > >> objeto da classes polígono espacial, e este tem > >> alguns slot's > >> como "labpt" e "area". Então venho pedir ajuta em duas questões. > >> > >> 1º Qual a função, e o pacote, que me fornecerá apenas o slot > "coords" > >> ? > >> > >> Sei existe uma forma de fazer com que o programa selecione > >> as coordenadas de polígono, atualmente utilizo a > >> função scan(), mas > >> gostaria de > >> utilizar uma outra. > >> > >> > >> 2º Caso venha a conseguir alguma estimativa, Kernel e.g., existirá > >> um viés > >> se os dados não tiverem projeção ? > >> > >> Na cartografia, mais especificamente no ArcMAP, utilizo o > sistema > >> de coordenadas adotado atualmente pelo IBGE > >> (SIRGAS 2000), no R > >> esse processo está em desenvolvimento. > >> > >> > >> Abraço a todos. > >> > >> > >> > >> > >> ______________________________**_________________ > >> 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. > >> > >> > >> > >> > > _______________________________________________ > > 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/20110825/6e9d121e/atta… > > > > ------------------------------ > > Message: 9 > Date: Thu, 25 Aug 2011 08:53:12 -0300 > From: Pedro Rafael <pedro.rafael.marinho(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] [OFF- Dúvida] Linguagem de programação para trabalhar > com estatística. > Message-ID: > <CAKwavr=E4ek4iJW-eObwQOvKWkGZx9DdvJNr+AcXxQGSZW5Ljw(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Pessoal qual seria uma outra boa linguagem de programação para se trabalhar > com estatística? Estudo C já faz algum tempo e vi ser bastante interessante > pelo fato de puder unir C com o R. Estava pensando em migrar para o C++, > será que é mais interessante? Existe outras linguagens de programação que > vocês indicam? > > -- > Saudações, > Pedro Rafael Diniz Marinho - Estatístico SES-PB > Currículo Lattes: > http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110825/7f0485af/atta… > > > > ------------------------------ > > Message: 10 > Date: Thu, 25 Aug 2011 08:53:28 -0300 > From: Pedro Rafael <pedro.rafael.marinho(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Pedro Rafael convidou você para abrir uma conta de > e-mail do Google > Message-ID: > <CAKwavrmoEjirnh+=uDN1Jx7ubXgYs_EAZxL9jWa4GASf++k9ng(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Estou usando o Gmail e acho que vale a pena conferir. Isto é um convite > para > você criar uma conta. > > > Você recebeu um convite para o Gmail! > > Pedro Rafael convidou você para abrir uma conta do Gmail. > > O Gmail é o serviço de e-mail gratuito do Google, desenvolvido com base na > ideia de que o e-mail pode ser intuitivo, eficiente e divertido. O Gmail > tem: > > *Menos spam* > Mantenha as mensagens indesejadas fora da sua caixa de entrada com a > tecnologia inovadora do Google. > > *Muito espaço* > Espaço de armazenamento suficiente para que você nunca precise excluir uma > mensagem. > > *Bate-papo integrado* > Converse por vídeo ou texto com Pedro Rafael e outros amigos em tempo real. > > *Acesso pelo celular* > Receba seus e-mails onde estiver com o Gmail no celular. > > Você pode até mesmo importar seus contatos e mensagens do Yahoo!, do > Hotmail, do AOL ou de qualquer outro webmail ou conta POP. > > Quando você criar sua conta, Pedro Rafael receberá uma notificação sobre > seu > novo endereço do Gmail para que vocês mantenham contato. Saiba > mais<http://mail.google.com/mail/help/intl/pt_BR/about.html>ou comece > agora< > http://mail.google.com/mail/a-3b4fd3d7b1-e4e8f2bf80-mBMo5xYWFYZNoKIRvN68yIo… > > > ! > Inscreva-se< > http://mail.google.com/mail/a-3b4fd3d7b1-e4e8f2bf80-mBMo5xYWFYZNoKIRvN68yIo… > > > > Google Inc. | 1600 Amphitheatre Parkway | Mountain View, California 94043 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110825/0878755e/atta… > > > > ------------------------------ > > Message: 11 > Date: Thu, 25 Aug 2011 08:02:45 -0400 > From: Gilson Sanchez <gilsonsch(a)gmail.com> > To: Grupo R-Br <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] [OFF TOPIC] - Learn R Webinar > Message-ID: > <CAJATvULeM74-TR9X3CyN1LJg3hi_Ky=1dVrned78Q9FG=HKbww(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Pode ser que alguem esteja precissando umas aulas complementares ou reforço > com o R, hoje terá mais um webinar, é minsitrado em inglês. > > Abç's > > ---------- Forwarded message ---------- > From: Paul H Geissler <phgeissler(a)usgs.gov> > Date: Wed, Aug 24, 2011 at 6:22 PM > Subject: Learn R Webinar Tomorrow > To: > > > > We now have our GoToWebinar license reinstated, but that meant that we have > TWO webinars scheduled for tomorrow. To avoid confusion, I canceled one of > them. However, that may add to the confusion, because you will get both a > reminder of the webinar and another saying that it is canceled. > > We will go with the webinar that Tom scheduled, rather than switching back > to the old one. If you have not registered, the link to register is * > https://www1.gotomeeting.com/register/511594025*< > https://www1.gotomeeting.com/register/511594025> > After you register, the system will send you a link to join. Tom, > again, > thank you very much for making the webinar available at such short notice. > > Sorry about all the confusion. > > On Thursday, August 25 we will cover Topic 13 - spatial statistics and > survival analysis, Crawley (2007) Chapters 24 and 25. The script is at > ftp://ftpext.usgs.gov/pub/cr/co/fort.collins/Geissler/LearnR/LearnR11a-13.R > > > Cheers, > Paul > ------------------------------------------------------------------- > Paul H. Geissler, Ph.D. > USGS Status & Trends of Biological Resources Program > Coordinator, National Park Monitoring Project > Assistant Program Coordinator > USGS Fort Collins Science Center > 2150 Centre Ave., Building C > Fort Collins, CO 80526-8118 > 970-226-9482, FAX 970-226-9452 > Paul_Geissler(a)usgs.gov > Do NOT send e-mail to PGeissler(a)usgs.gov (NOT me) > Home: 720-890-8785 Cell:303-815-9425 > Personal email: PaulGeisslerCo(a)gmail.com > > It is easy to lie with statistics. > It is hard to tell the truth without statistics. > Andrejs Dunkels, quoted by Maindonald & Braun > > > > > -- > Eng. Agr.Gilson Sánchez Chia, MSc. > > Laboratório de Fisiologia Vegetal > > Embrapa Amazônia Ocidental > > Fone: (92) 3303-7841 > --- > Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente. > Please consider the environment before printing this email. > Renew, Reduce and Recycle. The Planet Knows. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110825/b0f8c5fc/atta… > > > > ------------------------------ > > Message: 12 > Date: Thu, 25 Aug 2011 09:17:58 -0300 > From: Leandro Marino <leandromarino(a)leandromarino.com.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Plotar dados em Mapa > Message-ID: > <CAKSaaFnX3faWvdzgt2heC0KH8EPKip=3BnaxYJr4Vqi2T2nR6w(a)mail.gmail.com > > > Content-Type: text/plain; charset="windows-1252" > > Caros, > > > Mais uma vez venho pedir um caminho... Tenho um indicador por município que > gostaria que fosse ?plotado? em um mapa do Brasil aonde tivesse uma escala > de cores para cada um dos níveis deste indicador. Sei que existe alguma > forma de fazer, no entanto, nunca trabalhei com mapas no R e não sei nem > por > onde começar (Pacote, Função e principalmente MAPA). > > > Agradeço qualquer ajuda!!! > > > Um abraço > > Leandro > > > Atenciosamente, > Leandro Marino > http://www.leandromarino.com.br (Fotógrafo) > http://est.leandromarino.com.br/Blog (Estatístico) > Cel.: + 55 21 9845-7707 > Cel.: + 55 21 8777-7907 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110825/449bf3d6/atta… > > > > ------------------------------ > > Message: 13 > Date: Thu, 25 Aug 2011 10:02:08 -0300 > From: Pedro Rafael <pedro.rafael.marinho(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] [Dúvida] Invocar funções C++ no R > Message-ID: > <CAKwavrm=aWvvtRMrHzSAjPAjr4F7=ZbXi7DEw=+Cc3dU2W-gTg(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Como se procede para importar uma função C++ ou até mesmo C para o R? > > -- > Saudações, > Pedro Rafael Diniz Marinho - Estatístico SES-PB > Currículo Lattes: > http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110825/1bf2ef37/atta… > > > > ------------------------------ > > Message: 14 > Date: Thu, 25 Aug 2011 14:19:01 +0100 > From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] [Dúvida] Invocar funções C++ no R > Message-ID: > <CAO-arWMq0Z8r=oUR++ahXt5cKWuW=2qvZ=g6qO2Z0H3uS2WM0g(a)mail.gmail.com > > > Content-Type: text/plain; charset=ISO-8859-1 > > Sua funcao C++ precisa estar contida numa chamada extern... > > extern "C" { > void suaFuncao () ...... > } > > Depois de compila-la, vc a chama com .C() > > Sua funcao em C nao precisa do wrapper e pode ser compilada > diretamente e chamada via .C(). > > Para detalhes e tambem exemplos, veja o capitulo 5 do manual Writing R > Extensions. > > Nao faz muito tempo, eu postei aqui um link para um material que > tambem inclui exemplos de uso de codigo C em R... e inclui outra opcao > de chamada, .Call(). > > b > > > ------------------------------ > > Message: 15 > Date: Thu, 25 Aug 2011 10:25:39 -0300 > From: Pedro Rafael <pedro.rafael.marinho(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] [Dúvida] Invocar funções C++ no R > Message-ID: > <CAKwavr=eL9S0EnPku4o5YppZ3mr6Dy2vCZJVXnkvoc1cwuojcQ(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Obrigado Benilton pela ajuda. Vou dar uma olhada no material do CRAN e > tentar encontrar o material que você postou aqui na comunidade. > > Saudações, > Pedro Rafael Diniz Marinho > > Em 25 de agosto de 2011 10:19, beniltoncarvalho [via R-br] < > ml-node+3768192-710127126-223914(a)n4.nabble.com> escreveu: > > > Sua funcao C++ precisa estar contida numa chamada extern... > > > > extern "C" { > > void suaFuncao () ...... > > } > > > > Depois de compila-la, vc a chama com .C() > > > > Sua funcao em C nao precisa do wrapper e pode ser compilada > > diretamente e chamada via .C(). > > > > Para detalhes e tambem exemplos, veja o capitulo 5 do manual Writing R > > Extensions. > > > > Nao faz muito tempo, eu postei aqui um link para um material que > > tambem inclui exemplos de uso de codigo C em R... e inclui outra opcao > > de chamada, .Call(). > > > > b > > _______________________________________________ > > R-br mailing list > > [hidden email] <http://user/SendEmail.jtp?type=node&node=3768192&i=0> > > 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-Duvida-Invocar-funcoes-C-no-R-tp3768… > > To unsubscribe from R-br, click here< > http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib… > >. > > > > > > > > -- > Saudações, > Pedro Rafael Diniz Marinho - Estatístico SES-PB > Currículo Lattes: > http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110825/e2784a70/atta… > > > > ------------------------------ > > Message: 16 > Date: Thu, 25 Aug 2011 14:34:26 +0100 > From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] [Dúvida] Invocar funções C++ no R > Message-ID: > <CAO-arWMc=estcahNyV7oM8NnHuK1TULeZSnnLde8A9j0eh5D0A(a)mail.gmail.com > > > Content-Type: text/plain; charset=ISO-8859-1 > > Procure pelo thread chamado: "Pacotes e C, Fortran ou outros". > > > ------------------------------ > > Message: 17 > Date: Thu, 25 Aug 2011 10:45:56 -0300 > From: Pedro Rafael <pedro.rafael.marinho(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] [Dúvida] Invocar funções C++ no R > Message-ID: > <CAKwavrnPzPTcZkFYjs-3qoRU4mqLKhUBMBsVxrH2LYXnN4L9NQ(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > To olhando o exemplo abaixo: > > void convolve(double *a, int *na, double *b, int *nb, double *ab) > { > R_len_t i, j, nab = *na + *nb - 1; > > for(i = 0; i < nab; i++) > ab[i] = 0.0; > for(i = 0; i < *na; i++) > for(j = 0; j < *nb; j++) > ab[i + j] += a[i] * b[j]; > } > > called from R by > > conv <- function(a, b) > .C("convolve", > as.double(a), > as.integer(length(a)), > as.double(b), > as.integer(length(b)), > ab = double(length(a) + length(b) - 1))$ab > > Não estou conseguindo rodar. Devo gerar um arquivo .o ou o R tem algum > compilador por trás para ler o código? Estou perdido :( > > Pedro > > Em 25 de agosto de 2011 10:39, beniltoncarvalho [via R-br] < > ml-node+3768249-113061051-223914(a)n4.nabble.com> escreveu: > > > Procure pelo thread chamado: "Pacotes e C, Fortran ou outros". > > _______________________________________________ > > R-br mailing list > > [hidden email] <http://user/SendEmail.jtp?type=node&node=3768249&i=0> > > 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-Duvida-Invocar-funcoes-C-no-R-tp3768… > > To unsubscribe from R-br, click here< > http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib… > >. > > > > > > > > -- > Saudações, > Pedro Rafael Diniz Marinho - Estatístico SES-PB > Currículo Lattes: > http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110825/97855dd2/atta… > > > > ------------------------------ > > Message: 18 > Date: Thu, 25 Aug 2011 10:56:22 -0300 > From: Jakson Alves de Aquino <jalvesaq(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Plotar dados em Mapa > Message-ID: > <CAGBu4CP3Zg=ezdPAiRNwxw-+SHP20N_TE409N3i-SEZCsy73sQ(a)mail.gmail.com > > > Content-Type: text/plain; charset=UTF-8 > > 2011/8/25 Leandro Marino <leandromarino(a)leandromarino.com.br>: > > Mais uma vez venho pedir um caminho... Tenho um indicador por município > que > > gostaria que fosse ?plotado? em um mapa do Brasil aonde tivesse uma > escala > > de cores para cada um dos níveis deste indicador. Sei que existe alguma > > forma de fazer, no entanto, nunca trabalhei com mapas no R e não sei nem > por > > onde começar (Pacote, Função e principalmente MAPA). > > Tem um exemplo de como fazer isso no capítulo sobre mapas do R para > Cientistas Sociais: > > http://www.lepem.ufc.br/aquino.php > > Espero que lhe ajude... > > Jakson > > > ------------------------------ > > Message: 19 > Date: Thu, 25 Aug 2011 15:05:58 +0100 > From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] [Dúvida] Invocar funções C++ no R > Message-ID: > <CAO-arWOJXrNeYg4U3ZiGsP2F7QU05yg97W5YvEuFrAx2czVYwg(a)mail.gmail.com > > > Content-Type: text/plain; charset=ISO-8859-1 > > A primeira parte (tudo dentro do void convolve) e' uma funcao em C. > > Vc precisa te-la num arquivo como: convolve.c > > Depois, vc precisa ter compiladores C instalados no seu computador. > > Dai' compilar o codigo usando: R CMD SHLIB convolve.c > > Isso vai gerar um .so, que vc carregara' dentro do R usando dyn.load() > > Uma vez carregado no R, vc pode usar o .C(). > > Note que isso sera' um pouco mais complexo em Windows... E vc vai > precisar instalar um bundle chamado RTools, possivelmente setar > algumas variaveis de ambiente e substituir "R CMD" por Rcmd... mas > isso, alguem q use Windows eh quem tem condicoes de explicar em > detalhes. > > b > > > ------------------------------ > > Message: 20 > Date: Thu, 25 Aug 2011 07:17:29 -0700 (PDT) > From: Emmanuel Arnhold <emmanuelarnhold(a)yahoo.com.br> > To: R-br <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Agrupamento de Tocher no R > Message-ID: > <1314281849.70363.YahooMailNeo(a)web43143.mail.sp1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Alguém conhece função no R para construir o agrupamento de Tocher? > Aproveitando ... alguma função do R obtem a distância euclidiana média? > Conheço algumas que fazem apenas a distância euclidiana, como a vegdist() e > dist(). > > Obrigado por qualquer dica. > > Emmanuel > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110825/93ede2a1/atta… > > > > ------------------------------ > > Message: 21 > Date: Thu, 25 Aug 2011 11:38:30 -0300 > From: Leandro Marino <leandromarino(a)leandromarino.com.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Plotar dados em Mapa > Message-ID: > <CAKSaaF=r61P9LWAf2P4iYyh9Zi5vVOYD-7_OgMnR-j=ciEn7yA(a)mail.gmail.com > > > Content-Type: text/plain; charset="windows-1252" > > *Jackson, > > já foi de grande ajuda! Muito obrigado! > > Aproveitando, tenho os shapefiles do BR por Municipio e tenho outro por UF. > Gostaria de saber se existe uma possibilidade de plotar os resultados que > tenho por municipio, mas aparecendo as linhas das divisões do estados > apenas. É possível?! > > Seria o equivalente no exemplo de sua apostila a plotar o mapa do brasil > com > as cores das dos estados nas regiões mas apenas com a linha de divisão > entre > as regiões, ou seja, sem as divisões estaduais. > > > *Atenciosamente, > Leandro Marino > http://www.leandromarino.com.br (Fotógrafo) > http://est.leandromarino.com.br/Blog (Estatístico) > Cel.: + 55 21 9845-7707 > Cel.: + 55 21 8777-7907 > > > > Em 25 de agosto de 2011 10:56, Jakson Alves de Aquino > <jalvesaq(a)gmail.com>escreveu: > > > 2011/8/25 Leandro Marino <leandromarino(a)leandromarino.com.br>: > > > Mais uma vez venho pedir um caminho... Tenho um indicador por município > > que > > > gostaria que fosse ?plotado? em um mapa do Brasil aonde tivesse uma > > escala > > > de cores para cada um dos níveis deste indicador. Sei que existe > alguma > > > forma de fazer, no entanto, nunca trabalhei com mapas no R e não sei > nem > > por > > > onde começar (Pacote, Função e principalmente MAPA). > > > > Tem um exemplo de como fazer isso no capítulo sobre mapas do R para > > Cientistas Sociais: > > > > http://www.lepem.ufc.br/aquino.php > > > > Espero que lhe ajude... > > > > Jakson > > _______________________________________________ > > 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/20110825/b794e9c7/atta… > > > > ------------------------------ > > Message: 22 > Date: Thu, 25 Aug 2011 11:50:57 -0300 > From: Rodrigo Plei <rodrigo.plei(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Plotar dados em Mapa > Message-ID: > <CACs_BFOdc9bKRRTUN2O2+VZsNrB4C0JH-UqxMpd+Lit-0XjBow(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Pessoal, > > Aproveitando o gancho... Onde eu posso arrumar dados de mapas do litoral > sudeste-sul do Brasil (com as batimetrias de 100, 200, 500 e 1000 m) para > poltar com o R? > > Cheers, > > R. > > -- > =8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8--> > > Rodrigo Silvestre Martins, PhD > Bolsista FAPESP Pós-Doutorado > > Laboratório de Ecossistemas Pesqueiros (LabPesq) > Universidade de São Paulo, Instituto Oceanográfico > Praça do Oceanográfico, 191. Cidade Universitaria (sala 107-A/B) > Butantã - São Paulo/SP, Brasil > 05508-900 > Tel: +55 11 3091 6549 > Email: rodrigo.plei(a)gmail.com ; ocersm(a)lycos.com; rsmartins(a)usp.br > > CV Lattes: http://lattes.cnpq.br/5350064124902777 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110825/addd4e64/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 6, assunto 41 > **************************************** > -- Júlio Sílvio de Sousa Bueno Filho Departamento de Ciências Exatas Universidade Federal de Lavras Caixa Postal 3037 Lavras, MG - Brasil - 37200 000 tel: (+55) 35 3829 1369 fax: (+55) 35 3829 1371 "Natureza da gente não cabe em nenhuma certeza." João G. Rosa - Riobaldo (Grande Sertão: Veredas) "Homo sum; humani nihil a me alienum puto." Terêncio - Chremes (Heautontimoroumenos)

1 0

Como ler arquivo .ssd (S-Plus) no R
by Mauricio Cardeal 25 Ago '11

25 Ago '11

Olá caros RMasters. Estou tentando ler um dataset criado no S-Plus (não sei a versão) no R, mas ocorre mensagem de erro. Segundo o manual do provedor (http://www.who.int/childgrowth/software/en/) do dataset bastaria fazer: data.restore("bmianthro.sdd") depois de setado o diretório com setwd, mas essa orientação é para ser apliacada no S-Plus e não no R. Quando reproduzo a orientação no R, o resultado que sai no R é: > data.restore("bmianthro.sdd")Erro em ReadSdump(TRUE, " ") : S mode ‘bmianthro’ (near byte offset 50) not supportedAlém disso: Mensagens de aviso perdidas:In ReadSdump(TRUE, " ") : NAs introduzidos por coerção Segue o link do arquivo que postei no Data File Host: http://www.datafilehost.com/download-77a85580.html e o link do Readme: http://www.datafilehost.com/download-9c22749f.html Estou usando R versão 2.13.1 no ubuntu 11.04 Agradeço a ajuda. Mauricio Cardeal UFBA

2 3

Como ler arquivo .ssd (S-Plus) no R
by Mauricio Cardeal 25 Ago '11

25 Ago '11

Olá caros RMasters. Estou tentando ler um dataset criado no S-Plus (não sei a versão) no R, mas ocorre mensagem de erro. Segundo o manual do provedor (http://www.who.int/childgrowth/software/en/) do dataset bastaria fazer: data.restore("bmianthro.sdd") depois de setado o diretório com setwd, mas essa orientação é para ser apliacada no S-Plus e não no R. Quando reproduzo a orientação no R, o resultado que sai no R é: > data.restore("bmianthro.sdd")Erro em ReadSdump(TRUE, " ") : S mode ‘bmianthro’ (near byte offset 50) not supportedAlém disso: Mensagens de aviso perdidas:In ReadSdump(TRUE, " ") : NAs introduzidos por coerção Segue o link do arquivo que postei no Data File Host: http://www.datafilehost.com/download-77a85580.html e o link do Readme: http://www.datafilehost.com/download-9c22749f.html Estou usando R versão 2.13.1 no ubuntu 11.04 Agradeço a ajuda. Mauricio Cardeal UFBA

1 0

Objetos da Classe SpatialPolygonsDataFrame
by Leonel Lucena Rodrigues 25 Ago '11

25 Ago '11

Boa Tarde, Estou com duas dúvidas sobre objetos da classe SpatialPolygonsDataFrame. Tenho arquivos em formato shapefile, que é considerado pelo R como um objeto da classes polígono espacial, e este tem alguns slot's como "labpt" e "area". Então venho pedir ajuta em duas questões. 1º Qual a função, e o pacote, que me fornecerá apenas o slot "coords" ? Sei existe uma forma de fazer com que o programa selecione as coordenadas de polígono, atualmente utilizo a função scan(), mas gostaria de utilizar uma outra. 2º Caso venha a conseguir alguma estimativa, Kernel e.g., existirá um viés se os dados não tiverem projeção ? Na cartografia, mais especificamente no ArcMAP, utilizo o sistema de coordenadas adotado atualmente pelo IBGE (SIRGAS 2000), no R esse processo está em desenvolvimento. Abraço a todos.

2 6

Join my network on LinkedIn
by Lucas Cusinato via LinkedIn 25 Ago '11

25 Ago '11

LinkedIn ------------ Lucas Cusinato requested to add you as a connection on LinkedIn: ------------------------------------------ I'd like to add you to my professional network on LinkedIn. Accept invitation from Lucas Cusinato http://www.linkedin.com/e/-33g73l-grrvtt86-6z/vhj1xG42iBdUFK5d_HMEGL4iFMC51… View invitation from Lucas Cusinato http://www.linkedin.com/e/-33g73l-grrvtt86-6z/vhj1xG42iBdUFK5d_HMEGL4iFMC51… ------------------------------------------ DID YOU KNOW you can be the first to know when a trusted member of your network changes jobs? With Network Updates on your LinkedIn home page, you'll be notified as members of your network change their current position. Be the first to know and reach out! http://www.linkedin.com/ -- (c) 2011, LinkedIn Corporation

1 0

Agrupamento de Tocher no R
by Emmanuel Arnhold 25 Ago '11

25 Ago '11

Alguém conhece função no R para construir o agrupamento de Tocher? Aproveitando ... alguma função do R obtem a distância euclidiana média? Conheço algumas que fazem apenas a distância euclidiana, como a vegdist() e dist(). Obrigado por qualquer dica. Emmanuel

1 0

[OFF TOPIC] - Learn R Webinar
by Gilson Sanchez 25 Ago '11

25 Ago '11

Pode ser que alguem esteja precissando umas aulas complementares ou reforço com o R, hoje terá mais um webinar, é minsitrado em inglês. Abç's ---------- Forwarded message ---------- From: Paul H Geissler <phgeissler(a)usgs.gov> Date: Wed, Aug 24, 2011 at 6:22 PM Subject: Learn R Webinar Tomorrow To: We now have our GoToWebinar license reinstated, but that meant that we have TWO webinars scheduled for tomorrow. To avoid confusion, I canceled one of them. However, that may add to the confusion, because you will get both a reminder of the webinar and another saying that it is canceled. We will go with the webinar that Tom scheduled, rather than switching back to the old one. If you have not registered, the link to register is * https://www1.gotomeeting.com/register/511594025*<https://www1.gotomeeting.com/register/511594025> After you register, the system will send you a link to join. Tom, again, thank you very much for making the webinar available at such short notice. Sorry about all the confusion. On Thursday, August 25 we will cover Topic 13 - spatial statistics and survival analysis, Crawley (2007) Chapters 24 and 25. The script is at ftp://ftpext.usgs.gov/pub/cr/co/fort.collins/Geissler/LearnR/LearnR11a-13.R Cheers, Paul ------------------------------------------------------------------- Paul H. Geissler, Ph.D. USGS Status & Trends of Biological Resources Program Coordinator, National Park Monitoring Project Assistant Program Coordinator USGS Fort Collins Science Center 2150 Centre Ave., Building C Fort Collins, CO 80526-8118 970-226-9482, FAX 970-226-9452 Paul_Geissler(a)usgs.gov Do NOT send e-mail to PGeissler(a)usgs.gov (NOT me) Home: 720-890-8785 Cell:303-815-9425 Personal email: PaulGeisslerCo(a)gmail.com It is easy to lie with statistics. It is hard to tell the truth without statistics. Andrejs Dunkels, quoted by Maindonald & Braun -- Eng. Agr.Gilson Sánchez Chia, MSc. Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 --- Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente. Please consider the environment before printing this email. Renew, Reduce and Recycle. The Planet Knows.

1 0

Pedro Rafael convidou você para abrir uma conta de e-mail do Google
by Pedro Rafael 25 Ago '11

25 Ago '11

Estou usando o Gmail e acho que vale a pena conferir. Isto é um convite para você criar uma conta. Você recebeu um convite para o Gmail! Pedro Rafael convidou você para abrir uma conta do Gmail. O Gmail é o serviço de e-mail gratuito do Google, desenvolvido com base na ideia de que o e-mail pode ser intuitivo, eficiente e divertido. O Gmail tem: *Menos spam* Mantenha as mensagens indesejadas fora da sua caixa de entrada com a tecnologia inovadora do Google. *Muito espaço* Espaço de armazenamento suficiente para que você nunca precise excluir uma mensagem. *Bate-papo integrado* Converse por vídeo ou texto com Pedro Rafael e outros amigos em tempo real. *Acesso pelo celular* Receba seus e-mails onde estiver com o Gmail no celular. Você pode até mesmo importar seus contatos e mensagens do Yahoo!, do Hotmail, do AOL ou de qualquer outro webmail ou conta POP. Quando você criar sua conta, Pedro Rafael receberá uma notificação sobre seu novo endereço do Gmail para que vocês mantenham contato. Saiba mais<http://mail.google.com/mail/help/intl/pt_BR/about.html>ou comece agora<http://mail.google.com/mail/a-3b4fd3d7b1-e4e8f2bf80-mBMo5xYWFYZNoKIRvN68yIo…> ! Inscreva-se<http://mail.google.com/mail/a-3b4fd3d7b1-e4e8f2bf80-mBMo5xYWFYZNoKIRvN68yIo…> Google Inc. | 1600 Amphitheatre Parkway | Mountain View, California 94043

1 0

Re: [R-br] Representação usando expression()
by Alexandre Santos 25 Ago '11

25 Ago '11

Obrigado Benilton, Resolveu o meu problema e ficaram muito mais elegante os argumentos também, Abraço, Alexandre dos Santos Engenheiro Florestal, MSc. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 3829-5122 Tomei a liberdade de simplificar um pouco sua expressao (e' equivalente ao modelo que vc descreve): txt = expression(frac(1, 1+e^{0.9588-0.0841*x+0.0010*x^2})) plot(1, type='n') legend(1, 1, txt, bty='n') b

1 0

Representação usando expression()
by Alexandre Santos 25 Ago '11

25 Ago '11

Bom Bom dia pessoal, Gostaria de saber s é possível fazer com que em: text(33,0.9,expression(italic("y")==italic("exp")(-"0.9588"+"0.0841"*italic("*x")-"0.0010"*italic("*x"))/(1+italic("exp")(-"0.9588"+"0.0841"*italic("*x")-"0.0010"*italic("*x")))),cex=0.8) A equação aparecesse no gráfico como: ("exp")(-"0.9588"+"0.0841"*italic("*x")-"0.0010"*italic("*x")) sobre (1+italic("exp")(-"0.9588"+"0.0841"*italic("*x")-"0.0010"*italic("*x") e não separado por / , isso é possível? Obrigado, Alexandre dos Santos Engenheiro Florestal, MSc. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 3829-5122

2 1

Otimizar tempo de processamento
by Tiago Vieira 25 Ago '11

25 Ago '11

Estou tentando rodar o programa de simulação de percolação (rede de bethe), mas esta demandando um tempo de processamento muito alto, pois há preciso de muitas iterações. O computador que estou utilizando tem processador i3 e 4gb de memória. Alguém pode me dar uma dica de como posso diminuir o tempo de processamento? Abaixo segue o código: programa <- function(f,n,m) { # f: número máximo de filhos # n: número de gerações # m: número de iterações S <- matrix(0,101,2) for (p in seq(0,1,by=0.01)) { pp <- 0 for (g in 1:m) { y <- 1 for (k in 1:n) { y <- sum(rbinom(y,f,p)) if (y == 0) break } if (y > 0) pp <- pp + 1 } S[p/0.01+1,1] <- p S[p/0.01+1,2] <- pp/m } return(S) } prog <- programa(2,100,1000) plot(prog, xlab = "p", ylab = "Probabilidade de Percolação") Desde já agradeço pela atenção! Tiago Vieira.

6 12

Problema de memória: MATRIZES GIGANTES
by Márcio Luís Moreira de Souza 24 Ago '11

24 Ago '11

Gente... Tô com outro problema. Agora é em relação à memória do R. Não tô conseguindo gerar matrizes grandes (no meu caso são gigantes, na casa dos milhões de linhas por dezenas de colunas). Alguma dica?? Dá pra aumentar a memória utilizada no R (tenho um processador I5)?? E em computadores mais "fraquinhos", como faço?? Abraços!! -- Prof. Márcio Souza, M.Sc. (DE/UFJF) www.profmsouza.blogspot.com http://lattes.cnpq.br/4578008002785684

5 10

Invitation to connect on LinkedIn
by Marcelo Costa via LinkedIn 24 Ago '11

24 Ago '11

LinkedIn ------------ Marcelo Costa requested to add you as a connection on LinkedIn: ------------------------------------------ Fábio, I'd like to add you to my professional network on LinkedIn. - Marcelo Accept invitation from Marcelo Costa http://www.linkedin.com/e/-33g73l-grqxst4r-3t/vhj1xG42iBdUFK5d_HMEGL4iFMC51… View invitation from Marcelo Costa http://www.linkedin.com/e/-33g73l-grqxst4r-3t/vhj1xG42iBdUFK5d_HMEGL4iFMC51… ------------------------------------------ DID YOU KNOW LinkedIn can help you find the right service providers using recommendations from your trusted network? Using LinkedIn Services, you can take the risky guesswork out of selecting service providers by reading the recommendations of credible, trustworthy members of your network. http://www.linkedin.com/e/-33g73l-grqxst4r-3t/svp/inv-25/ -- (c) 2011, LinkedIn Corporation

1 0

Re: [R-br] [OFF] Web Framework usando R
by Pedro Rafael 24 Ago '11

24 Ago '11

Aproveitando o mesmo tópico, qual seria uma boa linguagem para a pessoa aprender fora o R para trabalhar com estatística? Em 24 de agosto de 2011 17:51, Pedro Rafael <pedro.rafael.marinho(a)gmail.com>escreveu: > Aproveitando o mesmo tópico, qual seria uma boa linguagem para a pessoa > aprender fora o R para trabalhar com estatística? > > > > 2011/8/24 beniltoncarvalho [via R-br] < > ml-node+3766456-521461530-223914(a)n4.nabble.com> > > algo alem do que ha' descrito na documentacao? >> >> http://cran.r-project.org/doc/FAQ/R-FAQ.html#R-Web-Interfaces >> >> E, aqui, algo que possa ser relevante (eh um bom comeco, IMHO): >> >> >> http://nsaunders.wordpress.com/2010/04/19/getting-your-web-application-and-… >> >> b >> _______________________________________________ >> R-br mailing list >> [hidden email] <http://user/SendEmail.jtp?type=node&node=3766456&i=0> >> 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-OFF-Web-Framework-usando-R-tp3766338… >> To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…>. >> >> > > > > -- > Saudações, > Pedro Rafael Diniz Marinho - Estatístico SES-PB > Currículo Lattes: > http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 > > -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

1 0

Factors no R
by Leandro Marino 24 Ago '11

24 Ago '11

Caros, Em algum ponto do tempo houve uma mudança na função factor. ### no passado: R.version # _ #platform i386-pc-mingw32 #arch i386 #os mingw32 #system i386, mingw32 #status #major 2 #minor 9.0 #year 2009 #month 04 #day 17 #svn rev 48333 #language R #version.string R version 2.9.0 (2009-04-17) set.seed(1) x <- round(runif(150,0,10),0) table(x) #x # 0 1 2 3 4 5 6 7 8 9 10 # 4 14 15 13 16 19 16 18 15 14 6 table(as.character(factor(x,levels=c(0:10),labels=c(rep('Tipo 1',3), rep('Tipo 2',4),rep('Tipo 3',2),'Tipo 4','Tipo 5')))) #Tipo 1 Tipo 2 Tipo 3 Tipo 4 Tipo 5 # 33 64 33 14 6 R.version # _ #platform i386-pc-mingw32 #arch i386 #os mingw32 #system i386, mingw32 #status #major 2 #minor 13.1 #year 2011 #month 07 #day 08 #svn rev 56322 #language R #version.string R version 2.13.1 (2011-07-08) set.seed(1) x <- round(runif(150,0,10),0) table(x) #x # 0 1 2 3 4 5 6 7 8 9 10 # 4 14 15 13 16 19 16 18 15 14 6 table(as.character(factor(x,levels=c(0:10),labels=c(rep('Tipo 1',3), rep('Tipo 2',4),rep('Tipo 3',2),'Tipo 4','Tipo 5')))) # #Tipo 1 Tipo 2 Tipo 3 Tipo 4 Tipo 5 # 33 64 33 14 6 #Mensagens de aviso perdidas: #In `levels<-`(`*tmp*`, value = c("Tipo 1", "Tipo 1", "Tipo 1", "Tipo 2", : # não se permite mais níveis duplicados em fatores Ele até faz mas com mensagem de erro. Como poderia fazer isto sem ter o erro!? Atenciosamente, Leandro Lins Marino Centro de Avaliação Fundação CESGRANRIO Rua Santa Alexandrina, 1011 - 2º andar Rio de Janeiro, RJ - CEP: 20261-903 ( (21) 2103-9600 R.:236 * <mailto:leandro@cesgranrio.org.br> leandro(a)cesgranrio.org.br P Antes de imprimir pense em sua responsabilidade e compromisso com o MEIO AMBIENTE Esta mensagem, incluindo seus anexos, pode conter informacoes privilegiadas e/ou de carater confidencial, nao podendo ser retransmitida sem autorizacao do remetente. Se voce nao e o destinatario ou pessoa autorizada a recebe-la, informamos que o seu uso, divulgacao, copia ou arquivamento sao proibidos. Portanto, se você recebeu esta mensagem por engano, por favor, nos informe respondendo imediatamente a este e-mail e em seguida apague-a.

3 6

[OFF] Web Framework usando R
by Pedro Rafael 24 Ago '11

24 Ago '11

Pessoal não sei se concordam mas a informática "atualmente" foca muito em Web. Muitas aplicações vem sendo feita para Web. Muitas empresas monitoram suas estatísticas de forma online via shmartphones, tablets via Web. Existe algum Framework para desenvolvimento de aplicações Web que utilize o R como linguagem principal? -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

3 2

(Resolvido) aov() e anova() no agricolae
by Alexandre Santos 24 Ago '11

24 Ago '11

Obrigado Ivan, Aumentou meu grau de confiança nas decisões que estou tomando, Abraço, Alexandre dos Santos Engenheiro Florestal, MSc. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 3829-5122 Alexandre, Até onde eu sei, o HSD.test é utilizado para erros normais. Mais você pode utilizar sem problemas o seu modelo glm com a função glht para comparar as médias. Não vejo problemas também, você utilizar tais análises para uma variável e tais análises para outra variável, desde que esteja justificada o uso de tais metodologias. Allaman (S,f,P) \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-5076 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. Alexandre dos Santos Engenheiro Florestal, MSc. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 3829-5122

1 0

ajuda para um scrip
by Mauro Sznelwar 24 Ago '11

24 Ago '11

Tenho este algoritmo p=1, N=0, c=e^-lambda Repita N=N+1 gere Ui atualize p<c.e^-lambda Faça X=N-1 Estou fazendo esta rotina abaixo, mas o resultado não está consistente! Qual o melhor jeito? p=1 n=0 lambda=2 c=exp(-lambda) for(n in 1:100){ u=runif(n) p=p*u p<c } x=n-1

4 10

Re: [R-br] aov() e anova() no agricolae
by Ivan Bezerra Allaman 24 Ago '11

24 Ago '11

Alexandre, Até onde eu sei, o HSD.test é utilizado para erros normais. Mais você pode utilizar sem problemas o seu modelo glm com a função glht para comparar as médias. Não vejo problemas também, você utilizar tais análises para uma variável e tais análises para outra variável, desde que esteja justificada o uso de tais metodologias. Allaman (S,f,P) \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-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

Re: [R-br] aov() e anova() no agricolae
by Ivan Bezerra Allaman 24 Ago '11

24 Ago '11

Alexandre, Primeiro nos tira uma dúvida. Porque você está avaliando n.ovulos por geração dentro de cada casal, sendo que você não avaliou a interação entre "tratamentos e geração"?. Não seria prudente você avaliar primeiro se há uma interação, para depois você avaliar o fator "geração" dentro de cada casal? Abraço! Allaman (S,f,P) \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-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

2 1

MANDANDO E-MAIL PELO R...
by Márcio Luís Moreira de Souza 24 Ago '11

24 Ago '11

Gente, Sei que já pode ter sido discutido antes, mas... Devo ter perdido o bonde. Alguém já conseguiu mandar e-mail pelo R?? Tipo: quando uma simulação acaba e tal, gostaria de ser notificado por e-mail. Já tentei a função *sendmail*... Sem sucesso. Eu uso GMAIL. Alguma especificação para esse tipo de e-mail? Abraços!! -- Prof. Márcio Souza, M.Sc. (DE/UFJF) www.profmsouza.blogspot.com http://lattes.cnpq.br/4578008002785684

4 5

aov() e anova() no agricolae
by Alexandre Santos 24 Ago '11

24 Ago '11

Bom Dia pessoal, Estou realizando alguns testes de Tukey com ajuda do pacote agricolae para algumas variáveis que possuem distribuição normal, então fiz para uma destas: # ##N.ovulos anova.o<-aov(n.ovulos~factor(tratamento)+factor(geracao)+factor(rep),data=dados6) summary(anova.o) # ##Tratamento 1 casal nas geracoes with(subset(dados6,tratamento == "1casal"),HSD.test(n.ovulos,geracao, DFerror=df.residual(anova.o),MSerro=deviance(anova.o)/df.residual(anova.o), group = TRUE)) ##Tratamento 10 casais nas geracoes with(subset(dados6,tratamento == "10casal"),HSD.test(n.ovulos,geracao, DFerror=df.residual(anova.o),MSerro=deviance(anova.o)/df.residual(anova.o), group = TRUE)) Ate ai tudo bem, porém ... tenho algumas outras variáveis explicativas dentro do mesmo experimento como razão sexual por exemplo, que é uma proporção (número de fêmeas/número de fêmeas+número de machos) que não tem distribuição normal e fiz o ajuste por GLM com distribuição de erros binomial. E agora gostaria de saber se não existe alguma forma de se utilizar HSD.test () do agricolae utilizando as informações da anova() do meu modelo binomial de maneira análoga ao aov() que utilizo para os dados com distribuição normal ou vou ter que tratar os dados com outras distribuições de erros com o glht() do multcomp? Ou deveria fazer toda a abordagem do problema pegando as variáveis que têm normalidade e utilizar distribuição gaussiana via glm() e tratar todo o experimento através de contrastes de modelo via multcomp(), visando a uniformizar a abordagem estatística dada ao mesmo experimento? Obrigado, Alexandre dos Santos Engenheiro Florestal, MSc. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 3829-5122

1 0

Re: [R-br] [OFF-TOPIC] Mudar default da extensão dos gráficos gerados pelo sweave
by Ivan Bezerra Allaman 23 Ago '11

23 Ago '11

Walmes, Vide pergunta novamente!! A solução ideal ficou a do Benilton. Faria, a sua solução é uma alternativa interessante, pois poderia gerar um gráfico png por exemplo para aquela(s) figura(s) em especial. Mais no caso de querermos gerar todas as figuras num formato específico, aí já ficaria complicado (pensemos nuns 40 gráficos por exemplo). Eu geralmente já utilizo um, \begin{figure}[H] <<fig=T>>= plot(1:10,1:10) @ \caption{} \end{figure} Com a sugestão do Benilton, ao compilar o arquivo, os gráficos gerados são em png ao invés do pdf (default), resolvendo o problema. O único inconveniente, é que teria que fazer isso direto do terminal. O ideal, é que o vim tivesse uma opção para mudar a extensão dos gráficos gerados. De qualquer maneira, obrigado a todos pelas sugestões! Agora, aproveitando o ensejo, outra pergunta. Tem algum comando específico, para que, ao compilar o arquivo sweave, os gráficos gerados vá para uma pasta específica diferente daquela onde está o arquivo.Rnw? Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

3 3

Declarar matriz Z modelos mistos
by Eder David Borges da Silva 23 Ago '11

23 Ago '11

Pessoal, dado um modelo y = Xb + Zu + e , modelo misto com Xb parte fixa e Zu parte aleatoria, estou tentando entender melhor a analise feita neste documento, http://www.infoteca.cnptia.embrapa.br/bitstream/doc/883425/1/Doc210.pdf , pagina 54, ao meu entendimento é um modelo misto, onde a mudança é a forma que é construída a matriz Z, que como o autor faz é produzido um efeito aleatório para cada marcador. desta forma gostaria de saber se alguém ja declarou um matriz Z em uma analise de modelo misto no R em especifico na lmer, essa analise é GWS (Seleção genômica ampla), caso alguém tenha familiaridade ficaria grato de umas dicas. dados <- data.frame(ind=c(1:5), diametro=c(9.87,14.48,8.91,14.64,9.55), M1=c(2,1,0,1,1), M2=c(0,1,2,0,0), M3=c(0,0,0,1,0), M4=c(0,0,0,0,1), M5=c(2,1,0,1,1), M6=c(0,1,0,0,1), M7=c(0,0,2,0,0)) dados require(lme4) Z <- model.matrix(~-1+M1+M2+M3+M4+M5+M6+M7,dados) Z ### Esta é a matriz que querro por no modelo m0 <- lmer(diametro~1+(1|?????),dados) summary(m0) ranef(m0) fixef(m0) # Meu interesse e ter os componentes de variância e o efeito fixo(Media geral), e os efeitos aleatórios (7, uma para cada marcador) OBS: Creio que não chegaremos ao valor exato do autor, visto que ele fixou (sigma.erro/(sigma.g/n)) =1. Obrigado

2 4

covariancia de beta
by Ricardo Pereira 23 Ago '11

23 Ago '11

Estou tentando fazer uns exercícios de outro programa no R. Como multiplico a matrix com o vetor? sendo que %*% multiplica matrizes............ sendo a matrix XPXINV: 6.867140 -2.001685 -2.001685 0.03283 proc IML MSE = 19.125957; /* Quadrado medio do erro */ varb = xpxinv*mse; /* 'Matriz de covariancia das estimativas de beta*/ Desde já agradeço ao pessoal da lista Ricardo Pereira

3 9

Mensagem de erro
by Carlos Mendonça 23 Ago '11

23 Ago '11

Caros, tenho um arquivo (arq1) com as variáveis codigo, valor, mes e ano. > str(arq1) 'data.frame': 1108685 obs. of 5 variables: $ codigo: int 10059020 10059020 10030000 10059020 10059020 10059020 10059020 10030000 10059020 10059020 ... $ valor : num 97159 39938 144 25002 800 ... $ mes : num 1 1 1 1 1 1 1 1 2 2 ... $ ano : chr "2009" "2009" "2009" "2009" ... $ trim : chr "1" "1" "1" "1" ... > head(arq1) codigo valor mes ano 1 10059020 97158.59 1 2009 2 10059020 39938.15 1 2009 3 10030000 143.95 1 2009 4 10059020 25002.10 1 2009 5 10059020 800.11 1 2009 6 10059020 263.32 1 2009 Desejo atribuir a variável trim o valor 1, 2, 3 ou 4 dependendo do mes. Só que o meu código está gerando a seguinte mensagem de erro: > if (arq1$mes<4) + { + arq1$trim = "1" + }else{ if(arq1$mes<7) { + arq1$trim = "2" + }else{ if(arq1$mes<10){ + arq1$trim = "3" + }else{ + arq1$trim = "4" + } + } + } Mensagens de aviso perdidas: In if (arq1$mes < 4) { : a condição tem comprimento > 1 e somente o primeiro elemento será usado Mensagens de aviso perdidas: In if (arq1$mes < 4) { : a condição tem comprimento > 1 e somente o primeiro elemento será usado E atribuindo a trim o valor 1 para todas as observações. > head(arq1) codigo valor mes ano trim 1 10059020 97158.59 1 2009 1 2 10059020 39938.15 1 2009 1 3 10030000 143.95 1 2009 1 4 10059020 25002.10 1 2009 1 5 10059020 800.11 1 2009 1 6 10059020 263.32 1 2009 1 O que fazer? Carlos Mendonça.

2 3

Atualizar biblioteca
by Andre Oliveira 23 Ago '11

23 Ago '11

Bom dia, eu gostaria de atualizar apenas uma biblioteca fiz o seguinte. update.packages(ask="brglm") Mas atualizou várias bibliotecas. Como devo proceder para atualizar uma apenas uma biblioteca? André Oliveira Souza

2 1

Curva da densidade truncada
by Mauro Sznelwar 23 Ago '11

23 Ago '11

Fiz esta rotina para o domínio 1,infinito. Porém como eu faço para a curva da densidade começar no 1 assim como no histograma? d2=runif(100)^-1/alfa d3=subset(d2,d2>=1) hist(d3,probability=T,nclass=60,col='green') dens=density(d3,n=200,na.rm=T) points(dens,type='l',col=2) abline(v=1,lty=4)

2 1

Re: [R-br] pacote instalado?
by Ivan Bezerra Allaman 22 Ago '11

22 Ago '11

Para verificar quais pacotes estão instalados na sua máquina, basta: help.start() (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-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

3 2

Legenda - expression
by Luís Gustavo 22 Ago '11

22 Ago '11

Boa noite a todos, tenho um gráfico abaixo e gostaria de saber se existe a possibilidade de colocar um subscrito no beta que está na legenda. Como por exemplo um beta_0. require(sn) x=seq(-5,5,by=0.01) a=dsn(x,0,1,0) b=dsn(x,0,1,-1) plot(x,a,type="l",ylim=c(0,0.8),xlim=c(-4,4),lwd=2.5,xlab="Z",ylab=expression(paste("Densidades NA(",lambda,")")),cex.lab=1.3) lines(x,b,lwd=2.5,lty=2,col=2) lines(x,c,lwd=2.5,lty=3,col=3) legend(1.5,0.79,expression(paste(beta,"=",0)),col=1,lty=1,box.lty=0) legend(1.5,0.75,expression(paste(beta,"=",-1)),col=2,lty=2,box.lty=0) Abraços -- Luís Gustavo Silva e Silva

2 2

pacote instalado?
by Luis Iván Ortiz Valencia 22 Ago '11

22 Ago '11

Existe um comando para verificar se um pacote ja esta instalado no R? atte IVAN

2 1

(sem assunto)
by Edson Lira 22 Ago '11

22 Ago '11

Bom dia a todos, tenho a estrutura abaixo: Descrição Quantidade Mês lab 1 17 ALFA HIDROXI PROGESTERONA 1 Jan H. Pardini 2 17 ALFA HIDROXI PROGESTERONA 1 Fev H. Pardini 3 17 ALFA HIDROXI PROGESTERONA 1 Fev H. Pardini 4 17 ALFA HIDROXI PROGESTERONA 2 Jan H. Pardini 5 ACIDO URICO 1 Jan Bioquímica 6 ACIDO URICO 1 Fev Bioquímica 7 ACIDO URICO 5 Jan Bioquímica 8 ACIDO URICO 1 Fev Bioquímica 9 ACIDO URICO 10 Jan Bioquímica 10 ACIDO URICO 1 Jan Bioquímica Com o comando: tapply(acj$Quantidade,acj$lab,sum,na.rm=T) Obtenho a soma das quantidades de exames por laboratório. Gostaria de buscar a soma das quantidades por laboratório e por Mês. O que devo acrescentar no comando, ou qual o comando para isso? Muito obrigado. Edson Lira Estatístico Manaus-Amazonas

4 3

Re: [R-br] Coeficiente de correlação em glm é correto?
by Alexandre Santos 22 Ago '11

22 Ago '11

Ok Benilton, Atendi a sua explicação e já corrigi. Agora estou utilizando os valores não diretamente em proporção, mais em duas coluna com número de insetos mortos e na outra o numero de insetos vivos. Mais agora existe uma outra questão que me inquieta e vou pedir sua ajuda. Existem alguns ajustes que faço e tenho problema de superdispersão, mais quando faço, apenas de maneira exploratória, alguma transformação dos dados como log(x+1), principamente quando utilizo a família Poisson eu tenho que a relação deviance residual/ df residual fica <1 e o Q-Q plot melhora também. Onde coloco a questão, o emprego de distribuições de erros mais adequadas aos dados não é justamente para fugir de transformações que tentavam fazer com que os dados mimetizassem normalidade e na verdade o que acontece com os dados quando faço esse tipo de transformação? Exemplificando: test<-rpois(100,30) m1<-glm(test~1,family="quasipoisson") summary(m1) plot(m1) m2<-glm((log(test+1))~1,family="quasipoisson") summary(m2) plot(m2) Obrigado, Alexandre dos Santos Engenheiro Florestal, MSc. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 3829-5122

2 1

Lista de discussão em inglês sobre R
by . . 22 Ago '11

22 Ago '11

Bom dia pessoal, Não faz muito que comecei a usar R e frequentar esta lista de discussão, que tem me ajudado a resolver várias dúvidas. Costumo acompanhar pelo nabble.com. As vantagens de um forum em pt_br são claras mas um inconveniente que vejo é que o número de usuários atingido é menor que em um forum internacional. Sobre Linux, um forum que gosto bastante é o linuxquestions.org e da mesma forma, sobre C/C++, gosto do cplusplus.com. Então, gostaria de perguntar quais os forums em inglês sobre R mais utilizados?

2 1

Coeficiente de correlação em glm é correto?
by Alexandre Santos 22 Ago '11

22 Ago '11

Boa noite pessoal, Venho colocar uma questão para os colegas estatísticos, até onde sei os modelos ajustados com GLM não são representados com o R^2. Porém, normalmente em revistas de entomologia (minha área) contam com editores da velha escola estatística que ficam batendo o pé para se colocar o R2 em todos os modelos ajustados, onde coloco a questão esta correto?Bom, tentei realizar o cálculo do R^2 usando aquelas abordagens que inúmeras vezes apareceram nas listasdo grupo do R sobre R^2 para modelos não lineares e fiz o seguinte cálculo: ##Ajustando o modelo com distribuição de erros binomial > m.model4<-glm((par/100)~densidade+I(densidade^2)+I(densidade^3), data=dados2,family="quasibinomial") > SQEm1<-summary(m.model4)$deviance^2*summary(m.model4)$df.residual#Soma de quadrados do erro > SQEm1 [1] 2387.366 > SQTm1<-var(par)*(length(par)-1)#Soma de quadrados total corrigida > R2<-1-SQEm1/SQTm1 > R2 [1] 0.8679119 Diante do problema apresentado isso esta correto? Obrigado, Alexandre dos Santos Engenheiro Florestal, MSc. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 3829-5122

3 4

comando para salvar o Editor do R
by Luis Iván Ortiz Valencia 22 Ago '11

22 Ago '11

Bom dia Existe algum comando no R para salvar o Editor do R a partir do console? atte IVAN

1 0

d, p, q , r - Quais as relações matemáticas entre estas funções
by . . 22 Ago '11

22 Ago '11

Olá a todos da lista, Gostaria de uma referência para estudar como são obtidas as funções d, p, q, r. A ideia é implementar as funções p, q, r à partir de uma função p particular. Obrigado

2 2

geeglm não reconhece os clusters...
by Paula Mendes Luz 21 Ago '11

21 Ago '11

oi pessoal alguém já usou a funcao geeglm do geepack? estou tentando usar e apesar de nao receber um erro, algo está errado pois os clusters nao sao reconhecidos: (o banco temp tem diversos clusters na var patient que tem tamanho 1-9) > geeglm(se ~ 1+offset(log(dfu/365.25)), data = temp, family = poisson(), id = patient, corstr = "independence") Call: geeglm(formula = se ~ 1 + offset(log(dfu/365.25)), family = poisson(), data = temp, id = patient, corstr = "independence") Coefficients: (Intercept) -1.671535 Degrees of Freedom: 100 Total (i.e. Null); 99 Residual Scale Link: identity Estimated Scale Parameters: [1] 4.447461 Correlation: Structure = independence Number of clusters: 1 Maximum cluster size: 100 Warning message: In diff(as.numeric(id)) : NAs introduced by coercion a maneira de especificar o modelo é igual ao glm e usando o glm o modelo estima corretamente: > fit2 <- glm(se ~ 1+offset(log(dfu/365.25)), data = bpal, family = poisson()) > summary(fit2) Call: glm(formula = se ~ 1 + offset(log(dfu/365.25)), family = poisson(), data = bpal) Deviance Residuals: Min 1Q Median 3Q Max -0.689 -0.688 -0.688 -0.688 13.198 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) -1.44036 0.01003 -143.6 <2e-16 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for poisson family taken to be 1) Null deviance: 61241 on 45028 degrees of freedom Residual deviance: 61241 on 45028 degrees of freedom AIC: 70587 Number of Fisher Scoring iterations: 7 entretanto, eu preciso corrigir pelas medidas repetidas nos pacientes... logo, preciso usar geeglm... se alguém tiver uma dica ou souber de alguma outra forma de rodar um GEE, favor orientar. obrigada, paula.

1 0

Pedindo ajuda para simulação em R
by Mauro Sznelwar 21 Ago '11

21 Ago '11

1) Este algoritmo da distribuição geométrica(p) Conhecendo a cdf: F(x)=1-(1-p)^x+1 F(x-1)=1-(1-p)^x Usando o método da transformada inversa F(x-1)<R<=F(x) 1-(1-p)^x<R<=1-(1-p)^x+1 (1-p)^x+1<R-1<=1-(1-p)^x ln(1-R)/ln(1-p)-1<=X<ln(1-R)/ln(1-p) X=[(ln(1-R)/ln(1-p))-1] 2)Gerar a distribuição de Poisson através do processo de Poisson que pode ser simulado gerando variáveis aleatórias exponenciais até sua soma exceder 1. Lembrar o método da inversa para distribuição exponencial X1= -1/A(ln(Ri)) Então A1+....+Ax<=1<A1+A2+....+Ax+Ax+1 torna-se -somatória de i=1 até x(1/A(ln(Ri)<=1< -somatória(i=1 até x+1(1/lambda(ln(Ri)) somatória de i=1 até x((ln(Ri)>=lambda>somatória(i=1 até x+1(ln(Ri)) ln(produtória(i=1 até x(Ri)>= -lambda>ln(produtória(i=1 até x(Ri) (produtória(i=1 até x(Ri)>=e^ -lambda>(produtória(i=1 até x(Ri) Isto leva ao algoritmo da rejeição Usando (produtória(i=1 até x(Ri)>=e^ -lambda>(produtória(i=1 até x(Ri) Passo 1: N=0 ; P=1 Passo 2: Gere número aleatório Rn-1 e substitua P por P*Rn-1 Passo 3: Se P<e^-lambda, então aceita N=n, senão rejeita o n, soma n por um e retorne ao passo 2. 3) Gerar a distribuição de Pareto que é definida por f(x|alfa)=alfa*x^-alfa-1 sobre (1,infinito). Mostre que pode ser gerada como potencia -1/alfa de uma variável uniforme. Eu tentei fazer a inversa mas não sei se é isto: a funçao de distribuição deve ser F(x)=1-x^-alfa e ao fazer a tranformada inversa: log(1-x)=log(u)/-alfa 4)Gerar a distribuição Beta onde U e V são i.i.d e a distribuição é (U^1/alfa)/((U^1/alfa)+V^(1/Beta)), condicional em U^1/alfa+V^1/B<=1.Tem que ser gerado com a potencia -1/alfa de uma uniforme. Desde já agradeço [edit] -------------------------------------------------------------------------------- _______________________________________________ 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.

2 2

[Dúvida] Salvar gráfico com boa qualidade
by Pedro Rafael 21 Ago '11

21 Ago '11

Pessoal como salvar um gráfico no R com uma ótima qualidade para que o mesmo possa ser ampliado sem prejudicar a qualidade da imagem. Existe formas de salvar a imagem e definir o tamanho da figura? Saudações, Pedro Rafael

4 6

OFF - Topic - procurando um artigo
by Pedro Emmanuel Alvarenga Americano do Brasil 20 Ago '11

20 Ago '11

Amigos de R, Será que alguns de voces poderia ter acesso ou tem o pdf desse artigo aqui? http://www.tandfonline.com/doi/abs/10.1080/03610910701212181#preview , Abraço forte e que a força esteja com você, Dr. Pedro Emmanuel A. A. do Brasil Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365 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

3 3

Função tapply considera comprimentos deiferentes com NA
by Alexandre Santos 20 Ago '11

20 Ago '11

Bom Dia Pessoal, Estava tentando tirar à média de 3 fatores para densidade de insetos, sendo 4, 16 e 64 insetos e cada fator com 10 repetições, porém em 4 insetos tenho 3 NA' s e quando tento usar tapply() a função dá erro obviamente, pois são comprimentos diferentes, sendo: > medias<-tapply(na.omit(n.provas),densidade,mean) Erro em tapply(na.omit(n.provas), densidade, mean) : argumentos devem ter o mesmo comprimento > medias 16prey 4prey 64prey 4.0 NA 12.4 Alguém saberia sugerir outra função, sendo que já tentei aggregate() também sem sucesso, Obrigado, Alexandre

3 3

R simulação: Entendendo as distribuições de probabilidade
by Jose Claudio Faria 20 Ago '11

20 Ago '11

Uma das coisas que considero mais interessantes ao adotar o R em um cursos introdutórios de Probabilidade e Estatística, é possibilitar que os estudante entendam (efetivamente) o que são as distribuições de probabilidades (inf. paramétrica) que eles estão começando a usar para tomar decisões em situação de incerteza. (Em outra palavras,: jogar a mágica na lata de lixo!) Uma das que mais gosto, que pode vir a ser útil para outras pessoas está abaixo: (Estou socializado: sujeito a críticas e sugestões...) # Simulação da origem da distribuição F (Snedecor) # Autor: J.C.Faria # Tamanho das amostras nn <- 05 # Tamanho da amostra do numerador nd <- 10 # Tamanho da amostra do denominador # Características da população mu <- 10 sigma <- 2 # Normal padrão # mu <- 0 # sigma <- 1 # Número de pares de amostras n <- 1e4 # Simulação s2n <- apply(matrix(rnorm(n*nn, mu, sigma), ncol=nn), 1, var) s2d <- apply(matrix(rnorm(n*nd, mu, sigma), ncol=nd), 1, var) require(fdth) plot(fdt(s2n/s2d, start=0, end=10, h=.05), type='d', xlim=c(0, 10), ylim=c(0, 0.8)) curve(df(x, nn-1, nd-1), col='darkgreen', add=TRUE, lwd=2) A conclusão óbvia a que se chega: A distribuição F descreve como se distribui a razão entre duas estimativas da variância de uma população normal! Portanto: Pode ser usada para decidir (sob as devidas condições) se duas estimativas da variância podem, ou não, serem consideradas como provenientes de uma mesma população normal. (Que é o uso mais básico da distribuição F, ou seja, teste de precisão). Gostaria de ver a galera desentocar esse tipo de script e postar aqui na lista! (Tenho mais uns poucos nessa linha) Abs, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

5 10

funcao para extrair o ano de uma variável data
by Paula Mendes Luz 19 Ago '11

19 Ago '11

oi pessoal acredito que deva existir uma funcao para isso... se alguem souber me avisa por favor! tenha uma variavel que é uma data (as.Date) no R e gostaria de extrair o ano. como posso fazer isso? obrigada, paula.

2 2

BATCH
by FHRB Toledo 19 Ago '11

19 Ago '11

Senhores, Gostaria de saber se existe uma forma elegante de rodar um script em BATCH passando argumentos como opção na linha de comando? Explico: tenho um script que faz uma tarefa "X", para tal no script tenho alguns argumentos especificados do tipo arg1 <- valorA. Entretanto queria rodar o mesmo escript variando o argumento parao valorB por exemplo. Na minha ignorância estou editando o script no valorB do arg e rodando denovo. Pois bem, existe uma forma melhor e menos tosca de se fazer isso? desde já grato, FH

2 2

Pedindo ajuda para simulação
by Mauro Sznelwar 19 Ago '11

19 Ago '11

1) Como para simular este algoritmo da distribuição geométrica(p)? Conhecendo a cdf: F(x)=1-(1-p)^x+1 F(x-1)=1-(1-p)^x Usando o método da transformada inversa F(x-1)<R<=F(x) 1-(1-p)^x<R<=1-(1-p)^x+1 (1-p)^x+1<R-1<=1-(1-p)^x ln(1-R)/ln(1-p)-1<=X<ln(1-R)/ln(1-p) X=[(ln(1-R)/ln(1-p))-1] 2)Gerar a distribuição de Poisson através do processo de Poisson que pode ser simulado gerando variáveis aleatórias exponenciais até sua soma exceder 1. Lembrar o método da inversa para distribuição exponencial X1= -1/A(ln(Ri)) Então A1+....+Ax<=1<A1+A2+....+Ax+Ax+1 torna-se -somatória de i=1 até x(1/A(ln(Ri)<=1< -somatória(i=1 até x+1(1/lambda(ln(Ri)) somatória de i=1 até x((ln(Ri)>=lambda>somatória(i=1 até x+1(ln(Ri)) ln(produtória(i=1 até x(Ri)>= -lambda>ln(produtória(i=1 até x(Ri) (produtória(i=1 até x(Ri)>=e^ -lambda>(produtória(i=1 até x(Ri) Isto leva ao algoritmo da rejeição Usando (produtória(i=1 até x(Ri)>=e^ -lambda>(produtória(i=1 até x(Ri) Passo 1: N=0 ; P=1 Passo 2: Gere número aleatório Rn-1 e substitua P por P*Rn-1 Passo 3: Se P<e^-lambda, então aceita N=n, senão rejeita o n, soma n por um e retorne ao passo 2. [edit]

1 0

Imprimir mais de um data.frame em um único arquivo
by Carolina Pena 19 Ago '11

19 Ago '11

Olá, Estou com o seguinte problema: gostaria de imprimir três data.frames em um mesmo arquivo csv. Eu gostaria de uma saída no seguinte formato: Titulo data.frame 1: Dados de Belo Horizontes Codigo Risco 1 2742 0.0004 2 2267 0.0006 3 2554 0.0007 4 2963 0.0008 Titulo data.frame 2: Dados de Minas Gerais Codigo Risco 1 4743 0.0004 2 4200 0.0003 3 4359 0.0009 4 4555 0.0001 Titulo data.frame 2: Dados do Brasil Codigo Risco 1 5743 0.0006 2 5200 0.0008 3 5359 0.0002 4 5555 0.0004 Eu pensei em colocar tudo em uma lista, mas até agora não consegui encontrar nenhum comando que cria arquivos a partir de uma lista. Eu costumo utilizar o write.table, mas nesse caso teria que criar um arquivo para cada data.frame. Obrigada! Abraços, Carolina

2 1

[OFF-Topic] Agenda de webnários sobre Mendeley
by Walmes Zeviani 18 Ago '11

18 Ago '11

Pessoal, Acredito que nem todos os interessados assistiram o webnário de terça-feira sobre Mendeley anunciado pelo Fernando Toledo. Bem, eu assisti e achei muito conveniente divulgar a agenda disponível no site para os próximos webnários. Recomendo que tomem conhecimento desde aplicativo. Para se cadastrar para os próximos webnários visite o link e escolha uma data. https://mendeley.webex.com/mw0306ld/mywebex/default.do?siteurl=mendeley 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