计算?2（为模型为直线的fine2定义一个函数）

我必须为给定数据计算卡方，其中包含三个变量：年，损失的分数和误差。我用过这个公式
 

其中y =损失的分数，x =年，而sigma是误差。 a和b常数。

fun.to.minimize <- function(params,x,y,z) {
p1 <- params[1]
p2 <- params[2]
output <- sum(((y- p1- x*p2)/z)^2)
return(output)
}
p1 = 20
p2 =100
solution <- optim(c(p1,p2),fun.to.minimize,x=frac$YEAR,y=frac$FRACLOST,z=frac$ERROR)
new.p1 <- solution$par[1]
print(solution)


输出为：

$par
[1] 38.52325439 -0.01525879
$value
[1] 607820.8
$counts
function gradient
41 NA
$convergence
[1] 0
$message
NULL


我是否正确实现了公式，是否有R内置函数可以交叉检查值以进行确认。