1
我有兩個data.frames,一個代表我的系統中兩個站點的定向連接,例如站點a1下面(看着行)連接到站點d1,但沒有另一種方式。而d1連接到自身和d2。r - 在主/明細層次結構中合併兩個數據幀
>connections=read.table("file1")
>connections
V1 V2 V3 V4 V5 Site
1 0 1 0 0 0 a1
2 0 1 1 0 0 d1
3 0 1 1 1 1 d2
4 1 0 0 0 0 d3
5 0 0 0 0 0 f1
什麼已經真的很難,我是寫在每個站點的亞單位而言,這連接,每個站點由5個亞基組成,下面我總結文件中看到,每個子單元ID都是唯一的。
> subunits=read.table("file2")
> subunits
Site minID maxId
1 a1 0 4
2 d1 5 9
3 d2 10 14
4 d3 15 19
5 f1 20 24
我的目標是建立由亞基連接的詳細信息表,結果應該是這個樣子
site subunit numconnections conectionids…
a1 0 5 5 6 7 8 9
a1 1 5 5 6 7 8 9
a1 2 5 5 6 7 8 9
a1 3 5 5 6 7 8 9
a1 4 5 5 6 7 8 9
d1 5 10 5 6 7 8 9 10 11 12 13 14
d1 6 10 5 6 7 8 9 10 11 12 13 14
d1 7 10 5 6 7 8 9 10 11 12 13 14
d1 8 10 5 6 7 8 9 10 11 12 13 14
d1 9 10 5 6 7 8 9 10 11 12 13 14
d2 10 20 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
d2 11 20 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
d2 12 20 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
d2 13 20 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
d2 14 20 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
d3 15 5 0 1 2 3 4
d3 16 5 0 1 2 3 4
d3 17 5 0 1 2 3 4
d3 18 5 0 1 2 3 4
d3 19 5 0 1 2 3 4
f1 20 0
f1 21 0
f1 22 0
f1 23 0
f1 24 0
…
這是file1和file2的
文件1
"V1" "V2" "V3" "V4" "V5" "Site"
"1" 0 1 0 0 0 "a1"
"2" 0 1 1 0 0 "d1"
"3" 0 1 1 1 1 "d2"
"4" 1 0 0 0 0 "d3"
"5" 0 0 0 0 0 "f1"
file2
"Site" "minID" "maxId"
"1" "a1" 0 4
"2" "d1" 5 9
"3" "d2" 10 14
"4" "d3" 15 19
"5" "f1" 20 24