- R-br - DINF UFPR

Diagrama de voronoi
by Fabrício Tavares 30 Mai '11

30 Mai '11

Boa tarde pssoal, Estou usando o pacote tripack para a construção de diagrama de voronoi. > data(tritest) > vm <- voronoi.mosaic(tritest) > plot.voronoi(vm, do.points=T) > points(tritest) Nesse exemplo, a "borda" do diagrama é um retângulo. É possível utilizar alguma função já existente para criar o diagrama de forma que a "borda" seja o polígono que eu especificar? (Algo como o que é feito no software MapInfo, para os que o conhecem) Fabrício. http://vizzual.net

2 2

Análise fatorial e regressão em painel
by Henrique Brilhante 30 Mai '11

30 Mai '11

Bom dia senhores, já acompanho a lista a algum tempo pois estou desenvolvendo um trabalho na graduação sobre o modelo APT. Dado meu conhecimento bem limitado do software gostaria de saber se alguém pode me auxiliar com o problema abaixo de forma passo-a-passo, mesmo que fora desta lista de discussão. VarAcao = PIBreal + Ptax + Ouro + IBOVESPA + Dow Jones + Nasdaq + Crise2008 (dummy) O método seria realizar a análise fatorial dos 6 fatores e salvar os betas para a regressão em painel considerando o efeito da crise de 2008. Abaixo um exemplo dos dados (3 das 40 empresas do estudo). Data Acao VarAcao PIBreal Ptax Ouro IBOVESPA Dow Jones Nasdaq Crise2008 1 T 2005 VALE - PNA 0,0949 -0,0623 -0,0552 -0,0757 -0,0439 -0,0838 -0,1420 0,0000 2 T 2005 VALE - PNA -0,1857 -0,0352 -0,1851 -0,1677 -0,1194 -0,0739 -0,0234 0,0000 3 T 2005 VALE - PNA 0,3734 -0,0651 -0,1112 -0,0357 0,1766 -0,0411 -0,0243 0,0000 4 T 2005 VALE - PNA -0,0332 -0,0499 -0,0078 0,0745 -0,0022 -0,0191 -0,0085 0,0000 1 T 2006 VALE - PNA 0,0960 -0,0374 -0,1294 -0,0117 0,0713 -0,0197 0,0036 0,0000 2 T 2006 VALE - PNA -0,0473 -0,0339 -0,0406 0,0280 -0,0723 -0,0476 -0,1256 0,0000 3 T 2006 VALE - PNA -0,0937 -0,0238 -0,0352 -0,0752 -0,0447 0,0111 0,0037 0,0000 4 T 2006 VALE - PNA 0,2958 -0,0309 -0,0593 0,0045 0,1564 0,0391 0,0412 0,0000 1 T 2007 VALE - PNA 0,1618 -0,0235 -0,0847 -0,0134 -0,0134 -0,0566 -0,0452 0,0000 2 T 2007 VALE - PNA 0,1198 -0,0211 -0,0996 -0,0899 0,1347 0,0383 0,0287 0,0000 3 T 2007 VALE - PNA 0,3559 -0,0257 -0,0832 -0,0022 0,0690 0,0071 0,0084 0,0000 4 T 2007 VALE - PNA -0,0349 -0,0232 -0,0780 0,0547 0,0145 -0,0803 -0,0522 0,0000 1 T 2008 VALE - PNA -0,0143 -0,0237 -0,0537 0,0210 -0,0878 -0,1210 -0,1941 0,0000 2 T 2008 VALE - PNA -0,0749 -0,0307 -0,1425 -0,1398 0,0159 -0,1294 -0,0460 0,0000 3 T 2008 VALE - PNA -0,3879 -0,0273 0,1415 0,0941 -0,3148 -0,0771 -0,1285 1 4 T 2008 VALE - PNA -0,2963 -0,0874 0,1551 0,1265 -0,3216 -0,2058 -0,2716 1 1 T 2009 VALE - PNA 0,1008 -0,0607 -0,0507 0,0053 0,0448 -0,1835 -0,0719 1 2 T 2009 VALE - PNA 0,1140 -0,0184 -0,2078 -0,1473 0,1922 0,0666 0,1449 1 3 T 2009 VALE - PNA 0,1976 -0,0028 -0,1213 -0,0717 0,1502 0,1164 0,1223 0,0000 4 T 2009 VALE - PNA 0,1440 -0,0073 -0,0525 0,0447 0,0773 0,0502 0,0460 0,0000 1 T 2010 VALE - PNA 0,1401 -0,0191 -0,0181 0,0517 -0,0150 0,0122 0,0272 0,0000 2 T 2010 VALE - PNA -0,2695 -0,0164 -0,0208 0,0537 -0,1762 -0,1288 -0,1521 0,0000 3 T 2010 VALE - PNA 0,1949 -0,0272 -0,0926 -0,0285 0,0993 0,0703 0,0876 0,0000 4 T 2010 VALE - PNA 0,0360 -0,0402 -0,0643 0,0509 -0,0495 0,0416 0,0843 0,0000 1 T 2005 PETROBRAS PN 0,0422 -0,0623 -0,0552 -0,0757 -0,0439 -0,0838 -0,1420 0,0000 2 T 2005 PETROBRAS PN 0,0424 -0,0352 -0,1851 -0,1677 -0,1194 -0,0739 -0,0234 0,0000 3 T 2005 PETROBRAS PN 0,2955 -0,0651 -0,1112 -0,0357 0,1766 -0,0411 -0,0243 0,0000 4 T 2005 PETROBRAS PN 0,0263 -0,0499 -0,0078 0,0745 -0,0022 -0,0191 -0,0085 0,0000 1 T 2006 PETROBRAS PN 0,1471 -0,0374 -0,1294 -0,0117 0,0713 -0,0197 0,0036 0,0000 2 T 2006 PETROBRAS PN 0,0141 -0,0339 -0,0406 0,0280 -0,0723 -0,0476 -0,1256 0,0000 3 T 2006 PETROBRAS PN -0,0664 -0,0238 -0,0352 -0,0752 -0,0447 0,0111 0,0037 0,0000 4 T 2006 PETROBRAS PN 0,2168 -0,0309 -0,0593 0,0045 0,1564 0,0391 0,0412 0,0000 1 T 2007 PETROBRAS PN -0,0797 -0,0235 -0,0847 -0,0134 -0,0134 -0,0566 -0,0452 0,0000 2 T 2007 PETROBRAS PN 0,1124 -0,0211 -0,0996 -0,0899 0,1347 0,0383 0,0287 0,0000 3 T 2007 PETROBRAS PN 0,1386 -0,0257 -0,0832 -0,0022 0,0690 0,0071 0,0084 0,0000 4 T 2007 PETROBRAS PN 0,3949 -0,0232 -0,0780 0,0547 0,0145 -0,0803 -0,0522 0,0000 1 T 2008 PETROBRAS PN -0,1893 -0,0237 -0,0537 0,0210 -0,0878 -0,1210 -0,1941 0,0000 2 T 2008 PETROBRAS PN 0,2043 -0,0307 -0,1425 -0,1398 0,0159 -0,1294 -0,0460 0,0000 3 T 2008 PETROBRAS PN -0,2857 -0,0273 0,1415 0,0941 -0,3148 -0,0771 -0,1285 1 4 T 2008 PETROBRAS PN -0,4048 -0,0874 0,1551 0,1265 -0,3216 -0,2058 -0,2716 1 1 T 2009 PETROBRAS PN 0,2109 -0,0607 -0,0507 0,0053 0,0448 -0,1835 -0,0719 1 2 T 2009 PETROBRAS PN 0,1264 -0,0184 -0,2078 -0,1473 0,1922 0,0666 0,1449 1 3 T 2009 PETROBRAS PN 0,0791 -0,0028 -0,1213 -0,0717 0,1502 0,1164 0,1223 0,0000 4 T 2009 PETROBRAS PN 0,0478 -0,0073 -0,0525 0,0447 0,0773 0,0502 0,0460 0,0000 1 T 2010 PETROBRAS PN -0,0565 -0,0191 -0,0181 0,0517 -0,0150 0,0122 0,0272 0,0000 2 T 2010 PETROBRAS PN -0,2710 -0,0164 -0,0208 0,0537 -0,1762 -0,1288 -0,1521 0,0000 3 T 2010 PETROBRAS PN 0,0181 -0,0272 -0,0926 -0,0285 0,0993 0,0703 0,0876 0,0000 4 T 2010 PETROBRAS PN -0,0089 -0,0402 -0,0643 0,0509 -0,0495 0,0416 0,0843 0,0000 1 T 2005 ITAUUNIBANCO PN 0,0771 -0,0623 -0,0552 -0,0757 -0,0439 -0,0838 -0,1420 0,0000 2 T 2005 ITAUUNIBANCO PN -0,0190 -0,0352 -0,1851 -0,1677 -0,1194 -0,0739 -0,0234 0,0000 3 T 2005 ITAUUNIBANCO PN 0,2088 -0,0651 -0,1112 -0,0357 0,1766 -0,0411 -0,0243 0,0000 4 T 2005 ITAUUNIBANCO PN 0,0466 -0,0499 -0,0078 0,0745 -0,0022 -0,0191 -0,0085 0,0000 1 T 2006 ITAUUNIBANCO PN 0,1268 -0,0374 -0,1294 -0,0117 0,0713 -0,0197 0,0036 0,0000 2 T 2006 ITAUUNIBANCO PN -0,0140 -0,0339 -0,0406 0,0280 -0,0723 -0,0476 -0,1256 0,0000 3 T 2006 ITAUUNIBANCO PN 0,0326 -0,0238 -0,0352 -0,0752 -0,0447 0,0111 0,0037 0,0000 4 T 2006 ITAUUNIBANCO PN 0,1629 -0,0309 -0,0593 0,0045 0,1564 0,0391 0,0412 0,0000 1 T 2007 ITAUUNIBANCO PN -0,0672 -0,0235 -0,0847 -0,0134 -0,0134 -0,0566 -0,0452 0,0000 2 T 2007 ITAUUNIBANCO PN 0,1719 -0,0211 -0,0996 -0,0899 0,1347 0,0383 0,0287 0,0000 3 T 2007 ITAUUNIBANCO PN 0,0653 -0,0257 -0,0832 -0,0022 0,0690 0,0071 0,0084 0,0000 4 T 2007 ITAUUNIBANCO PN -0,0243 -0,0232 -0,0780 0,0547 0,0145 -0,0803 -0,0522 0,0000 1 T 2008 ITAUUNIBANCO PN -0,1238 -0,0237 -0,0537 0,0210 -0,0878 -0,1210 -0,1941 0,0000 2 T 2008 ITAUUNIBANCO PN -0,0016 -0,0307 -0,1425 -0,1398 0,0159 -0,1294 -0,0460 0,0000 3 T 2008 ITAUUNIBANCO PN -0,0265 -0,0273 0,1415 0,0941 -0,3148 -0,0771 -0,1285 1 4 T 2008 ITAUUNIBANCO PN -0,1913 -0,0874 0,1551 0,1265 -0,3216 -0,2058 -0,2716 1 1 T 2009 ITAUUNIBANCO PN -0,0153 -0,0607 -0,0507 0,0053 0,0448 -0,1835 -0,0719 1 2 T 2009 ITAUUNIBANCO PN 0,1808 -0,0184 -0,2078 -0,1473 0,1922 0,0666 0,1449 1 3 T 2009 ITAUUNIBANCO PN 0,2339 -0,0028 -0,1213 -0,0717 0,1502 0,1164 0,1223 0,0000 4 T 2009 ITAUUNIBANCO PN 0,0709 -0,0073 -0,0525 0,0447 0,0773 0,0502 0,0460 0,0000 1 T 2010 ITAUUNIBANCO PN 0,0035 -0,0191 -0,0181 0,0517 -0,0150 0,0122 0,0272 0,0000 2 T 2010 ITAUUNIBANCO PN -0,1913 -0,0164 -0,0208 0,0537 -0,1762 -0,1288 -0,1521 0,0000 3 T 2010 ITAUUNIBANCO PN 0,2218 -0,0272 -0,0926 -0,0285 0,0993 0,0703 0,0876 0,0000 4 T 2010 ITAUUNIBANCO PN -0,0381 -0,0402 -0,0643 0,0509 -0,0495 0,0416 0,0843 0,0000 Agradeço a gentileza, Muito obrigado. Henrique. brilhante.henrique(a)gmail.com

2 1

[Dúvida] Plotando pontos centroides no MAPA
by Pedro Rafael 30 Mai '11

30 Mai '11

Pessoal como faço para plotar no map as coordenadas dos centróides dos municípios. Importei o mapa com os comandos que seguem require(SpatialEpi) require(maptools) dados = read.table("C:\\Users\\GORR-PEDRO\\Desktop\\Malha Municipal PB\\dados.txt", header = TRUE) attach(dados) mapa = readShapePoly("C:\\Users\\GORR-PEDRO\\Desktop\\Malha Municipal PB\\PB.shp") geo <- coordinates(mapa) geo <- latlong2grid(geo) Gostaria de plotar no mapa o objeto "geo" Pedro

1 0

Leitura de arquivo
by Emerson Cotta Bodevan 30 Mai '11

30 Mai '11

Prezados, como faço para ler meus dados no formato que está abaixo para um data frame? Meus dados devem ser compostos apenas por três colunas (coordenadas e pesos). O número que representa a quantidade de pontos não deve ser lido. 3 6.164639 6.164639 2 9.881965 3.717327 3 22.054118 3.217595 8 5 11.905148 11.905148 2 14.344441 2.439293 3 16.187096 1.842655 4 16.985997 0.798901 5 16.989722 0.003725 6 7 4.755562 4.755562 2 8.314285 3.558723 3 10.263021 1.948736 4 11.699839 1.432584 6 12.455548 0.755709 7 12.551566 0.096018 8 12.552208 0.000642 9 Agradeço qualquer ajuda, -- Emerson

3 6

Leitura de arquivo 2
by Emerson Cotta Bodevan 30 Mai '11

30 Mai '11

Prezados, Meus dados estao no formato abaixo, onde o numero sozinho representa o numero de informacoes que vem logo a seguir. As informacoes estao em tres colunas (coordenadas e peso). Preciso tanto do numero de pontos quanto das informacoes relativas a eles. O que e melhor de trabalhar nesse caso? Com lista ou com data frame? E como ler esse arquivo automaticamente? 3 6.164639 6.164639 2 9.881965 3.717327 3 22.054118 3.217595 8 5 11.905148 11.905148 2 14.344441 2.439293 3 16.187096 1.842655 4 16.985997 0.798901 5 16.989722 0.003725 6 7 4.755562 4.755562 2 8.314285 3.558723 3 10.263021 1.948736 4 11.699839 1.432584 6 12.455548 0.755709 7 12.551566 0.096018 8 12.552208 0.000642 9 Agradeço qualquer ajuda, -- Emerson

2 2

[Dúvida] Centroides dos Poligonos
by Pedro Rafael 30 Mai '11

30 Mai '11

Pessoal, como faço para gerar os centroides de cada um dos municípios do Estado da Paraíba, existe alguma forma de gerar apartir do arquivo .shp? Baxei a malha no site do IBGE mas dentro do arquivo dBase (.dbf) não tem as coordenadas dos controides de cada município. Obrigado desde já. Pedro

2 3

Ajuda em modelo misto
by Emmanuel Arnhold 28 Mai '11

28 Mai '11

Ajustei o seguinte modelo misto: ga=interaction(genótipo,ambiente) modelo=lmer(resposta~bloco/ambiente+ambiente+(1|genótipo)+(1|ga)) É um modelo de avaliação de genótipos em vários ambientes. Neste modelo tenho as seguintes variâncias: ga genótipo residual Gostaria que alguém me ajuda-se a ajustar o modelo acima para obter as seguintes variâncias, supondo 3 ambientes: ga genótipo residual (ambiente 1) residual(ambiente 2) residual (ambiente 3)

3 5

Package 'dicionariosIBGE'
by Arnaldo Gomes de Brito Junior 28 Mai '11

28 Mai '11

Prezados; Estou com dificuldades de utilizar o Package 'dicionariosIBGE'. Consigo instalr e visualizo os dicionários, mas como faço para ligara minha base da POF(por exemplo) com os dicionários do "dicionariosIBGE". Grato desde já. Arnaldo Gomes de Brito Junior - Estatístico (21) - 8133-4310 (21) - 3734-1043

3 3

Re: [R-br] Ajuda em plotar gráfico com mais de uma série.
by Marcos Batista 27 Mai '11

27 Mai '11

Prezado, Walmes, muito obrigado pela ajuda. Assim que eu tentar lhe digo se funcionou ou não. Cordialmente, Marcos Batista -- Marcos Roberto Dias Batista Departamento de Genética, Evolução e Bioagentes Instituto de Biologia, Unicamp Caixa Postal 6109, Campinas, São Paulo CEP 13083-970, Brasil Telefones: (19) 3521 1151 Message: 7 Date: Fri, 27 May 2011 11:32:24 -0300 From: Marcos Batista <mrdbatista(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: [R-br] Ajuda em plotar gráfico com mais de uma série. Message-ID: <BANLkTimUo1mq5SOXEpika+hv9mMSpKJErw(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Bom dia a todos da lista, sou iniciante no uso do programa e não achei em nenhum guia do programa como posso adicionar mais uma série de dados para fazer um gráfico. Vocês tem uma dica para me dar? Agradeço desde já, Cordialmente, Marcos Batista -- Marcos Roberto Dias Batista Departamento de Genética, Evolução e Bioagentes Instituto de Biologia, Unicamp Caixa Postal 6109, Campinas, São Paulo CEP 13083-970, Brasil Telefones: (19) 3521 1151 -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: < http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110527/86b9991b/atta… > ------------------------------ Message: 8 Date: Fri, 27 May 2011 11:45:05 -0300 From: Walmes Zeviani <walmeszeviani(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Ajuda em plotar gráfico com mais de uma série. Message-ID: <BANLkTikr1-3f=-xRNk2sq6R0=zdnZt6tvg(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Marcos, Você pode adicionar a nova série ao gráfico da primeira com points() e/ou lines(). Você pode também usar a lattice. x <- 1:100 y1 <- rnorm(x) y2 <- rnorm(y) plot(y1~x, type="l") lines(y2~x) require(lattice) xyplot(y1+y2~x) Código não testado. Verifique a documentação das funções e forneça CMR no futuro. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218

1 0

Ajuda em plotar gráfico com mais de uma série.
by Marcos Batista 27 Mai '11

27 Mai '11

Bom dia a todos da lista, sou iniciante no uso do programa e não achei em nenhum guia do programa como posso adicionar mais uma série de dados para fazer um gráfico. Vocês tem uma dica para me dar? Agradeço desde já, Cordialmente, Marcos Batista -- Marcos Roberto Dias Batista Departamento de Genética, Evolução e Bioagentes Instituto de Biologia, Unicamp Caixa Postal 6109, Campinas, São Paulo CEP 13083-970, Brasil Telefones: (19) 3521 1151

2 1