如何在R语言中提取满足≥50%列值>0的数据框行？

Hey there! Let's work through this problem together. Your current code uses all(x>0) to grab rows where every value is positive, but we need to tweak it to check the proportion of non-zero values instead—exactly what you're asking for.

First, let's align on the logic: we want rows where at least 50% of the sample columns (sample1 to sample4) have values greater than 0. Here's how to implement that:

Step-by-Step Breakdown #

1. Isolate numeric sample columns: We don't want to include the genename column in our calculation, so we'll target just the sample columns with df[, -1].
2. Calculate non-zero proportion per row: For each row, count how many sample columns have values >0, then divide by the total number of sample columns to get the proportion.
3. Filter rows: Keep only rows where this proportion meets or exceeds 50%.

Full Working Code #

# Define your original data frame
df <- data.frame(
  genename = c("A","B","C","D"),
  sample1 = c(10,0,50,0),
  sample2 = c(0,30,0,70),
  sample3 = c(50,0,0,30),
  sample4 = c(0,0,0,10)
)

# Filter rows where at least 50% of sample columns are >0
df_filtered <- df[apply(df[, -1], MARGIN = 1, function(x) {
  # Calculate the proportion of non-zero values in the row
  prop_non_zero <- sum(x > 0) / length(x)
  # Return TRUE if the proportion is 50% or higher
  prop_non_zero >= 0.5
}), ]

# Check the result
df_filtered

What This Code Does #

• df[, -1] selects all columns except the first (genename), so we only work with the numeric sample data.
• sum(x > 0) counts how many values in the row are greater than 0.
• length(x) gives the total number of sample columns (4 in this case).
• The function returns TRUE for rows that meet the 50% threshold, which we use to subset the original data frame.

Expected Result #

When you run this, you'll get exactly the rows you want—genename A and D:

genename sample1 sample2 sample3 sample4
1        A      10       0      50       0
4        D       0      70      30      10

If you ever need to adjust the threshold (like filtering for 70% non-zero columns instead), just change the 0.5 in the condition to your desired proportion.

内容的提问来源于stack exchange，提问作者zoe