- R-br - DINF UFPR

Teste de Aderência
by ana paula coelho madeira 22 Ago '16

22 Ago '16

Boa tarde. Tenho um conjunto de dados nos quais estou trabalhando com diferentes fdp's. Gostaria de aplicar diferentes testes de aderência, entre eles o Anderson Darling. No entanto, consegui apenas o ks.test. No pacote Nortest encontrei o ad.test para testar a normalidade, que não é o caso. Alguém tem algum material para me indicar? Desde já agradeço. Ana Paula

2 1

Teste T-pareado
by Simone D. Sartorio 22 Ago '16

22 Ago '16

Pessoal, boa noite! Surgiu uma dúvida conceitual, que não estou conseguindo resolver. Estava preparando minhas aulas e ao conferir um exercício no R que inicialmente tinha resolvido com esses comandos (Solução 1) observei que pelo p-value aceitaria a H0: não há diferença entre as médias. Contudo, pela minha resolução manual, deveria rejeitar H0. Logo, propus a solução 2, exatamente como realizei manualmente e dai, os resultados foram iguais. Gostaria de saber, qual é a diferença da solução 1 para a solução 2? Para mim, a sintaxe de ambas soluções deveria retornar a mesma conclusão. Problema:Uma distribuidora de combustíveis deseja verificar se um novo tipo de gasolina é eficaz na revitalização de motores velhos.Selecionou-se 12 automóveis de um mesmo modelo com mais de 8 anos de uso e, após regulagem dos motores, verifica-se o consumo de combustível. Em seguida, o carro é abastecido com o novo tipo de combustível durante 15 semanas e uma nova aferição é feita (km/litro). Como o desempenho dos automóveis foi medido antes e depois das 15 semanas, é razoável assumir que exista alguma dependência entre as variáveis. Gostaria de testar de realizar o teste de hipótese ao nível de 5%:H0: mD = 0 (O novo combustível não aumenta o rendimento) Ha: mD > 0 (O novo combustível aumenta o rendimento), onde mD é a verdade média da diferença do rendimento apos e antes. Aqueles que puderem me ajudar a entender, agradeço desde já!!!att.Simone ### Solução 1 consumo<- c(8.1, 7.9, 6.8, 7.8, 7.6, 7.9, 5.7, 8.4, 8.0, 9.5, 8.0, 6.8, 11.6, 8.8, 9.9, 9.5, 11.6, 9.1, 10.6, 10.8, 13.4, 10.6, 10.5, 11.4)momento<- factor(c(rep(‘Antes’,12), rep(‘Apos’,12))); momento t.test(consumo ~ momento, paired = T, conf.level = 0.95, alternative=‘greater’) # Paired t-test##data: consumo by momento#t = -6.5396, df = 11, p-value = 1#alternative hypothesis: true difference in means is greater than 0#95 percent confidence interval:# -3.7495 Inf#sample estimates:#mean of the differences # -2.941667 ### Solução 2X <- c(8.1, 7.9, 6.8, 7.8, 7.6, 7.9, 5.7, 8.4, 8.0, 9.5, 8.0, 6.8)# antesY <- c(11.6, 8.8, 9.9, 9.5, 11.6, 9.1, 10.6, 10.8, 13.4, 10.6, 10.5, 11.4)# aposD <- Y-X # Diferença > t.test(D, paired = F, conf.level = 0.95, alternative='greater')## One Sample t-test##data: D#t = 6.5396, df = 11, p-value = 2.097e-05#alternative hypothesis: true mean is greater than 0#95 percent confidence interval:# 2.133833 Inf#sample estimates:#mean of x # 2.941667 ***---***---***---***---***---***---***---***---***---***---***---***---**** Simone Daniela Sartorio* Professora Adjunta II no CCA/UFSCar, Araras/SP.* Doutora e Mestre em Estatística e Experimentação Agronômica - ESALQ/USP.* Licenciada em Matemática - UNESP/Rio Claro.* Contatos: * tel: (19)3543-2664* cel: (19)98182-0586 - Tim***---***---***---***---***---***---***---***---***---***---***---***---*** Tenha um bom dia! ;)

4 3

Re: [R-br] Teste T-pareado
by Ivan Bezerra Allaman 22 Ago '16

22 Ago '16

A resposta é simples! Você está testando a hipótese erra na solução 1. Perceba no output, que a diferença está negativa, ou seja, o teste está fazendo ANTES - DEPOIS. Logo, se nós esperamos que DEPOIS melhore o desempenho, fazendo mais km/litro, então a DIFERENÇA é negativa e a hipótese a ser testada é: H0: md >= 0 Ha: md < 0 That is it!

2 1

erro no loop
by Edimeire Alexandra Pinto 22 Ago '16

22 Ago '16

Oi Pessoal. Como não achei um pacote que trabalhasse com previsões dinâmicas,estou eu mesma tendo de montar isso. Já dei alguns passos, agora resta entender“for” de quando coloco for(j in 1:20), tendo em vista que o R diz que os índicesestão fora do limite, sendo que na realidade dava para realizar as previsões.Além disso, o tamanho das previsões, em predict, era para ser do tamanho doloop de 20+j+1 até 84 (total de observações) e isso não ocorre... Suponha que eu tenha a base de dados Canada (é só copiar ecolar os comandos), meu objetivo rodarvários modelos dynlm para cada amostra que começa em 1 e vai de 20+j. Asprevisões de cada modelo rodado devem ser de 20+j+1 até 76, pois a janela ésempre fixa de tamanho 8.Gente, qualquer ajuda é muito bem vinda. require(vars)require(dynlm)data(Canada)colnames(Canada)summary(Canada)for(j in 1:20){ deflator_prev<-predict(dynlm(e-prod~ I(L(e,1)-rw)+ I(U-rw) +L(U,1),data=ts(Canada[(1:20+j),],start=c(1980,01),freq=4)), newdata=ts(Canada[(20+j+1:20+j+8),],start=c(1985,02),freq=4 )) } deflator_prev

1 0

Gerar series sintéticas de precipitação no R
by Alisson Lucrécio 22 Ago '16

22 Ago '16

Caro Colegas, Bom dia. Algum colega do tem um script para geração de séries sintéticas de precipitação no R. Baseado nessa metodologia. Segue o link. https://www.google.com.br/url?sa=t&rct=j&q=&esrc=s&source= web&cd=1&cad=rja&uact=8&ved=0ahUKEwjTy-XKgtXOAhWIgZAKHWO3BA8QFggcMAA& url=ftp%3A%2F%2Fftp.ufv.br%2Fdea%2FGPRH%2Fteses%2Fds_ oliveira%2FVicente%2520de%2520Paulo%2520Santos%2520de% 2520Oliveira%2520-%2520Tese%2520DS.doc&usg=AFQjCNEWA- 2QE2NG05QQWr0sUUbG5Csw4g&bvm=bv.129759880,d.Y2I Muito obrigado. -- Alisson Lucrecio da Costa

1 0

Geração de series sintéticas de precipitação no R
by Alisson Lucrécio 22 Ago '16

22 Ago '16

Caro Colegas, Bom dia. Algum colega do tem um script para geração de séries sintéticas de precipitação no R. Baseado nessa metodologia. Segue o link. https://www.google.com.br/url?sa=t&rct=j&q=&esrc=s&source=web&cd=1&cad=rja&… Muito obrigado. -- Alisson Lucrecio da Costa

1 0

previsão dinâmica
by Edimeire Alexandra Pinto 19 Ago '16

19 Ago '16

OI gente. Preciso de uma ajudinha. Sei que parece bobo, mas nãoconsigo estender a janela de previsões sequenciais com largura de , com largura fixa de tamanho 8.Criar um modelo com sample de 1:70 e prever 71 até 78,depois reestimar a equação com sample de 1:72 e previsão de 73:80, reestimar omodelo de 1:73 e previsão de 74:83 eassim por adiante até fechar o range em 84, última observação.Vejamos um data frame já existente no R como nome de Canadaque são séries que vão do primeiro trimestre de 1980 ao último trimestre de2004:require("dynlm") library("vars") data("Canada") summary(Canada) class(Canada)colnames(Canada) suponha que eu definaque minha amostra para estimar o modelo com sequencia em j que vai até 10,assim teríamos um modelo estimado de 01/1980 a 3/1996, depois outros modelos estimado de 01/1980a 4/1996, depoismais outro modelo estimado de 01/1980 a 1/1997, assim por adiante até 4/1998, ou seja, paraem 4/1998. No R seria mais ou menos assim:eq<-list() for(j in1:10){ eq(i)<-dynlm(e~L(e,1)+prod+rw+u,data=window(Canada,start=c(1980,1),end=c(1996,2),freq=4,extend=j)) o problema é que a previsão tem de ser sequencial tanto nostart quanto no end e com tamanho de fixo de 8, e também com uma janela que vaiindo sequencialmente em j de 1 a 18. Por exemplo, a previsão teria de serassim: a primeira de 04/1996 a 4/1998 (observe que o primeiro modelo termina em3/1996, conforme falei antes), depois 01/1997 a 2/1999 (observe que o primeiromodelo termina em 1/1997, conforme falei antes), depois faria outra previsão de2/1997 a 4/1999 e assim por adiante até chegar na última observação que é04/2000. Note que aumentamos sempre de 8+j em que j vai de 1 a 18. No R, não consigo desenvolver, pois seria algo mais ou menosassim: ## sei que está errado, mas é só para ter ideiafor( i in1:18){ predict(eq(i), newdata=window(Canada,start=c(1996,4+i),end=c(1998,4+i+8),freq=4,extend=i) ## 8 é o tamanho fixo da janela de previsão Se tiver ficado confuso, desculpe-me, mas quis apenasmostrar que quero um previsão dinâmica e não estática, tendo e, vista que aequação que sempre estimada à medida que o tamanho da amostra vai variando.Qualquer dica é válida!Obrigada, gente.

3 3

Re: [R-br] Digest R-br, volume 68, assunto 17
by ana paula coelho madeira 19 Ago '16

19 Ago '16

Prezado Alexandro, Os comandos em questão tentam apenas a normalidade. Tenho interesse em ver ficar a aderência de uma distribuição qualquer. Agradeço > From: r-br-request(a)listas.c3sl.ufpr.br > Subject: Digest R-br, volume 68, assunto 17 > To: r-br(a)listas.c3sl.ufpr.br > Date: Wed, 17 Aug 2016 12:00:02 -0300 > > 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. previsão dinâmica (Edimeire Alexandra Pinto) > 2. Re: Anderson Darling (Alexandro (Yahoo)) > 3. OFF TOPIC - Banco de dados (Andre Oliveira) > 4. Re: OFF TOPIC - Banco de dados (luc.souza(a)usp.br) > 5. RES: OFF TOPIC - Banco de dados (Leonard Mendonça de Assis) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Tue, 16 Aug 2016 21:53:38 +0000 (UTC) > From: Edimeire Alexandra Pinto <economatistica(a)yahoo.com.br> > To: a lista Brasileira oficial de discussão do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] previsão dinâmica > Message-ID: > <228175548.15741202.1471384418129.JavaMail.yahoo(a)mail.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > OI gente. > > Preciso de uma ajudinha. Sei que parece bobo, mas nãoconsigo estender a janela de previsões sequenciais com largura de , com largura fixa de tamanho 8.Criar um modelo com sample de 1:70 e prever 71 até 78,depois reestimar a equação com sample de 1:72 e previsão de 73:80, reestimar omodelo de 1:73 e previsão de 74:83 eassim por adiante até fechar o range em 84, última observação.Vejamos um data frame já existente no R como nome de Canadaque são séries que vão do primeiro trimestre de 1980 ao último trimestre de2004:require("dynlm") > > library("vars") > > data("Canada") > > summary(Canada) > > class(Canada)colnames(Canada) > > > > suponha que eu definaque minha amostra para estimar o modelo com sequencia em j que vai até 10,assim teríamos um modelo estimado de 01/1980 a 3/1996, depois outros modelos estimado de 01/1980a 4/1996, depoismais outro modelo estimado de 01/1980 a 1/1997, assim por adiante até 4/1998, ou seja, paraem 4/1998. No R seria mais ou menos assim:eq<-list() > > for(j in1:10){ > > eq(i)<-dynlm(e~L(e,1)+prod+rw+u,data=window(Canada,start=c(1980,1),end=c(1996,2),freq=4,extend=j)) > > o problema é que a previsão tem de ser sequencial tanto nostart quanto no end e com tamanho de fixo de 8, e também com uma janela que vaiindo sequencialmente em j de 1 a 18. Por exemplo, a previsão teria de serassim: a primeira de 04/1996 a 4/1998 (observe que o primeiro modelo termina em3/1996, conforme falei antes), depois 01/1997 a 2/1999 (observe que o primeiromodelo termina em 1/1997, conforme falei antes), depois faria outra previsão de2/1997 a 4/1999 e assim por adiante até chegar na última observação que é04/2000. Note que aumentamos sempre de 8+j em que j vai de 1 a 18. No R, não consigo desenvolver, pois seria algo mais ou menosassim: > > ## sei que está errado, mas é só para ter ideiafor( i in1:18){ > > predict(eq(i), newdata=window(Canada,start=c(1996,4+i),end=c(1998,4+i+8),freq=4,extend=i) ## 8 é o tamanho fixo da janela de previsão > > Se tiver ficado confuso, desculpe-me, mas quis apenasmostrar que quero um previsão dinâmica e não estática, tendo e, vista que aequação que sempre estimada à medida que o tamanho da amostra vai variando.Qualquer dica é válida!Obrigada, gente. > > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160816/39b43423/atta…> > > ------------------------------ > > Message: 2 > Date: Tue, 16 Aug 2016 22:04:57 +0000 (UTC) > From: "Alexandro (Yahoo)" <vl.alexandro(a)yahoo.com.br> > To: ana paula coelho madeira <apcmadeira(a)hotmail.com>, a lista > Brasileira oficial de discussão do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Anderson Darling > Message-ID: > <1102947158.17317609.1471385097903.JavaMail.yahoo(a)mail.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > #Segue um exemplo!require(ADGofTest)#?ad.test > set.seed( 123 )x <- runif( 100 ) > ad.test( x )$p.valuead.test( x, pnorm, 0, 1 )$p.value #Espero ter ajudado! > Alexandro Vieira Lopes > > > > De: ana paula coelho madeira via R-br <r-br(a)listas.c3sl.ufpr.br> > Para: lista R <r-br(a)listas.c3sl.ufpr.br> > Enviadas: Segunda-feira, 15 de Agosto de 2016 15:28 > Assunto: [R-br] Anderson Darling > > Boa tarde. > Tenho um conjunto de dados nos quais estou trabalhando com diferentes fdp's. Gostaria de aplicar diferentes testes de aderência, entre eles o Anderson Darling. No entanto, consegui apenas o ks.test. No pacote Nortest encontrei o ad.test para testar a normalidade, que não é o caso. Alguém tem algum material para me indicar? > Desde já agradeço. > Ana Paula > _______________________________________________ > 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� c�igo m�imo reproduz�el. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160816/268f052c/atta…> > > ------------------------------ > > Message: 3 > Date: Wed, 17 Aug 2016 03:21:29 +0000 (UTC) > From: Andre Oliveira <andreolsouza(a)yahoo.com.br> > To: a lista Brasileira oficial de discussão do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] OFF TOPIC - Banco de dados > Message-ID: > <811409596.15775329.1471404089455.JavaMail.yahoo(a)mail.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > Boa noite,tenho um projeto com aluno de ensino médio sobre produzir interface em shiny para visualizar dados. É uma feira onde o tema é alimentação e pensei em explorar banco de dados de produção de milho, soja ou outro produtos por hectare ao longo de anos no Brasil, via shiny. Alguém saberia dizer onde posso encontrar banco de dados disponíveis com varias culturas e vários anos? Este aluno é do técnico em informática e vem trabalhando com estatística onde há uma integração entre as áreas. > 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 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160817/1a282242/atta…> > > ------------------------------ > > Message: 4 > Date: Wed, 17 Aug 2016 06:47:50 -0300 (BRT) > From: luc.souza(a)usp.br > To: Andre Oliveira <andreolsouza(a)yahoo.com.br>, "a lista Brasileira > oficial de discussão do programa R." <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] OFF TOPIC - Banco de dados > Message-ID: <3E69A4AD-23EB-49AC-B77D-2914AD7D210C(a)usp.br> > Content-Type: text/plain; charset="utf-8" > > Na conab tem uma série histórica de vários produtos e no ibge Tb. > > Enviado do meu iPhone > > > Em 17 de ago de 2016, às 00:21, Andre Oliveira via R-br <r-br(a)listas.c3sl.ufpr.br> escreveu: > > > > Boa noite, > > tenho um projeto com aluno de ensino médio sobre produzir interface em shiny para visualizar dados. É uma feira onde o tema é alimentação e pensei em explorar banco de dados de produção de milho, soja ou outro produtos por hectare ao longo de anos no Brasil, via shiny. Alguém saberia dizer onde posso encontrar banco de dados disponíveis com varias culturas e vários anos? Este aluno é do técnico em informática e vem trabalhando com estatística onde há uma integração entre as áreas. > > 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 > > _______________________________________________ > > 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/20160817/6a558ec2/atta…> > > ------------------------------ > > Message: 5 > Date: Wed, 17 Aug 2016 09:37:57 -0300 > From: Leonard Mendonça de Assis <assis.leonard(a)gmail.com> > To: 'a lista Brasileira oficial de discussão do programa R.' > <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] RES: OFF TOPIC - Banco de dados > Message-ID: <00ad01d1f884$30927090$91b751b0$(a)gmail.com> > Content-Type: text/plain; charset="utf-8" > > André, eu sei que o CEPEA/USP disponibiliza séries temporais com a cotação de diversos produtos agrícolas. Área plantada eu desconheço algum órgão que possa ter isso. Uma tentativa seria ministério da Agricultura ou censo agropecuário do ibge > > > > De: R-br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de Andre Oliveira via R-br > Enviada em: quarta-feira, 17 de agosto de 2016 00:21 > Para: a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> > Assunto: [R-br] OFF TOPIC - Banco de dados > > > > Boa noite, > > tenho um projeto com aluno de ensino médio sobre produzir interface em shiny para visualizar dados. É uma feira onde o tema é alimentação e pensei em explorar banco de dados de produção de milho, soja ou outro produtos por hectare ao longo de anos no Brasil, via shiny. Alguém saberia dizer onde posso encontrar banco de dados disponíveis com varias culturas e vários anos? Este aluno é do técnico em informática e vem trabalhando com estatística onde há uma integração entre as áreas. > > > > 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 > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160817/bcb2784f/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 68, assunto 17 > *****************************************

2 1

Re: [R-br] O R não abre no linux mint
by Sergio Pinto 18 Ago '16

18 Ago '16

Prezado Marcelo Obrigado pela ajuda. Infelizmente o problema persiste. Continua aparecendo a mesma saída no terminal. Att Sergio Medeiros Pinto Doutorando em Ciências Médicas UERJ Prof. de Educação Fisica Em 17/08/2016 21:50, "Marcelo Laia" <marcelolaia(a)gmail.com> escreveu: http://stackoverflow.com/questions/20967256/error-on- package-and-r-paths-and-installations On 17/08/16 at 08:53pm, Sergio Pinto via R-br wrote: > Prezados > tenho um computador 32 bits com linux mint e instalei o pacote r-base. > O ícone do R foi criado mas não abre o programa. > No terminal aparece a seguinte mensagem: > > sergio@sergio-Infoway ~ $ R > /usr/local/bin/R: line 236: /usr/local/lib/R/etc/ldpaths: No such file or > directory > /usr/local/bin/R: line 259: /usr/local/lib/R/bin/exec/R: No such file or > directory > > > O meu sistema é > Linux sergio-Infoway 3.13.0-39-generic #66-Ubuntu SMP Tue Oct 28 13:31:23 > UTC 2014 i686 i686 i686 GNU/Linux > > Poderiam me ajudar? > > Desde já agradeço a atenção > > *Sergio Medeiros Pinto* > Doutorando em Ciências Médicas UERJ > Prof. de Educação Física > _______________________________________________ > 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. -- Marcelo

2 1

Girar rótulo da lattice
by Andre Oliveira 18 Ago '16

18 Ago '16

Olá,como posso girar os rótulos da lattice? par(las=1) não deu certo! require(lattice) par(las=1) xyplot(GP~factor(Ano)|S,pch=19, ylab="Ganho de peso (g)",xlab="", main="Ganho de peso por sexo") 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

1 0