- R-br - DINF UFPR

geoR
by Fernanda De Bastiani 14 Dez '15

14 Dez '15

Caros, por favor, alguém poderia me ajudar com a melhor forma de ajusta um modelo sem intercepto no geoR, usando a função likfit? Ou alguma outra função que desempenhe o mesmo papel. Muito obrigada,Fernanda

2 2

Recorte de mapa de MG
by Andre Oliveira 14 Dez '15

14 Dez '15

Pessoal plotei o mapa de MG, gostaria de selecionar e colocar em cima só as cidades da zona da mata de MG! São em torno de 75 cidades das 853. Estou usando o comando subset. Tem algo mais simples de selecionar estas 75 cidades? require(rgdal) mg=readOGR("31MUE250GC_SIR.shp",layer="31MUE250GC_SIR",encoding = "latin1") plot(mg,axes=T) zmata=subset(mg, NM_MUNICIP%in%c("ABRE CAMPO","ACAIACA","ALVINOPOLIS","AMPARO DO SERRA")) plot(zmata,col="blue",add=TRUE) André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

3 2

Deviance e Scaled Deviance
by Felinto COSTA 11 Dez '15

11 Dez '15

Prezados. Tenho uma dúvida a respeito da estatística /de//v//iance (D) /que a função /glm{stats} /retorna ao estimar os parâmetros para um modelo Gama com ligação log. Devido à estimação do parâmetro (φˆ ) desconhecido da dispersão (φ) , essa /deviance/ é a /scaled deviance (D*=D///φˆ ). / A dúvida é se essa /scaled deviance /pode ser comparada a um valor da distribuição Chi-quadrado (χ2 (n−p,1−α) como uma (das primeiras) análises sobre a bondade do ajuste desse modelo. Há pesquisadores que afirmam ser isso válido apenas quando o parâmetro (φ) é conhecido (=1), como nas distribuições Binomial e de Poissone por essa razão apenas poderia se ter uma comparação grosseira de D* com os graus de liberdade (n-p) do modelo . Antecipadamente agradeço a ajuda e me desculpo se conceitualmente estou a cometer algum erro crasso . (estatística não é minha formação inicial). FCosta

1 0

Dúvidas gráficos Histograma
by Pedro Campelo 07 Dez '15

07 Dez '15

Pessoal, Existe algum comando no qual posso criar um histograma que o eixo do Y dê o volume (porcentagem) dos dados? Obrigado *----* *Pedro Henrique Campelo Felix* *Bacharel em Engenharia de Alimentos* *Mestre em Engenharia Química* *Doutorando em Ciências de Alimentos* *(31) 87427656 (OI)* *(31) 75601545 (TIM)* *"**Nós somos aquilo que fazemos repetidamente. Excelência, então, não é um modo de agir, mas um hábito." Aristóteles*

5 6

adicionar e configurar celula em uma tabela de Tktable
by Cleber N.Borges 06 Dez '15

06 Dez '15

olá a todos, ainda que seja muito especifico, pergunto caso alguém saiba dar uma dica estou tentando usar uma tabela do tk e tentando aplicar uma tag a uma determinada célula (que é o modo de configurar cada célula individualmente, até onde entendi ) mas não entendi como esse processo funciona... caso alguém saiba, agradeço por antecipação a ajuda. att., cleber ############################# library( tcltk ) tclRequire('Tktable') tt <- tktoplevel() tt.f <- ttkframe( tt, relief='sunken', borderwidth=10 ) tkpack( tt.f, expand=T, fill='both' ) tt.f.l <- ttklabel( tt.f, text='test' ) tkpack( tt.f.l, expand=T, fill='both' ) vartab <- tclArray() tt.f.table <- tkwidget( tt.f, 'table', titlecol=1, titlerow=1, background='white', multiline=0, variable=vartab ) tkpack( tt.f.table, expand=T, fill='both' ) tktag.add( tt.f.table, "celltag", "ZeroZero", "0,0" ) ### erro!!!!! --- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus

1 0

Covariância
by Rodrigo Lima 05 Dez '15

05 Dez '15

Prezados listeiros Bom dia. Eu tenho dados de produtividade de grãos em milho como variável resposta e dados de estande (número de plantas por parcela) como outra variável. Como nós tivemos problemas com irregularidade de estande, eu quero usar o número de plantas como covariável no no meu modelo. Abaixo o modelo: exp01.anova.2<-lm(gy ~ prog+rep+bloc%in%rep, data = exp01[exp01$local==2,]) Como eu coloco essa variável como covariável? O número de plantas ideal é 25 plantas por parcela, mas tenho parcelas com 16, 18 etc. Como coloco isso no aquivo? Eu devo subtrair o número de plantas ideal do real e, aquelas parcelas que tem o número ideal ficariam com zero, aquelas com mais de 25 plantas com número positivo e aquelas com menos de 25 plantas com número negativo? Isso mesmo? Ele vai para o modelo como fator ou numero? Muito obrigado. Att, Rodrigo. -- Prof. Rodrigo Oliveira de Lima Federal University of Viçosa Department of Plant Sciences Viçosa, MG, Brazil, 36570-900, Phone: +55(31)3899-2781 Fax: +55(31)3899-2614 Email: rodrigoodelima(a)ufv.br CV: http://lattes.cnpq.br/2761299616351806

3 3

[LEMBRETE] Meetup R São Paulo Users Group (IME-USP) - 4 de dezembro do 2015
by Luis Benites Sánchez 04 Dez '15

04 Dez '15

Caros, Hoje 4 de dezembro teremos nossa última reunião (numero 8) deste ano, do Meetup de usuários de R de São Paulo, no auditório Jacy Monteiro no Bloco B. Segue a programação: Programação de hoje, 4 de dezembro (14:00h-15:00h) - [14:00 - 14:55] Distribuições Assimétricas e Aplicações com R - Prof. Clécio da Silva Ferreira (UFJF) - [Online] - [14:55 - 15:00] Encerramento As palestras serão transmitidas ao vivo pela conexão IPTV ( http://iptv.usp.br/portal/transmission.action?idItem=31314) da USP e YouTube (https://plus.google.com/u/0/events/cgdvprrjafd27vba5pskbu0qghs) Para poder participar, pode confirmar sua participação gratuitamente no site do Meetup :) Para mais informações: Site: www.meetup.com/pt/useR-SP Canal de YouTube: www.youtube.com/channel/UCb9PdJYYaElE6Uqb1KxhE8w e-mail: meetuprsp(a)gmail.com Um abraço, Luis Benites Sánchez Rocío Paola Maehara Organizadores do São Paulo R Users Group

1 0

set environment
by Cleber Borges 03 Dez '15

03 Dez '15

olá a todos, Gostaria de saber se há um modo mais direto de ajustar um determinado environment e definir variáveis dentro dele e sub ambientes. agradeço por qualquer ajuda por antecedência cleber ############################## #Algo do tipo: Xenv <- new.env() * setCurrentEnv( Xenv ) * # Definir procedimentos e muitas variáveis * setCurrentEnv( baseenv() ) * --- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus

1 0

R jobs
by Elias Teixeira Krainski 01 Dez '15

01 Dez '15

Bom dia caros(as), É fato que "data science" virou uma profissão. E R é uma das ferramentas dessa profissão. Chamou minha atenção à diversidade de oportunidades de trabalho para quem é dessa área. Quando restringimos para R, surgem oportunidades, por exemplo, para desenvolver uma Shiny App. Segue uma lista de 16 destas oportunidades publicadas no r-bloggers: http://www.r-bloggers.com/14-new-r-jobs-from-around-the-world-for-2015-11-3… Elias

1 0

Rosa dos ventos e escala em mapas
by Andre Oliveira 28 Nov '15

28 Nov '15

Boa tarde,alguém saberia dizer como posso colocar Rosa dos ventos e escala nos mapas. Segue um exemplo. Mapas do IBGE no R | | | | | | | | | | | Mapas do IBGE no RVou postar uma série de comandos para fazer mapas coloridos no R usando os shape files disponíveis no site do IBGE. No R, vc vai precisar do pacote rgdal (Eu postei... | | | | Visualizar em leobastos.word... | Visualizado por Yahoo | | | | | André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

2 1

Sugestão sobre etiqueta na lista
by Cesar Rabak 28 Nov '15

28 Nov '15

Prez@dos, Fiz uma rápida consulta ao "Guia de postagem" desta lista, e como ele está em construção, gostaria de aproveitar o ensejo de uma postagem recente nesta lista, que ficará evidente em seguida, evito citá-la diretamente. Essa postagem, que insisto apenas me chamou a atenção de um número similar em teor dela, tem como estrutura geral o seguinte: OP: > Estou usando o programa/SO/pacote X e estou tendo as dificuldades/problemas/ dúvidas Y. Resposta que coloco para análise dos colegas: > O "programa/SO/pacote X" 'não funciona'/'é ruim'/'deve ser evitado' use o "programa/SO/pacote Z" que alem de resolver "Y" vai <várias expressões de julgamento de valor [que não ajudam em nada o OP]>. Minha proposta é que o tipo de resposta acima seja desaconselhado, especialmente se for baseado apenas em opinião "vazia", preconceito, ou preferência que só pode ser definida como "pessoal" e não objetiva. Por exemplo, se um certo pacote for "deprecado" tando pela equipa do CRAN como eventualmente tenha sido descontinuado pelo desenvolvedor, uma orientação sobre o novo pacote ou solução de contorno me parece OK, mas dizer que "A" é melhor que "B" quando ambos os pacotes estão à disposição me parece arbitrário, a menos que o proponente ofereça uma comparação objetiva sobre sob quais aspectos "A" seria superior a "B". Assim, críticas a determinada distribuição do Linux, uso do Windows em vez de outro SO, e até mesmo a pejorativa denominação do sistema operacional comercial da Microsoft de "ruindows" e equivalentes deveria ser desencorajada na mesma alínea dessa orientação no Guia. Deixo a sugestão para a lista para termos uma democrática e produtiva discussão e inclusive para poder melhorar minha ideia e posterior inclusão no Guia. []s -- Cesar Rabak

2 1

Simulação - Geo Estatística
by Andre Oliveira 28 Nov '15

28 Nov '15

Boa noite,como simular dados com (altitude e longitude + Y) para formar um banco de dados? Algo como o CMR abaixo, só que para Geo Estatística! da <- expand.grid(Tratamentos=gl(4,1),Repetição=gl(4,1)) da$Resposta<- rnorm(nrow(da),40,3) da obrigado André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

4 5

[LEMBRETE] Meetup R São Paulo Users Group (IME-USP) - 27 de novembro do 2015
by Luis Benites Sánchez 26 Nov '15

26 Nov '15

Caros, Amanha 27 de novembro teremos nossa reunião numero 7 do Meetup de usuários de R de São Paulo, no auditório Jacy Monteiro no Bloco B. Segue a programação: Programação da sexta 27 de novembro (14:00h-16:00h) - [14:00 - 14:55] Uma introdução aos modelos mistos em R - Ana Paula Zerbeto (IME-USP) - [Presencial] - [14:55 - 15:05] Coffee Break - [15:05 - 16:00] Web scraping usando o R - Julio Trecenti (IME-USP) - [Presencial] As palestras serão transmitidas ao vivo pela conexão IPTV da USP e YouTube ( https://plus.google.com/u/0/events/cpjg070m652k6su1gkeuivrr6bc) Para poder participar, pode confirmar sua participação gratuitamente no site do Meetup :) Para mais informações: Site: www.meetup.com/pt/useR-SP Canal de YouTube: www.youtube.com/channel/UCb9PdJYYaElE6Uqb1KxhE8w e-mail: meetuprsp(a)gmail.com Um abraço, Luis Benites Sánchez Rocío Paola Maehara Organizadores do São Paulo R Users Group --

1 0

Aplicar função em parte de um vetor
by joão Rodrigo Castro 25 Nov '15

25 Nov '15

Prezados, Estou com uma dificuldade de atribuir uma função à parte de um vetor. Suponhamos o seguinte data.frame: ano <-c(2008,2009,2010,2011,2008,2009,2010,2011,2008,2009,2010,2011,2008,2009,2010,2011) cultivar <-c("c1","c1","c1","c1","c2","c2","c2","c2","c3","c3","c3","c3","c4","c4","c4","c4") produtividade <- c(3000,4023,5432,3212,4567,6543,5467,2345,3245,99,5467,7654,3214,1234,5456,-99) exemplo <- data.frame(ano,cultivar, produtividade) Gostaria de calcular separadamente as médias das produtividades de cada uma das cultivares, para o intervalo de anos agrícolas (2008-2011). Desde já agradeço a colaboração, -- João Rodrigo de Castro Programa de Pós-Graduação em Meteorologia - Universidade Federal de Pelotas Bolsista Laboratório de Agrometeorologia - Embrapa Clima Temperado

3 5

Re: [R-br] Resíduos e MCMCglmm
by Filipe Cristovão 23 Nov '15

23 Nov '15

Olá, Muito obrigado pelas mensagens. Vou dar continuidade a discussão com alguns esclarecimentos e um pouco de código: 1- Elias pergunta como calculo os resíduos nesse modelo: >residuals_model<-(data$response_variable-predict(model)) 2- Wagner, sua pergunta é, de alguma maneira, minha grande dúvida: "Se vc já sabe que sua variável resposta é gamma pq não ajusta um modelo gamma? Ou mesmo qq outro modelo para resposta contínua como a Gaussian Inverse?”. Contudo, qual tipo de análise que não MCMCglmm me permitiria controlar meus dados quanto à filogenia? 3. E finalmente sobre os motivos de se usar MCMCglmm: recaem novamente na não independência dos dados e na possibilidade de se controlar filogeneticamente. Mas, a mesma pergunta que foi feita, é parte da minha questão: "Por que assumindo distribuição de probabilidade aos parâmetros, a interpretação dos resultados é baseada no fato de que também se obtem distribuições de probabilidade à posteriori para interpretar?” Em outras palavras, por que os resíduos devem seguir uma distribuição normal mesmo se os dados não são normais? Ou essa última frase está equivocada? O modelo segue: #Phylo é minha árvore filogenética. >MCMC1<-MCMCglmm(Response_var~var1*var2*+var3, + random = ~Phylo, family=“?”, + ginverse=list(Phylo=inv.phylo$Ainv), prior = prior, + nitt = 10000, + burnin = 1000, thin = 500, data = mob, + verbose = F ) # Os resíduos estão acima. E para checar a normalidade dos resíduos usei: > qqnorm(residuals_model) > shapiro.test(residuals_model) Enfim, após as mensagens que me fizeram estudar novos artigos =) refaço as perguntas: A) Meus resíduos precisam ser normais quando uso um modelo MCMCglmm? B) Que outro pacote me permite controlar a filogenia em um modelo? Abraços, Filipe Cristovão R. da Cunha Anthropological Institute & Museum University of Zurich - Campus Irchel Winterthurerstrasse 190, CH-8057 Zürich Switzerland www.aim.uzh.ch/Research/birdfamilies/filipe.html <http://www.aim.uzh.ch/Research/birdfamilies/filipe.html> <http://www.flickr.com/photos/filipe_cristovao>www.flickr.com/photos/filipe_cristovao <http://www.flickr.com/photos/filipe_cristovao>

3 3

Resíduos e MCMCglmm
by Filipe Cristovão 23 Nov '15

23 Nov '15

Olá, Bom dia. Talvez o assunto aqui abordado não seja tanto sobre o uso geral do software, mas uma questão teórica sobre estatística. Espero estar adequado às “guidelines” do grupo. Recentemente tenho optado por analisar meus dados usando MCMCglmm, isso se deu por várias questões, mas principalmente pela liberdade que o pacote me dá de criar meus “priors" e de poder controlar meus dados filogenéticamente. Contudo, sempre ouvi e li que “os resíduos de um modelo seguindo MCMCglmm devem ter distribuição normal”. Ai sequem algumas questões sobre isso: I) Isso é mesmo verdade absoluta? Pergunto isso porque tenho uma variável que não consigo normalizar. Essa variável apresenta uma distribuição gamma com uma obliqüidade positiva (ou "positive skewed”). Já tentei de tudo inclusive o “milagroso” boxcox. É uma variável que não se normaliza - se é que posso dizer isso =). II) Mesmo com meus resíduos não obedecendo uma distribuição normal, levei em consideração alguns dados importantes do “output” do meu modelo: i) não existe colinearidade entre as variáveis, ii) a influência da filogenia da minha amostra (ou melhor da diversidade filogenetica da minha amostra) é “irrelevante” - pois apresenta um intervalo de confiança exageradamente pequeno, iii) e os resultados são os mesmos (diferentes valores mesma direção) que quando conduzo um glm por exemplo. As questões são: a) poderia ignorar o controle filogenético? b) Ou existe outra maneira de normalizar meus dados? e por fim c) Em algum outro pacote posso controlar meu modelo quanto à filogenia e ainda assim executar meu modelo usando estatística "bayseana"? Obrigado pela paciência ao ler o email =) Atenciosamente, Filipe Cristovão R. da Cunha Anthropological Institute & Museum University of Zurich - Campus Irchel Winterthurerstrasse 190, CH-8057 Zürich Switzerland www.aim.uzh.ch/Research/birdfamilies/filipe.html <http://www.aim.uzh.ch/Research/birdfamilies/filipe.html> <http://www.flickr.com/photos/filipe_cristovao>www.flickr.com/photos/filipe_cristovao <http://www.flickr.com/photos/filipe_cristovao>

3 5

Re: [R-br] Digest R-br, volume 59, assunto 21
by Wilson Cohen 22 Nov '15

22 Nov '15

Uso R e RStudio atualizados em 3 máquinas: note com Ubuntu Trusty, note com Linux Mint 17 e PC com LInux Mint 17......tudo rodando redondinho!!.....com certeza problemas no Rwindows!! Seja Livre, Use Linux!! Em 22 de novembro de 2015 11: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. Problema na ultima versão do RStudio (Leonardo Monteiro) > 2. Re: Problema na ultima versão do RStudio (Jobenil - Gmail) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sat, 21 Nov 2015 15:28:43 +0000 > From: Leonardo Monteiro <monteiroleonardo6(a)gmail.com> > To: Lista R <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Problema na ultima versão do RStudio > Message-ID: > < > CAH4vHCHQjekqJVcS07WU+wa4rnBJBRHbARsFyhJXXP2Ub7-AqA(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Prezados, Venho rodando o R em Windows 10.. > > Fiz o update do RStudio dia 13/novembro e nao estou conseguindo a mesma > performance, ou melhor, nem estou conseguindo abrir, e quando o programa > abre, fica tudo travado e nao consigo rodar os meus codigos. > > Para desligar, eu estou conseguindo apenas pelo Gerenciador de tarefas > (ctrl+alt+del). > > Alguem sugere alguma coisa para resolver este problema? (que nao seja > trocar o OS kkk) > > obrigado. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20151121/5b72af39/atta… > > > > ------------------------------ > > Message: 2 > Date: Sat, 21 Nov 2015 13:41:05 -0200 > From: Jobenil - Gmail <pjobenil(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Problema na ultima versão do RStudio > Message-ID: <FDEE088F-7613-45F5-80E8-3E4F429CB5A9(a)gmail.com> > Content-Type: text/plain; charset=utf-8 > > Em três máquinas diferentes, todas com win 10, instalei e uso sem > problemas o RStudio e R. > Acho que você deve verificar a instalação do Windows. Quando ocorre o > travamento ele envia uma mensagem de log para Microsoft. Tente verificar > qual é o tipo de erro e consultar no próprio suporte da Microsoft como > reparar o erro. > At. > Jobenil > > Enviado do meu iPhone > > > Em 21 de nov de 2015, às 13:28, Leonardo Monteiro < > monteiroleonardo6(a)gmail.com> escreveu: > > > > Prezados, Venho rodando o R em Windows 10.. > > > > Fiz o update do RStudio dia 13/novembro e nao estou conseguindo a mesma > performance, ou melhor, nem estou conseguindo abrir, e quando o programa > abre, fica tudo travado e nao consigo rodar os meus codigos. > > > > Para desligar, eu estou conseguindo apenas pelo Gerenciador de tarefas > (ctrl+alt+del). > > > > Alguem sugere alguma coisa para resolver este problema? (que nao seja > trocar o OS kkk) > > > > obrigado. > > > > > > > > _______________________________________________ > > 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. > > > ------------------------------ > > Subject: Legenda do Digest > > _______________________________________________ > 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 59, assunto 21 > ***************************************** > -- *Prof. MSc. Wilson Martins da Silva* *Mestre em Ecologia Aquática e Pesca - PPGEAP/ICB/UFPA* *Especialista em Biodiversidade Amazônica - UFPA/ATMDepto de Morfologia e Ciências Fisiológicas* *Centro de Ciências Biológicas e da Saúde* *Universidade do Estado do Pará* *Campus IX - Altamira - Pará* http://lattes.cnpq.br/2440617911389759 wilsoncelula(a)gmail.com, wilsonms(a)ufpa.br

1 0

Problema na ultima versão do RStudio
by Leonardo Monteiro 21 Nov '15

21 Nov '15

Prezados, Venho rodando o R em Windows 10.. Fiz o update do RStudio dia 13/novembro e nao estou conseguindo a mesma performance, ou melhor, nem estou conseguindo abrir, e quando o programa abre, fica tudo travado e nao consigo rodar os meus codigos. Para desligar, eu estou conseguindo apenas pelo Gerenciador de tarefas (ctrl+alt+del). Alguem sugere alguma coisa para resolver este problema? (que nao seja trocar o OS kkk) obrigado.

2 1

Re: [R-br] : Extensão .DBC
by sznelwar 19 Nov '15

19 Nov '15

1 0

📊 Curso Básico de R + RStudio 📈
by Leonardo Aguirre 18 Nov '15

18 Nov '15

> De: Leonardo Aguirre <leoaguirre32(a)gmail.com> > Data: 18 de novembro de 2015 10:35:48 BRST > Assunto: 📊 Curso Básico de R + RStudio 📈 > > Pessoal, > > Segue divulgação de um curso de R em Brasília, caso alguém tenha interesse. > As inscrições podem ser feitas em http://oddsactions.com.br/cursos/inscricao/ > > Abraços, > > Enviado do meu iPhone > > Início da mensagem encaminhada > >> De: Odds&Actions - Inteligência Analítica <contato(a)oddsactions.com.br> >> Data: 18 de novembro de 2015 09:59:55 BRST >> Para: <leoaguirre32(a)gmail.com> >> Assunto: 📊 Curso Básico de R + RStudio 📈 >> Responder A: Odds&Actions - Inteligência Analítica <contato(a)oddsactions.com.br> >> >> >> Inscrições abertas para o Curso de R e RStudio Veja esse e-mail em seu navegador. >> >> Curso Básico de R + RStudio >> Aprenda a programar na língua universal de análise de dados. >> >> Quem deveria aprender R? >> >> Se você trabalha com análises de dados, sejam eles de planilhas de Excel, dados de pesquisas, bancos de dados da sua empresa ou qualquer outra fonte de informação, o R é para você. E a melhor parte: você não precisa saber programar para aprender R. >> >> O que o R é capaz de fazer? >> >> Manusear bases de dados, gerar relatórios estatísticos com tabelas, gráficos e mapas profissionais, criar aplicações estilo dashboards, criar modelos preditivos, fazer análises de redes sociais, trabalhar com Big Data e muito mais. >> >> Como funciona o curso? >> >> Nossa próxima turma será dos dias 7 a 11 de dezembro, no período da manhã (9:00-12:00), na SCLRN 710 Bloco A. Nesses dias, ensinaremos você: a estrutura básica da linguagem R; como importar, manipular e exportar dados; como gerar estatísticas básicas; como gerar seus primeiros gráficos; e como escrever suas próprias funções. Faremos vários exercícios ao longo do curso e espero que você se divirta com nossos desafios! >> >> Investimento: R$ 700 (setecentos reais) por pessoa. >> (à vista ou em 2x de R$ 350) >> Não é necessário levar notebook. >> >> Que tal aprender a programar em R? Inscreva-se agora! >> Essa é uma turma especial e estamos abrindo apenas 5 vagas! >> Inscreva-se agora >> Odds&Actions © 2015 Odds&Actions - Inteligência Analítica, Todos os direitos reservados. >> Você se inscreveu em nosso mailing pelo site http://oddsactions.com.br/cursos >> >> Nosso endereço: >> Odds&Actions - Inteligência Analítica >> SCN Quadra 2 Bloco D Entrada B Sala 823 Parte A >> Brasília, DF 70712-903 >> Brazil >> >> Add us to your address book >> me remova dessa lista editar preferências de recebimento >> >>

1 0

ajuste com glmmADMB
by Iábita Fabiana 18 Nov '15

18 Nov '15

Boa noite! Estou tentando analisar um conjunto de dados de germinação, utilizando o pacote glmmADMB, o meu conjunto de dados é como o que se segue: dose dias Rep germ Incrustrado 0 1 1 0 Sim 0 2 1 0 Sim 0 3 1 0 Sim 0 4 1 0 Sim 0 5 1 12 Sim 0 6 1 13 Sim 0 7 1 0 Sim 0 8 1 4 Sim 0 9 1 0 Sim 0 10 1 1 Sim 0 11 1 3 Sim 0 12 1 0 Sim 0 1 2 0 Sim (...) (...) (...) (...) 0 12 6 0 Sim (...) (...) (...) (...) 16 6 6 2 Sim 16 7 6 0 Sim 16 11 6 0 Sim 16 12 6 0 Sim 0 1 1 8 Não 0 2 1 23 Não (...) (...) (...) (...) 0 11 6 0 Não 0 12 6 0 Não (...) (...) (...) (...) (...) 8 12 4 0 Não 16 10 6 0 Não 16 11 6 0 Não 16 12 6 0 Não Em que Dose= são 6 diferentes doses (0,1,2,4,8,16) de composto para acelerar a germinação Rep= 6 repetições (gerbox com 50 sementes cada) Germ= quantidade de sementes germinadas no dia Dias= dias em que foi feita a contagem Incrustrado= Informação se era incrustrada ou não A cada dia, durante 12 dias era observado a quantidade de sementes germinadas. Os dados são inflacionados de zeros. Então estou tentando analisar usando uma distribuição binomial negativa com inflação de zeros. Com o seguinte código para ajustar, porém estou com dúvida quanto se essa seria a especificação correta do modelo, pois obtive umas estimativas negativas que não fazem sentido. Se alguém tiver mais familiaridade nesse tipo de análise qualquer ajuda seria bem vinda. M1<-glmmadmb(germ~(dose+dias)*incrustrado,random=~(1|rep),data=da,family="nbinom1",zeroInflation=TRUE) Desde já agradeço a atenção, Fabiana

3 4

Extensão .DBC
by Rafael Costa 17 Nov '15

17 Nov '15

Prezados, Existe alguma forma de ler o formato .dbc diretamente no R? A forma que faço geralmente é utilizar o programa tabwin para converter .dbc para .dbf. O problema todo reside no fatos que os dados ficam muito pesados. Certo de sua colaboração, agradeço desde já! Respeitosamente, -- *Rafael Carneiro da CostaDoutorando em EconomiaCentro de Pós-Graduação em Economia - CAENUniversidade Federal do Ceará - UFC*

5 13

Excluir linhas com NA de um dataframe
by Yury Duarte 17 Nov '15

17 Nov '15

Boa tarde! Estou tentando excluir as linhas de uma coluna em especifico, na qual ocorrem algumas falhas de dados, representadas por NA. Tentei o comando: is.na(dados$coluna) Com ele, obtive uma saída no formato T/F. Gostaria de saber como aplico ao meu objeto "dados" apenas os valores de TRUE, onde não há falha de dados. Desde já, agradeço pela ajuda! Yury Duarte Engenheiro Agrônomo - ESALQ/USP

5 6

Laço para gerar Subconjuntos
by Michelle Bau Graczyk 16 Nov '15

16 Nov '15

Caros, bom dia… Estou com uma dúvida bem boba mas não consigo resolver! Eu tenho um data frame que a primeira coluna é a data, a segunda é o horário e as outras seis são dados como preço e volume. Os dias estão em forma numérica e eu os deixei assim porque achei que seria mais fácil fazer o laço. O objetivo é que para cada daia eu ache o valor inicial e final do preço assim como o seu máximo e seu mínimo. Estou tendo problemas somente em trocar de dia. Abaixo segue um pedaço reduzido do data frame e o código. Alguém poderia me dar uma dica? Muito obrigada! AA<-read.table("/Users/bau/Dropbox/AA.N_ExemploReduzido.csv",header=TRUE + ,dec=".",as.is=TRUE, sep=",") > #AA > dput(AA) structure(list(X37988 = c(37988L, 37988L, 37988L, 37988L, 37988L, 37988L, 37991L, 37991L, 37991L, 37991L, 37992L, 37992L, 37992L, 37992L), X34244.059 = c(34251.038, 34263.36, 34296.73, 34301.369, 34318.91, 34324.707, 34398.295, 34400.585, 34416.511, 34425.892, 34250.443, 34252.127, 34266.757, 34273.928), X.1 = c(-1, -1, 38.05, 38.06, 38.06, 38.04, 37.8, 37.8, 37.85, 37.85, -1, -1, 38.85, 38.85), X.1.1 = c(-1L, -1L, 32L, 2L, 4L, 1L, 50L, 50L, 1L, 1L, -1L, -1L, 1L, 1L), X.1.2 = c(-1, -1, 38.08, 38.08, 38.07, 38.07, 38.3, 38.3, 37.95, 37.95, -1, -1, 38.95, 38.95), X.1.3 = c(-1L, -1L, 1L, 1L, 1L, 1L, 203L, 203L, 1L, 1L, -1L, -1L, 1L, 1L), X38 = c(38, 38, 38.05, 38.06, 38.06, 38.05, 37.9, 37.9, 37.9, 37.9, 38.9, 38.9, 38.9, 38.9), X55900 = c(600L, 5000L, 6000L, 1500L, 500L, 4500L, 161800L, 500L, 3700L, 1100L, 125000L, 1100L, 400L, 1000L )), .Names = c("X37988", "X34244.059", "X.1", "X.1.1", "X.1.2", "X.1.3", "X38", "X55900"), class = "data.frame", row.names = c(NA, -14L)) > > matriz<-matrix(NA,230,4)## numero de linhas é o numero de dias. > colnames(matriz)<-c("open","close","low","higth") > i<-1 > for( x in 37988:38000){ + + + #x<-37988 + + if(AA$data==x){ + AA1<-subset(AA,AA$data==x) + matriz[i,1]<-AA1[1,7] + matriz[i,2]<-AA1[nrow(AA1),7] + matriz[i,3]<-min(AA1$last) + matriz[i,4]<-max(AA1$last) + + i=i+1 + } + AA<-read.table("/Users/bau/Dropbox/AA.N_ExemploReduzido.csv",header=TRUE + ,dec=".",as.is=TRUE, sep=",") + #AA + dput(AA) + + matriz<-matrix(NA,3,4)## numero de linhas é o numero de dias. + colnames(matriz)<-c("open","close","low","higth") + i<-1 + for( x in 37988:37992){ + + + + + if(AA$data==x){ + AA1<-subset(AA,AA$data==x) + matriz[i,1]<-AA1[1,7] + matriz[i,2]<-AA1[nrow(AA1),7] + matriz[i,3]<-min(AA1$last) + matriz[i,4]<-max(AA1$last) + + i=i+1 + } + } + matriz + write.table(matriz,"Matriz4ParametrosPreço.txt") +

5 11

Reorganizar um data.frame
by Michelle Bau Graczyk 15 Nov '15

15 Nov '15

Caros, boa noite, Eu usei a função dcast() para reorganizar meu data frame pois meu objetivo era calcular a media do volume por hora. Eu tinha: AA1 data hora preco volume horaid 1 01.07.2004 09:31 51.37 27900 1 2 01.07.2004 09:32 51.32 100 2 3 01.07.2004 09:33 51.45 700 3 e fiz ficar: require(reshape2) AA2 <- dcast(AA, data~hora, value.var='volume') head(AA2[,1:8], 3) # data 09:31 09:32 09:33 09:34 09:35 09:36 09:37 # 1 2004-07-01 27900 100 700 1000 800 3000 1600 # 2 2004-07-02 145400 500 NA 500 6900 NA 4400 # 3 2004-07-06 3200 2300 3800 1500 NA 1100 3700 Porém, quero dividir o meu volume de AA1 pela media por hora de AA2, ou seja, para todos os volumes de 09:31 de AA1, quero dividir pela média de AA2[,2](media do volume da hora 09:31). Eu consigo fazer isso no formato de AA2: seria AA2[,2]/mean(AA2[,2]).Porém, eu quero que fique no formato de AA1. Alguém saberia me dar uma dica como eu devo proceder? Como a quantidade de horários e cada dia são diferentes, não consigo fazer de uma forma mais simples. Muito obrigada, Michelle

1 0

Re: [R-br] RES: Laço para gerar Subconjuntos
by sznelwar 12 Nov '15

12 Nov '15

1 0

[OFF-TOPIC] Teste sobre funcionamento da lista
by Walmes Zeviani 11 Nov '15

11 Nov '15

Mensagem de teste para avaliar se voltou a funcionar o serviço de email que hospeda a lista de discussão. W. ======================================================================== Walmes Marques Zeviani LEG - Laboratório de Estatística e Geoinformação [25.450418S 49.231759W] Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 linux user number: 531218 pages & profiles: http://www.leg.ufpr.br/~walmes | http://lattes.cnpq.br/4410617539281650 http://git.leg.ufpr.br/walmes | https://github.com/walmes https://www.mendeley.com/profiles/walmes-zeviani1/ ========================================================================

10 9

Re: [R-br] Laço para gerar Subconjuntos
by Robert Iquiapaza 11 Nov '15

11 Nov '15

Pode fazer com laço ou com funções da família apply: matriz<-matrix(NA,3,4)## numero de linhas é o numero de dias. colnames(matriz)<-c("open","close","low","higth") i<-1 for( x in unique(AA$X37988)){# x=37988 AA1<-subset(AA,AA$X37988==x) matriz[i,1]<-AA1[1,7] matriz[i,2]<-AA1[nrow(AA1),7] matriz[i,3]<-min(AA1[,7]) matriz[i,4]<-max(AA1[,7]) i=i+1 } matriz #ou t(sapply(split(AA,AA$X37988),function(AA1)c(AA1[1,7], AA1[nrow(AA1),7], min(AA1[,7]), max(AA1[,7])))) [,1] [,2] [,3] [,4] 37988 38.0 38.05 38.0 38.06 37991 37.9 37.90 37.9 37.90 37992 38.9 38.90 38.9 38.90 Att Em 11 de novembro de 2015 13:00, Michelle Bau Graczyk <mbgraczyk(a)gmail.com> escreveu: > Caros, bom dia… > > Estou com uma dúvida bem boba mas não consigo resolver! > > Eu tenho um data frame que a primeira coluna é a data, a segunda é o > horário e as outras seis são dados como preço e volume. > Os dias estão em forma numérica e eu os deixei assim porque achei que > seria mais fácil fazer o laço. O objetivo é que para cada daia eu ache o > valor inicial e final do preço assim como o seu máximo e seu mínimo. Estou > tendo problemas somente em trocar de dia. > Abaixo segue um pedaço reduzido do data frame e o código. > Alguém poderia me dar uma dica? > Muito obrigada! > AA<-read.table("/Users/bau/Dropbox/AA.N_ExemploReduzido.csv",header=TRUE > + ,dec=".",as.is=TRUE, sep=",") > > #AA > > dput(AA) > structure(list(X37988 = c(37988L, 37988L, 37988L, 37988L, 37988L, > 37988L, 37991L, 37991L, 37991L, 37991L, 37992L, 37992L, 37992L, > 37992L), X34244.059 = c(34251.038, 34263.36, 34296.73, 34301.369, > 34318.91, 34324.707, 34398.295, 34400.585, 34416.511, 34425.892, > 34250.443, 34252.127, 34266.757, 34273.928), X.1 = c(-1, -1, > 38.05, 38.06, 38.06, 38.04, 37.8, 37.8, 37.85, 37.85, -1, -1, > 38.85, 38.85), X.1.1 = c(-1L, -1L, 32L, 2L, 4L, 1L, 50L, 50L, > 1L, 1L, -1L, -1L, 1L, 1L), X.1.2 = c(-1, -1, 38.08, 38.08, 38.07, > 38.07, 38.3, 38.3, 37.95, 37.95, -1, -1, 38.95, 38.95), X.1.3 = c(-1L, > -1L, 1L, 1L, 1L, 1L, 203L, 203L, 1L, 1L, -1L, -1L, 1L, 1L), X38 = c(38, > 38, 38.05, 38.06, 38.06, 38.05, 37.9, 37.9, 37.9, 37.9, 38.9, > 38.9, 38.9, 38.9), X55900 = c(600L, 5000L, 6000L, 1500L, 500L, > 4500L, 161800L, 500L, 3700L, 1100L, 125000L, 1100L, 400L, 1000L > )), .Names = c("X37988", "X34244.059", "X.1", "X.1.1", "X.1.2", > "X.1.3", "X38", "X55900"), class = "data.frame", row.names = c(NA, > -14L)) > > > > matriz<-matrix(NA,230,4)## numero de linhas é o numero de dias. > > colnames(matriz)<-c("open","close","low","higth") > > i<-1 > > for( x in 37988:38000){ > + > + > + #x<-37988 > + > + if(AA$data==x){ > + AA1<-subset(AA,AA$data==x) > + matriz[i,1]<-AA1[1,7] > + matriz[i,2]<-AA1[nrow(AA1),7] > + matriz[i,3]<-min(AA1$last) > + matriz[i,4]<-max(AA1$last) > + > + i=i+1 > + } > + AA<-read.table("/Users/bau/Dropbox/AA.N_ExemploReduzido.csv",header=TRUE > + ,dec=".",as.is=TRUE, sep=",") > + #AA > + dput(AA) > + > + matriz<-matrix(NA,3,4)## numero de linhas é o numero de dias. > + colnames(matriz)<-c("open","close","low","higth") > + i<-1 > + for( x in 37988:37992){ > + > + > + > + > + if(AA$data==x){ > + AA1<-subset(AA,AA$data==x) > + matriz[i,1]<-AA1[1,7] > + matriz[i,2]<-AA1[nrow(AA1),7] > + matriz[i,3]<-min(AA1$last) > + matriz[i,4]<-max(AA1$last) > + > + i=i+1 > + } > + } > + matriz > + write.table(matriz,"Matriz4ParametrosPreço.txt") > + > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. >

1 0

Intervalo de confiança das estimativas dos parâmetros das EDO geradas pela função modFit, do pacote FME
by Adriele Giaretta Biase 11 Nov '15

11 Nov '15

Como pode ser calculado o intervalo de confiança das estimativas dos parâmetros das EDO (ou modelos não lineares) geradas pela função modFit, do pacote FME? Abaixo segue o summay obtido Estimate Std. Error t value Pr(>|t|) Beta1 4.285 0.0041 0.411 0.00687 Beta2 0.082 0.0769 0.530 0.00605 Beta3 2.041 0.0900 0.411 0.00068 Eu devo usar a distribuição t-student retornado pela função, fazendo: Estimate ± Std. Error * t value Em que CI_95%(Beta1)= (4.285- 0.0041* 0.411 ; 4.285+ 0.0041* 0.411 ) CI_95%(Beta2)= (0.082- 0.0769* 0.530 ; 0.082+ 0.0769* 0.530) CI_95%(Beta3)= (2.041- 0.0900* 0.411; 2.041+ 0.0900* 0.411) Por que os valores da distribuição t estão variando, os graus de liberdade não deveriam ser os mesmos (n-p)? Eem que o n é o tamanho da amostra e p o número de parâmetros? Ou devo usar a distribuição normal Z(α/2): CI_95%(Beta1)= (4.285- 0.0041* 1.96 ; 4.285+ 0.0041* 1.96 ) Onde eu posso conseguir um material falando sobre o cálculo dos Intervalos de Confiança produzidos pela função modFit do pacote FME. O manual do pacote não deixa claro esse cálculo. O cálculo dos intervalos sempre serão baseados em distribuições simétricas? Adriele -- Adriele Giaretta Biase. Mestre em Estatística e Experimentação Agropecuária - UFLA. Doutoranda em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 8861-0619.

1 0

Meetup R São Paulo Users Group (IME-USP) - 27 de novembro do 2015
by Luis Benites Sánchez 11 Nov '15

11 Nov '15

Caros, A sexta 27 novembro teremos nossa reunião número 7 do Meetup de usuários de R de São Paulo, no auditório Jacy Monteiro no Bloco B. Segue a programação desse dia: Programação da sexta *27 de novembro* (14h-16:00h) - [14:00 - 14:55] Uma introdução aos modelos mistos em R - Ana Paula Zerbeto (IME-USP) - [Presencial] - [14:55 - 15:05] Coffee Break - [15:05 - 16:00] Web scraping usando o R - Julio Trecenti (IME-USP) - [Presencial] As palestras serão transmitidas ao vivo pela conexão IPTV da USP e YouTube (*https://plus.google.com/u/0/events/cpjg070m652k6su1gkeuivrr6bc <https://plus.google.com/u/0/events/cpjg070m652k6su1gkeuivrr6bc>*) . Para poder participar, pode confirmar sua participação gratuitamente no site do Meetup :) Para mais informações: Site: www.meetup.com/pt/useR-SP Canal de YouTube: www.youtube.com/channel/UCb9PdJYYaElE6Uqb1KxhE8w e-mail: meetuprsp(a)gmail.com Um abraço, Luis Benites Sánchez Rocío Paola Maehara Organizadores do São Paulo R Users Group --

1 0

Trabalhar com subsets de data frame
by Michelle Bau Graczyk 11 Nov '15

11 Nov '15

Caros, bom dia, Eu estou com um data frame da forma : "data" "hora" "bid" "bidQ" "ask" "askQ" "last" "volume" "1" 2004-01-04 "09:30" -1 -1 -1 -1 38 55900 "2" 2004-01-04 "09:30" -1 -1 -1 -1 38 600 "3" 2004-01-04 "09:31" -1 -1 -1 -1 38 5000 "4" 2004-01-04 "09:31" 38.05 32 38.08 1 38.05 6000 "5" 2004-01-04 "09:31" 38.06 2 38.08 1 38.06 1500 "6" 2004-01-04 "09:31" 38.06 4 38.07 1 38.06 500 "7" 2004-01-04 "09:32" 38.04 1 38.07 1 38.05 4500 "8" 2004-01-04 "09:32" 38 1 38.2 1 38.03 1600 "9" 2004-01-04 "09:32" 38 1 38.2 1 38.03 3800 "10" 2004-01-04 "09:32" 38 1 38.04 1 38.04 500 "11" 2004-01-04 "09:32" 38 1 38.07 1 38.03 200 "12" 2004-01-04 "09:32" 38 1 38.08 2 38.07 100 "13" 2004-01-04 "09:32" 38 1 38.08 2 38.01 5000 "14" 2004-01-04 "09:32" 38 21 38.08 2 38.05 200 "15" 2004-01-04 "09:33" 38.03 9 38.13 8 38.08 700 "16" 2004-01-04 "09:33" 38.08 15 38.13 8 38.08 100 "17" 2004-01-04 "09:33" 38.11 16 38.13 10 38.12 900 onde as datas vão de 2004-01-04 até 2013-06-30 e os horários de 09:30 até 16:10. Eu tenho que achar o volume médio para cada minuto e depois normalizar esses volumes para cada minuto pelo volume médio. Eu tento montar um programa cheio de for() para data e for() para horário mas não está dando certo… Acho que deve ter uma lógica mais simples. Alguém saberia me dar uma dica? Muito obrigada! Michelle

2 2

equacao de S.E. para regressao
by Cleber Borges 10 Nov '15

10 Nov '15

Prezados UseRs, Usando o livro do Faraway (practical regression and anova using r) adaptei um trecho de código (abaixo) e consigo traçar um gráfico de modelo de regressão (reta) e com os intervalos de confiança (IC) da regressão. Com uso da função *predicit* eu obtenho os IC... Minha perguntinha é: Tem alguma formulazinha bonitinha e elegante para se obter esses *s.e.* da regressão tal como são bonitinhas e elegantes as formulazinhas de estimativas de coeficientes dos parâmetros e erros segundo OLS?? tipo: 1) b = ( X^t * X )^1 * X^t * y 2) Var( b ) = ( X^t * X )^1 * Var_residual 3) Algo do gênero de 1) e 2) ??? Agradeço pelo help com antecedência! Cleber ps.: vi uma formula mas é muita feia! ############ x <- rnorm( 100 ) y <- x + rnorm( 100, s=5 ) reg <- lm( y ~ x ) grid <- seq( -3, 3, 0.1 ) p <- predict( reg, data.frame( x=grid ), se=T ) IC <- cbind( p$fit - 1.96*p$se, p$fit, p$fit + 1.96*p$se ) matplot( grid, IC, lty=1, col=c(2,1,2), type="l", lwd=3 ) --- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus

2 2

Índice de concordância Kappa
by André Lucas de Oliveira Moreira 09 Nov '15

09 Nov '15

Olá pessoal, tudo bem? Aqui no trabalho, pela primeira vez me deparei com uma situação em que é necessário comparar a eficiência de técnicas que detectam que uma mesma espécie de parasita. A principio não sabia qual seria o melhor teste para verificar isso, mas me falaram pela primeira vez do kappa. Então tenho buscado sobre o assunto. Porém, nos livros que eu tive acesso fala disso muito rapidamente e não me deram segurança suficiente para executar no R e interpretar os resultados. Pensando que tenho resultados categóricos (presença e ausência dos parasitos para ambas técnicas), qual seria a melhor forma de analisar isso, melhor pacote e função? Também preciso do valor de p. Seria muito útil se alguém indicasse um tutorial e também algum material teórico que me dê mais segurança no uso desta técnica. Agradecido, André -- *MSc. André Lucas de O. Moreira* http://lattes.cnpq.br/7258065668864153 <http://www.wikiaves.com.br/perfil_andrelukinhas> http://www.wikiaves.com.br/perfil_andrelukinhas 79 8837-3562 79 9132-9093

5 10

Como aumentar o espaço de strip, para caber um nome de 2 linhas. bwplot (lattice)
by Augusto Ribas 09 Nov '15

09 Nov '15

Ola, estou com dificuldades de escrever o nome o titulo dos strips numa figura usando o lattice. Meu problema é o seguinte, eu tenho os seguintes dados, por exemplo resposta<-rnorm(90) preditor1<-rep(rep(c("a","b"),each=15),3) preditor2<-rep(c("sp1","sp2","sp3"),each=30) Dai preciso fazer uma figura, mas os nomes nos strips são espécies biologicas, então precisam estar em italico, resolvi da seguinte forma. library(lattice) bwplot(resposta~preditor1|preditor2, strip = strip.custom(factor.levels = c( expression(italic("P. paradoxa")), expression(italic("H. raniceps")), expression(italic("P. azurea"))) ) ) Mas além disso, gostaria de pular uma linha, e escrever uma segunda informação, mas essa segunda informação não é em italico, mas a primeira precisa. Pesquisando uma pouco, vi a possibilidade de usar o atop do expression, mas ai o texto não cabe no strip, como eu posso aumentar o tamanho do strip, para caber tudo? bwplot(resposta~preditor1|preditor2, strip = strip.custom(factor.levels = c( expression(atop(italic("P. paradoxa"),"outra info")), expression(atop(italic("H. raniceps"),"outra info")), expression(atop(italic("P. azurea"),"outra info"))) ) ) Bem é isso, agradeço desde ja a atenção. -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.com.br/ <http://augustoribas.heliohost.org> Github: https://github.com/Squiercg Lattes: http://lattes.cnpq.br/7355685961127056

2 2

Auxilio com Funções de Estatistica
by Yury Duarte 09 Nov '15

09 Nov '15

Bom dia! Meu nome é Yury Duarte e sou pós graduando em Engenharia de Biossistemas. Estou com uma necessidade de aplicar alguns índices estatísticos a uma serie de dados, entretanto, não consegui encontrar uma função no R que as faça. Os índices são: - R² - r - c - D - RMSE - Erro médio - Erro médio absoluto Dessa forma, gostaria de saber se já existem comandos para calcular esses índices ou se terei que desenvolver o código para cada um deles. Desde já, agradeço pela colaboração! Att Yury Duarte Engenheiro Agrônomo - ESALQ/USP

3 3

Re: [R-br] Chi-quadrado demonstra associação, mas e daí?
by andremoreirazoo 08 Nov '15

08 Nov '15

Oi pessoal, agradeço pela atenção e ajuda. Francisco, possuo dados categóricos apenas, presença ou ausência de parasitas. Felipe, acho que o caminho é esse mesmo, mas você poderia me ajudar com comentários nas linhas dos comandos? Apesar de gostar de estudar estatística a matemática da coisa não é o meu forte. Então, uma ajuda para sabe do que se tratacada comando/resultado me faria entender melhor. Obrigado a todos, André Enviado do tablet Samsung. -------- Mensagem original -------- De : putini(a)outlook.com Data:03/11/2015 19:45 (GMT-03:00) Para: r-br@listas.c3sl.ufpr.br,r-br@listas.c3sl.ufpr.br Assunto: Re: [R-br] Chi-quadrado demonstra associação, mas e daí? Não seria o caso de fazer uma matriz de correlação, para estabelecer a intensidade da relação entre dois fatores? De qualquer maneira esse é um passo para fazer a regressão, estou certo? Att. Francisco Putini Enviado pelo Outlook para Android On Tue, Nov 3, 2015 at 1:44 PM -0800, "Leonardo Ferreira Fontenelle" <leonardof(a)leonardof.med.br> wrote: Em Ter 3 nov. 2015, às 17:21, Leonardo Ferreira Fontenelle escreveu: [...] além do prop.diff(), que já vem de fábrica. [...] Ops! A função é prop.test(), não prop.diff(). Leonardo Ferreira Fontenelle

6 7

Re: [R-br] Média ponderada de salário para cada idade
by Thaís Brenda Martins 08 Nov '15

08 Nov '15

Boa noite pessoal, quando escrevi fiquei com medo de vocês terem essa dúvida e realmente aconteceu. Mas o que eu quero, conforme o Leandro e o Cesar falaram, é calcular a média de salário ponderada (pelo peso) para cada idade (que vai de 12 a 100 anos). Tem alguma função que faça isso? Cesar vou dar uma olhada nessa documentação da família de funções apply do R.Muito obrigada pessoal, espero ter explicado melhor.Abraços. Att.,Thaís Brenda Martins

2 1

Média ponderada de salário para cada idade
by Thaís Brenda Martins 06 Nov '15

06 Nov '15

Boa tarde pessoal, Tenho três vetores no R: idade, salario e peso.Preciso fazer uma média ponderada desses salários por idade. Por exemplo:Idade Média Salarial16 20017 35018 40019 800...Se fosse apenas a média salarial, desconsiderando a idade da pessoa usaria o comando weighted.mean(salario,peso). Mas preciso da média por idade. (os vetores são correspondentes, por exemplo, a posição 1 de cada vetor corresponde a mesma pessoa). O vetor de idade não está ordenado. Como posso fazer isso no R? Obrigada, Att.,Thaís Brenda Martins

3 2

Teste de distribuição de erros tabajara
by ASANTOS 06 Nov '15

06 Nov '15

Prezados listeiros, Estou tentando criar uma função para testar a distribuição de erros de normal, poisson e binomial negativa por enquanto, porém não sei o que estou fazendo de errado na minha função e solicito a ajuda de vocês para conseguir achar o erro, segue CRM: ## Teste de distribuição de erros tabajara require(MASS) BD<-rpois(1000,10) ### Valores a serem testados Dist<-c("negative binomial","normal","Poisson") for(typeD in 1:Dist){ k <- fitdistr(BD,typeD[Dist]) par <- k$estimate size <- par[1]#Parametro mu <- par[2]#Média SD<-sd(BD) N <- length(BD) if (typeD[Dist]=='negative binomial'){ ### Para binomial negativa est <-N*dnbinom(BD,size=size,mu=mu) ## Estimativas fecdf <- ecdf(BD) ###ecdf- Função de distribuição cumulativa empírica knotsX <- knots(fecdf) emp <- fecdf(c(knotsX,Inf)) # Empírica chisq.test(table(emp),table(est),correct=TRUE) ##Teste Chi - Quadrado para valiar a significância da distribuição } if (typeD[Dist]=='normal'){ ### Para normal est <-N*dnorm(BD,mean=mu, sd=SD) ## Estimativas fecdf <- ecdf(BD) ###ecdf- Função de distribuição cumulativa empírica knotsX <- knots(fecdf) emp <- fecdf(c(knotsX,Inf)) # Empírica chisq.test(table(emp),table(est),correct=TRUE) ##Teste Chi - Quadrado para valiar a significância da distribuição } if (typeD[Dist]=='Poisson'){ ### Para Poisson est <-N*dpois(BD,lambda=mu) ## Estimativas fecdf <- ecdf(BD) ###ecdf- Função de distribuição cumulativa empírica knotsX <- knots(fecdf) emp <- fecdf(c(knotsX,Inf)) # Empírica chisq.test(table(emp),table(est),correct=TRUE) ##Teste Chi - Quadrado para valiar a significância da distribuição } } # Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

2 2

Re: [R-br] OFF TOPIC - Criar função no emacs para trocar # por ###
by Rafael Tieppo 05 Nov '15

05 Nov '15

Prezado, espero colaborar,Favor add no seu '.emacs' ;;----------------------------------------------------------------------------- ;; Define C-y para fazer linha com 3 sinais de # (global-set-key [?\C-y] (kbd "C-u 0 3 #")) ;;----------------------------------------------------------------------------- Vc pode trocar o y pelo que for mais conveniente para vc. obs: adaptado do Walmes Saudações Rafael Tieppo State University of Mato Grosso - Department of Agricultural Engineering site: http://www2.unemat.br/rafaeltieppo blog: https://sistemasagricolas.wordpress.com "Evite o desperdício: antes de imprimir pense na sua responsabilidade com o ambiente". On Thursday, November 5, 2015 11:00 AM, "r-br-request(a)listas.c3sl.ufpr.br" <r-br-request(a)listas.c3sl.ufpr.br> wrote: 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: Chi-quadrado demonstra associação, mas e daí? (André Figueiras Rutz) 2. Re: Chi-quadrado demonstra associação, mas e daí? (Felipe) 3. Re: quem é o autor de tanks.R ??? (Marcos Vital) 4. Re: Campanha de crowdfunding envolvendo o R e Bioestatística (Marcos Vital) 5. OFF TOPIC - Criar função no emacs para trocar # por ### (Augusto Ribas) ---------------------------------------------------------------------- Message: 1 Date: Wed, 4 Nov 2015 16:31:08 -0200 From: André Figueiras Rutz <rutzaf(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Chi-quadrado demonstra associação, mas e daí? Message-ID: <CAE7MPajKj=PL_67o9n-m7bh=8NJxENkErQhUpS-h7my1HWyy1Q(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Caro Leonardo, Acompanhando as perguntas e respostas aprendi coisas interessantes como referência de início de estudo. Sua vontade em ajudar foi sim bastante útil. Obrigado pela sua contribuição. André Rutz Psicólogo Em 04/11/2015 11:53, "Leonardo Ferreira Fontenelle" < leonardof(a)leonardof.med.br> escreveu: > Em Ter 3 nov. 2015, às 22:33, andremoreirazoo escreveu: > > [...] > > Felipe, acho que o caminho é esse mesmo, mas você poderia me ajudar com > comentários nas linhas dos comandos? > > Apesar de gostar de estudar estatística a matemática da coisa não é o > meu forte. Então, uma ajuda para sabe do que se tratacada > comando/resultado me faria entender melhor. > > > André, > > Partindo do princípio de que você saiba os fundamentos de como usar o R, > você pode completar seu conhecimento lendo a ajuda para cada comando que > Felipe usou. Mesmo assim, eu não recomendo que você escreva código assim. > Como você disse, a parte matemática não é seu forte; se você tiver > conhecimento suficiente para escrever alguma coisa, ela provavelmente já > foi escrita. Pior ainda, você pode acabar deixando de usar uma alternativa > melhor, pois existe mais de uma forma de calcular intervalos de confiança. > > Por exemplo, o código que Felipe propões calcula o intervalo de confiança > da diferença de proporções usando o que se costuma chamar método de Wald. > Só que o R já dispõe de uma função para calcular o intervalo de confiança > para a diferença de proporções, que é o prop.test(). Essa função utiliza o > método dos escores de Wilson, que é uma alternativa superior ao método de > Wald > <http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.408.7107&rep=rep1&…>, > especialmente quando a amostra é pequena o suficiente para que isso faça > alguma diferença. > > Repito que existem pacotes como o "PropCIs" e o "binom" com funções que > provavelmente vão atender às suas necessidades, além dos pacotes de > epidemiologia. > > Se você precisa de ajuda para interpretar as estimativas fornecidas pelo > R, você precisa de um livro de estatística (ou epidemiologia). Se você > precisa ajuda para chegar a essas estimativas com o R, a gente pode lhe > ajudar. Mas, agora me ocorreu, talvez você nem precise do R. Dependendo da > complexidade das análises que você pretende realizar, você pode usar o > www.openepi.com. > > Espero ter ajudado mais do que confundido :) > > Leonardo Ferreira Fontenelle <http://lattes.cnpq.br/9234772336296638> > > > _______________________________________________ > 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/20151104/439fc86a/atta…> ------------------------------ Message: 2 Date: Wed, 4 Nov 2015 18:33:28 -0200 From: Felipe <felipe.e.barletta(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Chi-quadrado demonstra associação, mas e daí? Message-ID: <563A6B98.5070403(a)gmail.com> Content-Type: text/plain; charset="utf-8"; Format="flowed" André, Como o Leonardo disse no e-mail anterior, há pacotes que já calculam medidas como diferença de proporção OR, seus respectivos IC e outras medidas que podem atender suas necessidades no seu estudo. Além dos pacotes que ele já sugeriu, outro que pode consultar é o epiR: https://cran.r-project.org/web/packages/epiR/epiR.pdf Outra sugestão de leitura que gostaria de é o material da professora Silvia Shimakura: http://leg.ufpr.br/~silvia/CE008/ http://leg.ufpr.br/~silvia/CE001/node68.html Veja qual forma se apresenta mais interessante para seu aprendizado, mas quando escrevo as funções no R como calculadora, acredito que os exemplos se tornam mais didáticos mesmo que já implementados em alguns pacotes do R. E como solicitou segue alguns comentários acerca dos comandos que enviei anteriormente: ## Carregando os dados da tabela que enviou no e-mail dados<-matrix(c(250,15,34,14),ncol=2,byrow=T) ## Verificando a existência de associação entre os parasitas através da Estatística Qui-quadrado ## Quando utilizamos o teste o argumento sim=500, há um alerta pois há casela com frequência logo um pressuposto de validade do teste não foi atendido. ## Uma alternativa então é calcular o p-valor através de simulação ou o teste exato de Fisher. Note que quando simulamos o p-valor não é necessário usar a correção de continuidade de Yates. Q<-chisq.test(dados,sim=500) Q Q$observed ### frequência observada Q$expected ### frequência esperada ##Há evidências de se rejeitar H0 # Comandos para obtenção da diferença entre proporções e seu IC(95%) ## Calculando as proporções entre Cryptosporidium negativo e Cryptosporidium positivo p11<-(dados[1,1]/(sum(dados[1,]))) p22<-(dados[1,1]/(sum(dados[1,]))) d<-p11-p21 # diferença entre as proporções vd<-((p11*(1-p11))/(sum(dados[1,])-1)) + ((p21*(1-p21))/(sum(dados[2,])-1)) ## Estimativa para a variância dvd<-sqrt(vd) ## raíz quadrada da variância z<-qnorm(0.975) #percentil da Normal padrão li<- d - (z*dvd) # Limite inferior ls<- d + (z*dvd) # Limite superior cbind(d,li,ls) # Intervalo de Confiança de 95%. Como o valor zero não está contido no IC a diferença é significativa ao nível de 95% de confiança. ##Razão de Chances ou Odds Ratio (OR) e IC95%(OR) OR<-(dados[1,1]*dados[2,2])/(dados[1,2]*dados[2,1]) ## Calculando a /*odds ratio (n11*n22/n12*n21) */## Quando OR=1 indica chances iguais. Se for OR>1, o grupo 1 apresenta maior chance que o grupo 2. ## Para o cálculo do IC para a OR, usamos o logaritmo da OR na base /e./ vf<-(1/dados[1,1])+(1/dados[1,2])+(1/dados[2,1]+(1/dados[2,2])) ##Estimativa para variância dpf<-sqrt(vf) ## raíz quadrada da variância z<-qnorm(0.975) #Percentil da Normal padrão liOR<-exp(log(OR)-z*dpf) #Limite inferior lsOR<-exp(log(OR)+z*dpf) # Limite Superior cbind(OR,liOR,lsOR) ## A chance de não haver Cryptosporidium e Giardia é 6,8 vezes maior que a presença podendo variar entre 3 e 15,4 vezes ao nível de confiança de 95%. -- Atenciosamente Felipe E. Barletta Mendes Estatístico - Conre3 9766-A +55 (41)-92077191 +55 (41)-33287216 -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20151104/ae036f7d/atta…> ------------------------------ Message: 3 Date: Wed, 4 Nov 2015 18:17:27 -0300 From: Marcos Vital <marcosvital(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] quem é o autor de tanks.R ??? Message-ID: <CAKLUGQ_d+SyHWigkh=byZMer0+=jiDdYyJt1WQB8Aczc1dLncQ(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Olá, Rodrigo. Fiz uma postagem falando da sua função, e mencionando alguns outros joguinhos no R. Está aqui, caso queira dar uma olhada: https://cantinhodor.wordpress.com/2015/10/25/tiros-explosoes-e-destruicao-a… Abraços! -- Marcos Vinícius Carneiro Vital Universidade Federal de Alagoas Instituto de Ciências Biológicas e da Saúde Setor de Biodiversidade e Ecologia -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20151104/62339867/atta…> ------------------------------ Message: 4 Date: Wed, 4 Nov 2015 18:20:38 -0300 From: Marcos Vital <marcosvital(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Campanha de crowdfunding envolvendo o R e Bioestatística Message-ID: <CAKLUGQ9W5A2kT0jqZexPAu4u=nOTpuBgHknAVqUhLDjarnubEQ(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Oi, pessoal Estou passando só para divulgar nossa campanha de financiamento uma última vez, pois estamos a uma semana do final. Relembrando: estamos usando uma campanha de financiamento coletivo para equipar o nosso laboratório e produzir material didático sobre diversos assuntos, incluindo um bocado de Bioestatística e suas aplicações no R. A campanha está aqui, para quem se interessar: http://www.kickante.com.br/campanhas/ciencia-livre-ecologia-e-bioestatistic… Abraços! -- Marcos Vinícius Carneiro Vital Universidade Federal de Alagoas Instituto de Ciências Biológicas e da Saúde Setor de Biodiversidade e Ecologia -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20151104/1d033226/atta…> ------------------------------ Message: 5 Date: Thu, 5 Nov 2015 10:38:13 -0300 From: Augusto Ribas <ribas.aca(a)gmail.com> To: R-br Lista <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] OFF TOPIC - Criar função no emacs para trocar # por ### Message-ID: <CACMkfRwN6HdTp_912PKHpLWYzJjLCaL+zAxkffuu9y26VTWGHQ(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Ola a todos. Bem eu tenho usado o emacs como ide para o R, mas sempre tive alguns problemas com identação. Pelo que entendo, no emacs, ele le comentário como três #, tipo ###. Mais ou menos esse problema, mas com hashtag http://stackoverflow.com/questions/26312317/wrong-indentation-of-comments-i… Assim quando eu ponho o comentário como um #, se dou um enter logo a frente desse #, o emacs identa e joga la pra frente. Pesquisando no google, eu achei gente reclamando da mesma coisa, e a solução que vi foi desligar a identação automática, o que é ruim, porque é legal o código identado bonitinho. Dai eu pensei, se teria como fazer uma função e deixar no .emacs , para que o emacs pega-se toda hora que eu digita-se #, ele transformasse em ###, porque ### funciona perfeito, mas é chato digitar três ###. Ai estou lendo esse tutorial https://www.masteringemacs.org/article/mastering-key-bindings-emacs E vendo que da para alterar o key-map do teclado, se entendi bem, e trocar teclas para outras funções, mas não estou conseguindo implementar isso pra deixar no .emacs para quando iniciar o emacs, toda hora que eu digitar #, na tela vira três ### Será que alguém poderia dar uma luz aqui? -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.com.br/ <http://augustoribas.heliohost.org> Github: https://github.com/Squiercg Lattes: http://lattes.cnpq.br/7355685961127056 -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20151105/2c14d559/atta…> ------------------------------ Subject: Legenda do Digest _______________________________________________ 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 59, assunto 5 ****************************************

2 1

OFF TOPIC - Criar função no emacs para trocar # por ###
by Augusto Ribas 05 Nov '15

05 Nov '15

Ola a todos. Bem eu tenho usado o emacs como ide para o R, mas sempre tive alguns problemas com identação. Pelo que entendo, no emacs, ele le comentário como três #, tipo ###. Mais ou menos esse problema, mas com hashtag http://stackoverflow.com/questions/26312317/wrong-indentation-of-comments-i… Assim quando eu ponho o comentário como um #, se dou um enter logo a frente desse #, o emacs identa e joga la pra frente. Pesquisando no google, eu achei gente reclamando da mesma coisa, e a solução que vi foi desligar a identação automática, o que é ruim, porque é legal o código identado bonitinho. Dai eu pensei, se teria como fazer uma função e deixar no .emacs , para que o emacs pega-se toda hora que eu digita-se #, ele transformasse em ###, porque ### funciona perfeito, mas é chato digitar três ###. Ai estou lendo esse tutorial https://www.masteringemacs.org/article/mastering-key-bindings-emacs E vendo que da para alterar o key-map do teclado, se entendi bem, e trocar teclas para outras funções, mas não estou conseguindo implementar isso pra deixar no .emacs para quando iniciar o emacs, toda hora que eu digitar #, na tela vira três ### Será que alguém poderia dar uma luz aqui? -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.com.br/ <http://augustoribas.heliohost.org> Github: https://github.com/Squiercg Lattes: http://lattes.cnpq.br/7355685961127056

1 0

Re: [R-br] Campanha de crowdfunding envolvendo o R e Bioestatística
by Marcos Vital 04 Nov '15

04 Nov '15

Oi, pessoal Estou passando só para divulgar nossa campanha de financiamento uma última vez, pois estamos a uma semana do final. Relembrando: estamos usando uma campanha de financiamento coletivo para equipar o nosso laboratório e produzir material didático sobre diversos assuntos, incluindo um bocado de Bioestatística e suas aplicações no R. A campanha está aqui, para quem se interessar: http://www.kickante.com.br/campanhas/ciencia-livre-ecologia-e-bioestatistic… Abraços! -- Marcos Vinícius Carneiro Vital Universidade Federal de Alagoas Instituto de Ciências Biológicas e da Saúde Setor de Biodiversidade e Ecologia

1 0

Re: [R-br] quem é o autor de tanks.R ???
by Marcos Vital 04 Nov '15

04 Nov '15

Olá, Rodrigo. Fiz uma postagem falando da sua função, e mencionando alguns outros joguinhos no R. Está aqui, caso queira dar uma olhada: https://cantinhodor.wordpress.com/2015/10/25/tiros-explosoes-e-destruicao-a… Abraços! -- Marcos Vinícius Carneiro Vital Universidade Federal de Alagoas Instituto de Ciências Biológicas e da Saúde Setor de Biodiversidade e Ecologia

1 0

Chi-quadrado demonstra associação, mas e daí?
by André Lucas de Oliveira Moreira 03 Nov '15

03 Nov '15

Prezados membros da lista R-br, Solicito a ajuda de vocês para resolver o seguinte problema: Tenho duas espécies de parasitas e alguns fatores que podem estar associados a incidência deles, como o período climático e a própria presença de um dos parasitas. Então executei o qui-quadrado para os seguintes dados: Giardia negativa Giardia positiva Cryptosporidium negativo 79,9% (250/313) 4,8% (15/313) Cryptosporidium positivo 10,9% (34/313) 4,5% (14/313) O resultado do chi-quadrado é siginificativo, mas como interpretar/extrair mais informações além do fato de os fatores estarem associados? O teste fala que há associação, mas apenas 4,5% das amostras continham os dois parasitas ao mesmo tempo. Pesquisei bastante sobre isso e o que achei foi algo muito sucinto sobre odds ratio, mas não entendi, nem saberia fazer no R. Agradecido, André -- *MSc. André Lucas de O. Moreira* http://lattes.cnpq.br/7258065668864153 <http://www.wikiaves.com.br/perfil_andrelukinhas> http://www.wikiaves.com.br/perfil_andrelukinhas 79 8837-3562 79 9132-9093

4 4

(sem assunto)
by Lester Cavero 03 Nov '15

03 Nov '15

1 0

Fwd: Gerar tabela pelo comando latex ( tabular ) a partir de 2 colunas de fatores pertencentes a um dataframe
by André Figueiras Rutz 03 Nov '15

03 Nov '15

Prezados, Vi a sugestão de gerar tabelas no markdown a partir do comando função Hmisc::latex() e li a documentação do pacote. Não localizei muitos exemplos para adaptar, apesar da documentação extensa ainda não consegui gerar a tabela da maneira que desejo que é mais ou menos assim: Segue abaixo CRM do que eu consegui fazer e os arquivos dos dados nos links --- title: "tabela" output: pdf_document --- ```{r, include=F} #--------------------------------------------------------- # Preparo da sessao pkg <- c( "Hmisc", "reshape", "plyr", "tables", "gmodels", "xtable") sapply(pkg, require, character.only = TRUE) #------------------------------------------------------------ # lendo os dados da <- read.csv("D:/data.csv", sep =";") #------------------------------------------------------------ # reordenando fatores da$escopai <- factor(da$escopai, levels = c("Primário", "Médio", "Graduação", "Pós-Graduação", "Mestrado", "Doutorado")) da$escomae <- factor(da$escomae, levels = c("Primário", "Médio", "Graduação", "Pós-Graduação", "Mestrado", "Doutorado")) ``` Tabela de escolaridade dos genitores ```{r, echo=FALSE, results='asis'} latex (tabular ( (Factor(escopai, "Níveis de escolaridade")) ~ (Factor(escopai, "Pai")) + (Factor(escomae, "Mãe") ), data=da )) ``` _________________________________________________ Gostaria também saber como colocar aquelas famosas notinhas de rodapé em letrinhas miúdas no final da tabela e se alguem conseguir um modelo de outra tabela que inclua uma coluna de porcentagem para cada coluna de dados links: https://dl.dropboxusercontent.com/u/64332051/data.csv https://dl.dropboxusercontent.com/u/64332051/tabela.Rmd Saudações cordiais, André Rutz

2 2

Gerar tabela pelo comando latex ( tabular ) a partir de 2 colunas de fatores de um dataframe
by André Figueiras Rutz 02 Nov '15

02 Nov '15

Prezados, Vi a sugestão de gerar tabelas no markdown a partir do comando função Hmisc::latex() e li a documentação do pacote. Não localizei muitos exemplos para adaptar, apesar da documentação extensa ainda não consegui gerar a tabela da maneira que desejo que é mais ou menos assim: [image: Inline image 1] Segue abaixo CRM do que eu consegui fazer e os arquivos dos dados nos links --- title: "tabela" output: pdf_document --- ```{r, include=F} #--------------------------------------------------------- # Preparo da sessao pkg <- c( "Hmisc", "reshape", "plyr", "tables", "gmodels", "xtable") sapply(pkg, require, character.only = TRUE) #------------------------------------------------------------ # lendo os dados da <- read.csv("D:/data.csv", sep =";") #------------------------------------------------------------ # reordenando fatores da$escopai <- factor(da$escopai, levels = c("Primário", "Médio", "Graduação", "Pós-Graduação", "Mestrado", "Doutorado")) da$escomae <- factor(da$escomae, levels = c("Primário", "Médio", "Graduação", "Pós-Graduação", "Mestrado", "Doutorado")) ``` Tabela de escolaridade dos genitores ```{r, echo=FALSE, results='asis'} latex (tabular ( (Factor(escopai, "Níveis de escolaridade")) ~ (Factor(escopai, "Pai")) + (Factor(escomae, "Mãe") ), data=da )) ``` _________________________________________________ Gostaria também saber como colocar aquelas famosas notinhas de rodapé em letrinhas miúdas no final da tabela e se alguem conseguir um modelo de outra tabela que inclua uma coluna de porcentagem para cada coluna de dados links: https://dl.dropboxusercontent.com/u/64332051/data.csv https://dl.dropboxusercontent.com/u/64332051/tabela.Rmd Saudações cordiais, André Rutz

1 0

Meetup R São Paulo Users Group (IME-USP) - 6 de novembro do 2015
by Luis Benites Sánchez 02 Nov '15

02 Nov '15

Caros, Neste mês de novembro teremos dois Meetups (6 e 27) de usuários de R de São Paulo, no auditório Jacy Monteiro no Bloco B. Segue a programação do mês de novembro. Programação da sexta *6 de novembro* (14h-16:30h) - [14:00 - 14:20] Como publicar seu livro em uma editora internacional - Robinson dos Santos (Springer) - [Presencial] - [14:20 - 15:10] Misturas Finitas de Distribuições Usando o R - Prof. Celso Rômulo (UFAM) - [Presencial] - [15:10 - 15:25] Coffee Break - [15:25 - 16:15] Regressão não linear com o R - Prof. Walmes Zeviani (LEG/DEST/UFPR) - [Online] - [16:15 - 16:25] Sorteio de livros As palestras serão transmitidas ao vivo pela conexão IPTV da USP e YouTube (https://plus.google.com/events/cc1iefrbck0f9m1offtub0fd8p4) . Para poder participar, pode confirmar sua participação gratuitamente no site do Meetup :) Para mais informações: Site: www.meetup.com/pt/useR-SP Canal de YouTube: www.youtube.com/channel/UCb9PdJYYaElE6Uqb1KxhE8w e-mail: meetuprsp(a)gmail.com Um abraço, Luis Benites Sánchez Rocío Paola Maehara Organizadores do São Paulo R Users Group

1 0

Modelo Regressão Beta Inflacionado
by Arlene Guimaraes Leite 31 Out '15

31 Out '15

Prezados, Ospina sugere a construção de bandas de confiança (envelopes) simuladas com o seguinte comando: #==== normal probability plots ==== envelope.BIc(fit, type="quantil", main="quantal residual") mas não estou conseguindo executá-lo no R. Para executar o modelo utilizei as seguintes bibliotecas: library(gamlss) library(gamlss.dist) library(gamlss.nl) library(gamlss.mx) library(gamlss.cens) library(gamlss.tr) library(foreign) library(sas7bdat) Alguém pode me ajudar? Att., Arlene

2 1

Erro em Acesso ao Servidor - Pacote RCurl
by Paulo Nogueira Starzynski 29 Out '15

29 Out '15

Pessoal, estou desenvolvendo um programa que realiza login em um servidor e depois faz uma requisição POST. estava tudo funcionando bem durante o desenvolvimento até que começou a retornar "Internal Server Error: 500" ao fazer o POST. Porém, quando eu acesso via browser, tudo funciona direito. O fato de funcionar no browser indica necessariamente que o erro é no meu código, certo? devem faltar opções de controle no objeto curl... Meu código está abaixo, mas infelizmente não é reproduzĩvel porque o acesso ao servido não é público. Mesmo assim, se alguém com familiaridade em curl puder olhar e dar alguma idéia... testei configurar várias opções do curl e nada contronou esse erro 500. ----------- require(RCurl) #------------------------------------------ curl <- getCurlHandle() curlSetOpt( verbose = T, cookiejar = 'cookies.txt' , followlocation = TRUE , ssl.verifypeer = FALSE, useragent = 'Mozilla/5.0 (Windows NT 5.1) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/44.0.2403.157 Safari/537.36', curl = curl ) #------------------------------------------ # list parameters to pass to the website (pulled from the source html) your.username <- 'usuario' your.password <- 'senha' params <- list( 'txt_email' = your.username, 'txt_senha' = your.password) # Realiza Login no Servidor webAddr <- "http://hostname.br/dirname/ct_login.php" login <- postForm(webAddr, .params = params, curl = curl, style = 'POST') login [1] "1{Sucesso controller{Login efetuado" attr(,"Content-Type") charset "text/html" "ISO-8859-1" #------------------------------------------ # Realiza Busca webBuscaId <- "http://hostname.br/dirname/executa_busca.php" busca <- postForm(webBuscaId, .opts = list(httpheader = c('Content-Type' = 'application/json', ssl.verifypeer = FALSE)), .params = list('idConsulta' = '13678616F', 'adicPar' = 'formatoNulo'), curl = curl, style="POST") * Hostname was found in DNS cache * Trying 10.65.0.169... * Connected to xxxxxx.xxxxxx.xxx.br (10.65.0.169) port 80 (#3) > POST /xxxxx/executa_busca.php HTTP/1.1 User-Agent: Mozilla/5.0 (Windows NT 5.1) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/44.0.2403.157 Safari/537.36 Host: xxxxx.xxxx.xxx.br Accept: */* Cookie: PHPSESSID=3dug68beqvmgpe0v19cc78a947 Content-Type: application/json ssl.verifypeer: FALSE Content-Length: 28 * upload completely sent off: 28 out of 28 bytes * HTTP 1.0, assume close after body < HTTP/1.0 500 Internal Server Error < Date: Thu, 29 Oct 2015 12:53:27 GMT < Server: Apache/2.2.15 (Red Hat) < X-Powered-By: PHP/5.3.3 < Expires: Thu, 19 Nov 1981 08:52:00 GMT < Cache-Control: no-store, no-cache, must-revalidate, post-check=0, pre-check=0 < Pragma: no-cache < Content-Length: 0 < Connection: close < Content-Type: text/html; charset=ISO-8859-1 < * Closing connection 3 Error: Internal Server Error > Abraços, Paulo

1 1