23

私がこの質問をした最初の人だとは想像できませんが、まだ解決策を見つけていません (ここでも他の場所でも)。

Rで平均化したい列がいくつかあります。最小限のトリッキーな側面は、一部の列にNAが含まれていることです。

例えば：

Trait Col1 Col2 Col3
DF    23   NA   23
DG    2    2    2
DH    NA   9    9

NA を無視して、最初の 3 列のエントリを平均する Col4 を作成したいと考えています。そう：

 Trait Col1 Col2 Col3 Col4
 DF    23   NA   23   23
 DG    2    2    2    2
 DH    NA   9    9    9

理想的には、次のようなものが機能します。

data$Col4 <- mean(data$Chr1, data$Chr2, data$Chr3, na.rm=TRUE)

しかし、そうではありません。

4

2 に答える 2

33

必要ですが、重要なのは、に設定したい引数があることrowMeans()に注意してください。例えば：na.rmTRUE

> mat <- matrix(c(23,2,NA,NA,2,9,23,2,9), ncol = 3)
> mat
     [,1] [,2] [,3]
[1,]   23   NA   23
[2,]    2    2    2
[3,]   NA    9    9
> rowMeans(mat)
[1] NA  2 NA
> rowMeans(mat, na.rm = TRUE)
[1] 23  2  9

あなたの例に合わせるには：

> dat <- data.frame(Trait = c("DF","DG","DH"), mat)
> names(dat) <- c("Trait", paste0("Col", 1:3))
> dat
  Trait Col1 Col2 Col3
1    DF   23   NA   23
2    DG    2    2    2
3    DH   NA    9    9
> dat <- transform(dat, Col4 = rowMeans(dat[,-1], na.rm = TRUE))
> dat
  Trait Col1 Col2 Col3 Col4
1    DF   23   NA   23   23
2    DG    2    2    2    2
3    DH   NA    9    9    9
于 2012-09-18T23:02:48.177 に答える