[R] QQ plot

[Peter Dalgaard]

Madan Sigdel wrote:
> Dear users
>  
> I have applied following for my works:
>  
>> library(car)
>> x<-scan()
> 1: 0.92545
> 2: 0.89321
> 3: 0.9846
> 4: 2.9
> 5: 0.85968
> 6: 5.2
> 7: 4.66
> 8: 1.18788
> 9: 1.07683
> 10: 1.07683
> 11: 8.38
> 12: 7.423
> 13: 0.972
> 14: 3.73
> 15: 1.06474
> 16: 1.48
> 17: 0.92876
> 18: 2.26493
> 19: 0.85696
> 20: 1.89313
> 21: 2.71
> 22: 5
> 23: 3.02
> 24: 0.90369
> 25: 8.81
> 26: 1.69466
> 27: 1.07055
> 28: 1.17077
> 29: 2.31647
> 30: 0.83481
> 31: 5.42
> 32: 9.68
> 33: 1.27294
> 34: 5.49
> 35: 2.48
> 36: 1.55876
> 37: 1.41419
> 38: 0.94503
> 39: 5.24
> 40: 
> Read 39 items
>> x<-rgamma(39,shape=1.7,scale=1.6)
>> qq.plot(x,dist="gamma",shape=1.7,scale=1.6)
>>
> 
> In this way can I get correct figure? My doubt is in the resulting
figure, gamma quantile are in X-axis and emperical are in Y-axis,
however the values in the Y-axis are up to 15. But my maximum point in
data set is 9.68. How can it be correct?

As far as I can see, it is not your data that you are plotting, but 
rgamma(......).


-- 
    O__  ---- Peter Dalgaard             Øster Farimagsgade 5, Entr.B
   c/ /'_ --- Dept. of Biostatistics     PO Box 2099, 1014 Cph. K
  (*) \(*) -- University of Copenhagen   Denmark      Ph:  (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)              FAX: (+45) 35327907