- 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

Duvida na organozacao de dados hdf5
by Humberto Hazin (hotmail) 26 Mai '11

26 Mai '11

Boa Noite Pessoal, Estou com um problema e pode ser que vcs consigam resolver! A dúvida é a seguinte: Baixei uma arquivo de formato hdf4 que contem valores de latitude, longitude e SST. Em seguida converti o arquivo para hdf5 (usando o programa hdf4toh5) quando tento extrair os valores por colunas com seus respectivos valores ex. Latitude longitude SST 24 30 26 24 30 22 24 30 25 Ele não vai! Ele se apresenta assim!!! Alguém tem alguma solução para esse problema que poderia compartilhar comigo? Humberto $longitude [1] -129.9711 -129.9272 -129.8832 -129.8393 -129.7953 -129.7514 -129.7074 -129.6635 -129.6196 [10] -129.5756 -129.5316 -129.4877 -129.4438 -129.3998 -129.3559 -129.3119 -129.2680 -129.2240 [19] -129.1801 -129.1361 -129.0922 -129.0482 -129.0043 -128.9604 -128.9164 -128.8725 -128.8285 [28] -128.7846 -128.7406 -128.6967 -128.6527 -128.6088 -128.5648 -128.5209 -128.4769 -128.4330..... $latitude [1] 20.01753 20.06148 20.10542 20.14937 20.19331 20.23726 20.28121 20.32515 20.36910 20.41305 20.45699 [12] 20.50094 20.54489 20.58883 20.63278 20.67672 20.72067 20.76462 20.80856 20.85251 20.89646 20.94040 [23] 20.98435 21.02829 21.07224 21.11619 21.16013 21.20408 21.24803 21.29197 21.33592 21.37986 21.42381 [34] 21.46776 21.51170 21.55565 21.59960 21.64354 21.68749 21.73143 21.77538 21.81933 21.86327 21.90722 [45] 21.95117 21.99511 22.03906 22.08300 22.12695 22.17090 22.21484 22.25879 22.30274 22.34668 22.39063 [56] 22.43457 22.47852 22.52247 22.56641 22.61036 22.65431 22.69825 22.74220 22.78614 22.83009 22.87404 [67] 22.91798 22.96193 23.00588 23.04982 23.09377 23.13771 23.18166 23.22561 23.26955 23.31350 23.35745 $param [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 21.150 21.150 21.150 21.150 21.150 21.075 21.075 21.075 21.075 21.075 21.075 21.075 21.075 [2,] 21.150 21.150 21.150 21.150 21.075 21.075 21.075 21.075 21.075 21.075 21.075 21.075 21.075 [3,] 21.225 21.225 21.075 21.075 21.075 21.075 21.075 21.075 21.075 21.075 21.075 21.000 21.000 [4,] 21.300 21.150 21.075 21.075 21.075 21.075 21.075 21.075 21.075 21.000 21.000 21.000 21.000 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [,26] [1,] 21.075 21.000 21.000 21.000 21.000 21.000 21.000 21.000 20.925 20.925 20.925 20.925 20.925 [2,] 21.000 21.000 21.000 21.000 21.000 21.000 20.925 20.925 20.925 20.925 20.925 20.925 20.925 [3,] 21.000 21.000 21.000 21.000 20.925 20.925 20.925 20.925 20.925 20.925 20.925 21.450 21.525 [4,] 21.000 21.000 20.925 20.925 20.925 20.925 20.925 20.925 20.925 20.925 20.925 20.925 21.450 [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [,38] [,39] [1,] 20.925 20.925 21.525 21.525 21.525 21.525 21.450 21.525 21.450 21.450 21.375 21.375 21.375 [2,] 20.925 20.925 20.925 20.850 20.625 21.450 21.450 21.450 21.450 21.375 21.375 21.375 21.375 [3,] 21.450 20.925 20.775 20.625 20.625 21.450 21.450 21.450 21.450 21.450 21.375 21.375 21.375 [4,] 21.450 21.450 20.850 20.850 21.450 21.450 21.450 21.450 21.450 21.375 21.375 21.375 21.375 [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [,50] [,51] [,52] [1,] 21.375 21.375 21.300 21.300 21.300 21.300 21.300 21.225 21.225 21.225 21.225 21.225 21.225 [2,] 21.375 21.300 21.300 21.300 21.300 21.300 21.225 21.225 21.225 21.225 21.225 21.150 21.225 [3,] 21.375 21.375 21.300 21.300 21.300 21.300 21.225 21.225 21.150 21.150 21.150 21.150 21.150 [4,] 21.375 21.300 21.300 21.300 21.225 21.225 21.225 21.225 21.225 21.150 21.150 21.150 21.075 [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [,62] [,63] [,64] [,65] [1,] 21.000 21.000 21.075 21.075 21.075 20.850 21.075 20.475 20.475 21.075 21.000 21.000 21.075 [2,] 21.150 21.150 21.150 21.150 20.850 21.075 21.075 21.075 20.475 21.000 21.075 21.075 21.000 [3,] 21.150 21.075 21.075 21.075 21.075 21.075 21.000 21.000 21.075 21.075 21.000 21.000 21.000 [4,] 21.150 21.075 21.075 21.075 21.075 21.000 21.000 21.075 21.075 21.000 21.000 21.000 21.000 [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [,74] [,75] [,76] [,77] [,78]

1 0

Aviso de erro
by Rodrigo Pereira 26 Mai '11

26 Mai '11

Caros colegas, bom dia! Estou fazendo a atualização dos pacotes que estão instalados no R, mas está aparecendo a seguinte mensagem de erro: Aviso em install.packages(update[instlib == l, "Package"], l, contriburl = contriburl, : 'lib = "C:/Program Files/R/R-2.13.0/library"' is not writable Erro em install.packages(update[instlib == l, "Package"], l, contriburl = contriburl, : unable to install packages Alguém pode me explicar o que é isso e como posso contornar caso isso seja um problema? Desde já obrigado! Abraço! Rodrigo Pereira

2 2

Ajuda com erro em Loop
by Victor Eduardo 26 Mai '11

26 Mai '11

Galera, estou com um banco de dados que tem várias estações, e preciso criar histograma para cada estação, então criei um loop para criar histograma para cada estacao, sendo que algumas estações que eu pego da o seguinte erro quando eu uso o comando abaixo: ifelse(length(na.omit(obs_clcl))>=1,hist(obs_clcl,main=tit.obs),message("sem grafico")) Erro em ans[!test & !nas] <- rep(no, length.out = length(ans))[!test & : substituto tem comprimento zero O grande problema é que o loop trava quando aparece esse erro e simplesmente acaba. Abraços, Victor Eduardo

2 5

Swave não gera figura
by Andre Oliveira 26 Mai '11

26 Mai '11

Bom dia, estou no ubuntu e tenho um arquivo Sweave muito grande que estava tudo funcionado mas, ontem fui processar o mesmo e deu pau. Agora vai saber onde é este erro dentro deste arquivo. Descobri, depois de muita tentativas, que em algum lugar por algum motivo este pedaço <<fig=T>>= boxplot(rnorm(999)) @ Não gera a imagem imagem apesar de inserir a mesma no arquivo LaTeX. Logo meu problema é com o arquivo .tex e não com os códigos R. Estou no RStudio e depois de processar o códigos R aparece este erro já na parte de processo do .tex este erro. . . . . Output written on macro_latex.pdf (5 pages, 221782 bytes). Transcript written on macro_latex.log. Erro em tools:::texi2dvi(file = fileName, pdf = TRUE, clean = TRUE, index = index) : unable to run 'pdflatex' on 'macro_latex.tex' Se alguém tiver alguma dica ficarei grato. ps.As configurações especificar de .rnw no preâmbulo são. \usepackage{Sweave} \SweaveOpts{eps=F} \SweaveOpts{keep.source=TRUE} AOS André Oliveira Souza

2 1

[Dúvida] Importar dados para um mapa
by Pedro Rafael 25 Mai '11

25 Mai '11

Pessoal estou querendo criar um mata de níveis de cores no R e baixei os arquivos para gerar o mapa e vi dentro do arquivo Dbase(.dbf) o geocodigo de cada município do Estado da Paraíba. Contudo as linhas da base de dados não estão na mesma sequência que está no arquivo dbf isso tem algum problema? Obrigado à todos desde já, Pedro

2 1

Re: [R-br] Outro problema com datas!!
by Ivan Bezerra Allaman 25 Mai '11

25 Mai '11

Olá Tura! Está acontecendo examente isso, as datas são interpretadas como strings. Elas são do tipo "Mês-Ano". Tive que fazer todo aquele procedimento que o Benilton passou porque estava com problemas ao plotar os gráficos. Na verdade ainda estou, mais contornei de uma outra forma menos elegante. Já tentei fazer um "as.Date" e outras coisas, mais nada disso faz com o R entenda que aqueles números são datas. É provável que haja algum padrão do tipo "AA - MM - DD", ou algo do gênero. Abraços. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

3 2

Outro problema com datas!!
by Ivan Bezerra Allaman 25 Mai '11

25 Mai '11

Novamente venho recorrer a ajuda dos senhores!! Toda vez que trabalho com datas a vontade que eu tenho é de matar o R. Os meus dados estão dispostos da seguinte maneira: dao[1:15,] trati dat vari 1 1 03-08 0.0000000 2 1 04-08 0.3472222 3 1 05-08 0.0000000 4 1 06-08 1.7361111 5 1 07-08 2.7777778 6 1 08-08 2.0833333 7 1 09-08 2.4305556 8 1 10-08 0.0000000 9 1 11-08 0.0000000 10 1 12-08 0.0000000 11 1 01-09 0.0000000 12 1 02-09 26.7361111 13 1 03-09 70.4861111 14 1 04-09 99.3055556 15 1 05-09 99.3055556 No entanto, a boneca do R ordena as datas da seguinte maneira: levels(dao$dat) [1] "01-09" "01-10" "02-09" "02-10" "03-08" "03-09" "03-10" "04-08" "04-09" [10] "04-10" "05-08" "05-09" "05-10" "06-08" "06-09" "06-10" "07-08" "07-09" [19] "07-10" "08-08" "08-09" "08-10" "09-08" "09-09" "09-10" "10-08" "10-09" [28] "10-10" "11-08" "11-09" "11-10" "12-08" "12-09" "12-10" Ou seja, ele ordena da maneira que ele acha mais conveniente e não da maneira que eu coloquei. Como faço para deixas os levels da maneira como eu VEJO os dados, ou seja, ordenado por ANO? Mais uma vez obrigado a todos e desculpem a ignorância. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

4 3

Dúvida Inferência Bayesiana
by Victor Eduardo 25 Mai '11

25 Mai '11

Galera, qual comando que me da a distribuição gama inversa? Como que calculo intervalo de credibilidade? Abraços, Victor Eduardo

5 6

Re: [R-br] Outro problema com datas!!
by Ivan Bezerra Allaman 24 Mai '11

24 Mai '11

Benilton, Nem sempre há a necessidade de colocar um exemplo reproduzível, tanto é verdade que a sua solução agora resolveu o problema. Obrigado!! Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

2 1

Re: [R-br] Outro problema com datas!!
by Ivan Bezerra Allaman 24 Mai '11

24 Mai '11

Não concordo contigo e nem como o R interpreta as coisas. Ele deveria interpretar da maneira como colocamos na base de dados raiz. Afinal de contas, se colocamos daquele jeito é pq é assim que deve ser. Quanto a sua sugestão eu já havia pensado antes, mais operacionalmente é desgastante, pois minha base tem 1 e poucas linhas com datas que vão de 2008 a 2010 e de maneira alternada. Ficaria até amanhã procedendo daquele modo. Outras sugestões? Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

2 1

ajuda teste hipotese duas variáveis
by carla eloize carducci 24 Mai '11

24 Mai '11

oi pessoal..Oi Daniel, obrigada pela ajuda, respondendo a pergunta.. sim a variável salario é numerica e a variavel instrução é fator..será isso o problema? como mudo, o arquivo é um data.frame. reposta Daniel: Não deu certo como? O que saiu no R? A variável salário esta adequadamente reconhecida como numérica? 2011/5/24 carla eloize carducci <elocarducci(a)hotmail.com> > OI PESSOAL > sou uma aprendiz bem iniciante do R, preciso de ajuda com um teste de > hipótese, alguém pode me ajudar? > bom, é o seguinte: tenho um conj de dados com duas variáveis, quero fazer > um t.test para verificar se as médias dos salarios (variável numerica) > de funcionários com diferente grau de instrução (2ºgrau e superior, > variável categorica) são diferentes ou não. > tentei fazer o seguinte: meu arquivo é da apostila do prof. justiniano o > velho arquivo 'milsa'.(bem de iniciante) > # teste de hipótese: H0:salario 2grau = superior# formulei essa hipotese > # H1:salario 2grau diferente superior > > milsa$instrucao <- factor(milsa$instrucao, label = c("2oGrau", "Superior"), > lev = 1:2, ord = T) > var.test(milsa$salario~milsa$instrucao) > alpha <- 0.05 # especifica valor da signific^ancia > if (var.test$p.value > alpha) > diff.test <- t.test(milsa$salario ~ milsa$instrucao,alternative = > "two.sided", var.equal = TRUE, conf.level = 1-alpha) else > diff.test <- t.test(milsa$salario ~ milsa$instrucao, alternative = > "two.sided", var.equal = FALSE,conf.level = 1-alpha) > diff.test > não deu certo, o que estou esquecendo? não acho o erro. Atenciosamente, Carla Eloize Carducci Engº Agronoma (UFGD) Msc. Ciência do Solo (UFLA) Doutoranda em Ciência do Solo (UFLA)

1 0

como dividir um mapa em áreas retangulares?
by Gustavo Henrique de Carvalho 24 Mai '11

24 Mai '11

Olá pessoal, Por favor, como eu faço para dividir um mapa em uma grade? Eu preciso que o tamanho da grade seja ajustável e que todos os retângulos tenham pela menos uma porção da sua área no continente: require('maps') require('spatial') require('maptools') brazil <- map("world", "brazil", fill = TRUE, plot = FALSE) ids <- sapply(strsplit(brazil$names, ":"), function(x) x[1]) llCRS <- CRS("+proj=longlat +ellps=WGS84") brazil.poly <- map2SpatialPolygons(brazil, IDs = ids, proj4string = llCRS) A representação gráfica aqui não é importante. Eu preciso apenas das coordenadas dos quatro cantos de cada retângulo passar para um webservice. Como há um limite de requisições que eu posso fazer, preciso dar um jeito de não requisitar dados para áreas no meio do oceano. Farei isso para todo o continente americano. Eu nunca trabalhei com mapas mas imaginei que, como eu não tenho um objeto da classe SpatialGrid, eu poderia criar um mapa de classe SpatialPolygon, usar a função spsample para conseguir pontos no continente e depois transformar esses pontos em uma grade. De qualquer forma, não funcionou: spsample(brazil.poly, 100, type = "regular") # Error in .local(x, n, type, ...) : cannot sample in zero-area polygons Realmente há polígonos com área 0 nesse mapa e eu não sei como resolver isso. Obrigado, Gustavo. > sessionInfo() R version 2.13.0 (2011-04-13) Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit) locale: [1] C/en_US.UTF-8/C/C/C/C attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] maptools_0.8-7 sp_0.9-81 spatial_7.3-3 maps_2.1-6 foreign_0.8-43 sos_1.3-0 [7] brew_1.0-6 igraph_0.5.5-2 Matrix_0.999375-50 lattice_0.19-23 loaded via a namespace (and not attached): [1] grid_2.13.0 raster_1.8-22 tools_2.13.0

1 1

OT: LaTeX - Nota de Rodape - Era: Re: OFF topic
by Benilton Carvalho 24 Mai '11

24 Mai '11

Assumindo que vc esteja se referindo a LaTeX, vc tbm deve prover um CMR... O que eu tentei aqui funciona normalmente, assumido que eu tenha entendido o que vc quer dizer: %% comeco \documentclass{article} \begin{document} afdkjfdalk fda\footnote{blah} lakjfsalk ajdjlfkdaj lf da\footnote{bleh}. \end{document} %% fim nao e' isso que vc esta' obtendo? 2011/5/24 Fernando Henrique Toledo <fernandohtoledo(a)gmail.com>: > Senhores, boa tarde... > > Totalmente fora de tópico, mas vocês são minha ultima alternativa a se > recorrer! > > Como fazer para produzir notas de rodapé com referência numérica! Já tentei > de todo jeito, e ainda continua com "fnsymbol"! Qualquer sugestão é bem > vinda! > > abraço, > FH > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > -- Successful people ask better questions, and as a result, they get better answers. (Tony Robbins)

2 6

De colunas para linhas
by Roberto Pinho 24 Mai '11

24 Mai '11

Caros, Uma dúvida básica: Eu tenho um data-frame no formato: $ SIGLA: chr $ X2000: num $ X2001: num $ X2002: num $ X2003: num $ X2004: num $ X2005: num $ X2006: num $ X2007: num $ X2008: num $ X2009: num Eu quero colocar os dados no formato : $ SIGLA: chr $ Ano: num ou chr $ Valor: num Desde já obrigado. -- Roberto Pinho robertodepinho(a)gmail.com msn: robertopinho(a)yahoo.com.br Skype: roberto.pinho http://www.ascoisas.com

3 4

OFF topic
by Fernando Henrique Toledo 24 Mai '11

24 Mai '11

Senhores, boa tarde... Totalmente fora de tópico, mas vocês são minha ultima alternativa a se recorrer! Como fazer para produzir notas de rodapé com referência numérica! Já tentei de todo jeito, e ainda continua com "fnsymbol"! Qualquer sugestão é bem vinda! abraço, FH

2 2

ajuda teste hipótese duas variáveis
by carla eloize carducci 24 Mai '11

24 Mai '11

OI PESSOAL sou uma aprendiz bem iniciante do R, preciso de ajuda com um teste de hipótese, alguém pode me ajudar? bom, é o seguinte: tenho um conj de dados com duas variáveis, quero fazer um t.test para verificar se as médias dos salarios (variável numerica) de funcionários com diferente grau de instrução (2ºgrau e superior, variável categorica) são diferentes ou não. tentei fazer o seguinte: meu arquivo é da apostila do prof. justiniano o velho arquivo 'milsa'.(bem de iniciante) # teste de hipótese: H0:salario 2grau = superior# formulei essa hipotese # H1:salario 2grau diferente superior milsa$instrucao <- factor(milsa$instrucao, label = c("2oGrau", "Superior"), lev = 1:2, ord = T) var.test(milsa$salario~milsa$instrucao) alpha <- 0.05 # especifica valor da signific^ancia if (var.test$p.value > alpha) diff.test <- t.test(milsa$salario ~ milsa$instrucao,alternative = "two.sided", var.equal = TRUE, conf.level = 1-alpha) else diff.test <- t.test(milsa$salario ~ milsa$instrucao, alternative = "two.sided", var.equal = FALSE,conf.level = 1-alpha) diff.test não deu certo, o que estou esquecendo? não acho o erro. muito obrigada... Atenciosamente, Carla Eloize Carducci Engº Agronoma (UFGD) Msc. Ciência do Solo (UFLA) Doutoranda em Ciência do Solo (UFLA)

3 2

Re: [R-br] Problemas ao formatar datas!!
by Ivan Bezerra Allaman 24 Mai '11

24 Mai '11

Valeu Dallas e Daniel! Fabrício, esta e várias dicas no:http://ridiculas.wordpress.com/ Esta foi obra do Walmes. Abraços. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

2 1

Re: [R-br] Problemas ao formatar datas!!
by Ivan Bezerra Allaman 24 Mai '11

24 Mai '11

Tudo bem jovem, vamos facilitar as coisas para ver se há solução. Inci <- read.table("dadosRegInc.txt",header=TRUE) # passamos os dados para o formato longo, criou uma coluna para mês. Inci <- melt(Inci, id=c("Trat", "bloco")) names(Inci) <- c("Trat","bloco","Data","Incidencia") O problema ainda persiste. Alguma solução? Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

3 3

Problemas ao formatar datas!!
by Ivan Bezerra Allaman 24 Mai '11

24 Mai '11

Bom dia senhores!! Tenho o seguinte cabeçalho na base de dados original. Tratbloco03-0804-0805-08 ... .... .... Quando chamo no console do R, este cabeçalho aparece da seguinte forma: Trat bloco X03.08 X04.08 X05.08 ... ... Como faço para que as datas fiquem no formato original da base de dados? Desde já agradeço a atenção. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

3 2

Re: [R-br] Instalação do pacote rparallel
by Ivan Bezerra Allaman 24 Mai '11

24 Mai '11

Me parece que este pacote foi deletado dos repositórios. Entre no cran para verificar que este pacote não está mais lá e para procurar por pacotes que podem lhe auxiliar. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Instalação do pacote rparallel
by Junior Beleti 24 Mai '11

24 Mai '11

Boa noite. Estou tentando utilizar o pacote rparallel, mas estou com o seguinte problema: Após baixar o pacote, realizar a instalação (R CMD INSTALL ...), ao tentar carregar o pacote "library(rparallel)", recebo a seguinte mensagem: Erro: pacote 'rparallel' foi compilado antes do R 2.10.0: por favor instale-o novamente Ou por meio de "require(rparallel)", recebo a mensagem: Carregando pacotes exigidos: rparallel Failed with error: ‘pacote 'rparallel' foi compilado antes do R 2.10.0: por favor instale-o novamente’ Já tentei re-instalá-lo mas ocorre o mesmo problema. Alguém poderia me ajudar? Obrigado. -- Carlos Roberto Beleti Junior Mestrado em Ciência da Computação Departamento de Informática Universidade Estadual de Maringá

2 1

Teste de Dunnet
by Eduardo Vargas 23 Mai '11

23 Mai '11

Será que alguem poderia me informar os comandos para eu usar o teste de Dunnet? obrigado galera -- José Eduardo Vargas Lopes de Aráujo Engenheiro Agrônomo Contato:(38)8811-6985 duvargas(a)gmail.com

2 4

pacote snow windows 64
by Fernando Henrique Toledo 23 Mai '11

23 Mai '11

Senhores, Alguém aí instalou ou sabe como fugir da instalação do pacote 'snow' no windows? Estava tentando instalar outro pacote 'Geneland' e ele acusa a dependência do snow! No Linux, isso foi facilmente resolvido baixando o 'snow' pelo aptitude, mas não consegui fazer isso no windows! grato pela ajuda, FH

2 5

Re: [R-br] Modelo Animal - Pedigreemm
by Ivan Bezerra Allaman 23 Mai '11

23 Mai '11

Olhe novamente para o seu exemplo. O fator 'Animal' está com seis observações!! Deveria estar como mesmo número de observações das demais fontes de variação e variáveis, ou seja, cinco observações. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

4 6

SensoMineR
by Peterson Gardin 20 Mai '11

20 Mai '11

Pessoal, comecei a trabalhar com análise sensorial de alimentos e vi que existe o pacote SensoMiner e que pode ser integrado ao excel, programa que uso muito. Gostaria de saber se me aconselham usar este pacote, se devo me dedicar ao estudo do mesmo. Sei que tem gente na lista que conhece, se puder me dar uma dica. Abraço João Peterson Pereira Gardin Pesquisador Epagri - Videira - SC Doutor em Fisiologia Vegetal Engenheiro Agrônomo Fones: (49) 3566-0054

2 2

[Dúvida] Arquivos .shp para construção de mapas
by Pedro Rafael 20 Mai '11

20 Mai '11

Pessoal onde consigo mapas .shp (Shapefile) por municípios do Estado da Paraíba? Existe malhas gratuitas por estados dividido pelos seus respectivos municípios? Saudações, Pedro Rafael Diniz Marinho

3 2

divisão de banco de dados
by Rodrigo Pereira 20 Mai '11

20 Mai '11

Caros, colegas, bom dia! Estou com a linha de programa abaixo para dividir um banco de dados de 720 linhas em dois de 360, porém ao pedir para armazenar os valores em cada um deles, apenas a última condição satisfeita é armazenada e não todas. Já tentei colocar outros contadores como j e k para os objetos criados, mas não estou conseguindo. Alguém poderia me ajudar? reftestR<-NULL reftestT<-NULL for (i in 1:length(medidas$formulation)) { if (medidas$formulation[i]=="R") { reftestR<-data.frame(medidas$concentration[i],medidas$formulation[i],medidas$subject[i]) } else if (medidas$formulation[i]=="T") { reftestT<-data.frame(medidas$concentration[i],medidas$formulation[i],medidas$subject[i]) } } Muito obrigado! Abraço! Rodrigo Pereira

3 4

[Dúvida] Conglomerado Espacial (Método Skater)
by Pedro Rafael 20 Mai '11

20 Mai '11

Caros amigos, gostaria de saber se posso utilizar o Método Skater para detecção de conglomerados espaciais via arvore geradora mínima em que tenho apenas uma variável de interessa, a saber taxa de mortalidade infantil. Vi que muitas pessoas utilizam essa metodologia com o uso de covariáveis contudo o interesse está na variável resposta. Fui ler mais sobre a metodologia e o método é bem determinístico e não faz uso de distribuição de probabilidade e na minha opinião acredito que o método pode ser utilizado com a penas uma variável. Alguém domina esta metodologia? Pedro

1 0

Modelo Animal - Pedigreemm
by Marcelo Cardoso mello 20 Mai '11

20 Mai '11

Não estou conseguindo reproduzir o pequeno exemplo abaixo, usando o modelo animal. Aparece o erro abaixo. Alguém poderia me dar uma luz?? É possível predizer valores ao nível de indivíduos (Animal) neste exemplo, usando o R (pedigreemm)?? O que esta errado?? Herd<-as.factor(c(1,2,1,2,1)) Animal<-as.factor(c(1,2,3,4,5,6)) Sire<-c(0,0,0,1,4) Dam<-c(0,0,0,2,3) Phenotyp<-c(100,130,120,110,140) require(pedigreemm) p<- pedigree(sire = Sire, dam = Dam, label= 1:5) model<-pedigreemm(Phenotyp~ (1|Animal)+ Herd, pedigree=list(Animal=p)) Erro em function (fr, FL, start, REML, verbose) : Number of levels of a grouping factor for the random effects must be less than the number of observations Herd Animal Sire Dam Phenotyp 1 1 0 0 100 2 2 0 0 130 1 3 0 0 120 2 4 1 2 110 1 5 4 3 140

1 0

equação hiperbole não rectangular
by Gilson Sanchez 19 Mai '11

19 Mai '11

Pessoal, especialmente se alguem trabalha com fisiologia vegetal, então coleite dados de fotossíntese, e agora estou querendo saber se da para estimar pelo R o ajuste de curva luz estimada pela equação da hiperbole não rectangular. Tenho intensidades de luz de 0, 25, 50, 75, 100, 500, 1000, 1500 e 2000 PPDF, então eu quero utilizar esta equação pra fazer o ajuste da curva de luz. O fato é que me falaram que só o software Statistca sabe estimar este ajuste de curva. Se alguem souver alguma coisa que me possa ajudar agradeceria muito. -- MSc. Gilson Sánchez Chia Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 --- Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente. Please consider the environment before printing this email. Renew, Reduce and Recycle. The Planet Knows.

4 3

Conflito de resultados + plotagem
by Rodrigo Plei 19 Mai '11

19 Mai '11

Prezados, Utilizei a função "nls" do R para estimar as constantes de uma função potencial ajustada aos meus dados de comprimento e peso de lulas utilizando o método dos mínimos quadrados: m <- nls(y~a*x^b, data = ds, start = c(a=.1, b=2)) Obtive as seguintes constantes: a = 0.0003967 , b = 2.3933797 e o R2 = 0.9681 Ocorre que o mesmo ajuste no Excel 2007 resulta em a = 0.0002 , b = 2.6097 e R2 = 0.9819 O mesmo ajuste no Statistica resulta em constantes e R2 idênticas aos encontrados no R. Não deveria ter dado o mesmo resultado nos três programas? Não sei em qual resultado confiar. Também gostaria de saber como faço para plotar a curva dessa função potencial sobreposta a minha dispersão de dados (peso vs. comprimeto). Agradeço desde logo. Rodrigo -- =8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8--> Rodrigo Silvestre Martins, PhD Bolsista FAPESP Pós-Doutorado Laboratório de Ecossistemas Pesqueiros (LabPesq) Universidade de São Paulo, Instituto Oceanográfico Praça do Oceanográfico, 191. Cidade Universitaria (sala 107-A/B) Butantã - São Paulo/SP, Brasil 05508-900 Tel: +55 11 3091 6549 Email: rodrigo.plei(a)gmail.com ; ocersm(a)lycos.com; rsmartins(a)usp.br CV Lattes: http://lattes.cnpq.br/5350064124902777

4 4

R studio - linux
by Marcelo Claro de Souza 19 Mai '11

19 Mai '11

Eu costumo abrir um terminal dentro da pasta onde estão meus arquivos e utilizo o comando write.csv2('nome do arquivo.csv') para salvar minhas planilhas dentro da própria pasta (não sei salvar os arquivos fora da pasta onde abro o terminal). Entretanto, ao atualizar o ubuntu para a versão 11.04 esse comando não está funcionando, o R roda o comando mas não salva o arquivo. Gostaria de saber como possa fazer isso pelo R studio como uma alternativa, ou alguma outra solução p o R clássico. Muito obrigado.

3 2

fatorial fracionado em R
by jonas bodini alonso 19 Mai '11

19 Mai '11

Bom dia aos membros da lista Alguem sabe uma função do R que trabalhe com experimentos fatoriais fracionados? Se sim como faço para fazer para escolher a fração de interesse? Cordialmente Jonas

2 2

Desdobramento da Interação! Função split()
by Thiago De paula protásio 18 Mai '11

18 Mai '11

Pessoal; Recorro novamente a vcs para pedir um auxílio. Estou querendo utilizar a função split para desdobrar o efeito da interação entre fatores.Contudo, não estou obtendo sucesso. Meu experimento é um fatorial completo 3 x 2 instalado em esquema fatorial, ou seja, tenho 3 clones plantados em 2 locais distintos. Alguém já trabalhou com essa função ou me aconselha a utilizar outra?Posso enviar alguns dados para elucidar o problema. Muito obrigado!! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

2 1

Erro ao executar arquivo Rpatch.sh!
by Ivan Bezerra Allaman 18 Mai '11

18 Mai '11

Bom dia senhores!! Estou tentando atualizar a versão do R por meio do arquivo Rpatch.sh que o professor PJ disponibilizou em um tópico anterior. Está ocorrendo o seguinte erro: ivan-Aspire-5920 Área de Trabalho # chmod 777 Rpatch.sh ivan-Aspire-5920 Área de Trabalho # ./Rpatch.sh ./Rpatch.sh: linha 1: -------------------------------------------------------: comando não encontrado ============================================= ============================================= ARQUIVOS FONTE R - VERSAO 2-13 ============================================= ============================================= ============================================= ============================================= ARQUIVOS FONTES R - VERSAO 2-13 ============================================= ============================================= Nao instalando dependencias, o R ou alguns pacotes podem nao compilar dependendo do que voce tem instalado no seu sistema. Para instalar dependencias automaticamente utilize opcao --deps Veja opcoes com --help ============================================= ============================================= ARQUIVOS FONTE R-VERSAO 2-13 ============================================= ============================================= ./Rpatch.sh: linha 70: svn: comando não encontrado R version 2.11.1 (2010-05-31) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R é um software livre e vem sem GARANTIA ALGUMA. Você pode redistribuí-lo sob certas circunstâncias. Digite 'license()' ou 'licence()' para detalhes de distribuição. R é um projeto colaborativo com muitos contribuidores. Digite 'contributors()' para obter mais informações e 'citation()' para saber como citar o R ou pacotes do R em publicações. Digite 'demo()' para demonstrações, 'help()' para o sistema on-line de ajuda, ou 'help.start()' para abrir o sistema de ajuda em HTML no seu navegador. Digite 'q()' para sair do R. Vejam que não atualizou! Eu copiei exatamente o que estava no tópico. Se alguém puder me ajudar eu ficarei grato. Abraços. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

3 3

Leitura de arquivos
by Emerson Cotta Bodevan 18 Mai '11

18 Mai '11

Prezados. Tenho um arquivo texto com 100 linhas, sendo que as 50 primeiras possuem três colunas e as 50 últimas, quatro colunas. Como faço para ler esse arquivo, separando-o em dois data-frames (Um com 3 colunas e o outro com 4 colunas)? Agradeço a ajuda, -- Emerson

3 4

combinar tipos diferentes de gráficos
by Rodrigo Plei 18 Mai '11

18 Mai '11

Prezados, Como combinar um gráfico de barras com um gráfico de dispersão? Ambas plotagens possuem os mesmos dados para o eixo x. Gostaria que as barras ficassem atreladas ao eixo y principal (esquerdo) e os pontos da dispersão atrelados ao eixo y secundário (direito). Agradeço desde logo. Rodrigo -- =8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8--> Rodrigo Silvestre Martins, PhD Bolsista FAPESP Pós-Doutorado Laboratório de Ecossistemas Pesqueiros (LabPesq) Universidade de São Paulo, Instituto Oceanográfico Praça do Oceanográfico, 191. Cidade Universitaria (sala 107-A/B) Butantã - São Paulo/SP, Brasil 05508-900 Tel: +55 11 3091 6549 Email: rodrigo.plei(a)gmail.com ; ocersm(a)lycos.com; rsmartins(a)usp.br CV Lattes: http://lattes.cnpq.br/5350064124902777

2 2

Re: [R-br] Erro ao executar arquivo Rpatch.sh!
by Ivan Bezerra Allaman 18 Mai '11

18 Mai '11

Resolvi do jeito tradicional. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Usando AIC para seleção de modelo com o comando drop1(). Alguém poderia dar uma olhadinha :)
by Augusto Ribas 18 Mai '11

18 Mai '11

Ola colegas de lista. Eu estava lendo sobre seleção de modelos e fiquei na duvida se estou indo pelo caminho certo no R. Outra coisa, eu não achei exatamente nos livros que li o que eu deveria reportar em um artigo, ao longo destes passos (se estes estiverem fazendo sentido). Eu tentei fazer um exemplo com uma com uma uma variavel resposta e duas variaveis explicativas. Usei o dataset do anscombe que vem no R, usando o primeiro exemplo dele mais uma variavel ao acaso para exemplificar. Aqui vai um codigo de como estava pensando: #Abrindo os dados anscombe x2<-runif(11) #fazendo um plot para olhar panel.hist <- function(x, ...) { usr <- par("usr"); on.exit(par(usr)) par(usr = c(usr[1:2], 0, 1.5) ) h <- hist(x, plot = FALSE) breaks <- h$breaks; nB <- length(breaks) y <- h$counts; y <- y/max(y) rect(breaks[-nB], 0, breaks[-1], y, col="cyan", ...) } panel.cor <- function(x, y, digits=2, prefix="", cex.cor, ...) { usr <- par("usr"); on.exit(par(usr)) par(usr = c(0, 1, 0, 1)) r <- abs(cor(x, y)) txt <- format(c(r, 0.123456789), digits=digits)[1] txt <- paste(prefix, txt, sep="") if(missing(cex.cor)) cex.cor <- 0.8/strwidth(txt) text(0.5, 0.5, txt, cex = cex.cor * r+1) } panel.lm = function(x, y, ...) { tmp <- lm(y ~ x, na.action = na.omit) abline(tmp) points(x, y, ...)} #legal, y1 é relacionado com x1 mas nao com x2 olhando esse grafico pairs(cbind(y1=anscombe$y1,x1=anscombe$x1,x2),diag.panel=panel.hist, lower.panel=panel.cor,upper.panel=panel.lm) #agora fazendo um modelo modelo1<-lm(anscombe$y1~anscombe$x1+x2) modelo1 #vendo premissas par(mfrow=c(2,2)) plot(modelo1) #vendo o quanto os parametros ajudam drop1(modelo1) #testando quem ajuda significativamente drop1(modelo1,test="F") #novo modelo que eu vou ficar ja que vi no drop1 que posso tirar o x2 modelo2<-lm(anscombe$y1~anscombe$x1) modelo2 #vendo premissas pro novo modelo par(mfrow=c(2,2)) plot(modelo2) #o que reportar num artigo de todo o processo? #todo esse processo pode ser inverso tb, começando de um modelo simples e indo pra um mais complexo add1(modelo2, ~anscombe$x1+x2,test="F") É isso. Agradeço a paciencia de te -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

2 1

coordenadas valor máximo num grafico xy
by jose coelho 18 Mai '11

18 Mai '11

Vivam , ainda sou aprendiz no R... mas gostaria de saber como calcular as coordenada do valor máximo de y em função de x para colocar num gráfico xy. Existe alguma função ? max(y) eu sei qual valor máximo do y mas também quero saber qual o valor de x no qual o y atinge o seu máximo. Muito obrigado jmmcoelho

3 3