@@ -1291,22 +1291,22 @@ output_table$LRT.p.value <- res_lrt@listData$padj
1291
1291
tcounts <- as.data.frame(t(normCounts ))
1292
1292
tcounts $ group <- sampleTable $ condition [match(rownames(tcounts ), rownames(sampleTable ))]
1293
1293
1294
- # Aggregate group means and standard deviations
1294
+ # ## Aggregate group means and standard deviations ###
1295
1295
agg_means <- aggregate(. ~ group , data = tcounts , FUN = mean , na.rm = TRUE )
1296
1296
agg_stdev <- aggregate(. ~ group , data = tcounts , FUN = sd , na.rm = TRUE )
1297
1297
1298
- # Save group names
1298
+ # ## Save group names ###
1299
1299
group_names <- agg_means $ group
1300
1300
1301
- # Remove the 'group' column and transpose to match expected structure
1301
+ # ## Remove the 'group' column and transpose to match expected structure ###
1302
1302
group_means <- as.data.frame(t(agg_means [- 1 ]))
1303
1303
group_stdev <- as.data.frame(t(agg_stdev [- 1 ]))
1304
1304
1305
1305
# Assign column names based on group names
1306
1306
colnames(group_means ) <- paste0(" Group.Mean_(" , group_names , " )" )
1307
1307
colnames(group_stdev ) <- paste0(" Group.Stdev_(" , group_names , " )" )
1308
1308
1309
- # Add computed group means and standard deviations to output_table
1309
+ # ## Add computed group means and standard deviations to output_table ###
1310
1310
output_table <- cbind(output_table , group_means , group_stdev )
1311
1311
1312
1312
# ## Read in GeneLab annotation table for the organism of interest ###
0 commit comments