- R-br - DINF UFPR

Inserir Equação e Ajuste de Dados
by Gustavo Dias Azevedo 03 Ago '12

03 Ago '12

Prezados, Estou trabalhando com ajuste de dados para uma equação e gostaria de saber se é possível, utilizando o R, saber se meus dados se ajustam a um modelo não-linear de equação e como eu faço isso. Obrigado,

1 0

Join my network on LinkedIn
by Julio Araujo via LinkedIn 03 Ago '12

03 Ago '12

LinkedIn ------------ Julio Araujo requested to add you as a connection on LinkedIn: ------------------------------------------ I'd like to add you to my professional network on LinkedIn. Accept invitation from Julio Araujo http://www.linkedin.com/e/-33g73l-h5ev5rfp-4y/vhj1xG42iBdUFK5d_HMEGL4iFMC51… View profile of Julio Araujo http://www.linkedin.com/e/-33g73l-h5ev5rfp-4y/rso/197235503/pw-P/name/82777… ------------------------------------------ You are receiving Invitation emails. This email was intended for Fábio Corrêa. Learn why this is included: http://www.linkedin.com/e/-33g73l-h5ev5rfp-4y/plh/http%3A%2F%2Fhelp%2Elinke… (c) 2012, LinkedIn Corporation. 2029 Stierlin Ct, Mountain View, CA 94043, USA.

1 0

Valores ajustados de um modelo ARIMA(0,2,1)
by Fatima do Nascimento Silva 02 Ago '12

02 Ago '12

Olá pessoal, Bom dia! Fiz um ajuste de um modelo ARIMA(0,2,1) para uma determinada série temporal. Agora estou precisando obter o vetor de valores que foram ajustados pelo modelo, para que eu possa criar um gráfico da série observada versus valores ajustados. O comando que utilizei para ajustar o modelo foi: (mdl=arima(x,order=c(0,2,1),include.mean=F)) Alguém saberia que comando utilizo no software R para obter o vetor de valores que foram ajustados pelo modelo? Obrigada. Atenciosamente, Fátima Nascimento Programa de Pós-Graduação em Ciência e Engª de Petróleo-PPGCEP/UFRN Fone: (84)9936-9725

4 5

Re: [R-br] barplot + lines com dois eixos Y
by Marcelo Claro de Souza 02 Ago '12

02 Ago '12

Pode crer..hehe Eu tb não me atentei a esse detalhe. Vou mandar um cc p a lista :-) Abraço Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil ________________________________ De: Augusto Ribas <ribas.aca(a)gmail.com> Para: Marcelo Claro de Souza <marcelo_claro(a)yahoo.com.br> Enviadas: Quarta-feira, 1 de Agosto de 2012 21:09 Assunto: Re: [R-br] barplot + lines com dois eixos Y Nossa agora que vi, a gente tava trocando e-mail pessoais. Tava respondendo e nem vi. Era bom deixa as resposta la na lista, pra aparecer no nabble, ai se alguem procurar via google ou sei la ja acha como fazer. Eu sempre acho muitas soluções assim :) Em 31 de julho de 2012 23:09, Marcelo Claro de Souza <marcelo_claro(a)yahoo.com.br> escreveu: > Olá Augusto. > Agradeço de coração todas as dicas que vc me deu. > Finalmente deu certo a figura que eu queria fazer. > > Abraço. > > Marcelo Claro de Souza > Biologist, PhD student in Plant Biology > Institute of Bioscience - UNESP, Brazil > > > ________________________________ > De: Augusto Ribas <ribas.aca(a)gmail.com> > Para: Marcelo Claro de Souza <marcelo_claro(a)yahoo.com.br> > Enviadas: Terça-feira, 31 de Julho de 2012 20:26 > > Assunto: Re: [R-br] barplot + lines com dois eixos Y > > x=c(seq(1:12)) > rainfall=c(17.54,7.68,9.87,5.45,0.15,1.61,0.01,0.33,1.53,4.85,3.77,9.35)##barras > tmax=c(30.21,30.71,27.23,28.23,25.19,24.39,26.59,28.39,29.77,28.49,28.61,29.79)##linhas > tmin=c(16.82,16.36,15.91,15.8,11.57,9.07,13.07,12.65,12.98,15.47,16.21,17.26)##linhas > tmed=c(23.51,23.53,21.57,22.02,18.38,16.73,19.83,20.52,21.38,21.98,22.41,23.53)##linhas > > > #o que vc vai ter que fazer é mudar as margens > #vc faz isso com um argumento dentro do par, o argumento mar > #basicamente é a distancia das coisas assim: c(bottom, left, top, right) > #os valores padrao são: c(5, 4, 4, 2) + 0.1., so subir o 2 pra mais > linhas pra caber o texto > #da ?par e da uma olhada na documentação > > par(mar=c(5, 4, 4, 5) + 0.1) > barras<-barplot(rainfall,axes=F,ylab="Nome do eixo esquerdo aqui") > axis(2) > axis(1,col="red") > > par(new=TRUE,mar=c(5, 4, 4, 5) + 0.1) > plot(tmed~barras,type="b",ylim=c(5,35),xlim=c(0.2,14.4),frame.plot=F,axes=F,ylab=NA) > #note que eu coloquei ylab=NA aqui em cima pq se vc escrever algo aqui > vai sobrepor > #com o que ja deixou de leganda no barplot() > points(tmin~barras,type="b",lty=2) > points(tmax~barras,type="b",lty=2) > axis(4) > #vc vai usar o comando mtext, que é pra escrever texto fora da area do plot. > #de ?mtext e veja as muitas possibilidades > #basicamente eu falei o lado no side, e a distancia do plot, line. > mtext(text="Nome do eixo direito aqui", side = 4,line = 3) > axis(1) > > #Mas lembre-se que existem mil jeitos de fazer as mesma coisa > #Talvez essa não seja a melhor forma, mas espero que possa resolver > seus problemas > > Em 31 de julho de 2012 08:50, Marcelo Claro de Souza > <marcelo_claro(a)yahoo.com.br> escreveu: >> Augusto, >> Muito obrigado. Vc me ajudou muito. >> Uma última dúvida. >> Para alterar a legenda do eixo Y do lado esquerdo eu uso o comando ylab, >> mas >> o mesmo não funciona para o segundo eixo Y. Qual comando devo usar para >> inserir a legenda do lado direito? >> Abraço. >> >> Marcelo Claro de Souza >> Biologist, PhD student in Plant Biology >> Institute of Bioscience - UNESP, Brazil >> >> >> ________________________________ >> De: Augusto Ribas <ribas.aca(a)gmail.com> >> Para: Marcelo Claro de Souza <marcelo_claro(a)yahoo.com.br> >> Enviadas: Segunda-feira, 30 de Julho de 2012 21:07 >> >> Assunto: Re: [R-br] barplot + lines com dois eixos Y >> >> So mais uma coisa, note que os eixos vão ficar bagunçados: >> >> barras<-barplot(rainfall,axes=F) >> axis(2) >> #vou colocar aqui em vermelho o eixo x do grafico de barras pra gente >> ver como fica em relação ao outro >> axis(1,col="red") >> par(new=TRUE) >> >> plot(tmed~barras,type="b",ylim=c(5,35),frame.plot=F,axes=F) >> points(tmin~barras,type="b",lty=2) >> points(tmax~barras,type="b",lty=2) >> axis(4) >> axis(1) >> >> A bolinha era pra ficar no centro da barra, mas não tava ficando. >> Isso por causa da escala de x, que do plot ficou diferente do barblot >> pra dar certo, bolinhas no meio da barrinha. >> dai tem que arrumar o limite do eixo x no plot: >> >> #denovo >> barras<-barplot(rainfall,axes=F) >> axis(2) >> #denovo vermelho pra gente ver >> axis(1,col="red") >> par(new=TRUE) >> >> >> >> plot(tmed~barras,type="b",ylim=c(5,35),xlim=c(0.2,14.4),frame.plot=F,axes=F) >> points(tmin~barras,type="b",lty=2) >> points(tmax~barras,type="b",lty=2) >> axis(4) >> #agora nao ficou bagunçado >> axis(1) >> >> #isso por causa do jeito que o barplot constroi o eixo X é diferente do >> plot >> #mas eu não entendo perfeitamente a diferença pra explicar >> #espero estar mais ajudando que atrapalhando agora. >> >> >> Em 30 de julho de 2012 19:28, Augusto Ribas <ribas.aca(a)gmail.com> >> escreveu: >>> Humm >>> >>> Veja se isso te ajuda: >>> >>> #Vai faze o grafico de barras >>> barplot(rainfall,axes=F) >>> #o eixo dele >>> axis(2) >>> >>> #vai recomeça o plot, na mesma janela >>> par(new=TRUE) >>> >>> #um plot normal de linha >>> plot(tmed~barras,type="b",ylim=c(5,35),frame.plot=F,axes=F) >>> #as outras 2 linhas >>> points(tmin~barras,type="b",lty=2) >>> points(tmax~barras,type="b",lty=2) >>> #vai fazer o eixo do outro lado (direito) >>> axis(4) >>> >>> #note que eu fiz os graficos sem eixos, e os adicionei a parte com o >>> comando axis, >>> #com o axis vc pode faze o eixo da forma que preferir. >>> #ai da pra controla cada grafico independente quanto a cores, tamanhos >>> eixos tudo >>> #se quiser mudar a escala, de uma olhada no comando ?scale(), mas acho >>> que talvez a solução acima seja mais >>> #o que vc deseja >>> >>> Em 30 de julho de 2012 17:27, Marcelo Claro de Souza >>> <marcelo_claro(a)yahoo.com.br> escreveu: >>>> Olá Augusto, >>>> Muito obrigado pela ajuda. >>>> Eu preciso que os dois eixos Y sejam independentes (escalas diferentes), >>>> de >>>> modo que as linhas fiquem sobre as barras. >>>> Tem como fazer isso? >>>> Muito obrigado. >>>> Abraço. >>>> >>>> >>>> Marcelo Claro de Souza >>>> Biologist, PhD student in Plant Biology >>>> Institute of Bioscience - UNESP, Brazil >>>> >>>> >>>> ________________________________ >>>> De: Augusto Ribas <ribas.aca(a)gmail.com> >>>> Para: r-br(a)listas.c3sl.ufpr.br; Marcelo Claro de Souza >>>> <marcelo_claro(a)yahoo.com.br> >>>> Enviadas: Segunda-feira, 30 de Julho de 2012 17:40 >>>> Assunto: Re: [R-br] barplot + lines com dois eixos Y >>>> >>>> Bom, de inicio acredito que vai ser algo assim: >>>> >>>> barras<-barplot(rainfall,ylim=c(0,35)) >>>> points(tmed~barras,type="l") >>>> points(tmin~barras,type="l",lty=2) >>>> points(tmax~barras,type="l",lty=2) >>>> axis(4) >>>> >>>> #vc vai salvar a posição das barras no vertor barras >>>> #ai depois note que vc esta fazendo as linhas >>>> # e por ultimo use o comando axes pra fazer o segundo ou mais eixos, >>>> se quiser algo no eixo X tb >>>> #leia a documentação pra deixa do jeito que gostar ?axis >>>> >>>> >>>> #agora não sei se entra nas suas necessidades mas >>>> #de uma olhada nesse pacote e nos graficos que da pra fazer com ele >>>> #é exatamente o tipo de dado que vc ilustrou >>>> #por exemplo: >>>> >>>> install.packages("climatol") >>>> require(climatol) >>>> >>>> data(datcli) >>>> #exemplo dele >>>> datcli >>>> #seus dados >>>> t(data.frame(rainfall,tmax,tmin,tmed)) >>>> >>>> #exemplo do comando >>>> diagwl(datcli,est="Example station",alt=100,per="1961-90",mlab="en") >>>> #seus dados de exemplo >>>> diagwl(t(data.frame(rainfall,tmax,tmin,tmed)),est="Example >>>> station",alt=100,per="1961-90",mlab="en") >>>> >>>> >>>> Em 30 de julho de 2012 15:18, Marcelo Claro de Souza >>>> <marcelo_claro(a)yahoo.com.br> escreveu: >>>>> Preciso plotar um gráfico de variáveis climáticas sendo o eixo Y >>>>> esquerdo >>>>> referente à precipitação e o eixo Y direito referente às temperaturas >>>>> min, >>>>> med e max. >>>>> A precipitação deve ser apresentada em barras e as temperaturas em >>>>> linhas. >>>>> Como posso fazer isso? >>>>> >>>>> x=c(seq(1:12)) >>>>> >>>>> >>>>> >>>>> rainfall=c(17.54,7.68,9.87,5.45,0.15,1.61,0.01,0.33,1.53,4.85,3.77,9.35)##barras >>>>> >>>>> >>>>> >>>>> tmax=c(30.21,30.71,27.23,28.23,25.19,24.39,26.59,28.39,29.77,28.49,28.61,29.79)##linhas >>>>> >>>>> >>>>> >>>>> tmin=c(16.82,16.36,15.91,15.8,11.57,9.07,13.07,12.65,12.98,15.47,16.21,17.26)##linhas >>>>> >>>>> >>>>> >>>>> tmed=c(23.51,23.53,21.57,22.02,18.38,16.73,19.83,20.52,21.38,21.98,22.41,23.53)##linhas >>>>> >>>>> Muito obrigado pela ajuda. >>>>> Abraço. >>>>> >>>>> Marcelo Claro de Souza >>>>> Biologist, PhD student in Plant Biology >>>>> Institute of Bioscience - UNESP, Brazil >>>>> >>>>> >>>>> _______________________________________________ >>>>> 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. >>>> >>>> >>>> >>>> -- >>>> Grato >>>> Augusto C. A. Ribas >>>> >>>> Site Pessoal: http://augustoribas.heliohost.org >>>> Lattes: http://lattes.cnpq.br/7355685961127056 >>>> >>>> >>> >>> >>> >>> -- >>> Grato >>> Augusto C. A. Ribas >>> >>> Site Pessoal: http://augustoribas.heliohost.org >>> Lattes: http://lattes.cnpq.br/7355685961127056 >> >> >> >> -- >> Grato >> Augusto C. A. Ribas >> >> Site Pessoal: http://augustoribas.heliohost.org >> Lattes: http://lattes.cnpq.br/7355685961127056 >> >> > > > > -- > Grato > Augusto C. A. Ribas > > Site Pessoal: http://augustoribas.heliohost.org > Lattes: http://lattes.cnpq.br/7355685961127056 > > -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

1 0

Agregar fatores
by Sérgio Henrique almeida da silva ju 01 Ago '12

01 Ago '12

Qual a melhor forma de agregar fatores? Eu tenho um banco com as variáveis cidade, fumo (1 = sim, 2 = não) e diversas outras, porém eu gostaria de criar um outro banco só com a cidade e a frequência de quem fuma ou não. Exemplo cidade sim não 1 aracaju 163 1843 2 belem 185 1824 3 belo horizonte 280 1727 4 boa vista 191 1828 5 brasilia 230 1780 6 campo grande 228 1782 Eu conseguir fazer de forma não muito elegante: tabela <- aggregate(dados$fumante, by = list(dados$CIDADE),FUN = function(x) c(table(x))) cidade <- data.frame(tabela) Como fazer de forma mais elegante? Outra pergunta, gostaria de criar dento desse banco (cidade) uma variável que é o % de sim, como fazer? Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

2 2

barplot + lines com dois eixos Y
by Marcelo Claro de Souza 30 Jul '12

30 Jul '12

Preciso plotar um gráfico de variáveis climáticas sendo o eixo Y esquerdo referente à precipitação e o eixo Y direito referente às temperaturas min, med e max. A precipitação deve ser apresentada em barras e as temperaturas em linhas. Como posso fazer isso? x=c(seq(1:12)) rainfall=c(17.54,7.68,9.87,5.45,0.15,1.61,0.01,0.33,1.53,4.85,3.77,9.35)##barras tmax=c(30.21,30.71,27.23,28.23,25.19,24.39,26.59,28.39,29.77,28.49,28.61,29.79)##linhas tmin=c(16.82,16.36,15.91,15.8,11.57,9.07,13.07,12.65,12.98,15.47,16.21,17.26)##linhas tmed=c(23.51,23.53,21.57,22.02,18.38,16.73,19.83,20.52,21.38,21.98,22.41,23.53)##linhas Muito obrigado pela ajuda. Abraço. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

2 1

Importar dados EpiData
by Pedro Emmanuel Alvarenga Americano do Brasil 30 Jul '12

30 Jul '12

Amigos de R, Há um software de entrada de dados que gosto muito chamado EpidData. ( www.epidata.dk) O Epidata utilizava até o ano passado o mesmo formato de dados que o antigo EpiInfo 6.4 em DOS. Os arquivos .rec. Então, para importar dados basta: library('foreign') dados <- read.epiinfo(...) Acontece que desde o ano passado o EpiData 3.1 esta sendo abandonado e substituído por um novo EpiData que está tentando se adequar a boas praticas clínicas com acesso hierarquizado aos dados entre outras coisas. O novo formato de arquivos são versoes XML do velho .rec... os epx e epz Alguem ja conseguiu importar dados nesse formato a partir do R? Isso está causando algumas dores de cabeça, pois não estou conseguindo migrar do velho para o novo EpiData devido a essa limitação! Cheguei mandar uma mensagem numa lista da turma do core do R mas o moderador rejeita as minhas mensagens. :-( Abraço, Dr. Pedro Emmanuel A. A. do Brasil Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365, CEP 21040-360, Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

2 2

Encontrar caminho de um arquivo no computador
by Luis Iván Ortiz Valencia 27 Jul '12

27 Jul '12

Prezados Gostaria de saber se é possível no R encontrar o caminho de um arquivo específico no computador. Por exemplo, existe um arquivo indicadores.pdf na pasta C:\ivan\R, existe uma função no R que responda com o caminho "C:\ivan\R" dando como entrada indicadores.pdf? atte IVAN

5 9

[Dúvida] Salvando data.frame em DB
by Pedro Rafael 27 Jul '12

27 Jul '12

Pessoal estou usando o comando dbWriteTable para salvar em um banco de dados PostgreSQL. O problema que mesmo colocando o argumento row.names = FALSE ele está salvando o banco com um ID. Só que eu gostaria que ele não criasse o ID pois ele já está criado nas tabelas do BD. Alguem sabe como resolver esse problema? [ ], Pedro Rafael Diniz Marinho.

1 0

regressão polinomial
by Tito Conte 26 Jul '12

26 Jul '12

Pessoal tenho o seginte conjunto de dados: f=c(-26.0000,-25.9862,-25.9343,-25.8822,-25.8433,-25.8054,-25.7948,-25.7872,-25.7668,-25.7284,-25.7015,-25.6282,-25.4612,-25.3016,-25.2564,-25.2412,-25.2412,-25.2232,-25.0869,-25.0000,-25.0000,-24.9856,-24.9397,-24.8976,-24.8533,-24.6587,-24.6373,-24.5740,-24.5406,-24.4379,-24.3934,-24.3628,-24.3277,-24.3042,-24.2972,-24.2973,-24.3020,-24.3068,-24.3068,-24.2956,-24.2920) w=c(-45.8529,-45.8302,-45.7575,-45.6680,-45.6107,-45.5050,-45.4574,-45.4058,-45.3538,-45.2913,-45.2634,-45.2264,-45.1667,-45.0719,-45.0213,-45.0000,-45.0000,-44.9748,-44.8726,-44.8405,-44.8405,-44.8351,-44.8259,-44.7972,-44.7764,-44.6867,-44.6674,-44.5912,-44.5587,-44.4927,-44.4515,-44.4115,-44.3314,-44.2664,-44.2171,-44.1729,-44.1288,-44.0820,-44.0509,-44.0118,-44.0000) se vocês plotarem eles, obteram uma curva, gostaria de obter a função desta curva ou uma regressão polinomial de grau elevado dela já tentei usar a função lm e não deu um bom resultado, quando chego em graus elevados e calculo utilizando os coeficientes obtidos os valores dão 10 vezes maiores (não sei bem porque) mas quando ploto ela aparece bonitinha no lugar que deveria aprecer em outro caso tentei na raça mesmo por multiplicação matricial mas também deu errado pois os valores da multiplicação da inversa são muito pequenos alguém pode me dar uma sugestão ou uma luz do que eu estou fazendo de errado? Tito Conte

4 14