gamma(lambda,k)
近似モデルの PDF をヒストグラムに重ね合わせたいと考えています。私は書きます :
hist(pressure)
curve(dgamma(x, lambda, k), add=TRUE, col="red")
しかし、値「x」が何であるかについて混乱しています。誰でも助けてください。
gamma(lambda,k)
近似モデルの PDF をヒストグラムに重ね合わせたいと考えています。私は書きます :
hist(pressure)
curve(dgamma(x, lambda, k), add=TRUE, col="red")
しかし、値「x」が何であるかについて混乱しています。誰でも助けてください。
x <- rgamma(100,2,1) #sample
h <- hist(x, plot=FALSE) #generate hist
plot(h, col="grey") #plot hist
xlines <-seq(min(h$breaks),max(h$breaks),length.out=100) #seq of x for pdf
lines(x = xlines,y=dgamma(xlines,2,1) *length(x)*diff(h$breaks)[1])
x
PDF の値を取得するサポートの値のベクトルです。比べてみてもいいかもしれません
plot(dgamma(1:20, shape=1))
http://en.wikipedia.org/wiki/Gamma_distribution (パラメーター 1)の最初のプロットで