所有列上的簡單z變換

Question

我剛剛爲所有列上的Matrix變換中的z變形寫了一個非常基本的函數。它看起來像這樣：

sapply(MyObject, function(x){(x-mean(x))/sd(x)})

我隨機抽查我的函數矩陣內的一些細胞，它似乎很好地工作。不過我想確認這個功能還行，因爲我對R非常陌生，在互聯網上找不到任何好的例子。

Answer 1

你應該看看scale()這是爲你做的。

您的功能接近正確;您應該將na.rm = TRUE添加到sd()和mean()函數調用中。

我會寫（如果不使用scale()）函數使用sweep()而不是sapply()。例如。

ztran <- function(x, na.rm = TRUE) { 
    mns <- colMeans(x, na.rm = na.rm) 
    sds <- apply(x, 2, sd, na.rm = na.rm) 
    x <- sweep(x, 2, mns, "-") 
    x <- sweep(x, 2, sds, "/") 
    x 
} 

在使用中，我們有

> df <- data.frame(matrix(1:9, ncol = 3)) 
> ztran(df) 
    X1 X2 X3 
1 -1 -1 -1 
2 0 0 0 
3 1 1 1 
> scale(df) 
    X1 X2 X3 
[1,] -1 -1 -1 
[2,] 0 0 0 
[3,] 1 1 1 
attr(,"scaled:center") 
X1 X2 X3 
2 5 8 
attr(,"scaled:scale") 
X1 X2 X3 
1 1 1 

sweep是這種操作的一個非常有用的向量化工具。還要注意，sapply()簡化爲矩陣，這可能不是你想要的。 sweep()不這樣做：

> class(ztran(df)) 
[1] "data.frame" 
> class(sapply(df, function(x){(x-mean(x))/sd(x)})) 
[1] "matrix"