Dezembro 2012 - R-br

Correlação Espacial
by Ricardo Dobrovolski 07 Dez '12

07 Dez '12

Caros colegas, Alguém conhece pacotes do R que permitam realizar correlações com correção dos graus de liberdade devido à autocorrelação espacial? Os métodos mais comuns para essa correção estão descritos em: Dutilleul, P. 1993. Modifying the t-test for assessing the correlation between two spatial processes. Biometrics 49: 305-314. Clifford, P., Richardson, S. & Hémon, D. (1989) Assessing the significance of the correlation between two spatial processes. Biometrics, 45, 123–134. Atenciosamente, Ricardo Dobrovolski

2 1

Semivariogramas experimentais com envelopes simulados
by Hélio Gallo Rocha 07 Dez '12

07 Dez '12

Prezados, boa tarde. Lendo uma dissertação sobre geoestatística aplicada ao levantamento de bicho-mineiro em café. Ela usa um recurso para verificar se o atributo possui dependência espacial, chamado de envelope simulado, interessante. Minha dúvida é sobre o parâmetro nsim, assim descrito: nsim - número de simulações utilizadas para calcular os envelopes. O padrão é 99. Seria o nível de significância? Se não for, teria como verificar e mudar esse valor? Exemplos s100.ml <- likfit (s100, ini c = (0,5, 0,5), fix.nugget = TRUE) s100.vario <- variog (s100, max.dist = 1) s100.env <- variog.model.env (s100, obj.v = s100.vario, model.pars = s100.ml) parcela (s100.vario, env = s100.env) Grato Hélio

2 2

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 07 Dez '12

07 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