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 4444 | E: louise.wilson@csiro.au

	[[alternative HTML version deleted]]


