Использование ave в R без значений NA?

Question

У меня есть обширная структура данных. Это выглядит примерно так:

> b 
     fips  vix 
1400 08005 18.58862 
1401 47155 10.93712 
1402 51191 10.93712 
1403 47059 10.93712 
1404 08005 10.93712 
1405 08059 10.93712 
1406 47063 10.93712 
1407 37021 10.93712 
1408 08031 10.93712 
1409 45083 10.93712 
1410 37089 10.93712 
1411 37113 10.93712 
1412 13207 10.93712 
1413 08041 10.93712 
1414 47093 21.50425 
1415 08031 21.50425 
1416 37009 21.50425 
1417 36103 21.50425 
1418 08035 21.50425 
1419 08031 53.58363 
1420 08035 53.58363 
1421 08013 53.58363 
1422 55105 21.17450 
1423 08001 21.17450 
1424 08031 21.17450 
1425 47179 21.17450 
1426 08059 21.17450 
1427 37009 17.35675 
1428 08041 17.35675 
1429 08031 17.35675 
1430 08005 17.35675 
1431 08001  NA 
1432 08031  NA 
1433 47059  NA 
1434 47145  NA 
1435 13207  NA 
1436 37021  NA 
1437 37113  NA 
1438 37089  NA 

я взял некоторые из столбцов для простоты и показывается только часть тех строк. Я пытаюсь изменить столбец Vix. То, что я пытаюсь сделать это:

b$vix <- b$vix - ave(b$vix,b$fips) 

Что, что следует сделать, это вычесть группа средства из каждого значения Vix. Например, для наблюдения 1400 я хочу взять среднее значение всех наблюдений с fips == 08005, а затем сделать 18.58862 минус это среднее значение. Однако проблема заключается в том, что существуют значения NA. Я хочу, чтобы средняя функция IGNORE значения NA. Вместо того, что происходит в том, что любая группа ФИПС коды, который имеет один NA поворачивает вверх, как NA:

> b$vix <- b$vix - ave(b$vix,b$fips) 
> b 
     fips  vix 
1400 08005 2.961125 
1401 47155 0.000000 
1402 51191 0.000000 
1403 47059   NA 
1404 08005 -4.690375 
1405 08059 -5.118688 
1406 47063 0.000000 
1407 37021   NA 
1408 08031   NA 
1409 45083 0.000000 
1410 37089   NA 
1411 37113   NA 
1412 13207   NA 
1413 08041 -3.209812 
1414 47093 0.000000 
1415 08031   NA 
1416 37009 2.073750 
1417 36103 0.000000 
1418 08035 -16.039688 
1419 08031   NA 
1420 08035 16.039688 
1421 08013 0.000000 
1422 55105 0.000000 
1423 08001   NA 
1424 08031   NA 
1425 47179 0.000000 
1426 08059 5.118688 
1427 37009 -2.073750 
1428 08041 3.209812 
1429 08031   NA 
1430 08005 1.729250 
1431 08001   NA 
1432 08031   NA 
1433 47059   NA 
1434 47145   NA 
1435 13207   NA 
1436 37021   NA 
1437 37113   NA 
1438 37089   NA 

Как вы можете видеть, любой из ФИПСА, которые Ч теперь дает числовую апертуру для всех остальных строк с теми же fips. Я попытался добавить в na.rm = TRUE, но это ничего не делает. Я также думал о добавлении в другую функцию, то есть ave (b $ vix, b $ fips, FUN = ...), но я не знал, что добавить. Возможно, есть еще один способ сделать это вообще. Надеюсь, я смог четко объяснить проблему. Любая помощь приветствуется!

Answer 1

Вы можете просто использовать пользовательскую функцию передать na.rm=TRUE флаг mean() с

b$vix <- b$vix - ave(b$vix,b$fips, FUN=function(x) mean(x, na.rm=T)) 

Испытано с

b<-read.table(text="  fips  vix 
08005 18 
08005 19 
08005 20 
47155 10 
47155 NA 
47155 20", header=T) 

b$vix <- b$vix - ave(b$vix,b$fips, FUN=function(x) mean(x, na.rm=T)) 
b 
# fips vix 
# 1 8005 -1 
# 2 8005 0 
# 3 8005 1 
# 4 47155 -5 
# 5 47155 NA 
# 6 47155 5