Olá Pedro, Rodei seu código,ok. Mas não estou conseguindo visualizar a Rosa. Não estou conseguindo coloca-la no mapa, ou até mesmo, plotá-la. -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pela comissão de Comunicação e do Ministério Jovem do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig http://jodavidferreira.com
Jodavid, Você utiliza assim:
RosaDosVentos (c(6.5,14),0.5)
Onde o primeiro argumento é a localização no mapa de onde deve estar a Rosa dos Ventos e o segundo informa o tamanho dela... Você precisa coloca-la em uma região visível do seu mapa. Grande abraço. Pedro De: r-br-bounces@listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de Jodavid Ferreira Enviada em: Monday, February 27, 2012 2:02 PM Para: r-br@listas.c3sl.ufpr.br Assunto: [R-br] RES: Rosa dos Ventos Olá Pedro, Rodei seu código,ok. Mas não estou conseguindo visualizar a Rosa. Não estou conseguindo coloca-la no mapa, ou até mesmo, plotá-la. -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pela comissão de Comunicação e do Ministério Jovem do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig http://jodavidferreira.com
Olá Pedro, Eu estou até entendendo o seu código. Mas ainda não consegui colocar a Rosa dos Ventos. O problema está, em como fazer ela aparecer na plotagem. se ela vai como legenda... se você puder me enviar algum exemplo dela funcionando agradeceria bastante. abraço, Em 27 de fevereiro de 2012 14:01, Jodavid Ferreira <jodavid.arts@gmail.com>escreveu:
Olá Pedro,
Rodei seu código,ok. Mas não estou conseguindo visualizar a Rosa.
Não estou conseguindo coloca-la no mapa, ou até mesmo, plotá-la.
-- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pela comissão de Comunicação e do Ministério Jovem do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig http://jodavidferreira.com
-- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pelo Ministério Jovem e de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com
Jodavid, Segue um exemplo: require(maptools) require(classInt) columbus <- readShapePoly(system.file("etc/shapes/columbus.shp", package="spdep")[1]) mypal <- c("wheat1", "red3") f5 <- classIntervals(columbus$CRIME, n=5, style="fisher") f5Colours <- findColours(f5, mypal) plot(columbus, col=f5Colours,border ="black", axes = TRUE) title(main="Criminalidade em Columbus.",xlab="Longitude",ylab="Latitude") RosaDosVentos <- function(loc,size,bearing=0,cols,cex=1,...) { if(missing(loc)) stop("loc is missing") if(missing(size)) stop("size is missing") if(missing(cols)) cols <- rep(c("white","black"),8) radii <- rep(size/c(1,4,2,4),4) x <- radii[(0:15)+1]*cos((0:15)*pi/8+bearing)+loc[1] y <- radii[(0:15)+1]*sin((0:15)*pi/8+bearing)+loc[2] for (i in 1:15) { x1 <- c(x[i],x[i+1],loc[1]) y1 <- c(y[i],y[i+1],loc[2]) polygon(x1,y1,col=cols[i]) } polygon(c(x[16],x[1],loc[1]),c(y[16],y[1],loc[2]),col=cols[16]) b <- c("L","N","O","S") for (i in 0:3) text((size+par("cxy")[1])*cos(bearing+i*pi/2)+loc[1], (size+par("cxy")[2])*sin(bearing+i*pi/2)+loc[2],b[i+1],cex=cex) } RosaDosVentos (c(5,13),0.4) Por favor, informe se funcionou. Grande abraço. Pedro De: r-br-bounces@listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de Jodavid Ferreira Enviada em: Friday, March 02, 2012 3:13 PM Para: r-br@listas.c3sl.ufpr.br Assunto: Re: [R-br] Rosa dos Ventos Olá Pedro, Eu estou até entendendo o seu código. Mas ainda não consegui colocar a Rosa dos Ventos. O problema está, em como fazer ela aparecer na plotagem. se ela vai como legenda... se você puder me enviar algum exemplo dela funcionando agradeceria bastante. abraço, Em 27 de fevereiro de 2012 14:01, Jodavid Ferreira <jodavid.arts@gmail.com> escreveu: Olá Pedro, Rodei seu código,ok. Mas não estou conseguindo visualizar a Rosa. Não estou conseguindo coloca-la no mapa, ou até mesmo, plotá-la. -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pela comissão de Comunicação e do Ministério Jovem do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig http://jodavidferreira.com -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pelo Ministério Jovem e de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com
Jodavid, Rodei o exemplo do Pedro no meu PC e funcionou. Observe que a função "RosaDosVentos (c(5,13),0.4)" substitua o valor 5 por 6.5 que aparecera a Rosa dos ventos Humberto Em 3/2/2012 6:30 PM, Pedro Henrique Melo Albuquerque escreveu:
Jodavid,
Segue um exemplo:
require(maptools)
require(classInt)
columbus <- readShapePoly(system.file("etc/shapes/columbus.shp", package="spdep")[1])
mypal <- c("wheat1", "red3")
f5 <- classIntervals(columbus$CRIME, n=5, style="fisher")
f5Colours <- findColours(f5, mypal)
plot(columbus, col=f5Colours,border ="black", axes = TRUE)
title(main="Criminalidade em Columbus.",xlab="Longitude",ylab="Latitude")
RosaDosVentos <- function(loc,size,bearing=0,cols,cex=1,...)
{
if(missing(loc)) stop("loc is missing")
if(missing(size)) stop("size is missing")
if(missing(cols)) cols <- rep(c("white","black"),8)
radii <- rep(size/c(1,4,2,4),4)
x <- radii[(0:15)+1]*cos((0:15)*pi/8+bearing)+loc[1]
y <- radii[(0:15)+1]*sin((0:15)*pi/8+bearing)+loc[2]
for (i in 1:15)
{
x1 <- c(x[i],x[i+1],loc[1])
y1 <- c(y[i],y[i+1],loc[2])
polygon(x1,y1,col=cols[i])
}
polygon(c(x[16],x[1],loc[1]),c(y[16],y[1],loc[2]),col=cols[16])
b <- c("L","N","O","S")
for (i in 0:3) text((size+par("cxy")[1])*cos(bearing+i*pi/2)+loc[1],
(size+par("cxy")[2])*sin(bearing+i*pi/2)+loc[2],b[i+1],cex=cex)
}
RosaDosVentos (c(5,13),0.4)
Por favor, informe se funcionou.
Grande abraço.
Pedro
*De:*r-br-bounces@listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] *Em nome de *Jodavid Ferreira *Enviada em:* Friday, March 02, 2012 3:13 PM *Para:* r-br@listas.c3sl.ufpr.br *Assunto:* Re: [R-br] Rosa dos Ventos
Olá Pedro, Eu estou até entendendo o seu código.
Mas ainda não consegui colocar a Rosa dos Ventos.
O problema está, em como fazer ela aparecer na plotagem.
se ela vai como legenda...
se você puder me enviar algum exemplo dela funcionando agradeceria bastante.
abraço,
Em 27 de fevereiro de 2012 14:01, Jodavid Ferreira <jodavid.arts@gmail.com <mailto:jodavid.arts@gmail.com>> escreveu:
Olá Pedro,
Rodei seu código,ok. Mas não estou conseguindo visualizar a Rosa.
Não estou conseguindo coloca-la no mapa, ou até mesmo, plotá-la.
-- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pela comissão de Comunicação e do Ministério Jovem do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig <http://de.ufpb.br/%7Eleapig> http://jodavidferreira.com
--
Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pelo Ministério Jovem e de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig <http://de.ufpb.br/%7Eleapig> Site Pessoal - http://jodavidferreira.com
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Eu tentei rodar este exemplo e não consegui! Já deu erro no terceiro comando: columbus <- readShapePoly(system.file("etc/shapes/columbus.shp", package="spdep")[1]) Erro em getinfo.shape(filen) : Error opening SHP file Jodavid, Rodei o exemplo do Pedro no meu PC e funcionou. Observe que a função "RosaDosVentos (c(5,13),0.4)" substitua o valor 5 por 6.5 que aparecera a Rosa dos ventos Humberto Em 3/2/2012 6:30 PM, Pedro Henrique Melo Albuquerque escreveu: Jodavid, Segue um exemplo: require(maptools) require(classInt) columbus <- readShapePoly(system.file("etc/shapes/columbus.shp", package="spdep")[1]) mypal <- c("wheat1", "red3") f5 <- classIntervals(columbus$CRIME, n=5, style="fisher") f5Colours <- findColours(f5, mypal) plot(columbus, col=f5Colours,border ="black", axes = TRUE) title(main="Criminalidade em Columbus.",xlab="Longitude",ylab="Latitude") RosaDosVentos <- function(loc,size,bearing=0,cols,cex=1,...) { if(missing(loc)) stop("loc is missing") if(missing(size)) stop("size is missing") if(missing(cols)) cols <- rep(c("white","black"),8) radii <- rep(size/c(1,4,2,4),4) x <- radii[(0:15)+1]*cos((0:15)*pi/8+bearing)+loc[1] y <- radii[(0:15)+1]*sin((0:15)*pi/8+bearing)+loc[2] for (i in 1:15) { x1 <- c(x[i],x[i+1],loc[1]) y1 <- c(y[i],y[i+1],loc[2]) polygon(x1,y1,col=cols[i]) } polygon(c(x[16],x[1],loc[1]),c(y[16],y[1],loc[2]),col=cols[16]) b <- c("L","N","O","S") for (i in 0:3) text((size+par("cxy")[1])*cos(bearing+i*pi/2)+loc[1], (size+par("cxy")[2])*sin(bearing+i*pi/2)+loc[2],b[i+1],cex=cex) } RosaDosVentos (c(5,13),0.4) Por favor, informe se funcionou. Grande abraço. Pedro De: r-br-bounces@listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de Jodavid Ferreira Enviada em: Friday, March 02, 2012 3:13 PM Para: r-br@listas.c3sl.ufpr.br Assunto: Re: [R-br] Rosa dos Ventos Olá Pedro, Eu estou até entendendo o seu código. Mas ainda não consegui colocar a Rosa dos Ventos. O problema está, em como fazer ela aparecer na plotagem. se ela vai como legenda... se você puder me enviar algum exemplo dela funcionando agradeceria bastante. abraço, Em 27 de fevereiro de 2012 14:01, Jodavid Ferreira <jodavid.arts@gmail.com> escreveu: Olá Pedro, Rodei seu código,ok. Mas não estou conseguindo visualizar a Rosa. Não estou conseguindo coloca-la no mapa, ou até mesmo, plotá-la. -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pela comissão de Comunicação e do Ministério Jovem do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig http://jodavidferreira.com -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pelo Ministério Jovem e de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com
Mauro, Instale o pacote spdep e tente novamente. Caso queira fazer de maneira direta você pode baixar o shapefile nesse link: http://geodacenter.org/downloads/data-files/columbus.zip De: r-br-bounces@listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de Mauro Sznelwar Enviada em: Friday, March 02, 2012 8:44 PM Para: r-br@listas.c3sl.ufpr.br Assunto: Re: [R-br] RES: Rosa dos Ventos Eu tentei rodar este exemplo e não consegui! Já deu erro no terceiro comando: columbus <- readShapePoly(system.file("etc/shapes/columbus.shp", package="spdep")[1]) Erro em getinfo.shape(filen) : Error opening SHP file Jodavid, Rodei o exemplo do Pedro no meu PC e funcionou. Observe que a função "RosaDosVentos (c(5,13),0.4)" substitua o valor 5 por 6.5 que aparecera a Rosa dos ventos Humberto Em 3/2/2012 6:30 PM, Pedro Henrique Melo Albuquerque escreveu: Jodavid, Segue um exemplo: require(maptools) require(classInt) columbus <- readShapePoly(system.file("etc/shapes/columbus.shp", package="spdep")[1]) mypal <- c("wheat1", "red3") f5 <- classIntervals(columbus$CRIME, n=5, style="fisher") f5Colours <- findColours(f5, mypal) plot(columbus, col=f5Colours,border ="black", axes = TRUE) title(main="Criminalidade em Columbus.",xlab="Longitude",ylab="Latitude") RosaDosVentos <- function(loc,size,bearing=0,cols,cex=1,...) { if(missing(loc)) stop("loc is missing") if(missing(size)) stop("size is missing") if(missing(cols)) cols <- rep(c("white","black"),8) radii <- rep(size/c(1,4,2,4),4) x <- radii[(0:15)+1]*cos((0:15)*pi/8+bearing)+loc[1] y <- radii[(0:15)+1]*sin((0:15)*pi/8+bearing)+loc[2] for (i in 1:15) { x1 <- c(x[i],x[i+1],loc[1]) y1 <- c(y[i],y[i+1],loc[2]) polygon(x1,y1,col=cols[i]) } polygon(c(x[16],x[1],loc[1]),c(y[16],y[1],loc[2]),col=cols[16]) b <- c("L","N","O","S") for (i in 0:3) text((size+par("cxy")[1])*cos(bearing+i*pi/2)+loc[1], (size+par("cxy")[2])*sin(bearing+i*pi/2)+loc[2],b[i+1],cex=cex) } RosaDosVentos (c(5,13),0.4) Por favor, informe se funcionou. Grande abraço. Pedro De: r-br-bounces@listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de Jodavid Ferreira Enviada em: Friday, March 02, 2012 3:13 PM Para: r-br@listas.c3sl.ufpr.br Assunto: Re: [R-br] Rosa dos Ventos Olá Pedro, Eu estou até entendendo o seu código. Mas ainda não consegui colocar a Rosa dos Ventos. O problema está, em como fazer ela aparecer na plotagem. se ela vai como legenda... se você puder me enviar algum exemplo dela funcionando agradeceria bastante. abraço, Em 27 de fevereiro de 2012 14:01, Jodavid Ferreira <jodavid.arts@gmail.com> escreveu: Olá Pedro, Rodei seu código,ok. Mas não estou conseguindo visualizar a Rosa. Não estou conseguindo coloca-la no mapa, ou até mesmo, plotá-la. -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pela comissão de Comunicação e do Ministério Jovem do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig <http://de.ufpb.br/%7Eleapig> http://jodavidferreira.com -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pelo Ministério Jovem e de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig <http://de.ufpb.br/%7Eleapig> Site Pessoal - http://jodavidferreira.com
O exemplo funcionou perfeitamente! A Rosa dos Ventos ficou demais! É interessante fazer isto para mapas do Brasil!
________________________________ De: Pedro Henrique Melo Albuquerque <pedro.melo.albuquerque@gmail.com> Para: r-br@listas.c3sl.ufpr.br Enviadas: Sexta-feira, 2 de Março de 2012 20:49 Assunto: [R-br] RES: RES: Rosa dos Ventos
Mauro, Instale o pacote spdep e tente novamente. Caso queira fazer de maneira direta você pode baixar o shapefile nesse link: http://geodacenter.org/downloads/data-files/columbus.zip De:r-br-bounces@listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de Mauro Sznelwar Enviada em: Friday, March 02, 2012 8:44 PM Para: r-br@listas.c3sl.ufpr.br Assunto: Re: [R-br] RES: Rosa dos Ventos Eu tentei rodar este exemplo e não consegui! Já deu erro no terceiro comando: columbus <- readShapePoly(system.file("etc/shapes/columbus.shp", package="spdep")[1]) Erro em getinfo.shape(filen) : Error opening SHP file
Jodavid, Rodei o exemplo do Pedro no meu PC e funcionou. Observe que a função "RosaDosVentos (c(5,13),0.4)" substitua o valor 5 por 6.5 que aparecera a Rosa dos ventos
Humberto
Em 3/2/2012 6:30 PM, Pedro Henrique Melo Albuquerque escreveu: Jodavid, Segue um exemplo: require(maptools) require(classInt) columbus <- readShapePoly(system.file("etc/shapes/columbus.shp", package="spdep")[1]) mypal <- c("wheat1", "red3") f5 <- classIntervals(columbus$CRIME, n=5, style="fisher") f5Colours <- findColours(f5, mypal) plot(columbus, col=f5Colours,border ="black", axes = TRUE) title(main="Criminalidade em Columbus.",xlab="Longitude",ylab="Latitude") RosaDosVentos <- function(loc,size,bearing=0,cols,cex=1,...) { if(missing(loc)) stop("loc is missing") if(missing(size)) stop("size is missing") if(missing(cols)) cols <- rep(c("white","black"),8) radii <- rep(size/c(1,4,2,4),4) x <- radii[(0:15)+1]*cos((0:15)*pi/8+bearing)+loc[1] y <- radii[(0:15)+1]*sin((0:15)*pi/8+bearing)+loc[2] for (i in 1:15) { x1 <- c(x[i],x[i+1],loc[1]) y1 <- c(y[i],y[i+1],loc[2]) polygon(x1,y1,col=cols[i]) } polygon(c(x[16],x[1],loc[1]),c(y[16],y[1],loc[2]),col=cols[16]) b <- c("L","N","O","S") for (i in 0:3) text((size+par("cxy")[1])*cos(bearing+i*pi/2)+loc[1], (size+par("cxy")[2])*sin(bearing+i*pi/2)+loc[2],b[i+1],cex=cex) } RosaDosVentos (c(5,13),0.4) Por favor, informe se funcionou. Grande abraço. Pedro De:r-br-bounces@listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de Jodavid Ferreira Enviada em: Friday, March 02, 2012 3:13 PM Para: r-br@listas.c3sl.ufpr.br Assunto: Re: [R-br] Rosa dos Ventos Olá Pedro, Eu estou até entendendo o seu código. Mas ainda não consegui colocar a Rosa dos Ventos. O problema está, em como fazer ela aparecer na plotagem. se ela vai como legenda... se você puder me enviar algum exemplo dela funcionando agradeceria bastante. abraço, Em 27 de fevereiro de 2012 14:01, Jodavid Ferreira <jodavid.arts@gmail.com> escreveu: Olá Pedro, Rodei seu código,ok. Mas não estou conseguindo visualizar a Rosa. Não estou conseguindo coloca-la no mapa, ou até mesmo, plotá-la.
-- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pela comissão de Comunicação e do Ministério Jovem do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig http://jodavidferreira.com
-- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pelo Ministério Jovem e de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com
R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Essa programação eu tirei daqui: <http://www.jstatsoft.org/v19/c01/paper> http://www.jstatsoft.org/v19/c01/paper De: r-br-bounces@listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de Alexandro (Yahoo) Enviada em: Friday, March 02, 2012 9:26 PM Para: r-br@listas.c3sl.ufpr.br Assunto: Re: [R-br] RES: RES: Rosa dos Ventos O exemplo funcionou perfeitamente! A Rosa dos Ventos ficou demais! É interessante fazer isto para mapas do Brasil! _____ De: Pedro Henrique Melo Albuquerque <pedro.melo.albuquerque@gmail.com> Para: r-br@listas.c3sl.ufpr.br Enviadas: Sexta-feira, 2 de Março de 2012 20:49 Assunto: [R-br] RES: RES: Rosa dos Ventos Mauro, Instale o pacote spdep e tente novamente. Caso queira fazer de maneira direta você pode baixar o shapefile nesse link: http://geodacenter.org/downloads/data-files/columbus.zip De: r-br-bounces@listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de Mauro Sznelwar Enviada em: Friday, March 02, 2012 8:44 PM Para: r-br@listas.c3sl.ufpr.br Assunto: Re: [R-br] RES: Rosa dos Ventos Eu tentei rodar este exemplo e não consegui! Já deu erro no terceiro comando: columbus <- readShapePoly(system.file("etc/shapes/columbus.shp", package="spdep")[1]) Erro em getinfo.shape(filen) : Error opening SHP file Jodavid, Rodei o exemplo do Pedro no meu PC e funcionou. Observe que a função "RosaDosVentos (c(5,13),0.4)" substitua o valor 5 por 6.5 que aparecera a Rosa dos ventos Humberto Em 3/2/2012 6:30 PM, Pedro Henrique Melo Albuquerque escreveu: Jodavid, Segue um exemplo: require(maptools) require(classInt) columbus <- readShapePoly(system.file("etc/shapes/columbus.shp", package="spdep")[1]) mypal <- c("wheat1", "red3") f5 <- classIntervals(columbus$CRIME, n=5, style="fisher") f5Colours <- findColours(f5, mypal) plot(columbus, col=f5Colours,border ="black", axes = TRUE) title(main="Criminalidade em Columbus.",xlab="Longitude",ylab="Latitude") RosaDosVentos <- function(loc,size,bearing=0,cols,cex=1,...) { if(missing(loc)) stop("loc is missing") if(missing(size)) stop("size is missing") if(missing(cols)) cols <- rep(c("white","black"),8) radii <- rep(size/c(1,4,2,4),4) x <- radii[(0:15)+1]*cos((0:15)*pi/8+bearing)+loc[1] y <- radii[(0:15)+1]*sin((0:15)*pi/8+bearing)+loc[2] for (i in 1:15) { x1 <- c(x[i],x[i+1],loc[1]) y1 <- c(y[i],y[i+1],loc[2]) polygon(x1,y1,col=cols[i]) } polygon(c(x[16],x[1],loc[1]),c(y[16],y[1],loc[2]),col=cols[16]) b <- c("L","N","O","S") for (i in 0:3) text((size+par("cxy")[1])*cos(bearing+i*pi/2)+loc[1], (size+par("cxy")[2])*sin(bearing+i*pi/2)+loc[2],b[i+1],cex=cex) } RosaDosVentos (c(5,13),0.4) Por favor, informe se funcionou. Grande abraço. Pedro De: r-br-bounces@listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de Jodavid Ferreira Enviada em: Friday, March 02, 2012 3:13 PM Para: r-br@listas.c3sl.ufpr.br Assunto: Re: [R-br] Rosa dos Ventos Olá Pedro, Eu estou até entendendo o seu código. Mas ainda não consegui colocar a Rosa dos Ventos. O problema está, em como fazer ela aparecer na plotagem. se ela vai como legenda... se você puder me enviar algum exemplo dela funcionando agradeceria bastante. abraço, Em 27 de fevereiro de 2012 14:01, Jodavid Ferreira <jodavid.arts@gmail.com> escreveu: Olá Pedro, Rodei seu código,ok. Mas não estou conseguindo visualizar a Rosa. Não estou conseguindo coloca-la no mapa, ou até mesmo, plotá-la. -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pela comissão de Comunicação e do Ministério Jovem do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig <http://de.ufpb.br/%7Eleapig> http://jodavidferreira.com -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pelo Ministério Jovem e de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig <http://de.ufpb.br/%7Eleapig> Site Pessoal - http://jodavidferreira.com _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Prezados, Estou com uma dúvida que pode parecer besteira, porém gostarai da ajuda de vocês! A dúvida é a seguinte como calcular o percentual de "Bias" (percentage prediction error) por linha de um modelo glm? Meus dados tem 500 linhas e 4 colunas (anexo). O modelo é: m1<- glm(lgcpue=fy+fFF+fq+fA, family=gaussian, data=teste) Meu objetivo é calcular o percentual de desvio (bias) de um modelo por ano via box plot para saber o quanto ele esta sobre ou sub estimando minhas predições com a variável fFF e sem a variável fFF, por exemplo! todos as variaveis são fatores Agradeço desde já Humberto Hazin
Consegui rodar e funcionou perfeito! Obrigado Pedro, e a você também Humberto pela dica da localização! Me ajudaram bastante Em 2 de março de 2012 15:13, Jodavid Ferreira <jodavid.arts@gmail.com>escreveu:
Olá Pedro, Eu estou até entendendo o seu código.
Mas ainda não consegui colocar a Rosa dos Ventos.
O problema está, em como fazer ela aparecer na plotagem.
se ela vai como legenda...
se você puder me enviar algum exemplo dela funcionando agradeceria bastante.
abraço,
Em 27 de fevereiro de 2012 14:01, Jodavid Ferreira <jodavid.arts@gmail.com
escreveu:
Olá Pedro,
Rodei seu código,ok. Mas não estou conseguindo visualizar a Rosa.
Não estou conseguindo coloca-la no mapa, ou até mesmo, plotá-la.
-- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pela comissão de Comunicação e do Ministério Jovem do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig http://jodavidferreira.com
-- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pelo Ministério Jovem e de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com
Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com
-- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pelo Ministério Jovem e de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com
participantes (5)
-
Alexandro (Yahoo) -
Humberto Hazin -
Jodavid Ferreira -
Mauro Sznelwar -
Pedro Henrique Melo Albuquerque