hlo1 <- c(beta1hat-sigma*solve(t(x1)%*%x1)*qnorm(0.025)) Error in sigma * solve(t(x1) %*% x1) : non-numeric argument to binary operator >

这个错误通常表示sigma的类型不是数值型，无法与其他数值型变量进行乘法运算。请检查sigma的类型，确保它是数值型变量，并且已经被成功计算出来。如果sigma是一个矩阵或数组，可以使用as.numeric()函数将其转换为数值型变量。例如，可以在计算hlo1之前先将sigma转换为数值型变量：

sigma <- as.numeric(sqrt(sum((y - x1 %*% beta1hat - x2 %*% beta2hat)^2) / (n - 2 - 1)))
hlo1 <- c(beta1hat - sigma * solve(t(x1) %*% x1) * qnorm(0.025))

这样就可以避免这个错误了。