A long dataframe can be reduced by mergeing certain rows into a
single one. These new variables are constructed as a SimpleList
containing all the original values. Invariant columns, i.e columns
that have the same value along all the rows that need to be
merged, can be shrunk into a new variables containing that
invariant value (rather than in list columns). The grouping of
rows, i.e. the rows that need to be shrunk together as one, is
defined by a vector.
The opposite operation is expand. But note that for a
DataFrame to be expanded back, it must not to be simplified.
reduceDataFrame(x, k, count = FALSE, simplify = TRUE, drop = FALSE)
expandDataFrame(x, k = NULL)The DataFrame to be reduced or expanded.
A ‘vector’ of length nrow(x) defining the grouping
based on which the DataFrame will be shrunk.
logical(1) specifying of an additional column
(called by default .n) with the tally of rows shrunk into on
new row should be added. Note that if already existing, .n
will be silently overwritten.
A logical(1) defining if invariant columns
should be converted to simple lists. Default is TRUE.
A logical(1) specifying whether the non-invariant
columns should be dropped altogether. Default is FALSE.
An expanded (reduced) DataFrame.
Missing values do have an important effect on reduce. Unless all
values to be reduces are missing, they will result in an
non-invariant column, and will be dropped with drop = TRUE. See
the example below.
The presence of missing values can have side effects in higher
level functions that rely on reduction of DataFrame objects.
library("IRanges")
k <- sample(100, 1e3, replace = TRUE)
df <- DataFrame(k = k,
x = round(rnorm(length(k)), 2),
y = seq_len(length(k)),
z = sample(LETTERS, length(k), replace = TRUE),
ir = IRanges(seq_along(k), width = 10),
r = Rle(sample(5, length(k), replace = TRUE)),
invar = k + 1)
df
#> DataFrame with 1000 rows and 7 columns
#> k x y z ir r invar
#> <integer> <numeric> <integer> <character> <IRanges> <Rle> <numeric>
#> 1 93 -2.94 1 D 1-10 2 94
#> 2 25 0.94 2 F 2-11 1 26
#> 3 48 0.31 3 E 3-12 5 49
#> 4 59 -0.57 4 K 4-13 5 60
#> 5 8 -0.78 5 Z 5-14 1 9
#> ... ... ... ... ... ... ... ...
#> 996 19 0.58 996 U 996-1005 5 20
#> 997 44 0.35 997 P 997-1006 1 45
#> 998 46 0.10 998 S 998-1007 2 47
#> 999 7 -0.15 999 D 999-1008 1 8
#> 1000 39 0.02 1000 Q 1000-1009 2 40
## Shinks the DataFrame
df2 <- reduceDataFrame(df, df$k)
df2
#> DataFrame with 100 rows and 7 columns
#> k x y z
#> <integer> <NumericList> <IntegerList> <CharacterList>
#> 1 1 0.05,-0.82,-0.32,... 98,136,156,... T,P,S,...
#> 2 2 -0.92,-1.81,-0.92,... 264,328,432,... F,Z,R,...
#> 3 3 -0.88, 0.80, 0.87,... 7,105,263,... N,O,T,...
#> 4 4 1.36,0.46,1.39,... 92,117,160,... A,W,V,...
#> 5 5 -0.99, 1.29, 1.64,... 60,75,103,... L,G,X,...
#> ... ... ... ... ...
#> 96 96 1.53,0.60,0.82,... 303,344,356,... T,Q,N,...
#> 97 97 -2.78, 1.32,-1.70,... 159,244,425,... A,Z,X,...
#> 98 98 -0.28, 0.92,-0.92,... 33,74,283,... G,L,V,...
#> 99 99 -0.66,-1.10,-0.80,... 179,267,309,... L,Y,O,...
#> 100 100 -0.97, 0.38, 0.57,... 37,64,304,... G,D,P,...
#> ir r invar
#> <IRangesList> <RleList> <numeric>
#> 1 98-107,136-145,156-165,... 5,5,4,... 2
#> 2 264-273,328-337,432-441,... 5,4,4,... 3
#> 3 7-16,105-114,263-272,... 2,4,3,... 4
#> 4 92-101,117-126,160-169,... 1,2,5,... 5
#> 5 60-69,75-84,103-112,... 5,3,5,... 6
#> ... ... ... ...
#> 96 303-312,344-353,356-365,... 1,5,4,... 97
#> 97 159-168,244-253,425-434,... 1,4,1,... 98
#> 98 33-42,74-83,283-292,... 1,3,1,... 99
#> 99 179-188,267-276,309-318,... 3,2,4,... 100
#> 100 37-46,64-73,304-313,... 5,2,1,... 101
## With a tally of the number of members in each group
reduceDataFrame(df, df$k, count = TRUE)
#> DataFrame with 100 rows and 8 columns
#> k x y z
#> <integer> <NumericList> <IntegerList> <CharacterList>
#> 1 1 0.05,-0.82,-0.32,... 98,136,156,... T,P,S,...
#> 2 2 -0.92,-1.81,-0.92,... 264,328,432,... F,Z,R,...
#> 3 3 -0.88, 0.80, 0.87,... 7,105,263,... N,O,T,...
#> 4 4 1.36,0.46,1.39,... 92,117,160,... A,W,V,...
#> 5 5 -0.99, 1.29, 1.64,... 60,75,103,... L,G,X,...
#> ... ... ... ... ...
#> 96 96 1.53,0.60,0.82,... 303,344,356,... T,Q,N,...
#> 97 97 -2.78, 1.32,-1.70,... 159,244,425,... A,Z,X,...
#> 98 98 -0.28, 0.92,-0.92,... 33,74,283,... G,L,V,...
#> 99 99 -0.66,-1.10,-0.80,... 179,267,309,... L,Y,O,...
#> 100 100 -0.97, 0.38, 0.57,... 37,64,304,... G,D,P,...
#> ir r invar .n
#> <IRangesList> <RleList> <numeric> <integer>
#> 1 98-107,136-145,156-165,... 5,5,4,... 2 12
#> 2 264-273,328-337,432-441,... 5,4,4,... 3 10
#> 3 7-16,105-114,263-272,... 2,4,3,... 4 11
#> 4 92-101,117-126,160-169,... 1,2,5,... 5 9
#> 5 60-69,75-84,103-112,... 5,3,5,... 6 15
#> ... ... ... ... ...
#> 96 303-312,344-353,356-365,... 1,5,4,... 97 10
#> 97 159-168,244-253,425-434,... 1,4,1,... 98 10
#> 98 33-42,74-83,283-292,... 1,3,1,... 99 10
#> 99 179-188,267-276,309-318,... 3,2,4,... 100 9
#> 100 37-46,64-73,304-313,... 5,2,1,... 101 9
## Much faster, but more crowded result
df3 <- reduceDataFrame(df, df$k, simplify = FALSE)
df3
#> DataFrame with 100 rows and 7 columns
#> k x y z
#> <IntegerList> <NumericList> <IntegerList> <CharacterList>
#> 1 1,1,1,... 0.05,-0.82,-0.32,... 98,136,156,... T,P,S,...
#> 2 2,2,2,... -0.92,-1.81,-0.92,... 264,328,432,... F,Z,R,...
#> 3 3,3,3,... -0.88, 0.80, 0.87,... 7,105,263,... N,O,T,...
#> 4 4,4,4,... 1.36,0.46,1.39,... 92,117,160,... A,W,V,...
#> 5 5,5,5,... -0.99, 1.29, 1.64,... 60,75,103,... L,G,X,...
#> ... ... ... ... ...
#> 96 96,96,96,... 1.53,0.60,0.82,... 303,344,356,... T,Q,N,...
#> 97 97,97,97,... -2.78, 1.32,-1.70,... 159,244,425,... A,Z,X,...
#> 98 98,98,98,... -0.28, 0.92,-0.92,... 33,74,283,... G,L,V,...
#> 99 99,99,99,... -0.66,-1.10,-0.80,... 179,267,309,... L,Y,O,...
#> 100 100,100,100,... -0.97, 0.38, 0.57,... 37,64,304,... G,D,P,...
#> ir r invar
#> <IRangesList> <RleList> <NumericList>
#> 1 98-107,136-145,156-165,... 5,5,4,... 2,2,2,...
#> 2 264-273,328-337,432-441,... 5,4,4,... 3,3,3,...
#> 3 7-16,105-114,263-272,... 2,4,3,... 4,4,4,...
#> 4 92-101,117-126,160-169,... 1,2,5,... 5,5,5,...
#> 5 60-69,75-84,103-112,... 5,3,5,... 6,6,6,...
#> ... ... ... ...
#> 96 303-312,344-353,356-365,... 1,5,4,... 97,97,97,...
#> 97 159-168,244-253,425-434,... 1,4,1,... 98,98,98,...
#> 98 33-42,74-83,283-292,... 1,3,1,... 99,99,99,...
#> 99 179-188,267-276,309-318,... 3,2,4,... 100,100,100,...
#> 100 37-46,64-73,304-313,... 5,2,1,... 101,101,101,...
## Drop all non-invariant columns
reduceDataFrame(df, df$k, drop = TRUE)
#> DataFrame with 100 rows and 2 columns
#> k invar
#> <integer> <numeric>
#> 1 1 2
#> 2 2 3
#> 3 3 4
#> 4 4 5
#> 5 5 6
#> ... ... ...
#> 96 96 97
#> 97 97 98
#> 98 98 99
#> 99 99 100
#> 100 100 101
## Missing values
d <- DataFrame(k = rep(1:3, each = 3),
x = letters[1:9],
y = rep(letters[1:3], each = 3),
y2 = rep(letters[1:3], each = 3))
d
#> DataFrame with 9 rows and 4 columns
#> k x y y2
#> <integer> <character> <character> <character>
#> 1 1 a a a
#> 2 1 b a a
#> 3 1 c a a
#> 4 2 d b b
#> 5 2 e b b
#> 6 2 f b b
#> 7 3 g c c
#> 8 3 h c c
#> 9 3 i c c
## y is invariant and can be simplified
reduceDataFrame(d, d$k)
#> DataFrame with 3 rows and 4 columns
#> k x y y2
#> <integer> <CharacterList> <character> <character>
#> 1 1 a,b,c a a
#> 2 2 d,e,f b b
#> 3 3 g,h,i c c
## y isn't not dropped
reduceDataFrame(d, d$k, drop = TRUE)
#> DataFrame with 3 rows and 3 columns
#> k y y2
#> <integer> <character> <character>
#> 1 1 a a
#> 2 2 b b
#> 3 3 c c
## BUT with a missing value
d[1, "y"] <- NA
d
#> DataFrame with 9 rows and 4 columns
#> k x y y2
#> <integer> <character> <character> <character>
#> 1 1 a NA a
#> 2 1 b a a
#> 3 1 c a a
#> 4 2 d b b
#> 5 2 e b b
#> 6 2 f b b
#> 7 3 g c c
#> 8 3 h c c
#> 9 3 i c c
## y isn't invariant/simplified anymore
reduceDataFrame(d, d$k)
#> DataFrame with 3 rows and 4 columns
#> k x y y2
#> <integer> <CharacterList> <CharacterList> <character>
#> 1 1 a,b,c NA,a,a a
#> 2 2 d,e,f b,b,b b
#> 3 3 g,h,i c,c,c c
## y now gets dropped
reduceDataFrame(d, d$k, drop = TRUE)
#> DataFrame with 3 rows and 2 columns
#> k y2
#> <integer> <character>
#> 1 1 a
#> 2 2 b
#> 3 3 c