Maio 2011 - R-br - DINF UFPR

Converter datas em nº de dias do ano
by Cristiano Santos 13 Jun '11

13 Jun '11

Caros, Sou novo no R e não estou conseguindo manipular datas. Tenho uma variável com data e hora, e quero transformar no número de dias do ano, por exemplo: "4/02/2011 03:24:13" deverá transformar-se em "35" Desde já agradeço pela ajuda. Cristiano Santos Economia-UFC

5 8

Desdobrar modelos de regresão em fatorial duplo (um fator quantitativo)
by Marcelo Cardoso mello 12 Jun '11

12 Jun '11

Estou com dificuldades em desdobrar a interação em uma regressão em fatorial. Segue um exemplo abaixo. Como obter um modelo de regressão de r em função de b para a1 e a2?????? a<-as.factor(c(rep(1,9),rep(2,9))) b<-c(rep(0,3),rep(5,3),rep(10,3), rep(0,3),rep(5,3),rep(10,3)) r<-c(1,1.2,1.3,2,2.3,2.25,4,3.9,3.8,5,4.8,4.7,3.5,3,3.4,1.7,1.5,1.2) b<-as.factor(b) m<-lm(r~a*b) anova(m) b<-as.numeric(b) regressão<-lm(r~b) regressão anova(regressão) Como desdobrar b dentro de a1 e a2????? Obrigado. Marcelo.

4 9

teste de Homogeneidade e Normalidade
by Eduardo Vargas 02 Jun '11

02 Jun '11

Boa noite pessoal,preciso de uma ajuda.Estou fazendo uma comparação de médias de germinação de sementes, e as varianveis são tratamentos(trat) fator e Germinação (germ) numerica.Preciso fazer o teste de homogeneidade de variancias com o Bartlett e normalidade com o teste de shapiro-Wilk, porém dá um mensagem de erro(Erro em bartlett.test(germ, trat) : objeto 'germ' não encontrado) e no Shapiro -wilk essa mensagem Erro em `[.data.frame`(x, complete.cases(x)) : undefined columns selected Será que alguem pode me ajudar? -- José Eduardo Vargas Lopes de Aráujo Engenheiro Agrônomo Contato:(38)8811-6985 duvargas(a)gmail.com

10 24

modifica dataframe
by Samuel Carvalho 01 Jun '11

01 Jun '11

Boa noite caros(a) senhores(as) membros Por gentileza alguém poderia me auxliar em um código que a princípio parece ser simples mas nao consegui uma solução Abaixo mostro o codigo e a seguir venho comentando ### Codigo R dados<-data.frame(idfuste=seq(1,5), dap=c(8.5, 11.0, 10.5, 17.5, 12.5), ht=c(10, 13, 12.5, 18, 14)) dados dados$hi<-seq(0.1, dados$ht, by=0.5) ### O que preciso é gerar uma sequencia para cada "idfuste" de 0.1 até a "ht" que corresponde aquele "idfuste", replicando também cada observação para cada "idfuste". A sequencia é de 0.5 em 0.5 até a altura total da árvore - 1 pois esta pode nao ser um número inteiro dividido por 0.5 sendo que a ultima observação de cada idfuste deve ser igual ao valor de ht Em resumo queria um data.frame de saída a partir dos dados de entrada, nesta estrutura idfuste dap ht hi 1 8,5 10 0 1 8,5 10 0,5 1 8,5 10 1 1 8,5 10 1,5 1 8,5 10 2 1 8,5 10 2,5 1 8,5 10 3 1 8,5 10 3,5 1 8,5 10 4 1 8,5 10 4,5 1 8,5 10 5 1 8,5 10 5,5 1 8,5 10 6 1 8,5 10 6,5 1 8,5 10 7 1 8,5 10 7,5 1 8,5 10 8 1 8,5 10 8,5 1 8,5 10 9 1 8,5 10 9,5 1 8,5 10 10 2 11 13,2 0 2 11 13,2 0,5 2 11 13,2 1 2 11 13,2 1,5 2 11 13,2 2 2 11 13,2 2,5 2 11 13,2 3 2 11 13,2 3,5 2 11 13,2 4 2 11 13,2 4,5 2 11 13,2 5 2 11 13,2 5,5 2 11 13,2 6 2 11 13,2 6,5 2 11 13,2 7 2 11 13,2 7,5 2 11 13,2 8 2 11 13,2 8,5 2 11 13,2 9 2 11 13,2 9,5 2 11 13,2 10 2 11 13,2 10,5 2 11 13,2 11 2 11 13,2 11,5 2 11 13,2 12 2 11 13,2 12,5 2 11 13,2 13 2 11 13,2 13,2 e assim por diante para todos idfustes Desde já agradeço Para quem é da area florestal estou tentando estruturar uma base de cubagem a partir de dados de parcelas Abraço e boa noite a todos Att, Samuel ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] =============================================

6 8

Ordenar gráfico bipartido
by Jéssica Caroline Falcão 01 Jun '11

01 Jun '11

Olá, Eu gostaria de saber como eu ordeno o gráfico bipartido gerado pelo pacote bipartite de acordo com o número de interações das espécies na rede, em ordem decrescente. Att. Jéssica Falcão -- ------------------------------------------------------------------ Jéssica Falcão Mestranda em Ecologia e Conservação da Biodiversidade Universidade Federal de Mato Grosso IB- Departamento de Ecologia e Botânica Lab. Interação Inseto-Planta / Mirmecologia Cuiabá, Mato Grosso, BRASIL Email: jecafalcao(a)gmail.com

2 1

R e LaTex
by Cassiano 31 Mai '11

31 Mai '11

Olá a todos, Gostaria de aproveitar a lista para difulgar um trabalho sobre o Sweave - uma interface entre o R e o LaTex. O trabalho está publicado na The PracTex Journal 2010-2 e pode ser acessado no endereço: http://tug.org/pracjourn/2010-2/toc.html Att, -- ======================================= MSc Cassiano S. Rosa Universidade Federal de Viçosa -UFV Dept. Biologia Animal Pós-graduação em Entomologia Lab. Termitologia Viçosa-MG, Brasil CEP 36570-000 Tel: (31) 38994017 =======================================

2 1

Distribuição Exponencial
by Fernanda Balbino 31 Mai '11

31 Mai '11

Alexandro, Muito obrigada, os comandos dão o resultado que eu precisava sobre a distribuição exponencial. Fernanda. ________________________________ De: "r-br-request(a)listas.c3sl.ufpr.br" <r-br-request(a)listas.c3sl.ufpr.br> Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Terça-feira, 31 de Maio de 2011 12:00:02 Assunto: Digest R-br, volume 3, assunto 36 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. [Dúvida] Plotando pontos centroides no MAPA (Pedro Rafael) 2. Re: Diagrama de voronoi (Wagner bonat) 3. Re: Análise fatorial e regressão em painel (Felipe Emanoel Barletta Mendes) 4. Re: Diagrama de voronoi (Fabrício Tavares) 5. Distribuição Exponencial (Fernanda Balbino) 6. Re: Distribuição Exponencial (Fernando Mayer) 7. Re: Distribuição Exponencial (Alexandro (Yahoo)) 8. nao consigo instalar o pacote "rimage" - Ubuntu (Guilherme Moraes Ferraudo) 9. Re: nao consigo instalar o pacote "rimage" - Ubuntu (Daniel C Bezerra) 10. Re: nao consigo instalar o pacote "rimage" - Ubuntu (Benilton Carvalho) 11. Livro (Peterson Gardin) ---------------------------------------------------------------------- Message: 1 Date: Mon, 30 May 2011 12:00:20 -0300 From: Pedro Rafael <pedro.rafael.marinho(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: [R-br] [Dúvida] Plotando pontos centroides no MAPA Message-ID: <BANLkTik2Nu2xDKh2BqDKUdYX6pa=j_W42Q(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" 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 -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110530/1eef92f9/atta…> ------------------------------ Message: 2 Date: Mon, 30 May 2011 12:18:55 -0300 From: Wagner bonat <wbonat(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Diagrama de voronoi Message-ID: <BANLkTinCo6d5NgvxV3VhbVQn+y7PMFNOBg(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Veja este exemplo, fiz isto a algum tempo deve estar desatualizado require(gpclib) require(tripack) require(sp) require(spdep) ## Lendo o arquivo das bordas bordas <- read.table(" http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/points_higi.csv ",header=TRUE,sep=";",dec=".") bordas = bordas[,-1] bordas[44,] <- bordas[1,] ## Lendo o arquivo que contem as coordenadas mosq = read.table(" http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/mosqhigi.csv ",header=TRUE,sep=";",dec=".") ## Funcoes extras source("http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/voronoi.R") source("http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/grp2sp.r") coords.x <- c() coords.y <- c() cod.mosq <- unique(mosq$mosquitrap) for(i in 1:40){ coords.x[i] = unique(mosq[which(mosq$mosquitrap == cod.mosq[i]),]$x) coords.y[i] = unique(mosq[which(mosq$mosquitrap == cod.mosq[i]),]$y)} mapa.mosq <- voronoi(coords.x1=coords.x,coords.x2=coords.y,poligono = bordas) mapa.mosq = grp2sp(mapa.mosq,ID = cod.mosq) mapa.mosq.nb <- poly2nb(mapa.mosq) plot(mapa.mosq) -- Wagner Hugo Bonat LEG - Laboratório de Estatística e Geoinformação UFPR - Universidade Federal do Paraná -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110530/13fe55d8/atta…> ------------------------------ Message: 3 Date: Mon, 30 May 2011 13:13:07 -0300 From: "Felipe Emanoel Barletta Mendes" <felipe(a)leg.ufpr.br> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Análise fatorial e regressão em painel Message-ID: <500a2d6b5bb74891c63e982af2c8d472.squirrel(a)leg.ufpr.br> Content-Type: text/plain;charset=iso-8859-1 Henrique, há o pacote plm para dados em painel. ?plm Henrique Brilhante > 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 > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Felipe E. Barletta Mendes (41)9189-5198 (41)3025-2150 (41)3328-7216 http://www.leg.ufpr.br/doku.php/pessoais:felipe ------------------------------ Message: 4 Date: Mon, 30 May 2011 17:28:00 -0300 From: Fabrício Tavares <fabriciotavareso(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Diagrama de voronoi Message-ID: <BANLkTin-nDpN30Z1vbkjZpW7VDhFfSk+aQ(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Wagner, funcionou perfeitamente! Muito obrigado! Fabrício. http://vizzual.net Em 30 de maio de 2011 12:18, Wagner bonat <wbonat(a)gmail.com> escreveu: > Veja este exemplo, fiz isto a algum tempo deve estar desatualizado > > require(gpclib) > require(tripack) > require(sp) > require(spdep) > ## Lendo o arquivo das bordas > bordas <- read.table(" > http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/points_higi.csv > ",header=TRUE,sep=";",dec=".") > bordas = bordas[,-1] > bordas[44,] <- bordas[1,] > > ## Lendo o arquivo que contem as coordenadas > mosq = read.table(" > http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/mosqhigi.csv > ",header=TRUE,sep=";",dec=".") > > ## Funcoes extras > source("http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/voronoi.R") > source("http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/grp2sp.r") > > coords.x <- c() > coords.y <- c() > cod.mosq <- unique(mosq$mosquitrap) > for(i in 1:40){ > coords.x[i] = unique(mosq[which(mosq$mosquitrap == cod.mosq[i]),]$x) > coords.y[i] = unique(mosq[which(mosq$mosquitrap == cod.mosq[i]),]$y)} > > mapa.mosq <- voronoi(coords.x1=coords.x,coords.x2=coords.y,poligono = > bordas) > mapa.mosq = grp2sp(mapa.mosq,ID = cod.mosq) > mapa.mosq.nb <- poly2nb(mapa.mosq) > plot(mapa.mosq) > > -- > Wagner Hugo Bonat > LEG - Laboratório de Estatística e Geoinformação > UFPR - Universidade Federal do Paraná > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110530/644eef83/atta…> ------------------------------ Message: 5 Date: Mon, 30 May 2011 13:51:16 -0700 (PDT) From: Fernanda Balbino <ferbalbino(a)yahoo.com.br> To: r-br(a)listas.c3sl.ufpr.br Subject: [R-br] Distribuição Exponencial Message-ID: <25518.42840.qm(a)web161501.mail.bf1.yahoo.com> Content-Type: text/plain; charset="iso-8859-1" Olá, boa tarde a todos, Tenho um conjunto de dados que devido a natureza e aos gráficos de análise exploratória acredito que sigam uma distribuição exponencial. Gostaria de saber se no R é possível testar os dados para confirmar se realmente eles se comportam conforme uma distribuição exponencial. Muito obrigada. Fernanda Balbino. -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110530/21cbd860/atta…> ------------------------------ Message: 6 Date: Mon, 30 May 2011 18:28:14 -0300 From: Fernando Mayer <fernandomayer(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Distribuição Exponencial Message-ID: <BANLkTi=ofAuy684SbjuU1doRt1z9nyY3Hw(a)mail.gmail.com> Content-Type: text/plain; charset=ISO-8859-1 Veja a função fitdistr() do pacote MASS [], --- Fernando Mayer URL: http://sites.google.com/site/fernandomayer e-mail: fernandomayer [@] gmail.com 2011/5/30 Fernanda Balbino <ferbalbino(a)yahoo.com.br>: > Olá, boa tarde a todos, > Tenho um conjunto de dados que devido a natureza e aos gráficos de análise > exploratória acredito que sigam uma distribuição exponencial. > Gostaria de saber se no R é possível testar os dados para confirmar se > realmente eles se comportam conforme uma distribuição exponencial. > Muito obrigada. > Fernanda Balbino. > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > ------------------------------ Message: 7 Date: Mon, 30 May 2011 16:11:06 -0700 (PDT) From: "Alexandro $Yahoo$" <alexandro2x(a)yahoo.com.br> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Distribuição Exponencial Message-ID: <499551.91899.qm(a)web130108.mail.mud.yahoo.com> Content-Type: text/plain; charset="iso-8859-1" Boa noite Fernanda, seja x, a variável que tenha os seus dados, tu podes testar o seguinte: #EXEMPLO DE NÃO TER DIST. EXP #x=rnorm(1000,10,2) #ks.test(x, pexp, rate=5) #pvalor=ks.test(x, pexp, rate=5)$p.value #EXEMPLO DE DIST. EXP x <- rexp(30, rate=5) ks.test(x, pexp, rate=5) pvalor=ks.test(x, pexp, rate=5)$p.value pvalor alfa=0.05 teste<-function(pvalor) { if(pvalor<alfa) { cat("rejeita Ho >>> pvalor é ",pvalor, " NÃO TEM DIST. EXP ","\n") } else if(pvalor>=alfa) { cat("NÃO rejeita Ho >>> pvalor é ",pvalor, " TEM DIST. EXP. ","\n") } } teste(pvalor) Espero ter ajudado ! Alexandro >________________________________ >De: Fernanda Balbino <ferbalbino(a)yahoo.com.br> >Para: r-br(a)listas.c3sl.ufpr.br >Enviadas: Segunda-feira, 30 de Maio de 2011 17:51 >Assunto: [R-br] Distribuição Exponencial > > >Olá, boa tarde a todos, > > >Tenho um conjunto de dados que devido a natureza e aos gráficos de análise >exploratória acredito que sigam uma distribuição exponencial. >Gostaria de saber se no R é possível testar os dados para confirmar se realmente >eles se comportam conforme uma distribuição exponencial. > > >Muito obrigada. >Fernanda Balbino. > > > > >_______________________________________________ >R-br mailing list >R-br(a)listas.c3sl.ufpr.br >https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110530/c7b1fdad/atta…> ------------------------------ Message: 8 Date: Mon, 30 May 2011 21:58:02 -0300 From: Guilherme Moraes Ferraudo <guiferraudo(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: [R-br] nao consigo instalar o pacote "rimage" - Ubuntu Message-ID: <BANLkTim9AKvn2zsX_t=oeBffhX=0k+beKQ(a)mail.gmail.com> Content-Type: text/plain; charset=windows-1252 Boa noite pessoal, estou no Ubuntu e não consigo instalar o pacote "rimage" - Ubuntu. Olha o erro: ERROR: configuration failed for package ?rimage? * removing ?/home/gferrau/R/i486-pc-linux-gnu-library/2.10/rimage? The downloaded packages are in ?/tmp/RtmpcEgQZS/downloaded_packages? Warning message: In install.packages("rimage", dep = T) : installation of package 'rimage' had non-zero exit status O que fazer?? Abs -- Guilherme Moraes Ferraudo Jaboticabal/Campinas - SP http://lattes.cnpq.br/2096118558794430 ------------------------------ Message: 9 Date: Mon, 30 May 2011 22:10:42 -0300 From: Daniel C Bezerra <danielcbezerra(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] nao consigo instalar o pacote "rimage" - Ubuntu Message-ID: <BANLkTikZ=cNKNEQN0AH_0WbXXCQgbWd1PA(a)mail.gmail.com> Content-Type: text/plain; charset="windows-1252" Parece que vc está usando a versão 2.10 do R. Tente atualizar primeiro. Abs, D 2011/5/30 Guilherme Moraes Ferraudo <guiferraudo(a)gmail.com> > Boa noite pessoal, > > estou no Ubuntu e não consigo instalar o pacote "rimage" - Ubuntu. > > Olha o erro: > > ERROR: configuration failed for package ?rimage? > * removing ?/home/gferrau/R/i486-pc-linux-gnu-library/2.10/rimage? > > The downloaded packages are in > ?/tmp/RtmpcEgQZS/downloaded_packages? > Warning message: > In install.packages("rimage", dep = T) : > installation of package 'rimage' had non-zero exit status > > > O que fazer?? > > Abs > > > -- > Guilherme Moraes Ferraudo > Jaboticabal/Campinas - SP > http://lattes.cnpq.br/2096118558794430 > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110530/51828ce7/atta…> ------------------------------ Message: 10 Date: Tue, 31 May 2011 10:43:51 +0100 From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] nao consigo instalar o pacote "rimage" - Ubuntu Message-ID: <BANLkTinyw7kdkhY-6ZSFNuJ8jk7cndNfLw(a)mail.gmail.com> Content-Type: text/plain; charset=windows-1252 Adicionalmente, a mensagem de erro parece incompleta. 2011/5/31 Daniel C Bezerra <danielcbezerra(a)gmail.com>: > Parece que vc está usando a versão 2.10 do R. Tente atualizar primeiro. > > Abs, > > D > > 2011/5/30 Guilherme Moraes Ferraudo <guiferraudo(a)gmail.com> >> >> Boa noite pessoal, >> >> estou no Ubuntu e não consigo instalar o pacote "rimage" - Ubuntu. >> >> Olha o erro: >> >> ERROR: configuration failed for package ?rimage? >> * removing ?/home/gferrau/R/i486-pc-linux-gnu-library/2.10/rimage? >> >> The downloaded packages are in >> ?/tmp/RtmpcEgQZS/downloaded_packages? >> Warning message: >> In install.packages("rimage", dep = T) : >> installation of package 'rimage' had non-zero exit status >> >> >> O que fazer?? >> >> Abs >> >> >> -- >> Guilherme Moraes Ferraudo >> Jaboticabal/Campinas - SP >> http://lattes.cnpq.br/2096118558794430 >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > _______________________________________________ > 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) ------------------------------ Message: 11 Date: Tue, 31 May 2011 07:48:21 -0700 (PDT) From: Peterson Gardin <petersongardin(a)yahoo.com.br> To: r-br(a)listas.c3sl.ufpr.br Subject: [R-br] Livro Message-ID: <596446.76610.qm(a)web125801.mail.ne1.yahoo.com> Content-Type: text/plain; charset="iso-8859-1" Precisava deste livro, alguém tem ? "R Through Excel" Obrigado João Peterson Pereira Gardin Pesquisador Epagri - Videira - SC Doutor em Fisiologia Vegetal Engenheiro Agrônomo Fones: (49) 3566-0054 -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/9f3ca333/atta…> ------------------------------ _______________________________________________ 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 3, assunto 36 ****************************************

1 0

Livro
by Peterson Gardin 31 Mai '11

31 Mai '11

Precisava deste livro, alguém tem ? "R Through Excel" Obrigado João Peterson Pereira Gardin Pesquisador Epagri - Videira - SC Doutor em Fisiologia Vegetal Engenheiro Agrônomo Fones: (49) 3566-0054

4 3

Alterar a função codebook
by Daniel Marcelino 31 Mai '11

31 Mai '11

Caros, esta função codebook tem uma saída para valores characters e factors diferente do que eu preciso para montar um livro de códigos. No código abaixo, quando a coluna é string a única informação que mostra é: "A character vector". O que eu preciso é que a saída me mostre pelo menos uns 3 exemplos de strings. Alguém sabe como alterar isso? Daniel function (dataFrame = .data) { cat("\n", attr(dataFrame, "datalabel"), "\n", "\n") x1 <- dataFrame[1, ] for (i in 1:ncol(dataFrame)) { cat(paste(names(dataFrame)[i], "\t", ":", "\t", attr(dataFrame, "var.labels")[i]), "\n") if (all(is.na(dataFrame[, i]))) { cat(paste("All elements of ", names(dataFrame)[i], " have a missing value", "\n")) } else { if (any(class(x1) == "data.frame")) { x2 <- x1[, i] } else { x2 <- x1 } if (any(class(x2) == "character") | any(class(x2) == "AsIs")) { cat("A character vector", "\n") } else { if (any(class(x2) == "difftime")) { print(summary(x2)) } else { if (is.logical(x2)) x2 <- as.factor(x2) if (any(class(x2) == "factor")) { table1 <- (t(t(table(dataFrame[, i])))) table1 <- cbind(table1, format(table1/sum(table1) * 100, digits = 3)) colnames(table1) <- c(.frequency1, .percent) if (is.null(attr(dataFrame, "val.labels")[i])) { print.noquote(table1, right = TRUE) } else { if (any(is.na(attr(dataFrame, "label.table")))) { print.noquote(table1, right = TRUE) } else { attr(dataFrame, "label.table")[which(is.na(attr(attr(dataFrame, "label.table"), "names")))] <- "" index <- attr(attr(dataFrame, "label.table"), "names") == attr(dataFrame, "val.labels")[i] index <- na.omit(index) if (suppressWarnings(!all(rownames(as.data.frame(attr(dataFrame, "label.table")[index])) == levels(x2)))) { print.noquote(table1, right = TRUE) } else { table2 <- data.frame(attr(dataFrame, "label.table")[index], table1) colnames(table2) <- c("code", colnames(table1)) cat("Label table:", attr(dataFrame, "val.labels")[i], "\n") print.noquote(table2, right = TRUE) } } } } else { print(summ(dataFrame[, i], graph = FALSE)) } } } } cat("\n", "==================", "\n") } } -- Daniel Marcelino http://danielmarcelino.zip.net Skype: d_marcelino

1 0

RCurl - Ubuntu 11.04
by Fabrício Tavares 31 Mai '11

31 Mai '11

Boa tarde pessoal. Estou com dificuldades pra instalar o pacote RCurl no Ubuntu 11.04 (uso há pouco tempo). > sessionInfo() R version 2.13.0 (2011-04-13) Platform: x86_64-pc-linux-gnu (64-bit) locale: [1] LC_CTYPE=pt_BR.UTF-8 LC_NUMERIC=C [3] LC_TIME=pt_BR.UTF-8 LC_COLLATE=pt_BR.UTF-8 [5] LC_MONETARY=C LC_MESSAGES=en_US.UTF-8 [7] LC_PAPER=pt_BR.UTF-8 LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C [11] LC_MEASUREMENT=pt_BR.UTF-8 LC_IDENTIFICATION=C attached base packages: [1] stats graphics grDevices utils datasets methods base loaded via a namespace (and not attached): [1] tcltk_2.13.0 tools_2.13.0 > install.packages('RCurl') Installing package(s) into ‘/home/fabricio/R/x86_64-pc-linux-gnu-library/2.13’ (as ‘lib’ is unspecified) --- Please select a CRAN mirror for use in this session --- Loading Tcl/Tk interface ... done trying URL 'http://cran-r.c3sl.ufpr.br/src/contrib/RCurl_1.6-4.tar.gz' Content type 'application/x-gzip' length 703548 bytes (687 Kb) opened URL ================================================== downloaded 687 Kb * installing *source* package ‘RCurl’ ... checking for curl-config... no Cannot find curl-config ERROR: configuration failed for package ‘RCurl’ * removing ‘/home/fabricio/R/x86_64-pc-linux-gnu-library/2.13/RCurl’ The downloaded packages are in ‘/tmp/Rtmp8Pq7iI/downloaded_packages’ Warning message: In install.packages("RCurl") : installation of package 'RCurl' had non-zero exit status --- O que é esse "curl-config"? Estou deixando de fazer algo? (o curl 7.21 está instalado) Abs, Fabrício. http://vizzual.net

2 2