【R】chi^2分布でフィットして、ヒストグラムを書きたい

#fitdiftrを使うためのパッケージ
library(MASS)

# tableでデータ読み込み、1列目のみ抜き出しているが必要かはわからん
data1 <- read.table("hoge.txt")
data2 <- data1[, 1]

# chi square分布でフィットする
# chi_kには推定された自由度を代入しておく
chi_df <- fitdistr(data2, "chi-squared", start=list(df=2), method="BFGS")
chi_k <- chi_df[[1]][1]

# histではなくtruehistで確率密度関数にする
chi_hist <- truehist(data2, h=0.1, prob = TRUE, col="gray", main="hogehoge")
curve(dchisq(x, df=chi_k),add=TRUE,col="red", lwd=3)