R在表达式中添加换行符答案

【问题标题】：R add line breaks to expressionsR在表达式中添加换行符
【发布时间】：2021-06-07 03:35:46
【问题描述】：

我在.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)

【问题讨论】：

标签： r line-breaks large-data

【解决方案1】：

this question的回答很好地描述了这种情况的原因

在您的情况下，由于您在文件中获得了表达式，因此您可以 source() R 中的文件，而不是复制和粘贴。您需要在文本文件中添加作业，或者您可以这样做：

source("filename.txt")$value

得到计算结果

您在 cmets 中提到您希望能够在优化中使用它 - 一个简单的一维示例以找到最小化表达式的 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"

【讨论】：

谢谢。如果您对最终值感兴趣，这似乎是一种合理的方法。如果您也对参数感兴趣，它似乎不是特别适合在优化的上下文中。
@salomon 希望编辑显示您如何在优化中使用它。也许如果需要更多细节，您可以提出一个新问题，因为这与原始问题的主题有点偏离