我需要在数据框中创建一个新变量,该变量是带有许多嵌套括号的方程式的输出。此等式的一部分以下面最后一行的形式]
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 ) )