Dear UseRs, I have a dataset in the following form.> dput(Da)structure(c(0.0238095238095238, 0.0476190476190476, 0.0714285714285714, 0.0952380952380952, 0.119047619047619, 0.142857142857143, 0.166666666666667, 0.19047619047619, 0.214285714285714, 0.238095238095238, 0.261904761904762, 0.285714285714286, 0.30952380952381, 0.333333333333333, 0.357142857142857, 0.380952380952381, 0.404761904761905, 0.428571428571429, 0.452380952380952, 0.476190476190476, 0.5, 0.523809523809524, 0.547619047619048, 0.571428571428571, 0.595238095238095, 0.619047619047619, 0.642857142857143, 0.666666666666667, 0.69047619047619, 0.714285714285714, 0.738095238095238, 0.761904761904762, 0.785714285714286, 0.80952380952381, 0.833333333333333, 0.857142857142857, 0.880952380952381, 0.904761904761905, 0.928571428571429, 0.952380952380952, 0.976190476190476, 100, 110, 120, 130, 140, 150, 160, 170, 180, 190, 200, 210, 220, 230, 240, 250, 260, 270, 280, 290, 300, 310, 320, 330, 340, 350, 360, 370, 380, 390, 400, 410, 420, 430, 440, 450, 460, 470, 480, 490, 500), .Dim = c(41L, 2L), .Dimnames = list(NULL, c("P", "RP"))) Column number 1 is the probability of non execeedance and column 2 is return period of the event. I plotted the dataset in the following way>plot(Da[,2],Da[,1],log="x")I now want to fit gumbel distribution, LN3 and GEV distributions on the data and plot them on the same window to see which distribution actually fits good to the data. I know there are other ways of determining best fitted distribution in R but I just wanted to know that if I can actually plot it on actual distribution or not. I want it to try for my thesis. Thankyou very much in advance Eliza [[alternative HTML version deleted]]