0

以下に示すマトリックス「結果」の各行について

            A   B   C   D   E   F   G   H   I   J      
       1    4   6   3   5   9   9   9   3   4   4
       2    5   7   5   5   8   8   8   7   4   5
       3    7   5   4   4   7   9   7   4   4   5
       4    6   6   6   6   8   9   8   6   3   6
       5    4   5   5   5   8   8   7   4   3   7
       6    7   9   7   6   7   8   8   5   7   6
       7    5   6   6   5   8   8   7   3   3   5
       8    6   7   4   5   8   9   8   4   6   5
       9    6   8   8   6   7   7   7   7   6   6

以下に示すように、3 つのビンを使用して各行のヒストグラムをプロットしたいと思います。

samp<-result[1,]
hist(samp, breaks = 3, col="lightblue", border="pink")

行 1 のヒストグラム (3 ビン)

ここで必要なのは、次のようにヒストグラム頻度カウントを配列に変換することです。たとえば、4 つのビンがあり、最初のビンのカウントが 5 で、2 番目のビンのカウントが 2 で、4 番目のビンが 3 であるとします。今、これらの各ビンのすべての値のベクトルが必要です。これは、出力としてベクトル内のデータ結果 (すべての行) から得られます。

       row1  5 2 0 3

何百もの行について、自動化された方法で実行したいので、この質問を投稿しました。

最終的に、マトリックスは次のようになります

             bin 2-4 bin 4-6 bin6-8 bin8-10
      row 1   5       2       0     3
      row 2
      row 3
      row 4
      row 5
      row 6
      row 7
      row 8
      row 9
4

2 に答える 2

2 
DF <- read.table(text="A   B   C   D   E   F   G   H   I   J      
1    4   6   3   5   9   9   9   3   4   4
2    5   7   5   5   8   8   8   7   4   5
3    7   5   4   4   7   9   7   4   4   5
4    6   6   6   6   8   9   8   6   3   6
5    4   5   5   5   8   8   7   4   3   7
6    7   9   7   6   7   8   8   5   7   6
7    5   6   6   5   8   8   7   3   3   5
8    6   7   4   5   8   9   8   4   6   5
9    6   8   8   6   7   7   7   7   6   6", header=TRUE)

m <- as.matrix(DF)

apply(m,1,function(x) hist(x,breaks = 3)$count)
# $`1`
# [1] 5 2 0 3
# 
# $`2`
# [1] 5 0 2 3
# 
# $`3`
# [1] 6 3 1
# 
# $`4`
# [1] 1 6 2 1
# 
# $`5`
# [1] 3 3 4
# 
# $`6`
# [1] 3 4 2 1
# 
# $`7`
# [1] 2 5 3
# 
# $`8`
# [1] 6 3 1
# 
# $`9`
# [1] 4 4 0 2

ドキュメントによると、休憩の数は単なる提案であることに注意してください。すべての行に同じ数のブレークを持たせたい場合は、次の範囲外でビニングを行う必要がありますhist。

breaks <- 1:5*2
t(apply(m,1,function(x) table(cut(x,breaks,include.lowest = TRUE))))
#   [2,4] (4,6] (6,8] (8,10]
# 1     5     2     0      3
# 2     1     4     5      0
# 3     4     2     3      1
# 4     1     6     2      1
# 5     3     3     4      0
# 6     0     3     6      1
# 7     2     5     3      0
# 8     2     4     3      1
# 9     0     4     6      0
于 2013-07-24T14:37:07.913 に答える
1

countsによって返されるベクトルにアクセスできますhist(詳細については、を参照?histしてください)。

counts <- hist(samp, breaks = 3, col="lightblue", border="pink")$counts
于 2013-07-24T14:31:58.883 に答える