如何使用For循环在R中插入这个sigma方程？

用于计算两个样本（每个样本有128个观测值）之间的系数的方程式：

我需要用这个公式计算128个观测值的两个像素之间的系数。我的想法是在R中使用for循环

方程式为：

Rsi，Sj=西格玛N（si-Sj）-西格玛N（si）西格玛N（Sj）/sqrt[（西格玛N（si^2-（西格玛N-si）^2][N西格玛Sj^2-（西格玛Sj）^2]

我对编码非常陌生，因此不知道如何解决这个问题。

---

提前谢谢。

我尝试添加方程式图像，但网站不允许我添加。我真的不知道如何在R中添加sigma符号，所以我有点被卡住了。我尝试添加方程式图像，但网站不允许我添加。我真的不知道如何在R中添加sigma符号，所以我有点被卡住了。

My dataset: 2 pixels x 128 observations (spectral bands) which i have in two lists
si = value of the reference pixel
sj = value of the target pixel
N = number of bands (128)