dcast with custom fun.aggregate
问题 I have data that looks like this: sample start end gene coverage X 1 10 A 5 X 11 20 A 10 Y 1 10 A 5 Y 11 20 A 10 X 1 10 B 5 X 11 20 B 10 Y 1 10 B 5 Y 11 20 B 10 I added additional columns: data$length <- (data$end - data$start + 1) data$ct_lt <- (data$length * data$coverage) I reformated my data using dcast: casted <- dcast(data, gene ~ sample, value.var = "coverage", fun.aggregate = mean) So my new data looks like this: gene X Y A 10.00000 10.00000 B 38.33333 38.33333 This is the correct