R 在表达式中添加换行符
R add line breaks to expressions
我在 .txt 文件中存储了很长的表达式。通过将此表达式复制到 R 中，它在一行中表示。当我尝试 运行 它时，它不起作用并且出现“意外符号”或“意外数字常量”之类的错误。但是，如果我在任意位置插入换行符，例如每3000个字符后，它工作顺利。
有谁知道如何处理这个问题？特别是，是否可以自动使这些换行符或者是否有其他可能的解决方案？
alpha <- 0.025; d1 <- 0.027; d2 <- 0.037; mu1 <- 0.024; mu2 <- 0.026; Z1 <- 0.6; Z2 <- 0.7; rho11 <- 0.076;
rho12 <- 0.074; rho21 <- 0.073; rho22 <- 0.072; S11 <- 0.086; S12 <- 0.087; S21 <- 0.084; S22 <- 0.083; u <- 2
exp(mu1*(mu1*u*d1^2 + d2*mu2*u*d1) - log(1 - 2*S22*d2^2*u - 2*S12*d1*d2*u - 2*S21*d1*d2*u - 2*S11*d1^2*u)/2 + mu2*(mu2*u*d2^2 + d1*mu1*u*d2) - (S21*((2*(S12*u*d2^2 + S11*d1*u*d2)*(d1^2*mu1*u + alpha*d1*u + d1*d2*mu2*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((d2^2*mu2*u + alpha*d2*u + d1*d2*mu1*u)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((2*(S21*u*d1^2 + S22*d2*u*d1)*(d2^2*mu2*u + alpha*d2*u + d1*d2*mu1*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((d1^2*mu1*u + alpha*d1*u + d1*d2*mu2*u)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))*(2*d1^2*mu1*u + 2*alpha*d1*u + 2*d1*d2*mu2*u) - (S22*((2*(S12*u*d2^2 + S11*d1*u*d2)*(d1^2*mu1*u + alpha*d1*u + d1*d2*mu2*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((d2^2*mu2*u + alpha*d2*u + d1*d2*mu1*u)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((2*(S21*u*d1^2 + S22*d2*u*d1)*(d2^2*mu2*u + alpha*d2*u + d1*d2*mu1*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((d1^2*mu1*u + alpha*d1*u + d1*d2*mu2*u)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))*(2*d2^2*mu2*u + 2*alpha*d2*u + 2*d1*d2*mu1*u) + Z1*(rho11*(2*d1^2*mu1*u + 2*alpha*d1*u + 2*d1*d2*mu2*u) - rho21*(d2^2*u*(S22*((2*(S12*u*d2^2 + S11*d1*u*d2)*(4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((2*(S21*u*d1^2 + S22*d2*u*d1)*(4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S21*((2*(S12*u*d2^2 + S11*d1*u*d2)*(4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((2*(S21*u*d1^2 + S22*d2*u*d1)*(4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) - rho11*(d1^2*u*(S21*((2*(S12*u*d2^2 + S11*d1*u*d2)*(4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((2*(S21*u*d1^2 + S22*d2*u*d1)*(4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S22*((2*(S12*u*d2^2 + S11*d1*u*d2)*(4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((2*(S21*u*d1^2 + S22*d2*u*d1)*(4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) + rho21*(2*d2^2*mu2*u + 2*alpha*d2*u + 2*d1*d2*mu1*u) + 2*alpha*d1*u) + Z2*(rho12*(2*d1^2*mu1*u + 2*alpha*d1*u + 2*d1*d2*mu2*u) - rho22*(d2^2*u*(S22*((2*(S12*u*d2^2 + S11*d1*u*d2)*(4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((2*(S21*u*d1^2 + S22*d2*u*d1)*(4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S21*((2*(S12*u*d2^2 + S11*d1*u*d2)*(4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((2*(S21*u*d1^2 + S22*d2*u*d1)*(4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) - rho12*(d1^2*u*(S21*((2*(S12*u*d2^2 + S11*d1*u*d2)*(4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((2*(S21*u*d1^2 + S22*d2*u*d1)*(4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S22*((2*(S12*u*d2^2 + S11*d1*u*d2)*(4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((2*(S21*u*d1^2 + S22*d2*u*d1)*(4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) + rho22*(2*d2^2*mu2*u + 2*alpha*d2*u + 2*d1*d2*mu1*u) + 2*alpha*d2*u) + 2*alpha^2*u + Z1*(Z1*(rho11*(rho11*u*d1^2 + d2*rho21*u*d1) + rho21*(rho21*u*d2^2 + d1*rho11*u*d2) - rho11*(d1^2*u*(S21*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho11*u*d1^2 + d2*rho21*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho21*u*d2^2 + d1*rho11*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho21*u*d2^2 + d1*rho11*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho11*u*d1^2 + d2*rho21*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S22*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho11*u*d1^2 + d2*rho21*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho21*u*d2^2 + d1*rho11*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho21*u*d2^2 + d1*rho11*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho11*u*d1^2 + d2*rho21*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) - rho21*(d2^2*u*(S22*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho11*u*d1^2 + d2*rho21*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho21*u*d2^2 + d1*rho11*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho21*u*d2^2 + d1*rho11*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho11*u*d1^2 + d2*rho21*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S21*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho11*u*d1^2 + d2*rho21*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho21*u*d2^2 + d1*rho11*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho21*u*d2^2 + d1*rho11*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho11*u*d1^2 + d2*rho21*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) + d1^2*u) + Z2*(rho11*(rho12*u*d1^2 + d2*rho22*u*d1) + rho21*(rho22*u*d2^2 + d1*rho12*u*d2) - rho11*(d1^2*u*(S21*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho12*u*d1^2 + d2*rho22*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho22*u*d2^2 + d1*rho12*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho22*u*d2^2 + d1*rho12*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho12*u*d1^2 + d2*rho22*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S22*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho12*u*d1^2 + d2*rho22*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho22*u*d2^2 + d1*rho12*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho22*u*d2^2 + d1*rho12*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho12*u*d1^2 + d2*rho22*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) - rho21*(d2^2*u*(S22*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho12*u*d1^2 + d2*rho22*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho22*u*d2^2 + d1*rho12*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho22*u*d2^2 + d1*rho12*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho12*u*d1^2 + d2*rho22*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S21*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho12*u*d1^2 + d2*rho22*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho22*u*d2^2 + d1*rho12*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho22*u*d2^2 + d1*rho12*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho12*u*d1^2 + d2*rho22*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) + d1*d2*u)) + Z2*(Z2*(rho12*(rho12*u*d1^2 + d2*rho22*u*d1) + rho22*(rho22*u*d2^2 + d1*rho12*u*d2) - rho12*(d1^2*u*(S21*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho12*u*d1^2 + d2*rho22*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho22*u*d2^2 + d1*rho12*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho22*u*d2^2 + d1*rho12*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho12*u*d1^2 + d2*rho22*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S22*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho12*u*d1^2 + d2*rho22*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho22*u*d2^2 + d1*rho12*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho22*u*d2^2 + d1*rho12*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho12*u*d1^2 + d2*rho22*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) - rho22*(d2^2*u*(S22*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho12*u*d1^2 + d2*rho22*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho22*u*d2^2 + d1*rho12*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho22*u*d2^2 + d1*rho12*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho12*u*d1^2 + d2*rho22*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S21*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho12*u*d1^2 + d2*rho22*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho22*u*d2^2 + d1*rho12*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho22*u*d2^2 + d1*rho12*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho12*u*d1^2 + d2*rho22*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) + d2^2*u) + Z1*(rho12*(rho11*u*d1^2 + d2*rho21*u*d1) + rho22*(rho21*u*d2^2 + d1*rho11*u*d2) - rho12*(d1^2*u*(S21*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho11*u*d1^2 + d2*rho21*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho21*u*d2^2 + d1*rho11*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho21*u*d2^2 + d1*rho11*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho11*u*d1^2 + d2*rho21*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S22*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho11*u*d1^2 + d2*rho21*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho21*u*d2^2 + d1*rho11*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho21*u*d2^2 + d1*rho11*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho11*u*d1^2 + d2*rho21*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) - rho22*(d2^2*u*(S22*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho11*u*d1^2 + d2*rho21*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho21*u*d2^2 + d1*rho11*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho21*u*d2^2 + d1*rho11*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho11*u*d1^2 + d2*rho21*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S21*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho11*u*d1^2 + d2*rho21*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho21*u*d2^2 + d1*rho11*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho21*u*d2^2 + d1*rho11*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho11*u*d1^2 + d2*rho21*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) + d1*d2*u)) + 2*alpha*d1*mu1*u + 2*alpha*d2*mu2*u)
原因在 this question
的答案中有很好的描述
在您的情况下，由于您在文件中获得了表达式，因此您可以 source() 来自 R 的文件，而不是复制和粘贴。您需要在文本文件中添加作业，或者您可以这样做：
source("filename.txt")$value
获取计算结果
您在评论中提到您希望能够在优化中使用它 - 一个在一维上找到最小化表达式的 alpha 的简单示例是：
f <- function(alpha){source("o:\test.txt")$value}
optim(1, f, lower=0.5, upper=1.5, method="L-BFGS-B")
#$par
#[1] 1
#
#$value
#[1] 1.011915
#
#$counts
#function gradient
# 1 1
#
#$convergence
#[1] 0
#
#$message
#[1] "CONVERGENCE: NORM OF PROJECTED GRADIENT <= PGTOL"
我在 .txt 文件中存储了很长的表达式。通过将此表达式复制到 R 中，它在一行中表示。当我尝试 运行 它时，它不起作用并且出现“意外符号”或“意外数字常量”之类的错误。但是，如果我在任意位置插入换行符，例如每3000个字符后，它工作顺利。
有谁知道如何处理这个问题？特别是，是否可以自动使这些换行符或者是否有其他可能的解决方案？
alpha <- 0.025; d1 <- 0.027; d2 <- 0.037; mu1 <- 0.024; mu2 <- 0.026; Z1 <- 0.6; Z2 <- 0.7; rho11 <- 0.076;
rho12 <- 0.074; rho21 <- 0.073; rho22 <- 0.072; S11 <- 0.086; S12 <- 0.087; S21 <- 0.084; S22 <- 0.083; u <- 2
exp(mu1*(mu1*u*d1^2 + d2*mu2*u*d1) - log(1 - 2*S22*d2^2*u - 2*S12*d1*d2*u - 2*S21*d1*d2*u - 2*S11*d1^2*u)/2 + mu2*(mu2*u*d2^2 + d1*mu1*u*d2) - (S21*((2*(S12*u*d2^2 + S11*d1*u*d2)*(d1^2*mu1*u + alpha*d1*u + d1*d2*mu2*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((d2^2*mu2*u + alpha*d2*u + d1*d2*mu1*u)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((2*(S21*u*d1^2 + S22*d2*u*d1)*(d2^2*mu2*u + alpha*d2*u + d1*d2*mu1*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((d1^2*mu1*u + alpha*d1*u + d1*d2*mu2*u)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))*(2*d1^2*mu1*u + 2*alpha*d1*u + 2*d1*d2*mu2*u) - (S22*((2*(S12*u*d2^2 + S11*d1*u*d2)*(d1^2*mu1*u + alpha*d1*u + d1*d2*mu2*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((d2^2*mu2*u + alpha*d2*u + d1*d2*mu1*u)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((2*(S21*u*d1^2 + S22*d2*u*d1)*(d2^2*mu2*u + alpha*d2*u + d1*d2*mu1*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((d1^2*mu1*u + alpha*d1*u + d1*d2*mu2*u)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))*(2*d2^2*mu2*u + 2*alpha*d2*u + 2*d1*d2*mu1*u) + Z1*(rho11*(2*d1^2*mu1*u + 2*alpha*d1*u + 2*d1*d2*mu2*u) - rho21*(d2^2*u*(S22*((2*(S12*u*d2^2 + S11*d1*u*d2)*(4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((2*(S21*u*d1^2 + S22*d2*u*d1)*(4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S21*((2*(S12*u*d2^2 + S11*d1*u*d2)*(4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((2*(S21*u*d1^2 + S22*d2*u*d1)*(4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) - rho11*(d1^2*u*(S21*((2*(S12*u*d2^2 + S11*d1*u*d2)*(4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((2*(S21*u*d1^2 + S22*d2*u*d1)*(4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S22*((2*(S12*u*d2^2 + S11*d1*u*d2)*(4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((2*(S21*u*d1^2 + S22*d2*u*d1)*(4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) + rho21*(2*d2^2*mu2*u + 2*alpha*d2*u + 2*d1*d2*mu1*u) + 2*alpha*d1*u) + Z2*(rho12*(2*d1^2*mu1*u + 2*alpha*d1*u + 2*d1*d2*mu2*u) - rho22*(d2^2*u*(S22*((2*(S12*u*d2^2 + S11*d1*u*d2)*(4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((2*(S21*u*d1^2 + S22*d2*u*d1)*(4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S21*((2*(S12*u*d2^2 + S11*d1*u*d2)*(4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((2*(S21*u*d1^2 + S22*d2*u*d1)*(4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) - rho12*(d1^2*u*(S21*((2*(S12*u*d2^2 + S11*d1*u*d2)*(4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((2*(S21*u*d1^2 + S22*d2*u*d1)*(4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S22*((2*(S12*u*d2^2 + S11*d1*u*d2)*(4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((2*(S21*u*d1^2 + S22*d2*u*d1)*(4*d2^2*mu2*u + 4*alpha*d2*u + 4*d1*d2*mu1*u))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - ((4*d1^2*mu1*u + 4*alpha*d1*u + 4*d1*d2*mu2*u)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) + rho22*(2*d2^2*mu2*u + 2*alpha*d2*u + 2*d1*d2*mu1*u) + 2*alpha*d2*u) + 2*alpha^2*u + Z1*(Z1*(rho11*(rho11*u*d1^2 + d2*rho21*u*d1) + rho21*(rho21*u*d2^2 + d1*rho11*u*d2) - rho11*(d1^2*u*(S21*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho11*u*d1^2 + d2*rho21*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho21*u*d2^2 + d1*rho11*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho21*u*d2^2 + d1*rho11*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho11*u*d1^2 + d2*rho21*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S22*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho11*u*d1^2 + d2*rho21*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho21*u*d2^2 + d1*rho11*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho21*u*d2^2 + d1*rho11*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho11*u*d1^2 + d2*rho21*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) - rho21*(d2^2*u*(S22*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho11*u*d1^2 + d2*rho21*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho21*u*d2^2 + d1*rho11*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho21*u*d2^2 + d1*rho11*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho11*u*d1^2 + d2*rho21*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S21*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho11*u*d1^2 + d2*rho21*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho21*u*d2^2 + d1*rho11*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho21*u*d2^2 + d1*rho11*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho11*u*d1^2 + d2*rho21*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) + d1^2*u) + Z2*(rho11*(rho12*u*d1^2 + d2*rho22*u*d1) + rho21*(rho22*u*d2^2 + d1*rho12*u*d2) - rho11*(d1^2*u*(S21*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho12*u*d1^2 + d2*rho22*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho22*u*d2^2 + d1*rho12*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho22*u*d2^2 + d1*rho12*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho12*u*d1^2 + d2*rho22*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S22*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho12*u*d1^2 + d2*rho22*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho22*u*d2^2 + d1*rho12*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho22*u*d2^2 + d1*rho12*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho12*u*d1^2 + d2*rho22*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) - rho21*(d2^2*u*(S22*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho12*u*d1^2 + d2*rho22*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho22*u*d2^2 + d1*rho12*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho22*u*d2^2 + d1*rho12*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho12*u*d1^2 + d2*rho22*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S21*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho12*u*d1^2 + d2*rho22*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho22*u*d2^2 + d1*rho12*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho22*u*d2^2 + d1*rho12*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho12*u*d1^2 + d2*rho22*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) + d1*d2*u)) + Z2*(Z2*(rho12*(rho12*u*d1^2 + d2*rho22*u*d1) + rho22*(rho22*u*d2^2 + d1*rho12*u*d2) - rho12*(d1^2*u*(S21*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho12*u*d1^2 + d2*rho22*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho22*u*d2^2 + d1*rho12*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho22*u*d2^2 + d1*rho12*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho12*u*d1^2 + d2*rho22*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S22*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho12*u*d1^2 + d2*rho22*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho22*u*d2^2 + d1*rho12*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho22*u*d2^2 + d1*rho12*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho12*u*d1^2 + d2*rho22*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) - rho22*(d2^2*u*(S22*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho12*u*d1^2 + d2*rho22*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho22*u*d2^2 + d1*rho12*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho22*u*d2^2 + d1*rho12*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho12*u*d1^2 + d2*rho22*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S21*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho12*u*d1^2 + d2*rho22*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho22*u*d2^2 + d1*rho12*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho22*u*d2^2 + d1*rho12*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho12*u*d1^2 + d2*rho22*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) + d2^2*u) + Z1*(rho12*(rho11*u*d1^2 + d2*rho21*u*d1) + rho22*(rho21*u*d2^2 + d1*rho11*u*d2) - rho12*(d1^2*u*(S21*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho11*u*d1^2 + d2*rho21*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho21*u*d2^2 + d1*rho11*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho21*u*d2^2 + d1*rho11*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho11*u*d1^2 + d2*rho21*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S22*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho11*u*d1^2 + d2*rho21*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho21*u*d2^2 + d1*rho11*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho21*u*d2^2 + d1*rho11*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho11*u*d1^2 + d2*rho21*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) - rho22*(d2^2*u*(S22*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho11*u*d1^2 + d2*rho21*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho21*u*d2^2 + d1*rho11*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S12*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho21*u*d2^2 + d1*rho11*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho11*u*d1^2 + d2*rho21*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1))) + d1*d2*u*(S21*((4*(S12*u*d2^2 + S11*d1*u*d2)*(rho11*u*d1^2 + d2*rho21*u*d1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho21*u*d2^2 + d1*rho11*u*d2)*(2*S11*u*d1^2 + 2*S12*d2*u*d1 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)) + S11*((4*(S21*u*d1^2 + S22*d2*u*d1)*(rho21*u*d2^2 + d1*rho11*u*d2))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1) - (2*(rho11*u*d1^2 + d2*rho21*u*d1)*(2*S22*u*d2^2 + 2*S21*d1*u*d2 - 1))/(2*S11*d1^2*u + 2*S22*d2^2*u + 2*S12*d1*d2*u + 2*S21*d1*d2*u - 1)))) + d1*d2*u)) + 2*alpha*d1*mu1*u + 2*alpha*d2*mu2*u)
原因在 this question
的答案中有很好的描述在您的情况下，由于您在文件中获得了表达式，因此您可以 source() 来自 R 的文件，而不是复制和粘贴。您需要在文本文件中添加作业，或者您可以这样做：
source("filename.txt")$value
获取计算结果
您在评论中提到您希望能够在优化中使用它 - 一个在一维上找到最小化表达式的 alpha 的简单示例是：
f <- function(alpha){source("o:\test.txt")$value}
optim(1, f, lower=0.5, upper=1.5, method="L-BFGS-B")
#$par
#[1] 1
#
#$value
#[1] 1.011915
#
#$counts
#function gradient
# 1 1
#
#$convergence
#[1] 0
#
#$message
#[1] "CONVERGENCE: NORM OF PROJECTED GRADIENT <= PGTOL"