- R-br - DINF UFPR

ajuda com gráfico
by D'Apolito 07 Dez '12

07 Dez '12

Desculpem não saber o nome desse tipo de gráfico para colocar no título do e-mail. Alguém sabe uma função/pacote que faça esse tipo de gráfico (os da parte superior)?? http://ars.els-cdn.com/content/image/1-s2.0-S0367253011000703-gr2.jpg grato Carlos

8 8

Erro: Intersecção entre dois shapefiles no R
by Ricardo Dobrovolski 07 Dez '12

07 Dez '12

Oi Pessoal, Tentei realizar a intersecção entre dois shapes por meio da função gIntersect do pacote rgeos. No entanto, ocorreu o seguinte erro, conforme mostrado abaixo: > wdpa<- readShapePoly(file.choose()) > grid_05<- readShapePoly(file.choose()) > inter_wdpa_grid<-gIntersection(wdpa,grid_05) Erro em createPolygonsComment(p) : rgeos_PolyCreateComment: orphaned hole, cannot find containing polygon for hole at index 2 Em 6 de dezembro de 2012 20:22, Ricardo Dobrovolski <rdobrovolski(a)gmail.com>escreveu: > Oi Pessoal, > > Eu gostaria de saber se vocês conhecem pacotes no R nos quais eu possa > fazer operações entre shapefiles como intersecções, clips, etc, de maneira > semelhante às operações que se podem realizar no ArcGis. > > Atenciosamente, > > Ricardo Dobrovolski >

2 1

Tinn-R 2.4.1.5
by Jose Claudio Faria 07 Dez '12

07 Dez '12

Caros usuários do Tinn-R, Uma nova versão estável do editor Tinn-R está disponível no servidor SourceForge: http://sourceforge.net/projects/tinn-r/ As novidades: http://sourceforge.net/p/tinn-r/news/?source=navbar Ab, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)9100.7351 - TIM 55(73)8817.6159 - OI ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

1 0

Derivar equações com R
by Fernando Antonio de souza 07 Dez '12

07 Dez '12

Caros colegas Gerei uma equação da forma = log10(PB)~a+b*log10(PCVZ) o qual relaciona o conteúdo de proteína corporal animal em função de seu peso corporal vazio. Segundo a metodologia que estou seguindo esta equação deve ser derivada para se estimar o conteúdo de PB por Kg de PCVZ. Há uma forma de fazer esta derivada? Alguém poderia me ensinar como fazer? grato pela atenção. modelo5<-lm(log10(PB)~log10(PCVZ),data=PB) summary(modelo5) dput(PB)structure(list(Animal = c(25L, 38L, 42L, 50L, 53L, 18L, 22L, 34L, 41L, 44L, 48L, 2L, 12L, 14L, 39L, 46L, 4L, 5L, 21L, 26L, 49L, 3L, 32L, 37L, 43L, 47L, 8L, 15L, 30L, 40L, 1L, 17L, 23L, 28L, 7L, 9L, 19L, 35L), Gest = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L), Manej = c(1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L), rep = c(4L, 5L, 3L, 1L, 2L, 1L, 2L, 6L, 4L, 3L, 5L, 3L, 2L, 5L, 1L, 4L, 2L, 4L, 1L, 3L, 5L, 1L, 5L, 2L, 3L, 4L, 1L, 4L, 3L, 2L, 2L, 3L, 4L, 1L, 3L, 2L, 4L, 1L), PV = c(40.5, 40.9, 35.5, 34.5, 36.7, 26.9, 35.5, 31.2, 31.5, 30, 35.5, 40.8, 36.6, 36.6, 36.5, 37.9, 41.9, 35, 35.5, 32.5, 42, 51.8, 50, 46.5, 38.7, 42.7, 35.4, 38.4, 36.5, 31.2, 49, 46.4, 48.4, 42.8, 40.5, 42.3, 50.2, 33.5), PCVZ = c(34.7811, 33.6882, 22.8156, 26.4014, 26.7112, 19.1745, 26.6094, 22.4639, 25.6466, 22.1737, 25.773, 35.4571, 30.8157, 28.1626, 31.8387, 30.3712, 32.9897, 26.5665, 27.9632, 23.4118, 38.1716, 42.3432, 38.2408, 36.4839, 31.3696, 36.7408, 25.9655, 28.1035, 28.9001, 26.0481, 39.8412, 37.9058, 40.8529, 43.4908, 30.3399, 28.4752, 34.3194, 24.6319), PB = c(9.18910573, 9.11864608, 5.653976965, 6.870682232, 6.838822605, 5.119182736, 6.932858931, 6.047701201, 5.593529478, 5.379896679, 6.611844309, 9.230423109, 8.184234256, 7.20898144, 8.703687287, 8.645435859, 7.017519902, 5.976945686, 7.366941986, 4.644961891, 9.085415538, 11.06095042, 9.461960295, 8.85340357, 8.008399724, 7.058114299, 5.333367984, 5.670736687, 7.363494857, 6.014652251, 7.077788459, 7.890982758, 8.935434043, 9.484189031, 6.804250745, 7.034024007, 7.628728687, 5.124003503)), .Names = c("Animal", "Gest", "Manej", "rep", "PV", "PCVZ", "PB"), class = "data.frame", row.names = c(NA, -38L))

3 6

Intersecção entre dois shapefiles no R
by Ricardo Dobrovolski 06 Dez '12

06 Dez '12

Oi Pessoal, Eu gostaria de saber se vocês conhecem pacotes no R nos quais eu possa fazer operações entre shapefiles como intersecções, clips, etc, de maneira semelhante às operações que se podem realizar no ArcGis. Atenciosamente, Ricardo Dobrovolski

3 2

Ler um arquivo grande
by Daniel Marcelino 06 Dez '12

06 Dez '12

Caros, Estou tentando ler um arquivo de texto com 930 mb, mas sem sucesso desde ontem. Meu computador é multicore Macbook 8 gb memória, mas não consegue finalizar de ler o arquivo e "congela". Alguém tem alguma ideia de como superar isso sem precisar abrir o arquivo e dividir manualmente o banco? Daniel

4 5

Problemas com Encoding
by Marcos Silva 06 Dez '12

06 Dez '12

Caros amigos, Estou tentando usar o pacote XML para ler alguns dados em tabelas existentes em páginas html. Especificamente, estou tentando ler páginas da previdência social contendo informações sobre o CRP (Certificado de Regularidade Previdenciária) O problema é que estou obtendo uns caracteres estranhos, coisas do tipo "SituaÃ§Ã£o". Estou fazendo da seguinte forma: library(XML) html <- readHTMLTable("F:\\Regime Próprio\\html\\ExtratoPrevidenciario.htm") html <- html[[3]][, c('V1', 'V2')] Já tentei o iconv() sem resultado. Já tentei ler após mudar a opção "encoding" em options(). Ah, estou usando o windows velho de guerra... sorry... :) A página pode ser vista no seguinte link: http://www1.previdencia.gov.br/sps/app/crp/ExtratoRegularidadeRegimes.asp?C… Se alguém tiver alguma idéia de corrigir... Desde já obrigado. Abs. -- Marcos F. Silva http://sites.google.com/site/marcosfs2006

2 3

Ajuda Mapa
by Sérgio Henrique almeida da silva ju 06 Dez '12

06 Dez '12

Amigos Estou plotando um mapa, no qual, eu quero destacar que em um município qualquer tem uma doença x. no meu banco a categoria da variável muniprio é "sim" e "não", "sim" significa que ali tem a doença x. Fiz dessa forma, porém não ficou legal, como posso fazer isso? brks <- c("não", "sim") cols <- c("blue", "white") plot(mapa, col=cols[findInterval(mapa.dt$muniprio, brks)], border="white") Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

2 2

Poisson multivariada
by Leonard de Assis 06 Dez '12

06 Dez '12

Existe algum pacote que lida com Poisson multivariada? Eu achei alguma coisa no rseek: http://www.rseek.org/?cx=010923144343702598753%3Aboaz1reyxd4&q=bivariate+po… mas o pacote mencionado por eles (bivpois) não está mais no CRAN e está dando erro de instalação quando baixa a ultima versão disponível. Ahhh, meu interesse é uma Poisson bivariada inflada por zeros. -- []s Leonard de Assis http://about.me/ldeassis

2 1

Re: [R-br] Como criar data.frame com linhas contendo NA's
by Eder Comunello 05 Dez '12

05 Dez '12

Taynãna, boa tarde! Poderia enviar a saída do console no email? Para mim a saída do último comando aparece assim: > dados[order(dados$local),] ### ordenando dia local dado1 dado2 1 1 1 0.10 15 3 2 1 0.15 16 5 3 1 NA NA 7 4 1 0.20 20 2 1 2 0.27 23 4 2 2 NA NA 6 3 2 NA NA 8 4 2 0.15 25 > E seu pedido era: "Preciso criar um data.frame full, completando os dados com NA´s, tipo: dia local dado1 dado2 1 1 0.10 15 2 1 0.15 16 3 1 NA NA 4 1 0.20 20 1 2 0.27 23 2 2 NA NA 3 2 NA NA 4 2 0.15 25" Como aparece no teu sistema? ================================================ Éder Comunello Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ) Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W] Researcher at Embrapa Western Region Agriculture Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] ================================================ UTC-03:00

8 11