Короткий скрипт из MATLAB в R, и один результат оказался матрицей

У вас наоборот. В MATLAB .* - это матричное произведение, а * поэлементное, как * в R, так что s должно совпадать как написано. Смотрите лучше на w. c() даёт простой вектор, а не столбец, и R повторяет его по кругу по всей матрице.

Наоборот. В MATLAB * между матрицами - матричное произведение, а .* поэлементное. В R ровно противоположно. * поэлементное, а %*% - произведение.

s <- A %*% w
s
#>      [,1]
#> [1,] 1.00
#> [2,] 1.25

w как простой вектор годится. %*% считает его столбцом, когда так сходятся размеры. B совпало только потому, что A .* A и A * A оба поэлементные. Убрать точку было правильно в той строке и неправильно в другой.

Да, у меня они были перепутаны местами. Точка в MATLAB означает поэлементную операцию, поэтому R она и не нужна. А про c() всё в силе. У простого вектора нет ориентации, и %*% выбирает её сам.

Строка в R считает вполне осмысленную вещь, просто не A * w из MATLAB. R повторяет w вдоль каждого столбца, так что строка i умножается на w[i]. Это A .* w из MATLAB, где w столбец. Результат 2 на 2 верен для этого выражения.

И делает это молча, потому что длины делятся. R предупреждает о повторении вектора, только когда большая длина не кратна меньшей, а четыре элемента на два веса делятся ровно. С w из трёх элементов было бы хотя бы предупреждение. Когда веса действительно построчные, лучше sweep(A, 1, w, "*"), чтобы намерение было видно в коде.

Та же ловушка, только символов больше. A^2 в MATLAB - матричная степень, то есть A * A. A^2 в R возводит в квадрат каждый элемент. A' в MATLAB транспонирует, в R для вещественной матрицы это t(A). A \ b в MATLAB решает систему, в R это solve(A, b), и / между матрицами в MATLAB тоже решает систему, а R делит поэлементно. Пройти по скрипту и проверить каждый *, ^ и /, где по обе стороны матрица.

Да, дело было в этом, спасибо. s <- A %*% w теперь даёт 1 и 1.25. Результат приходит матрицей 2 на 1, поэтому перед таблицей он проходит через drop(). Ниже в полном скрипте нашёлся ещё A^2, и список pgr поймал его раньше, чем он стал следующей загадкой.