如何编写x^.3*(1-x)^.7函数？R中x>1返回NaN问题求解

首先得搞清楚为什么testFunc(2)会返回NaN：当x>1时，1-x变成了负数，而负数的非整数次幂在实数范围内是没有定义的——比如你没法算出一个实数的10次方等于-1（因为实数的偶次幂都是非负的），而0.7是7/10，相当于先开10次方再取7次幂，这在实数域里做不到，所以R返回NaN表示“不是一个实数”。

根据你的需求，有两种解决方向：

1. 只需要实数结果（推荐，若为统计场景） #

原函数其实是Beta分布的核函数，它的自然有效定义域就是[0,1]。你可以在函数里加入输入检查，当x超出这个范围时给出明确提示，而不是返回模糊的NaN：

testFunc <- function(x){
  # 检查x是否在有效实数定义域内
  if (x < 0 || x > 1) {
    warning("x must be between 0 and 1 to get a real-valued result!")
    return(NA_real_)
  }
  x^0.3 * (1 - x)^0.7
}

调用testFunc(2)时，会弹出警告并返回NA，比NaN更直观；而testFunc(0.5)会正常返回实数结果：[1] 0.4781493。

2. 需要复数结果 #

如果你确实需要x>1时的计算结果（比如复数域的数学运算），可以把输入强制转换为复数类型，这样R会计算复数域的结果：

testFunc <- function(x){
  # 转为复数后计算
  as.complex(x)^0.3 * as.complex(1 - x)^0.7
}

这时候调用testFunc(2)会得到复数结果：[1] 0.6178161+0.7864059i，这个结果在复数域是合法的。

内容的提问来源于stack exchange，提问作者Iamsomean