我试图获得总数量的变量的汇总数据帧prop.damage和crop.damage通过STATE变量使用aggregate()函数在用下面的代码R：
stormdata$prop.damage <- with(stormdata, ifelse(PROPDMGEXP == 'K', (PROPDMG * 10^3), ifelse(PROPDMGEXP == 'M', (PROPDMG * 10^6), ifelse(PROPDMGEXP == 'B', (PROPDMG * 10^9), NA))))
stormdata$crop.damage <- with(stormdata, ifelse(CROPDMGEXP == 'K', (CROPDMG * 10^3), ifelse(CROPDMGEXP == 'M', (CROPDMG * 10^6), ifelse(CROPDMGEXP == 'B', (CROPDMG * 10^9), NA))))
damagecost <- with(stormdata, aggregate(x = prop.damage + crop.damage, by = list(STATE), FUN = sum, na.rm = TRUE))
damagecost <- damagecost[order(damagecost$x, decreasing = TRUE), ]
这里的PROPDMGEXP和CROPDMGEXP变量用作乘法器PROPDMG和CROPDMG数值变量。 我的主要数据集是stormdata 。
而我得到以下几点：
> head(damagecost)
Group.1 x
8 CA 120211639720
13 FL 27302948100
38 MS 14804212820
63 TX 12550131850
20 IL 11655920860
2 AL 9505473250
但是，比如，如果我做加法“手动”为加州（“CA”）我得到这样的：
> sum(stormdata$prop.damage[stormdata$STATE == 'CA'], na.rm = TRUE) + sum(stormdata$crop.damage[stormdata$STATE == 'CA'], na.rm = TRUE)
[1] 127115859410
我不明白为什么我得到不同的结果。
