Question

使用泊松分布，对数链接函数和等式将广义线性模型拟合到这些数据：

g(λi) = log λi = β1 + β2xi

其中xi = log i。首先，从第一原则开始，计算出权重矩阵W的表达式以及迭代所需的其他项等式：

XTWXb（m）= XTWz

while (int <= 10)
{
    j11 <-  sum(exp(A%*%B(b1,b2)))

    j12 <- sum(A[,2]*exp(A%*%B(b1,b2)))
    j21 <- j12
    j22 <- sum(A[,2]^2*exp(A%*%B(b1,b2)))

    XTW1 <- sum((A%*%B(b1,b2))*exp(A%*%B(b1,b2))+Y-exp(A%*%B(b1,b2)))
    XTW2 <- sum(A[,2]*((A%*%B(b1,b2))*exp(A%*%B(b1,b2))+Y-exp(A%*%B(b1,b2))))

    XTWX <- matrix(c(j11,j12,j21,j22),nrow=2,byrow=TRUE)
    XTWz <- matrix(c(XTW1,XTW2),byrow = TRUE)
    BT = solve(XTWX)%*%XTWz

    if (abs(b1-BT[1,1]) > 0.00001 |abs(b2-BT[2,1]) > 0.00001)
    {
        b1 <- BT[1,1];
        b2 <- BT[2,1];
        nit <- nit+1
    }
    else
    {
        break
    }
}

b1
b2
solve(XTWX)

solve.default（XTWX）中的错误：系统在计算上是奇异的：倒数条件数= 1.42746e-18

我该怎么办？请帮助

当我运行此行BT = solve（XTWX）％*％XTWz我收到此错误

0 个答案: