Добавить средства столбцов и строк означает матрицу в г

Question

У меня есть новичка вопрос:

Example of a dataset: 
    dat = matrix(1:6,3,2) #create a dataset 
#  [,1] [,2] 
#[1,] 1 4 
#[2,] 2 5 
#[3,] 3 6 

    rmean = rowMeans(dat) # calculate the row means of dat 
#[1] 2.5 3.5 4.5 
    cmean = colMeans(dat) # calculate the colum means of dat 
#[1] 2 5 

Я хотел бы матрицу, как следующие:

 [,1] [,2] 
[1,] 5.5 11.5 
[2,] 7.5 13.5 
[3,] 9.5 15.5 

Первая строка рассчитывается как следующие:

1+rowmeans(2.5) + colmeans(2) = 5.5 
    4+rowmeans(2.5) + colmeans(5) = 11.5 

Это то, что я пробовал:

newdat = matrix(NA, 3,2) 
    for (row in 1:nrow(dat)) 
    for (col in 1:ncol(dat)) { 
    rmean = rowMeans(dat) 
    cmean = colMeans(dat) 
    newdat[row,col] = dat[row,] + rmean[row] + cmean[col] 
    } 

Любая помощь будет оценена и, пожалуйста, исправьте мой цикл.

Answer 1

Существует множество функций для работы с матрицами в R; обычно это вопрос поиска того, кого вы хотите. В этом случае sweep является то, что вы хотите:

sweep(sweep(dat, 1, rowMeans(dat), `+`), 2, colMeans(dat), `+`) 
#  [,1] [,2] 
# [1,] 5.5 11.5 
# [2,] 7.5 13.5 
# [3,] 9.5 15.5 

• Второй аргумент наценка, как apply, где 1 строки и 2 является колонки;
• Третий - это вектор, который должен управляться поперек строки/столбца; и
• Четвертый - операция, которая над ними работает. Вам нужно обратные ссылки около +, чтобы передать его как имя.

Выяснить, как организовать функции, может быть боль, но как только вы это сделаете, они весьма полезны. К сожалению, справочная страница для некоторых более сложных из них трудно понять, но here's some useful reading объясняет это немного лучше.

---

В качестве альтернативы, вы можете получить нехитрый и воспользоваться вектором переработки:

t(t(dat + rowMeans(dat)) + colMeans(dat)) 
#  [,1] [,2] 
# [1,] 5.5 11.5 
# [2,] 7.5 13.5 
# [3,] 9.5 15.5 

Answer 2

Нет необходимости для цикла. Преимущество R заключается в использовании операций Matrix/Vector, чтобы сделать то, что вы пытаетесь достичь.

dat = matrix(1:6,3,2) #create a dataset 

rmean = rowMeans(dat) 
cmean = colMeans(dat) 

newdat = dat + rmean 

newdat = t(apply(newdat, 1, function(x) x + cmean)) 

Answer 3

Это может быть сделан просто воспользовавшись тем, что + действует на матрицы элемента с помощью элемента, без необходимости использовать for петли или sweep. Опирается на рециркуляцию, чтобы заполнить матрицы с помощью столбцов и строк.

dat + matrix(rowMeans(dat),3,2) + matrix(colMeans(dat),3,2, byrow = T) 
##  [,1] [,2] 
##[1,] 5.5 11.5 
##[2,] 7.5 13.5 
##[3,] 9.5 15.5 

Answer 4

Мы можем использовать функцию row и col сделать репликацию

dat+rowMeans(dat)[row(dat)] + colMeans(dat)[col(dat)] 
#  [,1] [,2] 
#[1,] 5.5 11.5 
#[2,] 7.5 13.5 
#[3,] 9.5 15.5