3

列ごとのスコアを含む非常に大きな data.frame があるとします。

例えば：

MA0001.1 AGL3 MA0003.1 TFAP2A MA0004.1 Arnt  MA0005.1 AG   MA0006.1 Arnt::Ahr
7.789524e-09  0.4012127249    3.771518e-03   1.892011e-06  0.002733200
5.032498e-07  0.0001873801    9.947449e-05   3.284222e-05  0.001367041
1.194487e-06  0.0009357406    6.943634e-05   1.589373e-05  0.002551519
4.833494e-06  0.0150703600    1.003488e-04   1.197928e-03  0.001431416
6.865040e-05  0.0000732607    3.857193e-04   5.388744e-03  0.001363706

R data.frame:

testfr<-structure(list(`MA0001.1 AGL3` = c(7.78952366977488e-09, 5.03249791215203e-07, 
1.19448739380034e-06, 4.83349413748598e-06, 6.86504034402563e-05
), `MA0003.1 TFAP2A` = c(0.401212724871542, 0.000187380067026448, 
0.000935740631438077, 0.0150703600158589, 7.32607018758816e-05
), `MA0004.1 Arnt` = c(0.00377151826447817, 9.94744903768433e-05, 
6.94363387424972e-05, 0.000100348764966112, 0.00038571926458373
), `MA0005.1 AG` = c(1.89201084302835e-06, 3.2842217133538e-05, 
1.58937284554136e-05, 0.00119792816070882, 0.00538874414923338
), `MA0006.1 Arnt::Ahr` = c(0.00273319966783363, 0.00136704060025893, 
0.00255151921946167, 0.00143141576426544, 0.00136370552325235
)), .Names = c("MA0001.1 AGL3", "MA0003.1 TFAP2A", "MA0004.1 Arnt", 
"MA0005.1 AG", "MA0006.1 Arnt::Ahr"), class = "data.frame", row.names = c(4L, 
2L, 5L, 1L, 3L))

ここで、値が最も高い列を選択し、その列を最初に配置します。したがって、1 つの列の値は同じ列名の下にとどまり、列全体がランクごとに移動する必要があります。

私は次のことを試しました：

ranked<-unlist(lapply(testfr,rank))
testranked<-testfr[ranked, ]

これにより、2259obs*459vars のデータ フレームが生成されますが、元は 5*459 でした。

testfr は、行列のリストにシーケンスをスコアリングする関数から派生した data.frame であることに注意してください! そして、そのスコアを data.frame に戻します。行はシーケンスで、列は行列です。

インデックス作成またはリストからの削除に何か問題があることは知っていますが、これを修正する方法についての手がかりがありません。どんな助けでも大歓迎です。

4

2 に答える 2

7

これはどう？

> testfr[rev(order(sapply(testfr, max, na.rm = TRUE)))]

壊す：

sapply(test.fr, max, na.rm = TRUE) # get max of each column (after removing NA)
order(.) # get the order of these values in increasing order
rev(.)   # get the reverse order so that highest value index stays first
testfr[.] # get the columns in this order back
于 2013-03-07T10:49:45.523 に答える
1

apply読みやすさのために使用します、

testfr[order(apply(testfr, 2, max, na.rm = TRUE),decreasing=T)]

ここで各マージン、列にmaxを適用し、次に列を降順で並べ替えます。

于 2013-03-07T12:15:54.617 に答える