Error in `$<-.data.frame`(`tmp`, group, value = integer(0))

Question

This code works fine:

d0 = DGEList(raw)
group = as.factor(metadata$Benefit)
d0$samples$group = group
keep.exprs = edgeR::filterByExpr(d0, group = group)
d0 = d0[keep.exprs,, keep.lib.sizes = FALSE]
d0 = calcNormFactors(d0 , method = 'TMM' )

But when I try to make a function for it, it doesn't work:

limmaNorm = function(meta_data, dataset, column){
  
  d = DGEList(dataset)
  group = as.factor(meta_data$column)
  d$samples$group = group
  keep.exprs = edgeR::filterByExpr(d, group = group)
  d = d[keep.exprs,, keep.lib.sizes = FALSE]
  d = calcNormFactors(d , method = 'TMM' )
  return(d)
  
}

Calling for the function:

d0 = limmaNorm(metadata, raw, 'Benefit')

Gives the following error:

Error in `$<-.data.frame`(`*tmp*`, group, value = integer(0)) : 
  replacement has 0 rows, data has 49

This is a sample of the raw datarame:

structure(list(Pt1 = c(436L, 0L, 565L, 67L, 6824L, 239L, 1045L
), Pt10 = c(306L, 0L, 980L, 338L, 21827L, 192L, 706L), Pt101 = c(279L, 
0L, 1630L, 228L, 2485L, 74L, 1189L), Pt103 = c(1124L, 0L, 1069L, 
23L, 1337L, 129L, 484L), Pt106 = c(290L, 3L, 800L, 3108L, 9424L, 
159L, 1184L), Pt11 = c(108L, 0L, 429L, 69L, 5105L, 182L, 907L
), Pt17 = c(178L, 0L, 244L, 36L, 1481L, 151L, 1049L)), row.names = c("A1BG", 
"NAT2", "ADA", "CDH2", "AKT3", "ZBTB11-AS1", "MED6"), class = "data.frame")

And this is the metadata dataframe:

structure(list(Cohort = c("NIV3-PROG", "NIV3-NAIVE", "NIV3-PROG", 
"NIV3-PROG", "NIV3-PROG", "NIV3-NAIVE", "NIV3-PROG"), Response = c("PD", 
"SD", "PR", "PD", "PD", "PD", "PD"), `Dead/Alive
(Dead = True)` = c(TRUE, 
TRUE, FALSE, FALSE, TRUE, TRUE, TRUE), `Time to Death
(weeks)` = c(22.85714286, 
36.57142857, 119.1428571, 69.14285714, 13, 119.5714286, 8.142857143
), Subtype = c("CUTANEOUS", "CUTANEOUS", "CUTANEOUS", "CUTANEOUS", 
"MUCOSAL", "CUTANEOUS", "CUTANEOUS"), `Mutational
 Subtype` = c("NA", 
"NF1", "TripleWt", "TripleWt", "BRAF", "BRAF", "TripleWt"), `M Stage` = c("M1C", 
"M1A", "M1A", "M1B", "M1C", "NA", "M1C"), `Mutation Load` = c("NA", 
"75", "10", "21", "700", "106", "20"), `Neo-antigen Load` = c("NA", 
"33", "5", "5", "219", "67", "13"), `Neo-peptide Load` = c("NA", 
"56", "6", "11", "273", "187", "32"), `Cytolytic Score` = c("977.86911190000001", 
"65.840716889999996", "1392.1422339999999", "1108.8620289999999", 
"645.54163300000005", "602.6740413", "20.904544959999999"), Benefit = c("NoResponse", 
"NoResponse", "Response", "NoResponse", "NoResponse", "NoResponse", 
"NoResponse")), row.names = c("Pt1", "Pt10", "Pt101", "Pt103", 
"Pt106", "Pt11", "Pt17"), class = "data.frame")

What can cause a code to fail when putting it in a function?

@Basti apologies, I didn't include all columns in the sample. I edited it. — Programming Noob, Aug 12 '22 at 11:19

score 0 · Accepted Answer · answered Aug 12 '22 at 11:31

The issue comes from the line group = as.factor(meta_data$column).

meta_data$column is searching for a column called column. Indeed $ does not allow computed indices and you need to switch to an alternative.

Instead you should do : group = as.factor(meta_data[[column]]) as suggested in the "Recursive (list-like) objects" in ?Extract documentation.

Error in `$<-.data.frame`(`*tmp*`, group, value = integer(0))

1 Answers1

Error in `$<-.data.frame`(`tmp`, group, value = integer(0))