[R] Howto Draw Bimodal Gamma Curve with User Supplied Parameters
Gundala Viswanath
gundalav at gmail.com
Tue Jul 29 16:03:20 CEST 2008
Hi,
Suppose I have the following vector (data points):
> x
[1] 36.0 57.3 73.3 92.0 300.4 80.9 19.8 31.4 85.8 44.9 24.6 48.0
[13] 28.0 38.3 85.2 103.6 154.4 128.5 38.3 72.4 122.7 123.1 41.8 21.7
[25] 143.6 120.2 46.6 29.2 44.8 25.0 57.3 96.4 29.4 62.9 66.4 30.0
[37] 24.1 14.8 56.6 102.4 117.5 90.4 37.2 79.6 27.8 17.1 26.6 16.3
[49] 41.4 48.9 24.1 23.3 9.9 11.5 15.0 23.6 29.3 27.0 19.2 18.7
[61] 4.1 13.0 3.3 5.5 38.2 8.5 39.6 39.2 16.1 35.3 23.3 31.5
[73] 38.8 51.5 28.4 18.8 24.1 25.4 28.8 32.8 31.0 28.8 33.3 55.5
[85] 39.2 21.0 43.7 16.3 50.6 34.6 66.3 50.5 59.4 46.7 51.9 125.6
[97] 69.8 43.7 86.8 50.6 132.4 56.0 6.1 4.9 7.1 7.1 12.8 12.1
[109] 164.2 69.3 15.6 11.4 34.3 9.2 17.6 21.7 19.2 30.7 61.1 35.8
[121] 185.8 118.4 13.0 9.6 19.1 45.2 94.5 248.0 56.3 24.4 13.8 12.8
[133] 35.0 31.6 22.5 50.1 18.7 22.1 28.3 39.5 48.2 33.1 43.5 35.1
[145] 37.4 30.3 15.8 13.9 15.3 16.1 12.7 11.4 13.0 13.8 31.5 25.3
[157] 65.2 39.5
And the following parameter set (2 component) for gamma function.
comp.1 comp.2
alpha (shape) 2.855444 2.152056
beta (scale) 10.418785 39.296224
These params are predefined/precalculated by user.
My question is how can I create a bimodal gamma curve - based on the
two parameter set -
on top of the histogram of data points above?
- Gundala Viswanath
Jakarta - Indonesia
More information about the R-help
mailing list