Re: [R] Masking oceans using polypath

2013-07-17 Thread Louise Wilson 
Many thanks, this works beautifully.

Louise


On Wed, Jul 17, 2013 at 7:42 AM, Paul Murrell wrote:

> Hi
>
> There are a couple of problems:
>
> 1.
> Your 'outline' is much bigger than it needs to be.  For example, the
> following produces just Australia ...
>
> outline <- map("worldHires", regions="Australia", exact=TRUE,
>
>plot=FALSE) # returns a list of x/y coords
>
> 2.
> The outline you get still has NA values in it (the coastline of Australia
> in this database is a series of distinct lines;  I don't know why that is).
>  Fortunately, you can stitch Australia back together again like this ...
>
> subset <- !is.na(outline$x)
> polypath(c(outline$x[subset], NA, c(xbox, rev(xbox))),
>  c(outline$y[subset], NA, rep(ybox, each=2)),
>
>  col="light blue", rule="evenodd")
>
> With those two changes, I get a much better result.  You may want to
> fiddle a bit more to add Tasmania and bits of PNG and Indonesia back in the
> picture OR you could replace these problems with a different approach and
> use a more sophisticated map outline via a "shapefile" (see the 'sp'
> package and the 'maptools' package and possibly the R-sig-geo mailing list).
>
> Hope that helps.
>
> Paul
>
>
> On 07/16/13 23:17, Louise Wilson wrote:
>
>> library(mapdata)
>>
>> image(x=110:155, y =-40:-10, z = outer(1:45, 1:30, "+"),
>>
>>xlab = "lon", ylab = "lat")
>>
>> outline <- map("worldHires", plot=FALSE) # returns a list of x/y coords
>>
>> xrange <- range(outline$x, na.rm=TRUE) # get bounding box
>>
>> yrange <- range(outline$y, na.rm=TRUE)
>>
>> xbox <- xrange + c(-2, 2)
>>
>> ybox <- yrange + c(-2, 2)
>>
>> # create the grid path in the current device
>>
>> polypath(c(outline$x, NA, c(xbox, rev(xbox))),
>>
>>   c(outline$y, NA, rep(ybox, each=2)),
>>
>>   col="light blue", rule="evenodd")
>>
>
> --
> Dr Paul Murrell
> Department of Statistics
> The University of Auckland
> Private Bag 92019
> Auckland
> New Zealand
> 64 9 3737599 x85392
> p...@stat.auckland.ac.nz
> http://www.stat.auckland.ac.**nz/~paul/<http://www.stat.auckland.ac.nz/~paul/>
>



-- 

*Louise Wilson*

Climate Projections Project Officer

CSIRO Climate Adaptation Flagship

National Resource Management (NRM)

Pacific-Australia Climate Change Science and Adaptation Planning Program
(PACCSAP)



CSIRO Marine and Atmospheric Research

Private Bag 1 (107-121 Station St)

Aspendale, Victoria 3195
P: +61 3 9239 4619 | F: +61 3 9239  | E: louise.wil...@csiro.au

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

[R] Masking oceans using polypath

2013-07-16 Thread Louise Wilson 
Hi R-help

I am trying to mask the ocean from an image plot I have made.


Here is some example code:

library(mapdata)

image(x=110:155, y =-40:-10, z = outer(1:45, 1:30, "+"),

  xlab = "lon", ylab = "lat")

outline <- map("worldHires", plot=FALSE) # returns a list of x/y coords

xrange <- range(outline$x, na.rm=TRUE) # get bounding box

yrange <- range(outline$y, na.rm=TRUE)

xbox <- xrange + c(-2, 2)

ybox <- yrange + c(-2, 2)

# create the grid path in the current device

polypath(c(outline$x, NA, c(xbox, rev(xbox))),

 c(outline$y, NA, rep(ybox, each=2)),

 col="light blue", rule="evenodd")


As you can see, the polypath is on both sides of the country outline. Any
help is much appreciated.


Louise

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

[R] (no subject)

2012-10-11 Thread louise wilson 
Hi All,



I am performing GEV analysis on temperature/precipitation data using L
moments:


dim(data)

[1] 145 192 156


Lmoments <- apply(data, 1:2, function(x) samlmu(x,nmom=4,sort.data=TRUE))

params_GEV <- apply(Lmoments,2:3,pelgev)


location <- params_GEV[1,,] (xi)

scale <- params_GEV[2,,] (alpha)

shape <- params_GEV[3,,] (k)


I want to use the ks.boot function but I am unsure of how to implement it
using the pgev distribution.



For example:



ks.test(data[1,1],pgev,shape[1,1],location[1,1],scale[1,1])


ks.boot(data,(distribution parameters?),alternative=”pgev”,nboots=1000)



Any advice?



Regards,

Louise Wilson

-- 

*Louise Wilson*

Climate Projections Project Officer

CSIRO Climate Adaptation Flagship

National Resource Management (NRM)

Pacific-Australia Climate Change Science and Adaptation Planning Program
(PACCSAP)



CSIRO Marine and Atmospheric Research

Private Bag 1 (107-121 Station St)

Aspendale, Victoria 3195
P: +61 3 9239 4619 | F: +61 3 9239  | E: louise.wil...@csiro.au

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.