偏相关_netcdf

问题描述

我正在尝试使用 r 中的 ppcor 包计算单个 NetCDF 中 6 个（观察到的和建模的数据）变量之间的偏相关，但我还没有成功。我还在学习循环。这是我到目前为止尝试过的：

library(ppcor)
mods <- nc_open("Obs_model_1981-2018.nc")
lon <- ncvar_get(mods,varid = "lon")
lat <- ncvar_get(mods,varid = "lat")

nx <- length(lon)
ny <- length(lat)
    obss <- ncvar_get(mods,"obs")
    mods1 <- ncvar_get(mods,"mod1")
    mods2 <- ncvar_get(mods,"mod2")
    mods3 <- ncvar_get(mods,"mod3")
    mods4 <- ncvar_get(mods,"mod4")
    mods5 <- ncvar_get(mods,"mod5")

rl=array(0,c(nx,ny))
err=array(0,ny))

for (x in 1:nx) {
  print(x)
  for (y in 1:ny) {
    pp =data.frame(obss,mods1,mods2,mods3,mods4,mods5)         
    ## pcor need a data frame
    rl[x,y] = pcor(pp,"pearson")
   }
}

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）