Agosto 2012 - R-br

Empilhar data-frames com colunas diferentes
by Paulo Nogueira 20 Ago '12

20 Ago '12

Boa tarde grupo! Estou com um problema que acredito ser bastante simples de resolver, porém, pesquisei e ainda não encontrei alternativa prática. A grosso modo, quero empilhar bases para formar uma única. O problema é que elas não possuem o mesmo layout. Todas possuem algumas colunas em comum e outras não. Quando uma base não possuir determinada coluna gostaria de deixar NA nos registros oriundos dessa base. Para quem conhece SAS, procuro o resultado de fazer simplesmente: DATA base_final; SET Base1 Base2; RUN; No R ainda não consigo utilizar a função rbind com o mesmo resultado. Segue exemplo para reprodução. # Exemplo d1<-data.frame(index=1,a=1,b=91)d2<-data.frame(index=1,b=92,a=2,c=200)rbind.data.frame(d1,d2)Erro em rbind.data.frame(d1, d2) : números de colunas dos argumentos não correspondem O resultado que desejo é index a b c 1 1 91 NA 1 2 92 200 O rbind troca a posição das colunas baseada nos nomes. Isso ajuda. Mas não sei como fazer funcionar quando algumas colunas estão apenas em algumas bases, o que no exemplo acontece com a coluna c da base d2. Como vocês fariam? Abraços, Paulo Nogueira Starzynski

3 4

Re: [R-br] Regressão Gama
by Ivan Bezerra Allaman 20 Ago '12

20 Ago '12

Neste site existem "n" gráficos com intervalos de confiança. Provavelmente seja possível adaptar para o seu caso. http://ridiculas.wordpress.com/ \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}

3 2

Re: [R-br] TABELAS POR GRUPO
by Ivan Bezerra Allaman 20 Ago '12

20 Ago '12

?fdth ?descr \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

TABELAS POR GRUPO
by Edson Lira 20 Ago '12

20 Ago '12

Bom dia,estou usando esta rotina para analisar o banco postado abaixo. (Veja no link) g<-read.csv2("C:\\Consulta\\HEMOAM\\Elene_fisio\\gest.csv") g[g==""]<-NA names(g) tabelas<-function(x){ qtde<-table(x) perc<-100*prop.table(table(x)) cbind(qtde,perc) } options(digits=3) #imc g$IMC_f<-rep(NA,length(g$IMC)) summary(g$IMC) g$IMC_f[g$IMC<18.5]<-'Baixo peso' g$IMC_f[18.5<=g$IMC &g$IMC<25]<-'Normal' g$IMC_f[25<=g$IMC &g$IMC<30]<-'Pré-obesidade' g$IMC_f[30<=g$IMC &g$IMC<35]<-'Obesidade - I' g$IMC_f[35<=g$IMC &g$IMC<40]<-'Obesidade - II' g$IMC_f[g$IMC>=40]<-'Obesidade - III' #Idade g$ida<-rep(NA,length(g$Idade)) summary(g$Idade) g$ida[18<=g$Idade &g$Idade<25]<-'1. 18 a 25' g$ida[25<=g$Idade &g$Idade<32]<-'2. 25 a 32' g$ida[32<=g$Idade &g$Idade<39]<-'3. 32 a 39' apply(g[,c(-1,-2,-3,-7,-11,-14,-15,-16)],2,tabelas) http://www.datafilehost.com/download-64b7d14d.html Gostaria de obter um resulado semelhante a rotina tabelas, só que estratificado por grupo(Controle, Tratado). Alguém tem uma sugestão? Edson Lira Estatístico Manaus-Amazonas

1 0

Re: [R-br] Somar vendas por semana e por mês
by Hugo Sena 18 Ago '12

18 Ago '12

Boas. Obrigado pela dica, penso que a função aggregate() resolve mas não consigo utilizar. Tenho os seguintes dados: *Data* *Qtd* 28-10-2003 1 28-10-2003 1 31-10-2003 1 05-11-2003 1 11-11-2003 1 11-11-2003 1 11-11-2003 1 11-11-2003 1 11-11-2003 1 11-11-2003 1 E o que pretendo é somar as quantidades em cada dia de forma a ficar assim: Data Qtd 28-10-2003 2 31-10-2003 1 05-11-2003 1 11-11-2003 6 Podem dar uma ajuda? E depois é suposto exportar esta agregação para um ficheiro? Como? Já procurei no google mas ainda não consegui fazer no meu caso em particular. Obrigado. Cumps.

1 0

Re: [R-br] Anova em faixas
by Juliano Farias 18 Ago '12

18 Ago '12

Srs, Grato pela ajuda. Éder o link que você repassou tem exemplos com dois fatores. O que eu tenho é com três fatores. Em Inglês o delineamento é chamado de strip split plot. Segue o quadro da análise da variância para o meu tipo de experimento. Não consegui montar esse quadro da ANOVA no R. CV GL QM F(SobHo) B b-1 V1 V1/V5 A a-1 V2 V2/V5 C c-1 V3 V3/V5 A*C (a-1)(c-1) V4 V4/V5 Erro(ac) (ac-1)(b-1) V5 V5/V11 D d-1 V6 V6/V7 Erro(d) (b-1)(d-1) V7 V7/V11 A*D (a-1)(d-1) V8 V8/V11 C*D (c-1)(d-1) V9 V9/V11 A*C*D (a-1)(c-1)(d-1) V10 V10/V11 Erro (ac-1)(d-1)(b-1) V11 ----- Atenciosamente, Juliano Em 16 de agosto de 2012 12:00, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > > Tópicos de Hoje: > > 1. Re: Regressão Gama (Ivan Bezerra Allaman) > 2. Re: Regressão Gama (Walmes Zeviani) > 3. Anova em faixas (Juliano Farias) > 4. Re: Anova em faixas (Eder David Borges da Silva) > 5. Re: Anova em faixas (Listeiro 037) > 6. Re: Anova em faixas (Eder David Borges da Silva) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 15 Aug 2012 08:20:07 -0700 (PDT) > From: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> > To: R Brasil <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Regressão Gama > Message-ID: > <1345044007.77714.YahooMailNeo(a)web161804.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Neste site existem "n" gráficos com intervalos de confiança. Provavelmente > seja possível adaptar para o seu caso. > > http://ridiculas.wordpress.com/ > > > > \begin{signature} > <<>>= > Prof. Dr. Ivan Bezerra Allaman > Universidade Estadual de Santa Cruz > Departamento de Ciências Exatas e Tecnológicas > Ilhéus/BA - Brasil > Fone: +55 73 3680-5596 > E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com > @ > \end{signature} > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120815/b1f6625c/atta… > > > > ------------------------------ > > Message: 2 > Date: Wed, 15 Aug 2012 13:31:28 -0300 > From: Walmes Zeviani <walmeszeviani(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, Ivan Bezerra Allaman > <ivanalaman(a)yahoo.com.br> > Subject: Re: [R-br] Regressão Gama > Message-ID: > <CAFU=EkYGjAuNnRnBD7fV_5+n4Yd0QS2NM3E+pZ6CyKX= > Ez+U1A(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Além da sugestão do Ivan, como trata-se de um GLM, recomendo a leitura > dessa matéria > > > http://freakonometrics.blog.free.fr/index.php?post/2011/11/01/Fisher,-delta… > > À disposição. > Walmes. > > ========================================================================== > Walmes Marques Zeviani > LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) > Departamento de Estatística - Universidade Federal do Paraná > fone: (+55) 41 3361 3573 > VoIP: (3361 3600) 1053 1173 > e-mail: walmes(a)ufpr.br > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120815/9188fb40/atta… > > > > ------------------------------ > > Message: 3 > Date: Wed, 15 Aug 2012 21:32:40 -0300 > From: Juliano Farias <julianofarias(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Anova em faixas > Message-ID: > < > CAEer2Pik5nco26+F88tkCr02nAx0-kH3geR7GfF01p33+Zwz7g(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Srs, > > Preciso de uma ajuda. > > Preciso fazer uma ANOVA de um experimento com três fatores, porém o > problema está no delineamento que foi em faixas. > Tenho o quadro da análise da variança para esse tipo de delineamento, porém > não consegui montar a programação e não achei um pacote para esse tipo de > delineamento. > O experimento foi realizado em blocos (B) ao acaso com quatro repetições em > esquema fatorial em faixas 2 x 2 x 2. Os dois níveis do fator A, com os > dois níveis do fator C foram casualizados nas parcelas principais, e os > dois níveis do fator D foram casualizados em faixas perpendiculares às > parcelas principais. > > Grato pela atenção. > > > -- > Atenciosamente, > > Juliano Ricardo Farias > Eng. Agrônomo > Doutorando em Entomologia > Lab. de Resistência de Artrópodes a Pesticidas > ESALQ/USP > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120815/1097b154/atta… > > > > ------------------------------ > > Message: 4 > Date: Thu, 16 Aug 2012 07:35:30 -0300 > From: Eder David Borges da Silva <eder(a)leg.ufpr.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Anova em faixas > Message-ID: > <CALKkXXru-3qZThcWiKNAeHAN= > TMYnzF2uyDMLE4Mj1RAyZk47w(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Juliano, > Acho que o script neste link pode te ajudar, > http://www.leg.ufpr.br/~eder/cursoR/Exercicios.R > Att > Éder > > Em 15 de agosto de 2012 21:32, Juliano Farias <julianofarias(a)gmail.com > >escreveu: > > > Srs, > > > > Preciso de uma ajuda. > > > > Preciso fazer uma ANOVA de um experimento com três fatores, porém o > > problema está no delineamento que foi em faixas. > > Tenho o quadro da análise da variança para esse tipo de delineamento, > > porém não consegui montar a programação e não achei um pacote para esse > > tipo de delineamento. > > O experimento foi realizado em blocos (B) ao acaso com quatro repetições > > em esquema fatorial em faixas 2 x 2 x 2. Os dois níveis do fator A, com > > os dois níveis do fator C foram casualizados nas parcelas principais, e > > os dois níveis do fator D foram casualizados em faixas perpendiculares às > > parcelas principais. > > > > Grato pela atenção. > > > > > > -- > > Atenciosamente, > > > > Juliano Ricardo Farias > > Eng. Agrônomo > > Doutorando em Entomologia > > Lab. de Resistência de Artrópodes a Pesticidas > > ESALQ/USP > > > > _______________________________________________ > > 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/20120816/3100052d/atta… > > > > ------------------------------ > > Message: 5 > Date: Thu, 16 Aug 2012 09:27:05 -0300 > From: Listeiro 037 <listeiro_037(a)yahoo.com.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Anova em faixas > Message-ID: <375172.55957.bm(a)smtp146.mail.mud.yahoo.com> > Content-Type: text/plain; charset=UTF-8 > > Olá. > > Há como colocar shebang nesse script de R abaixo? > > > Em Thu, 16 Aug 2012 07:35:30 -0300 > Eder David Borges da Silva <eder(a)leg.ufpr.br> escreveu: > > > Juliano, > > Acho que o script neste link pode te ajudar, > > http://www.leg.ufpr.br/~eder/cursoR/Exercicios.R > > Att > > Éder > > > > Em 15 de agosto de 2012 21:32, Juliano Farias > > <julianofarias(a)gmail.com>escreveu: > > > > -- > "There are three kinds of lies: lies, damn lies, and statistics." > Benjamin Disraeli > > > ------------------------------ > > Message: 6 > Date: Thu, 16 Aug 2012 10:44:36 -0300 > From: Eder David Borges da Silva <eder(a)leg.ufpr.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Anova em faixas > Message-ID: > < > CALKkXXqS7WV+sdjDoWbdMyZgHZsOf_qapeZbeXL+1qMmEnpirQ(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > ?Shebang > > Em 16 de agosto de 2012 09:27, Listeiro 037 <listeiro_037(a)yahoo.com.br > >escreveu: > > > Olá. > > > > Há como colocar shebang nesse script de R abaixo? > > > > > > Em Thu, 16 Aug 2012 07:35:30 -0300 > > Eder David Borges da Silva <eder(a)leg.ufpr.br> escreveu: > > > > > Juliano, > > > Acho que o script neste link pode te ajudar, > > > http://www.leg.ufpr.br/~eder/cursoR/Exercicios.R > > > Att > > > Éder > > > > > > Em 15 de agosto de 2012 21:32, Juliano Farias > > > <julianofarias(a)gmail.com>escreveu: > > > > > > > -- > > "There are three kinds of lies: lies, damn lies, and statistics." > > Benjamin Disraeli > > _______________________________________________ > > 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/20120816/c74a1216/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 18, assunto 21 > ***************************************** > -- Atenciosamente, Juliano Ricardo Farias Eng. Agrônomo Doutorando em Entomologia Lab. de Resistência de Artrópodes a Pesticidas ESALQ/USP

2 1

Re: [R-br] Digest R-br, volume 18, assunto 21
by Juliano Farias 18 Ago '12

18 Ago '12

Listeiro 037, Pelo que entendi do shebang, sim pode colocar no script do R. Obrigado, Juliano. Em 16 de agosto de 2012 12:00, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > > Tópicos de Hoje: > > 1. Re: Regressão Gama (Ivan Bezerra Allaman) > 2. Re: Regressão Gama (Walmes Zeviani) > 3. Anova em faixas (Juliano Farias) > 4. Re: Anova em faixas (Eder David Borges da Silva) > 5. Re: Anova em faixas (Listeiro 037) > 6. Re: Anova em faixas (Eder David Borges da Silva) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 15 Aug 2012 08:20:07 -0700 (PDT) > From: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> > To: R Brasil <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Regressão Gama > Message-ID: > <1345044007.77714.YahooMailNeo(a)web161804.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Neste site existem "n" gráficos com intervalos de confiança. Provavelmente > seja possível adaptar para o seu caso. > > http://ridiculas.wordpress.com/ > > > > \begin{signature} > <<>>= > Prof. Dr. Ivan Bezerra Allaman > Universidade Estadual de Santa Cruz > Departamento de Ciências Exatas e Tecnológicas > Ilhéus/BA - Brasil > Fone: +55 73 3680-5596 > E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com > @ > \end{signature} > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120815/b1f6625c/atta… > > > > ------------------------------ > > Message: 2 > Date: Wed, 15 Aug 2012 13:31:28 -0300 > From: Walmes Zeviani <walmeszeviani(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, Ivan Bezerra Allaman > <ivanalaman(a)yahoo.com.br> > Subject: Re: [R-br] Regressão Gama > Message-ID: > <CAFU=EkYGjAuNnRnBD7fV_5+n4Yd0QS2NM3E+pZ6CyKX= > Ez+U1A(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Além da sugestão do Ivan, como trata-se de um GLM, recomendo a leitura > dessa matéria > > > http://freakonometrics.blog.free.fr/index.php?post/2011/11/01/Fisher,-delta… > > À disposição. > Walmes. > > ========================================================================== > Walmes Marques Zeviani > LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) > Departamento de Estatística - Universidade Federal do Paraná > fone: (+55) 41 3361 3573 > VoIP: (3361 3600) 1053 1173 > e-mail: walmes(a)ufpr.br > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120815/9188fb40/atta… > > > > ------------------------------ > > Message: 3 > Date: Wed, 15 Aug 2012 21:32:40 -0300 > From: Juliano Farias <julianofarias(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Anova em faixas > Message-ID: > < > CAEer2Pik5nco26+F88tkCr02nAx0-kH3geR7GfF01p33+Zwz7g(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Srs, > > Preciso de uma ajuda. > > Preciso fazer uma ANOVA de um experimento com três fatores, porém o > problema está no delineamento que foi em faixas. > Tenho o quadro da análise da variança para esse tipo de delineamento, porém > não consegui montar a programação e não achei um pacote para esse tipo de > delineamento. > O experimento foi realizado em blocos (B) ao acaso com quatro repetições em > esquema fatorial em faixas 2 x 2 x 2. Os dois níveis do fator A, com os > dois níveis do fator C foram casualizados nas parcelas principais, e os > dois níveis do fator D foram casualizados em faixas perpendiculares às > parcelas principais. > > Grato pela atenção. > > > -- > Atenciosamente, > > Juliano Ricardo Farias > Eng. Agrônomo > Doutorando em Entomologia > Lab. de Resistência de Artrópodes a Pesticidas > ESALQ/USP > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120815/1097b154/atta… > > > > ------------------------------ > > Message: 4 > Date: Thu, 16 Aug 2012 07:35:30 -0300 > From: Eder David Borges da Silva <eder(a)leg.ufpr.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Anova em faixas > Message-ID: > <CALKkXXru-3qZThcWiKNAeHAN= > TMYnzF2uyDMLE4Mj1RAyZk47w(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Juliano, > Acho que o script neste link pode te ajudar, > http://www.leg.ufpr.br/~eder/cursoR/Exercicios.R > Att > Éder > > Em 15 de agosto de 2012 21:32, Juliano Farias <julianofarias(a)gmail.com > >escreveu: > > > Srs, > > > > Preciso de uma ajuda. > > > > Preciso fazer uma ANOVA de um experimento com três fatores, porém o > > problema está no delineamento que foi em faixas. > > Tenho o quadro da análise da variança para esse tipo de delineamento, > > porém não consegui montar a programação e não achei um pacote para esse > > tipo de delineamento. > > O experimento foi realizado em blocos (B) ao acaso com quatro repetições > > em esquema fatorial em faixas 2 x 2 x 2. Os dois níveis do fator A, com > > os dois níveis do fator C foram casualizados nas parcelas principais, e > > os dois níveis do fator D foram casualizados em faixas perpendiculares às > > parcelas principais. > > > > Grato pela atenção. > > > > > > -- > > Atenciosamente, > > > > Juliano Ricardo Farias > > Eng. Agrônomo > > Doutorando em Entomologia > > Lab. de Resistência de Artrópodes a Pesticidas > > ESALQ/USP > > > > _______________________________________________ > > 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/20120816/3100052d/atta… > > > > ------------------------------ > > Message: 5 > Date: Thu, 16 Aug 2012 09:27:05 -0300 > From: Listeiro 037 <listeiro_037(a)yahoo.com.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Anova em faixas > Message-ID: <375172.55957.bm(a)smtp146.mail.mud.yahoo.com> > Content-Type: text/plain; charset=UTF-8 > > Olá. > > Há como colocar shebang nesse script de R abaixo? > > > Em Thu, 16 Aug 2012 07:35:30 -0300 > Eder David Borges da Silva <eder(a)leg.ufpr.br> escreveu: > > > Juliano, > > Acho que o script neste link pode te ajudar, > > http://www.leg.ufpr.br/~eder/cursoR/Exercicios.R > > Att > > Éder > > > > Em 15 de agosto de 2012 21:32, Juliano Farias > > <julianofarias(a)gmail.com>escreveu: > > > > -- > "There are three kinds of lies: lies, damn lies, and statistics." > Benjamin Disraeli > > > ------------------------------ > > Message: 6 > Date: Thu, 16 Aug 2012 10:44:36 -0300 > From: Eder David Borges da Silva <eder(a)leg.ufpr.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Anova em faixas > Message-ID: > < > CALKkXXqS7WV+sdjDoWbdMyZgHZsOf_qapeZbeXL+1qMmEnpirQ(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > ?Shebang > > Em 16 de agosto de 2012 09:27, Listeiro 037 <listeiro_037(a)yahoo.com.br > >escreveu: > > > Olá. > > > > Há como colocar shebang nesse script de R abaixo? > > > > > > Em Thu, 16 Aug 2012 07:35:30 -0300 > > Eder David Borges da Silva <eder(a)leg.ufpr.br> escreveu: > > > > > Juliano, > > > Acho que o script neste link pode te ajudar, > > > http://www.leg.ufpr.br/~eder/cursoR/Exercicios.R > > > Att > > > Éder > > > > > > Em 15 de agosto de 2012 21:32, Juliano Farias > > > <julianofarias(a)gmail.com>escreveu: > > > > > > > -- > > "There are three kinds of lies: lies, damn lies, and statistics." > > Benjamin Disraeli > > _______________________________________________ > > 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/20120816/c74a1216/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 18, assunto 21 > ***************************************** > -- Atenciosamente, Juliano Ricardo Farias Eng. Agrônomo Doutorando em Entomologia Lab. de Resistência de Artrópodes a Pesticidas ESALQ/USP

2 1

Re: [R-br] dúvida - modelo [OFF-TOPIC]
by Fernando Colugnati 17 Ago '12

17 Ago '12

Só um pitaco na modelagem do ponto de vista fisiológico...o fato do cara não fazer atividade física não implica ele em ter Gasto Energético 0 !!!!!! Se vc utilizar esta abordagem vai ter uma subestimação não real das médias ajustadas e uma superestimação de efeitos. Neste caso vcs deveriam ter no mínimo o gasto energético basal! Se vcs estiverem medindo o gasto energético pelo IPAQ (um questionário, que na verdade fornece METs, uma medida proxy de gasto energético) sugiro ver outros trabalhos que adequam este problema. Estes 0s devem ser pessoas que não responderam ao questionário por não fazerem atividades, mas isso não quer dizer que não gastam energia....neste caso, estes caras não deveriam entrar na modelagem. Outro ponto, mesmo medido como valores inteiros, também não faz sentido os dados seguirem uma distribuição Geométrica...na verdade vcs devem ter uma discretização deste dado, mas veja a definição da Geométrica e verá que não faz sentido. Nâo assumo uma distribuição pq os dados parecem com aquela distribuição, mas sim pelo processo que pode ter gerado aqueles dados, mas isso é pessoal...há escolas que fazem desta forma. Desculpem se dei pitacos desnecessários, mas pela breve descrição do problema isso me chamou a atenção! Abs Em 16 de agosto de 2012 22:47, Sérgio Henrique <sergio.edfisica(a)gmail.com>escreveu: > Oi > > Obrigado por sempre ajudar Benilton. > > Deixa eu tentar explicar o objetivo do trabalho. > > Na verdade que vai trabalhar com esses dados é um amigo no doutorado. > > A variável resposta é Atividade Física e ele quer modelar, na verdade, > seria gasto energético por AF, os 0 são indivíduos que não praticam AF, ou > seja, tem gasto 0, acima desse valor, são gastos energéticos de atividades > físicas quaisquer... > Os dados não são normais e quando tento normalizar com boxcox não roda por > conta dos 0. > Quando fiz um histograma dessa variável, me pareceu que os mesmos > apresentam distribuição geométrica, parei por aí... > > Obrigado mais uma vez. > > Enviado via iPod > > Em 16/08/2012, às 22:11, Benilton Carvalho <beniltoncarvalho(a)gmail.com> > escreveu: > > > modelagem envolve muito mais que os dados com os quais vc esta' > > trabalhando... na verdade, o processo gerador dos dados e' q e' > > relevante... vc precisa entender a "historia" dos dados... > > > > sobre os seus dados, pergunte-se: 1) eles sao sempre inteiros? 2) se > > sim, eles representam o numero de tentativas (de alguma coisa) ate' o > > primeiro sucesso (ou falha)? 3) se sim, cada uma dessas tentativas tem > > a mesma chance de sucesso (ou falha)? > > > > se sua resposta for "nao" para qq uma das perguntas acima, e' pouco > > provavel que uma geometrica seja a resposta q vc procura... > > > > b > > > > 2012/8/17 Sérgio Henrique almeida da silva ju <sergio.edfisica(a)gmail.com > >: > >> Olá > >> > >> Estou modelando um conjunto de dados e a variável resposta parece ter > >> distribuição geométrica, qual é a melhor forma de modelar esses dados? > >> > >> Abraços > >> > >> > >> dput(dados) > >> structure(list(AFtotal = structure(c(103L, 76L, 58L, 105L, 125L, > >> 163L, 106L, 146L, 227L, 187L, 5L, 117L, 147L, 39L, 232L, 2L, > >> 67L, 150L, 46L, 118L, 18L, 171L, 157L, 226L, 80L, 96L, 142L, > >> 162L, 128L, 152L, 109L, 10L, 108L, 102L, 47L, 208L, 122L, 78L, > >> 176L, 3L, 196L, 40L, 209L, 119L, 179L, 29L, 1L, 81L, 190L, 24L, > >> 204L, 28L, 92L, 136L, 104L, 85L, 91L, 15L, 141L, 13L, 11L, 121L, > >> 112L, 210L, 31L, 6L, 70L, 156L, 114L, 68L, 107L, 112L, 211L, > >> 167L, 1L, 215L, 88L, 213L, 194L, 136L, 220L, 199L, 136L, 9L, > >> 136L, 118L, 38L, 192L, 42L, 62L, 155L, 19L, 79L, 48L, 138L, 172L, > >> 134L, 231L, 231L, 216L, 136L, 195L, 130L, 101L, 217L, 205L, 174L, > >> 36L, 111L, 140L, 106L, 99L, 22L, 218L, 170L, 7L, 8L, 136L, 19L, > >> 62L, 25L, 123L, 1L, 137L, 220L, 97L, 52L, 169L, 190L, 216L, 1L, > >> 53L, 172L, 100L, 79L, 86L, 201L, 61L, 198L, 17L, 218L, 67L, 106L, > >> 193L, 1L, 138L, 164L, 200L, 190L, 133L, 129L, 74L, 186L, 72L, > >> 41L, 177L, 48L, 113L, 50L, 229L, 178L, 106L, 124L, 75L, 14L, > >> 148L, 218L, 191L, 19L, 219L, 30L, 1L, 71L, 112L, 172L, 21L, 221L, > >> 16L, 203L, 27L, 69L, 143L, 82L, 207L, 223L, 32L, 54L, 60L, 59L, > >> 1L, 97L, 149L, 63L, 132L, 212L, 214L, 1L, 1L, 165L, 77L, 20L, > >> 73L, 26L, 168L, 160L, 233L, 139L, 1L, 188L, 95L, 175L, 206L, > >> 93L, 144L, 135L, 56L, 6L, 189L, 197L, 45L, 30L, 35L, 1L, 4L, > >> 136L, 154L, 89L, 65L, 184L, 225L, 230L, 151L, 1L, 173L, 43L, > >> 44L, 33L, 182L, 34L, 153L, 1L, 1L, 90L, 106L, 110L, 98L, 1L, > >> 219L, 1L, 1L, 116L, 49L, 64L, 57L, 131L, 185L, 145L, 17L, 94L, > >> 1L, 180L, 1L, 120L, 159L, 51L, 172L, 181L, 183L, 7L, 87L, 106L, > >> 228L, 55L, 126L, 157L, 218L, 158L, 115L, 133L, 222L, 84L, 166L, > >> 1L, 37L, 1L, 1L, 202L, 66L, 12L, 1L, 83L, 23L, 110L, 233L, 231L, > >> 161L, 106L, 224L, 127L), .Label = c("0,00", "1.031,00", "1.035,00", > >> "1.042,50", "1.064,00", "1.080,00", "1.120,00", "1.155,00", "1.171,50", > >> "1.177,50", "1.179,00", "1.200,00", "1.204,50", "1.210,00", "1.218,00", > >> "1.220,00", "1.230,00", "1.233,00", "1.260,00", "1.275,00", "1.287,00", > >> "1.290,00", "1.314,00", "1.320,00", "1.356,00", "1.359,50", "1.400,00", > >> "1.410,00", "1.418,00", "1.440,00", "1.466,00", "1.480,00", "1.500,00", > >> "1.566,00", "1.572,00", "1.590,00", "1.611,00", "1.652,00", "1.680,00", > >> "1.722,00", "1.734,00", "1.740,00", "1.752,00", "1.800,00", "1.818,00", > >> "1.821,00", "1.833,00", "1.836,00", "1.854,00", "1.872,00", "1.900,00", > >> "1.920,00", "1.938,00", "10.080,00", "10.230,00", "10.557,00", > >> "10.680,00", "11.238,00", "12.000,00", "12.270,00", "12.420,00", > >> "120,00", "13.644,00", "13.860,00", "14.590,00", "150,00", "165,00", > >> "180,00", "2.002,50", "2.060,00", "2.078,00", "2.100,00", "2.120,00", > >> "2.157,00", "2.320,00", "2.332,50", "2.376,00", "2.382,00", "2.430,00", > >> "2.457,00", "2.475,00", "2.504,00", "2.520,00", "2.526,00", "2.532,00", > >> "2.556,00", "2.572,00", "2.580,00", "2.616,00", "2.640,00", "2.645,00", > >> "2.655,00", "2.706,00", "2.739,00", "2.745,00", "2.767,50", "2.784,00", > >> "2.798,00", "2.820,00", "2.841,00", "2.895,00", "2.937,00", "2.940,00", > >> "2.965,00", "21.515,00", "240,00", "247,50", "259,00", "26.862,00", > >> "264,00", "273,00", "297,00", "3.019,00", "3.035,00", "3.036,00", > >> "3.282,00", "3.306,00", "3.360,00", "3.497,00", "3.504,00", "3.510,00", > >> "3.540,00", "3.600,00", "3.615,00", "3.660,00", "3.710,00", "3.840,00", > >> "3.892,50", "3.900,00", "3.960,00", "3.977,00", "3.996,00", "300,00", > >> "318,00", "330,00", "360,00", "372,00", "396,00", "4.020,00", > >> "4.086,00", "4.092,00", "4.100,00", "4.140,00", "4.158,00", "4.200,00", > >> "4.219,50", "4.331,00", "4.359,00", "4.554,00", "4.618,00", "4.680,00", > >> "4.770,00", "407,50", "410,00", "420,00", "424,00", "480,00", > >> "495,00", "5.052,00", "5.160,00", "5.340,00", "5.622,00", "5.688,00", > >> "5.715,00", "5.745,00", "5.832,00", "5.940,00", "5.958,00", "510,00", > >> "516,00", "537,00", "540,00", "565,50", "570,00", "575,00", "576,00", > >> "585,00", "6.097,50", "6.285,00", "6.495,00", "6.556,00", "6.576,00", > >> "6.720,00", "6.746,00", "6.828,00", "639,00", "645,00", "655,00", > >> "66,00", "660,00", "678,00", "690,00", "693,00", "696,00", "7.056,00", > >> "7.074,00", "7.200,00", "7.236,00", "7.245,00", "7.500,00", "7.600,00", > >> "7.680,00", "7.860,00", "7.920,00", "720,00", "738,00", "754,00", > >> "754,50", "780,00", "795,00", "8.268,00", "8.467,50", "8.592,00", > >> "8.755,00", "8.820,00", "80,00", "801,00", "810,00", "812,00", > >> "840,00", "876,00", "882,00", "895,50", "9.120,00", "9.150,00", > >> "9.310,00", "9.508,00", "924,00", "930,00", "933,00", "960,00", > >> "99,00", "990,00"), class = "factor")), .Names = "AFtotal", class = > >> "data.frame", row.names = c(NA, > >> -299L)) > >> > >> > >> -- > >> Sérgio Henrique Almeida da Silva Junior > >> Doutorando em Epidemiologia em Saúde Pública > >> Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ > >> http://lattes.cnpq.br/1611345552843383 > >> Tel: (21) 94429486/78101651 id: 123*20942 > >> > >> > >> > >> _______________________________________________ > >> R-br mailing list > >> R-br(a)listas.c3sl.ufpr.br > >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código > >> mínimo reproduzível. > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br > > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > _______________________________________________ > 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. > -- Fernando A.B. Colugnati

3 3

Re: [R-br] Utilizando DEA com R
by Roney Fraga 17 Ago '12

17 Ago '12

Cleyton Fábio, Já utilizei o pacote Benchmarking, segue link com alguns exemplos: http://r-br.2285057.n4.nabble.com/R-br-Analise-Envoltoria-de-Dados-DEA-td37… utilizei os dados do livro: EVALUACION DE LA EFICIENCIA MEDIANTE EL ANALISIS ENVOLVENTE DE DATOS dos autores Vicente Coll Serrano e Olga M Blasco Blasco de 2006. Esse livro esta disponível no google livros. Todos os cálculos tem exemplos no excel/calc, isso ajuda a entender o método e ainda permite comparar os resultados apresentados pelo livro com os realizados por você no excel/calc e no R. Gosto de utilizar o pacote Benchmarking porque tenho o livro Bogetoft e Otto (2011) Benchmarking with DEA, SFA, and R. Como você pode notar, esse livro é dos mesmos autores do pacote, acho o texto didático. Outro livro que gosto é Coelli, Rao, O'Donnell, Battese (2005) Introduction to Efficiency and Productivity Analysis 2ed. Atenciosamente Roney

3 2

Matriz de proximidade
by Gilbert Queiroz 17 Ago '12

17 Ago '12

Pessoal, Alguém sabe se é possível obter a matriz de proximidade pelo R e em qual pacote se faz isto? Abs. Gilbert.

1 0