用R中的data.frame的平均值計算

Question

我是R的初學者我想用下面的代碼計算組間差異。

calcBetweenGroupsVariance <- function(variable,groupvariable) 
{ 
# find out how many values the group variable can take 
groupvariable2 <- as.factor(groupvariable[[1]]) 
levels <- levels(groupvariable2) 
numlevels <- length(levels) 
# calculate the overall grand mean: 
grandmean <- mean(variable) 
# get the mean and standard deviation for each group: 
numtotal <- 0 
denomtotal <- 0 
for (i in 1:numlevels) 
{ 
    leveli <- levels[i] 
    levelidata <- variable[groupvariable==leveli,] 
    levelilength <- length(levelidata) 
    # get the mean and standard deviation for group i: 
    meani <- mean(levelidata) 
    sdi <- sd(levelidata) 
    numi <- levelilength * ((meani - grandmean)^2) 
    denomi <- levelilength 
    numtotal <- numtotal + numi 
    denomtotal <- denomtotal + denomi 
} 
# calculate the between-groups variance 
Vb <- numtotal/(numlevels - 1) 
Vb <- Vb[[1]] 
return(Vb) 
} 

但是，我發現了以下錯誤，同時使用該功能，

calcBetweenGroupsVariance (data[3],data[2]) 

警告消息：在mean.default（變量）：參數不是數字 或邏輯：返回NA

我明白使用均值函數時會出錯。

這裏是str(data)

'data.frame': 45 obs. of 11 variables: 
$ V1 : int 2 3 3 2 3 2 2 2 3 2 ... 
$ V2 : num 1.3243 -2.4546 0.1352 0.0676 -1.1901 ... 
$ V3 : num 0.913 -2.644 0.663 1.217 -0.409 ... 
$ V4 : num -1.863 1.965 -0.698 -0.945 0.617 ... 
$ V5 : num -0.574 1.031 -0.308 -0.574 0.354 ... 
$ V6 : num -0.8963 2.5702 0.0736 -1.3671 0.9045 ... 
$ V7 : num 0.2276 0.0624 0.5945 0.6194 0.5473 ... 
$ V8 : num 1.304 -1.624 0.408 0.368 -0.559 ... 
$ V9 : num -0.1827 -0.9748 -0.5158 -0.0191 -0.3053 ... 
$ V10: num -0.964 0.67 -0.12 0.789 0.711 ... 
$ V11: num -0.833 -0.833 -0.833 -0.0539 -0.0539 ... 

輸出請建議如何擺脫這種錯誤的。

感謝和問候

Answer 1

有在你的腳本複式錯誤有關陣列的尺寸和bewteen一個向量和列表的區別

假設參數variable, groupvariable你的功能應該是向量/ 1D -arrays。

• 線groupvariable2 <- as.factor(groupvariable[[1]])應groupvariable2 <- as.factor(groupvariable)因爲groupvariable不是一個列表，並有啊再沒有在第一個元素，但都只是感興趣。

• levelidata <- variable[groupvariable==leveli,]應該levelidata <- variable[groupvariable==leveli]因爲variable只有一個尺寸線（它不是一個矩陣）

• 給你的函數調用應calcBetweenGroupsVariance(data[[3]], data[[2]])（帶雙括號[[]]）），或者calcBetweenGroupsVariance(data[, 3],data[, 2])否則你會傳遞一個列表而不是一個向量到函數。