したがって、整数のベクトル がありquotes
、データ ポイントの頻度をプロットし、x 軸と y 軸の両方を対数にすることによって、べき乗則分布を観察するかどうかを確認したいと考えています。ただし、Rでこれを達成する方法がよくわかりません。現在、次を使用してヒストグラムを作成できます
hist(quotes, breaks = max(quotes))
しかし、軸はすべて直線的です。
おそらくこれを行うためのより良い方法がありますが、これは(基本的に)機能します:
data = rnorm(1000,0,1)
r <- hist(log(data))
plot(r$breaks[-1],log(r$counts))
編集:より良い解決策:
r <- hist(data)
plot(r$breaks[-1], r$counts, log='xy', type='h')
# or alternatively:
barplot(r$counts, log="y", col="white", names.arg=r$breaks[-1])
棒グラフ バージョンには、変換された x 軸で試してみると明らかになる理由から、変換された x 軸がありません。