データフレーム内の値の範囲の IQR を見つけたいと思っています。これらの値もグループ化されているため、データフレーム内の各グループの IQR を見つける必要があります。次の表があります。
Block DNAname Spot_Size Molarity Cy3_Fluorescence
1 DNA 01 100pl 100 14266
1 DNA 01 100pl 100 16020
1 DNA 01 100pl 100 15705
1 DNA 01 100pl 100 15783
1 DNA 01 100pl 100 15834
1 DNA 01 100pl 50 12248
1 DNA 01 100pl 50 12209
1 DNA 01 100pl 50 12511
1 DNA 01 100pl 50 12316
1 DNA 01 100pl 50 12469
1 DNA 01 100pl 25 9626
1 DNA 01 100pl 25 9804
1 DNA 01 100pl 25 9794
1 DNA 01 100pl 25 10020
1 DNA 01 100pl 25 9739
1 DNA 01 100pl 10 7158
1 DNA 01 100pl 10 6802
1 DNA 01 100pl 10 7378
1 DNA 01 100pl 10 5949
1 DNA 01 100pl 10 7484
1 DNA 01 100pl 5 5257
1 DNA 01 100pl 5 5560
1 DNA 01 100pl 5 6076
1 DNA 01 100pl 5 5925
次のコードを実行して、IQR を見つけます。
aggregate(Cy3.DNA1.100pl.1uM$Cy3_Fluorescence, list(Molarity=
Cy3.DNA1.100pl.1uM$Molarity, Spot_Size=Cy3.DNA1.100pl.1uM$Spot_Size ), IQR)
これは私に出力を与えます：
Molarity Spot_Size x
5 100pl 384
10 100pl 576
25 100pl 65
50 100pl 221
100 100pl 129
この出力はすべてのモル濃度を正しくグループ化しますが、IQR は正しくありません。上記のコードが IQR ではなく関数として平均を持っている場合、x の値 (関数値) は次のように正しいです。
Molarity Spot_Size x
5 100pl 5752.4
10 100pl 6954.2
25 100pl 9796.6
50 100pl 12350.6
100 100pl 15521.6
予想される IQRS は次のようになります。
Molarity IQR
100 324.25
50 258
25 363
10 519.5
5 400
どんな助けでも大歓迎です。モル濃度カテゴリを含むスポット サイズのグループ (スポット サイズの範囲が 100pl から 400pl の範囲) がある場合に、IQR でこの機能を実行する方法を誰かが知っている場合は、それらを聞きたいと思います。
SOさんよろしくお願いします。