Question

我正在使用以下代码检查线性趋势的P值，但由于我看不到P值的二维图，只能看到一行，因此循环似乎无法正常工作

library(chron)
library(RColorBrewer)
library(lattice)
library(ncdf4)
#-------------------------------------------------------------------------------------------
options(warn=-1)

ncin <- nc_open("MOD04_10K_Winter.nc", readunlim=FALSE)
#print(ncin)
lon <- ncvar_get(ncin, varid="Longitude", start=NA, count=NA, verbose=FALSE,
          signedbyte=TRUE, collapse_degen=TRUE, raw_datavals=FALSE )
lat <- ncvar_get(ncin, varid="Latitude", start=NA, count=NA, verbose=FALSE,
                signedbyte=TRUE, collapse_degen=TRUE, raw_datavals=FALSE )
aod <- ncvar_get(ncin, varid="AOD", start=NA, count=NA, verbose=FALSE,
                signedbyte=TRUE, collapse_degen=TRUE, raw_datavals=FALSE )

px           <- matrix(nrow = 1:length(lon), ncol = 1:length(lat))
is.matrix(px)

for (lo in 1:length(lon)) {
  for (la in 1:length(lat)) {
   int1a    = aod[lo, la,]

   # if mean of int is finite then proceed else fill NA to all arrays
   mn = mean(int1a, trim = 0, na.rm = FALSE)
   if (is.finite(mn)) 
    {
    print("---------------- Reading Finite data -------------")
    xs       = 1:30                                     
    fn1a     = lm(int1a~xs)       # Function_NCP
    p_val    = summary(fn1a)$coefficients[2, 4]          # Saving p-value
    if (p_val < 0.05) {print("statisticlly significant")} else {print("statisticlly in-significant")}
    print(p_val)
    print(lo)
    print(la)
    px[lo][la] = p_val                                   # variables in [] only (?)
   } 

  }                                                      # latitude dimension
}

如果我使用的是[lo, la]而不是[lo][la]，则出现以下错误

[<-（*tmp*，lo，la，值= 0.0543481042240582）中的错误：
下标超出范围

抱歉，如果解决方案非常琐碎，我才刚开始在R中工作。

Answer 1

您只需对矩阵ssh git@gitlab.com声明进行一些小修改。现在，您将行和列的数量设置为向量：px和nrow = 1:length(lon)。 R静默地仅获取这些向量的第一个元素，并生成一对一的矩阵。（实际上，它会生成警告，因为警告被禁止了！）

所以，解决方法是

nrow = 1:length(lon)

无法将数据存储在矩阵中

1 个答案: