Dear useRs,
I have the following dataset.
c(77.724, 76.708, 84.836, 85.09, 118.11, 65.024, 121.412, 63.5, 102.87, 81.3, 108.7, 110.7, 71.9, 42.2, 101, 151.4, 94, 112, 48, 73.4, 76.6, 62.2, 59.4, 114.3, 214.3, 110.5, 46, 84.7, 128.1, 45.2, 109.5, 102.3, 77.5, 61, 97.3, 78, 142, 88.2, 54, 91.4, 54.1, 96, 143.3, 153.7, 101.5, 95.8, 101, 131, 140, 189.4)
I fitted generalized Extreme Value distribution on it by using following codes
library(nsRFA)
q=EB
lmom=Lmoments(q)
pr = par.GEV (lambda1=lmom["l1"], lambda2=lmom["l2"], tau3=lmom["lca"])
RP = c(1.01,2, 10, 20, 50, 100, 200, 500)
quant = invF.GEV (1-1/RP, pr$xi, pr$alfa, pr$k)
qs = sort(q)
pp = 1:length(qs)/(length(qs)+1)
RPpp = 1/(1-pp)
plot(RP, quant, type="l", log="x",col="black",ylim=c(0,500),xlim=c(0.1,500))
points(RPpp, qs)
What I want to do now is to calculate the distance of all the points from the "fitted line" and ultimately calculating RMSE of the data.
Is there a way of doing it?
Thankyou very much in advance
Eliza