Junho 2012 - R-br - DINF UFPR

Colorir alguns municípios em um mapa estadual com cores distintas
by Wirton Macedo Coutinho 27 Jun '12

27 Jun '12

Boa tarde pessoal, Gostaria de saber se alguém tem alguma dica de como colorir alguns municipios em um mapa estadual (extensão shp) com cores distintas. Preciso colorir 7 municipios no mapa do Estado de Goiás com cores distintas (cada município com uma cor diferente), mas não estou conseguindo. Tenho tentado, sem sucesso, da seguinte forma: plot(mapa, axes=TRUE, col=mapa$Nome_Munic=="Santa Helena de Goiás") plot(mapa, col=mapa$Nome_Munic=="Chapadão do Céu", add=TRUE) plot(mapa, col=mapa$Nome_Munic=="Montividiu", add=TRUE) plot(mapa, col=mapa$Nome_Munic=="Acreúna", add=TRUE) plot(mapa, col=mapa$Nome_Munic=="Perolândia", add=TRUE) plot(mapa, col=mapa$Nome_Munic=="Palmeiras de Goiás", add=TRUE) plot(mapa, col=mapa$Nome_Munic=="Cristalina", add=TRUE) Apesar de conseguir colorir todos os 7 municipios, com os comando acima, todos os 7 ficam com a mesma cor. Certo de contar com a boa vontade de todos, agradeço antecipadamente. Att., -- Wirton Macedo Coutinho Pesquisador Fitopatologia Embrapa Algodão Rua Oswaldo Cruz, 1143, Centenário Campina Grande PB CEP 28428-095

3 4

nls comparação de curvas
by Alisson Lucrecio 27 Jun '12

27 Jun '12

Eu estou com dificuldades para comparar curvas de sorção alguem pode me ajudar? Segue minha workspece e o scrip. Obrigado. Alisson Lucrécio da Costa

1 0

Enviar email pelo R
by Leonardo Aguirre 27 Jun '12

27 Jun '12

Pessoal, Alguém poderia me indicar uma referência, ou um exemplo de programa que envie email pelo R? Preciso de um exemplo que aceite autenticação de login e senha, e também envie anexos. Procurei no google em vários sites mas não consegui nenhum exemplo que funcionasse corretamente. Preciso enviar 300 emails e quero fazê-lo de forma automática. Obrigado, Leonardo.

2 1

Re: [R-br] Solução snow+Rmpi
by Junior Beleti 26 Jun '12

26 Jun '12

Entendi Benilton. Sabe me dizer quais ferramentas precisam estar instaladas no cluster? Por exemplo, o MPI? Ou apenas os pacotes são suficientes? Att, -- Junior Beleti

2 1

nls
by Alisson Lucrecio 26 Jun '12

26 Jun '12

Olá, Eu ajustei os resultados de sorção de 3 solos para um pestida pelo modelo de freudlich no R com a função nls, entretanto preciso comparar esses resultados para ver se extiste diferença entre a sorção dos solos e qual apresenta maior sorção. Alguem sabe algum teste que posso fazer para comparar modelos nls? Obrigado. Alisson Lucrécio da Costa

2 1

(sem assunto)
by Alisson Lucrecio 26 Jun '12

26 Jun '12

Olá, Eu ajustei os resultados de sorção de 3 solos para um pestida pelo modelo de freudlich no R com a função nls, entretanto preciso comparar esses resultados para ver se existe diferença entre a sorção dos solos e qual apresenta maior sorção. Eu preciso de uma função para comparar resultado de modelos nls. Obrigado. Alisson Lucrécio da Costa

1 0

nls
by Alisson Lucrecio 26 Jun '12

26 Jun '12

Olá, Eu ajustei os resultados de sorção de 3 solos para um pestida pelo modelo de freudlich no R com a função nls, entretanto preciso comparar esses resultadas para ver se extiste diferença entre a sorção dos solos e qual apresenta maior sorção. Obrigado. Alisson Lucrécio da Costa

1 0

legenda em gráficos
by Edson Ferreira 26 Jun '12

26 Jun '12

olá amigos do R estou tentando fazer uma legenda e não estou conseguindo segue o código que estou usando e as modificações que quero fazer f<-function(x,alfa,beta){(beta/(sqrt(2*pi)*x^2*pnorm(alfa)))*exp(-0.5*(alfa-beta/x)^2)} BA<-function(x,alfa,beta,a,b){(beta/x^2)*(1/beta(a,b))*(dnorm(alfa-beta/x)/pnorm(alfa))*(pnorm(alfa-beta/x)/pnorm(alfa))^(a-1)*(1-pnorm(alfa-beta/x)/pnorm(alfa))^(b-1)} ##gráfico 1 curve(f(x,alfa=1,beta=1),ylab = "f(x)",lty=1,ylim=c(0,2.5),xlim=c(0,2),main="alfa=1 e beta=1",lwd=2,col=1) curve(BA(x,alfa=1,beta=1,a=0.5,b=0.5),add = T,lty=2,ylim=c(0,2.5),xlim=c(0,2),lwd=2,col=2) curve(BA(x,alfa=1,beta=1,a=0.5,b=2),add = T,lty=3,ylim=c(0,2.5),xlim=c(0,2),lwd=2,col=3) curve(BA(x,alfa=1,beta=1,a=2,b=0.5),add = T,lty=4,ylim=c(0,2.5),xlim=c(0,2),lwd=2,col=4) legend("topright",cex=1,legend = c("alfa","a=0.5,b=0.5","a=0.5,b=2","a=2 ,b=0.5"),lty= 1:4,lwd=2,col=1:4) modificações que estou tentando fazer: Estou precisando que no título do gráfico apareça as letras alpha e beta e não os nomes das letras. O mesmo com a legenda, onde aparece o nome alfa deve permanecer o nome alfa, mas onde aparece as letras a e b, quero que apareça as letras lambda e sigma. Att Edson

2 1

Ajuda modelo hierárquico de poisson com variância robusta
by Sérgio Henrique almeida da silva ju 26 Jun '12

26 Jun '12

Olá Estou enviando uma dúvida de um amigo: Ele está rodando um modelo hierárquico de poisson com variância robusta structure(list(AFDIC = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L), F6rec = structure(c(1L, 2L, 1L, 5L, 4L, 5L, 5L, 5L, 2L, 5L, 1L, 4L, 5L, 5L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 4L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 4L, 2L, 2L, 2L, 3L, 2L, 3L, 2L, 2L, 5L, 1L, 1L, 2L, 3L, 2L, 2L, 1L, 2L, 4L, 3L, 1L, 3L, 1L, 3L, 1L, 2L, 1L, 1L, 2L, 2L, 3L, 1L, 5L, 1L, 4L, 2L, 1L, 2L, 2L, 1L, 2L, 4L, 2L, 1L, 3L, 5L, 2L, 1L, 3L, 2L, 4L, 1L, 2L, 4L, 4L, 2L, 2L, 1L, 3L, 2L, 4L, 2L, 2L), .Label = c("1", "2", "3", "4", "5"), class = "factor")), .Names = c("AFDIC", "F6rec"), row.names = c(NA, -99L), class = "data.frame") # Quando ele roda o modelo o resultado do p-valor é para cada categoria da variável, porém ele gostaria também de obter o valor da variável como um todo: banco$F6rec=as.factor(banco$F6rec) library(epicalc) library(lmtest) library(sandwich) model2=glm(AFDIC~F6rec,family=poisson,data=banco) coeftest(model2, vcov=sandwich) z test of coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) -2.77259 0.68465 -4.0496 5.13e-05 *** F6rec2 1.49165 0.73551 2.0281 0.04255 * F6rec3 1.16315 0.93207 1.2479 0.21206 F6rec4 -15.53000 0.74810 -20.7592 < 2.2e-16 *** F6rec5 0.47000 1.16994 0.4017 0.68788 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 #Observe que o resultado é para cada categoria de F6rec (F6rec2, F6rec3, F6rec4, F6rec,5), porém ele gostaria de obter também o p-valor para o F6rec como um todo, quando ele rodar o modelo multivariado ele que decidir se mantem essa variável ou não. Vocês têm alguma ideia? 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

4 7

Solução snow+Rmpi
by Junior Beleti 26 Jun '12

26 Jun '12

Boa noite Benilton. Buscando alternativas para a solução do problema apresentado no post: "[R-br] Alteração em função de pacote (novo tópico)", vi que em um dos tópicos de discussão você mencionou que utilizava os pacotes snow + Rmpi. Para a resolução deste meu "problema" esses pacotes poderiam ser utilizados? Att, -- Junior Beleti

2 1