如何按类别变量计算连续变量均值？R语言绘图报错求助

Hey there, let's figure out why you're getting that frustrating error and get your barplot working. That error message is telling us the aggregated mean values you're trying to plot don't contain any finite numbers—so either all your calculated means are NA/NaN, or there's a problem with how your duration data is structured. Here's how to troubleshoot and fix it step by step:

Step 1: Validate your duration variable's type and missing values #

First, let's make sure your duration..seconds. column is actually a numeric type (you can't calculate a mean on character data!). Run this check:

# Check variable class
class(alien$duration..seconds.)

# Count missing values
sum(is.na(alien$duration..seconds.))

• If the class is character, you'll need to convert it to numeric first. Note that any non-numeric characters in the column will turn into NA, so you'll get a warning about coercion—this is normal, but it tells you which entries are invalid:

  alien$duration..seconds. <- as.numeric(alien$duration..seconds.)
  

• If the missing value count is extremely high (close to the total number of rows), you might need to check your data source for incomplete records.

Step 2: Inspect your aggregated results #

After running your aggregate() call, print the result to see what you're working with:

a <- aggregate(duration..seconds. ~ shape, data=alien, FUN=mean, na.rm=TRUE)
print(a)

Look for rows where duration..seconds. is NaN—this happens when a shape category has no non-missing duration values (even with na.rm=TRUE, you can't compute a mean from zero valid points).

To confirm this, calculate the number of valid duration entries per shape:

aggregate(duration..seconds. ~ shape, data=alien, FUN=function(x) sum(!is.na(x)))

Step 3: Filter out invalid groups and re-plot #

Once you identify shape categories with no valid mean values, filter them out from your aggregated data frame:

# Keep only rows with finite mean values
a_clean <- a[is.finite(a$duration..seconds.), ]

Now you can safely create your barplot. If you have lots of shape categories, adding las=2 will rotate the x-axis labels so they're readable:

barplot(a_clean$duration..seconds., names.arg=a_clean$shape, las=2, main="Mean UFO Sighting Duration by Shape", ylab="Average Duration (Seconds)")

That should resolve the error and give you a usable barplot of average sighting durations per UFO shape!

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