膨大な数の度数分布があります。中央値と四分位数を計算したいのですが、R が文句を言います。少数の場合に機能するのは次のとおりです。
> TABLE <- data.frame(DATA = c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19), F = c(48,0,192,1152,5664,23040,77952,214272,423984,558720,267840,0,0,0,0,0,0,0,0))
> summary(rep(TABLE$DAT,TABLE$F))
Min. 1st Qu. Median Mean 3rd Qu. Max.
1.000 9.000 10.000 9.397 10.000 11.000
これが、膨大な数に対して得られるものです。
> TABLE <- data.frame(DATA = c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19), F = c(240,0,1200,9600,69600,470400,2992800,17859840,98312880,489292800,2164619760,8325820800,26865302400,68711068800,128967422400,153763315200,96770419200,26824089600,2395008000))
> summary(rep(TABLE$DAT,TABLE$F))
Error in rep(TABLE$DAT, TABLE$F) : invalid 'times' argument
In addition: Warning message:
In summary(rep(TABLE$DAT, TABLE$F)) :
NAs introduced by coercion to integer range
「rep」を使用して巨大なベクトルを作成したかったので、このエラーは私を驚かせません。しかし、これを回避して中央値と四分位数を計算する方法がわかりません。