Я пытаюсь взять производную функции
C (β) = 1/2 (Y-X * β) ^ 2, в которой Y - это матрица 5 * 2, Y - это 5* 1 матрица с известными значениями. β - вектор 2 значений;Я хочу найти производную от C (β).
Я пробовал что-то вроде следующего:
X & Y
X = as.matrix(data.frame(x1=c(1,3,5,6,9),x2=c(2,4,4,10,6)))
Y = as.matrix(c(0,2,6,2,12))
войти в функцию
formula2 <- function(z) 1/2 * ((Y-(X%*%z))^2)
взять производную функции
deriv <- Deriv(formula2, 'z')
Это даст мнеОшибка Не удалось получить тело «% % ()», но я узнал, что% % предназначен для умножения матриц, поэтому я думаю, что мне нужно это использовать.