r - 如何将多个嵌套括号与指数一起使用？

Question

我需要在数据框中创建一个新变量，它是带有许多嵌套括号的方程的输出。该方程的一部分采用下面最后一行的形式

temp=36  
Tc = 647.097    
( ( 1-273.15+temp )/Tc )^1.5

其中 temp 是一个变量，Tc 是一个常量。但是，当我运行代码时，结果始终为 NA。

但是，如果我将代码分解为我知道结果的数字

( 1-273.15+temp )/Tc

然后像这样添加指数

-0.3649376^1.5

然后代码就可以正常工作了。

为什么R无法正确输出计算结果( ( 1-273.15+temp )/Tc )^1.5 ？

更重要的是，如何让 R 给出 ( ( 1-273.15+temp )/Tc )^1.5 的结果，同时保留对常量和变量的对象的使用？

我需要解决这个问题，因为完整的方程甚至更糟糕，我上面描述的问题嵌套在其自身中:

e_sat_test <- Pc^( ( Tc/(273.15+temp ) ) * 
  ( a1*( (1-273.15+temp)/Tc ) + a2*( (1-273.15+temp)/Tc )^1.5 +
     a3*( (1-273.15+temp)/Tc )^3 + a4* ( (1-273.15+temp)/Tc )^3.5 +
    a5*( (1-273.15+temp)/Tc)^4 + a6*( (1-273.15+temp)/Tc )^7.5  ) )

Answer 1

问题在于

-0.3649376^1.5 

被解释为

-(0.3649376^1.5)

不是

(-0.3649376)^1.5

因为指数运算符具有更高的优先级。当你将某个值取为 0.5 指数时，就像取平方根一样，而这些指数并不是为 R 中的简单数值向量定义的(除非你想使用虚数)。您的数据值的计算结果只是 NaN，因为您的结果不是真实的。您可能想再次检查您的公式，